Network crosstalk dynamically changes during neutrophil polarization

Uncovering the role of cytoskeleton proteins in the formation of neutrophil extracellular traps

Neutrophils are a type of lymphocyte involved in innate immune defense. In response to specific stimuli, these phagocytic cells undergo a unique form of cell death, NETosis, during which they release neutrophil extracellular traps (NETs) composed of modified chromatin structures decorated with cytoplasmic and granular proteins. Multiple proteins and pathways have been implicated in the formation of NETs. The cytoskeleton, an interconnected network of filamentous polymers and regulatory proteins, plays a crucial role in resisting deformation, transporting intracellular cargo, and changing shape during movement of eukaryotic cells. It may also have evolved to defend eukaryotic organisms against infection. Recent research focuses on understanding the mechanisms underlying NETs formation and how cytoskeletal networks contribute to this process, by identifying enzymes that trigger NETosis or interact with NETs and influence cellular behavior through cytoskeletal dynamics. An enhanced understanding of the complex relationship between the cytoskeleton and NET formation will provide a framework for future research and the development of targeted therapeutic strategies, and supports the notion that the long-lived cytoskeleton structures may have a lasting impact on this area of research.

Local Myo9b RhoGAP activity regulates cell motility

To migrate, cells assume a polarized morphology, extending forward with a leading edge with their trailing edge retracting back toward the cell body. Both cell extension and retraction critically depend on the organization and dynamics of the actin cytoskeleton, and the small, monomeric GTPases Rac and Rho are important regulators of actin. Activation of Rac induces actin polymerization and cell extension, whereas activation of Rho enhances acto-myosin II contractility and cell retraction. To coordinate migration, these processes must be carefully regulated. The myosin Myo9b, a Rho GTPase-activating protein (GAP), negatively regulates Rho activity and deletion of Myo9b in leukocytes impairs cell migration through increased Rho activity. However, it is not known whether cell motility is regulated by global or local inhibition of Rho activity by Myo9b. Here, we addressed this question by using Myo9b-deficient macrophage-like cells that expressed different recombinant Myo9b constructs. We found that Myo9b accumulates in lamellipodial extensions generated by Rac-induced actin polymerization as a function of its motor activity. Deletion of Myo9b in HL-60-derived macrophages altered cell morphology and impaired cell migration. Reintroduction of Myo9b or Myo9b motor and GAP mutants revealed that local GAP activity rescues cell morphology and migration. In summary, Rac activation leads to actin polymerization and recruitment of Myo9b, which locally inhibits Rho activity to enhance directional cell migration.

A multi-modal data resource for investigating topographic heterogeneity in patient-derived xenograft tumors

Patient-derived xenografts (PDXs) are an essential pre-clinical resource for investigating tumor biology. However, cellular heterogeneity within and across PDX tumors can strongly impact the interpretation of PDX studies. Here, we generated a multi-modal, large-scale dataset to investigate PDX heterogeneity in metastatic colorectal cancer (CRC) across tumor models, spatial scales and genomic, transcriptomic, proteomic and imaging assay modalities. To showcase this dataset, we present analysis to assess sources of PDX variation, including anatomical orientation within the implanted tumor, mouse contribution, and differences between replicate PDX tumors. A unique aspect of our dataset is deep characterization of intra-tumor heterogeneity via immunofluorescence imaging, which enables investigation of variation across multiple spatial scales, from subcellular to whole tumor levels. Our study provides a benchmark data resource to investigate PDX models of metastatic CRC and serves as a template for future, quantitative investigations of spatial heterogeneity within and across PDX tumor models.

Chemokine-triggered microtubule polymerization promotes neutrophil chemotaxis and invasion but not transendothelial migration

Microtubules (MTs) are critically involved in the transport of material within cells, but their roles in chemotactic leukocyte motility and effector functions are still obscure. Resting neutrophils contain few MTs assembled in an MT organizing center (MTOC) behind their multilobular nuclei. Using a probe of real-time tubulin polymerization, SiR-tubulin, we found that neutrophils elongated their MTs within minutes in response to signals from the two prototypic chemotactic peptides, CXCL1 and fMLP. Taxol, a beta-tubulin binding and MT stabilizing drug, was found to abolish this CXCL1- and fMLP-stimulated MT polymerization. Nevertheless, taxol treatment as well as disruption of existing and de novo generated MTs did not impair neutrophil protrusion and squeezing through IL-1β-stimulated endothelial monolayers mediated by endothelial deposited CXCL1 and neutrophil CXCR2. Notably, CXCL1-dependent neutrophil TEM was not associated with neutrophil MT polymerization. Chemokinetic neutrophil motility on immobilized CXCL1 was also not associated with MT polymerization, and taxol treatment did not interfere with this motility. Nevertheless, and consistent with its ability to suppress MT polymerization induced by soluble CXCL1 and fMLP, taxol treatment inhibited neutrophil chemotaxis toward both chemotactic peptides. Taxol treatment also suppressed CXCL1- and fMLP-triggered elastase-dependent neutrophil invasion through collagen I barriers. Collectively, our results highlight de novo chemoattractant-triggered MT polymerization as key for neutrophil chemotaxis and elastase-dependent invasion but not for chemotactic neutrophil crossing of inflamed endothelial barriers.

Different Isoforms of the Neuronal Guidance Molecule Slit2 Directly Cause Chemoattraction or Chemorepulsion of Human Neutrophils

The movement of neutrophils between blood and tissues appears to be regulated by chemoattractants and chemorepellents. Compared with neutrophil chemoattractants, relatively little is known about neutrophil chemorepellents. Slit proteins are endogenously cleaved into a variety of N- and C-terminal fragments, and these fragments are neuronal chemorepellents and inhibit chemoattraction of many cell types, including neutrophils. In this report, we show that the ∼140-kDa N-terminal Slit2 fragment (Slit2-N) is a chemoattractant and the ∼110-kDa N-terminal Slit2 fragment (Slit2-S) is a chemorepellent for human neutrophils. The effects of both Slit2 fragments were blocked by Abs to the Slit2 receptor Roundabout homolog 1 or the Slit2 coreceptor Syndecan-4. Slit2-N did not appear to activate Ras but increased phosphatidylinositol 3,4,5-triphosphate levels. Slit2-N-induced chemoattraction was unaffected by Ras inhibitors, reversed by PI3K inhibitors, and blocked by Cdc42 and Rac inhibitors. In contrast, Slit2-S activated Ras but did not increase phosphatidylinositol 3,4,5-triphosphate levels. Slit2-S-induced chemorepulsion was blocked by Ras and Rac inhibitors, not affected by PI3K inhibitors, and reversed by Cdc42 inhibitors. Slit2-N, but not Slit2-S, increased neutrophil adhesion, myosin L chain 2 phosphorylation, and polarized actin formation and single pseudopods at the leading edge of cells. Slit2-S induced multiple pseudopods. These data suggest that Slit2 isoforms use similar receptors but different intracellular signaling pathways and have different effects on the cytoskeleton and pseudopods to induce neutrophil chemoattraction or chemorepulsion.

Clathrin-mediated endocytosis regulates fMLP-mediated neutrophil polarization

A cell's ability to establish polarization is one of the key steps in directional migration. Upon the addition of a chemoattractant, N-formylmethionyl-leucyl-phenylalanine (fMLP), neutrophils rapidly develop a front end marked by a wide and dense actin network which is a feature of cell polarization. Despite a general understanding of bi-directional crosstalk between endocytosis and polarization, it remains unclear how clathrin-mediated endocytosis (CME) induced by chemoattractant binding to formyl peptide receptor (FPR) affects neutrophil polarization. In this work, we characterized the spatial organization of FPR and clathrin-coated pits (CCPs), the functional unit of CME, with and without fMLP and found that fMLP induced different distributions of FPR and CCPs. We further found that cells had impaired polarization induced by fMLP when CME is inhibited by small molecule inhibitors. Under these conditions, pERK, pAkt₃₀₈, and pAkt₄₇₃ were all severely blocked or had altered dynamics. The spatial organization between actin and two major clathrin-mediated endocytic proteins, clathrin and β-arrestin, were distinct and supported clathrin and β-arrestin's functional roles in mediating neutrophil polarization. Together these results suggest that CME plays a pivotal role in a complex process such as cell polarization.

Chemoattractant concentration-dependent tuning of ERK signaling dynamics in migrating neutrophils

The directed migration (chemotaxis) of neutrophils toward the bacterial peptide N-formyl-Met-Leu-Phe (fMLP) is a crucial process in immune defense against invading bacteria. While navigating through a gradient of increasing concentrations of fMLP, neutrophils and neutrophil-like HL-60 cells switch from exhibiting directional migration at low fMLP concentrations to exhibiting circuitous migration at high fMLP concentrations. The extracellular signal-regulated kinase (ERK) pathway is implicated in balancing this fMLP concentration-dependent switch in migration modes. We investigated the role and regulation of ERK signaling through single-cell analysis of neutrophil migration in response to different fMLP concentrations over time. We found that ERK exhibited gradated, rather than all-or-none, responses to fMLP concentration. Maximal ERK activation occurred in response to about 100 nM fMLP, and ERK inactivation was promoted by p38. Furthermore, we found that directional migration of neutrophils reached a maximal extent at about 100 nM fMLP and that ERK, but not p38, was required for neutrophil migration. Thus, our data suggest that, in chemotactic neutrophils responding to fMLP, ERK displays gradated activation and p38-dependent inhibition and that these ERK dynamics promote neutrophil migration.

Examining Crosstalk among Transforming Growth Factor β, Bone Morphogenetic Protein, and Wnt Pathways

The integration of morphogenic signals by cells is not well understood. A growing body of literature suggests increasingly complex coupling among classically defined pathways. Given this apparent complexity, it is difficult to predict where, when, or even whether crosstalk occurs. Here, we investigated pairs of morphogenic pathways, previously reported to have multiple points of crosstalk, which either do not share (TGFβ and Wnt/β-catenin) or share (TGFβ and bone morphogenetic protein (BMP)) core signaling components. Crosstalk was measured by the ability of one morphogenic pathway to cross-activate core transcription factors and/or target genes of another morphogenic pathway. In contrast to previous studies, we found a surprising absence of crosstalk between TGFβ and Wnt/β-catenin. Further, we did not observe expected cross-pathway inhibition in between TGFβ and BMP, despite the fact that both use (or could compete) for the shared component SMAD4. Critical to our assays was a separation of timescales, which helped separate crosstalk due to initial signal transduction from subsequent post-transcriptional feedback events. Our study revealed fewer (and different) inter-morphogenic pathway crosstalk connections than expected; even pathways that share components can be insulated from one another.

Transcriptional Control of Developmental Cell Behaviors

Tissue-specific transcription regulators emerged as key developmental control genes, which operate in the context of complex gene regulatory networks (GRNs) to coordinate progressive cell fate specification and tissue morphogenesis. We discuss how GRNs control the individual cell behaviors underlying complex morphogenetic events. Cell behaviors classically range from mesenchymal cell motility to cell shape changes in epithelial sheets. These behaviors emerge from the tissue-specific, multiscale integration of the local activities of universal and pleiotropic effectors, which underlie modular subcellular processes including cytoskeletal dynamics, cell-cell and cell-matrix adhesion, signaling, polarity, and vesicle trafficking. Extrinsic cues and intrinsic cell competence determine the subcellular spatiotemporal patterns of effector activities. GRNs influence most subcellular activities by controlling only a fraction of the effector-coding genes, which we argue is enriched in effectors involved in reading and processing the extrinsic cues to contextualize intrinsic subcellular processes and canalize developmental cell behaviors. The properties of the transcription-cell behavior interface have profound implications for evolution and disease.

Deficiency in Toll-interacting protein (Tollip) skews inflamed yet incompetent innate leukocytes in vivo during DSS-induced septic colitis

Functionally compromised neutrophils contribute to adverse clinical outcomes in patients with severe inflammation and injury such as colitis and sepsis. However, the ontogeny of dysfunctional neutrophil during septic colitis remain poorly understood. We report that the dysfunctional neutrophil may be derived by the suppression of Toll-interacting-protein (Tollip). We observed that Tollip deficient neutrophils had compromised migratory capacity toward bacterial product fMLF due to reduced activity of AKT and reduction of FPR2, reduced potential to generate bacterial-killing neutrophil extra-cellular trap (NET), and compromised bacterial killing activity. On the other hand, Tollip deficient neutrophils had elevated levels of CCR5, responsible for their homing to sterile inflamed tissues. The inflamed and incompetent neutrophil phenotype was also observed in vivo in Tollip deficient mice subjected to DSS-induced colitis. We observed that TUDCA, a compound capable of restoring Tollip cellular function, can potently alleviate the severity of DSS-induced colitis. In humans, we observed significantly reduced Tollip levels in peripheral blood collected from human colitis patients as compared to blood samples from healthy donors. Collectively, our data reveal a novel mechanism in Tollip alteration that underlies the inflamed and incompetent polarization of neutrophils leading to severe outcomes of colitis.

PHOCOS: inferring multi-feature phenotypic crosstalk networks

Motivation: Quantification of cellular changes to perturbations can provide a powerful approach to infer crosstalk among molecular components in biological networks. Existing crosstalk inference methods conduct network-structure learning based on a single phenotypic feature (e.g. abundance) of a biomarker. These approaches are insufficient for analyzing perturbation data that can contain information about multiple features (e.g. abundance, activity or localization) of each biomarker.

Results: We propose a computational framework for inferring phenotypic crosstalk (PHOCOS) that is suitable for high-content microscopy or other modalities that capture multiple phenotypes per biomarker. PHOCOS uses a robust graph-learning paradigm to predict direct effects from potential indirect effects and identify errors owing to noise or missing links. The result is a multi-feature, sparse network that parsimoniously captures direct and strong interactions across phenotypic attributes of multiple biomarkers. We use simulated and biological data to demonstrate the ability of PHOCOS to recover multi-attribute crosstalk networks from cellular perturbation assays.

Availability and implementation: PHOCOS is available in open source at https://github.com/AltschulerWu-Lab/PHOCOS

Contact:steven.altschuler@ucsf.edu or lani.wu@ucsf.edu

An analysis toolbox to explore mesenchymal migration heterogeneity reveals adaptive switching between distinct modes

Mesenchymal (lamellipodial) migration is heterogeneous, although whether this reflects progressive variability or discrete, 'switchable' migration modalities, remains unclear. We present an analytical toolbox, based on quantitative single-cell imaging data, to interrogate this heterogeneity. Integrating supervised behavioral classification with multivariate analyses of cell motion, membrane dynamics, cell-matrix adhesion status and F-actin organization, this toolbox here enables the detection and characterization of two quantitatively distinct mesenchymal migration modes, termed 'Continuous' and 'Discontinuous'. Quantitative mode comparisons reveal differences in cell motion, spatiotemporal coordination of membrane protrusion/retraction, and how cells within each mode reorganize with changed cell speed. These modes thus represent distinctive migratory strategies. Additional analyses illuminate the macromolecular- and cellular-scale effects of molecular targeting (fibronectin, talin, ROCK), including 'adaptive switching' between Continuous (favored at high adhesion/full contraction) and Discontinuous (low adhesion/inhibited contraction) modes. Overall, this analytical toolbox now facilitates the exploration of both spontaneous and adaptive heterogeneity in mesenchymal migration.

Directional memory arises from long-lived cytoskeletal asymmetries in polarized chemotactic cells

Chemotaxis, the directional migration of cells in a chemical gradient, is robust to fluctuations associated with low chemical concentrations and dynamically changing gradients as well as high saturating chemical concentrations. Although a number of reports have identified cellular behavior consistent with a directional memory that could account for behavior in these complex environments, the quantitative and molecular details of such a memory process remain unknown. Using microfluidics to confine cellular motion to a 1D channel and control chemoattractant exposure, we observed directional memory in chemotactic neutrophil-like cells. We modeled this directional memory as a long-lived intracellular asymmetry that decays slower than observed membrane phospholipid signaling. Measurements of intracellular dynamics revealed that moesin at the cell rear is a long-lived element that when inhibited, results in a reduction of memory. Inhibition of ROCK (Rho-associated protein kinase), downstream of RhoA (Ras homolog gene family, member A), stabilized moesin and directional memory while depolymerization of microtubules (MTs) disoriented moesin deposition and also reduced directional memory. Our study reveals that long-lived polarized cytoskeletal structures, specifically moesin, actomyosin, and MTs, provide a directional memory in neutrophil-like cells even as they respond on short time scales to external chemical cues.

Time series modeling of live-cell shape dynamics for image-based phenotypic profiling

Live-cell imaging can be used to capture spatio-temporal aspects of cellular responses that are not accessible to fixed-cell imaging. As the use of live-cell imaging continues to increase, new computational procedures are needed to characterize and classify the temporal dynamics of individual cells. For this purpose, here we present the general experimental-computational framework SAPHIRE (Stochastic Annotation of Phenotypic Individual-cell Responses) to characterize phenotypic cellular responses from time series imaging datasets. Hidden Markov modeling is used to infer and annotate morphological state and state-switching properties from image-derived cell shape measurements. Time series modeling is performed on each cell individually, making the approach broadly useful for analyzing asynchronous cell populations. Two-color fluorescent cells simultaneously expressing actin and nuclear reporters enabled us to profile temporal changes in cell shape following pharmacological inhibition of cytoskeleton-regulatory signaling pathways. Results are compared with existing approaches conventionally applied to fixed-cell imaging datasets, and indicate that time series modeling captures heterogeneous dynamic cellular responses that can improve drug classification and offer additional important insight into mechanisms of drug action. The software is available at http://saphire-hcs.org.

mTOR and differential activation of mitochondria orchestrate neutrophil chemotaxis

Neutrophil chemotaxis is regulated by opposing autocrine purinergic signaling mechanisms, which are stimulated by mitochondrial ATP formation that is up-regulated via mTOR and P2Y2 receptors at the front and down-regulated via A2a receptors and cAMP/PKA signaling at the back of cells.

Neutrophils use chemotaxis to locate invading bacteria. Adenosine triphosphate (ATP) release and autocrine purinergic signaling via P2Y2 receptors at the front and A2a receptors at the back of cells regulate chemotaxis. Here, we examined the intracellular mechanisms that control these opposing signaling mechanisms. We found that mitochondria deliver ATP that stimulates P2Y2 receptors in response to chemotactic cues, and that P2Y2 receptors promote mTOR signaling, which augments mitochondrial activity near the front of cells. Blocking mTOR signaling with rapamycin or PP242 or mitochondrial ATP production (e.g., with CCCP) reduced mitochondrial Ca²⁺ uptake and membrane potential, and impaired cellular ATP release and neutrophil chemotaxis. Autocrine stimulation of A2a receptors causes cyclic adenosine monophosphate accumulation at the back of cells, which inhibits mTOR signaling and mitochondrial activity, resulting in uropod retraction. We conclude that mitochondrial, purinergic, and mTOR signaling regulates neutrophil chemotaxis and may be a pharmacological target in inflammatory diseases.

Intracellular Information Processing through Encoding and Decoding of Dynamic Signaling Features

Cell signaling dynamics and transcriptional regulatory activities are variable within specific cell types responding to an identical stimulus. In addition to studying the network interactions, there is much interest in utilizing single cell scale data to elucidate the non-random aspects of the variability involved in cellular decision making. Previous studies have considered the information transfer between the signaling and transcriptional domains based on an instantaneous relationship between the molecular activities. These studies predict a limited binary on/off encoding mechanism which underestimates the complexity of biological information processing, and hence the utility of single cell resolution data. Here we pursue a novel strategy that reformulates the information transfer problem as involving dynamic features of signaling rather than molecular abundances. We pursue a computational approach to test if and how the transcriptional regulatory activity patterns can be informative of the temporal history of signaling. Our analysis reveals (1) the dynamic features of signaling that significantly alter transcriptional regulatory patterns (encoding), and (2) the temporal history of signaling that can be inferred from single cell scale snapshots of transcriptional activity (decoding). Immediate early gene expression patterns were informative of signaling peak retention kinetics, whereas transcription factor activity patterns were informative of activation and deactivation kinetics of signaling. Moreover, the information processing aspects varied across the network, with each component encoding a selective subset of the dynamic signaling features. We developed novel sensitivity and information transfer maps to unravel the dynamic multiplexing of signaling features at each of these network components. Unsupervised clustering of the maps revealed two groups that aligned with network motifs distinguished by transcriptional feedforward vs feedback interactions. Our new computational methodology impacts the single cell scale experiments by identifying downstream snapshot measures required for inferring specific dynamical features of upstream signals involved in the regulation of cellular responses.

Single cell studies have shown that differential patterns in the dynamics of signaling proteins, transcription factor activity, gene expression, etc. produce distinct downstream outcomes. The opposite also holds true where particular cellular outcomes have been found to be associated with the dynamical pattern of one or more signaling molecules. Signaling pathways, therefore, serve as signal processing units to inform specific downstream regulation. However, the functional capabilities of the dynamic aspects of signaling are not well understood. To address this issue, we developed a new approach that evaluates information processing between dynamic features in signaling patterns and transcriptional regulatory activity. Our work demonstrates that the information transfer occur through decoding of temporal history of signals rather than only through instantaneous correlations. Moreover, our results identify regulatory network motifs as the critical components in the information processing and filtering of variability in signaling dynamics to produce distinct patterns of downstream transcriptional responses. Our methodology can be broadly applied to single cell scale data on experimentally accessible downstream measures to infer dynamic aspects of upstream signaling.

Disentangling Membrane Dynamics and Cell Migration; Differential Influences of F-actin and Cell-Matrix Adhesions

Cell migration is heavily interconnected with plasma membrane protrusion and retraction (collectively termed “membrane dynamics”). This makes it difficult to distinguish regulatory mechanisms that differentially influence migration and membrane dynamics. Yet such distinctions may be valuable given evidence that cancer cell invasion in 3D may be better predicted by 2D membrane dynamics than by 2D cell migration, implying a degree of functional independence between these processes. Here, we applied multi-scale single cell imaging and a systematic statistical approach to disentangle regulatory associations underlying either migration or membrane dynamics. This revealed preferential correlations between membrane dynamics and F-actin features, contrasting with an enrichment of links between cell migration and adhesion complex properties. These correlative linkages were often non-linear and therefore context-dependent, strengthening or weakening with spontaneous heterogeneity in cell behavior. More broadly, we observed that slow moving cells tend to increase in area, while fast moving cells tend to shrink, and that the size of dynamic membrane domains is independent of cell area. Overall, we define macromolecular features preferentially associated with either cell migration or membrane dynamics, enabling more specific interrogation and targeting of these processes in future.

Critical behavior of subcellular density organization during neutrophil activation and migration

Physical theories of active matter continue to provide a quantitative understanding of dynamic cellular phenomena, including cell locomotion. Although various investigations of the rheology of cells have identified important viscoelastic and traction force parameters for use in these theoretical approaches, a key variable has remained elusive both in theoretical and experimental approaches: the spatiotemporal behavior of the subcellular density. The evolution of the subcellular density has been qualitatively observed for decades as it provides the source of image contrast in label-free imaging modalities (e.g., differential interference contrast, phase contrast) used to investigate cellular specimens. While these modalities directly visualize cell structure, they do not provide quantitative access to the structures being visualized. We present an established quantitative imaging approach, non-interferometric quantitative phase microscopy, to elucidate the subcellular density dynamics in neutrophils undergoing chemokinesis following uniform bacterial peptide stimulation. Through this approach, we identify a power law dependence of the neutrophil mean density on time with a critical point, suggesting a critical density is required for motility on 2D substrates. Next we elucidate a continuum law relating mean cell density, area, and total mass that is conserved during neutrophil polarization and migration. Together, our approach and quantitative findings will enable investigators to define the physics coupling cytoskeletal dynamics with subcellular density dynamics during cell migration.

Non-monotonic cellular responses to heterogeneity in talin protein expression-level

Talin is a key cell-matrix adhesion component with a central role in regulating adhesion complex maturation, and thereby various cellular properties including adhesion and migration. However, knockdown studies have produced inconsistent findings regarding the functional influence of talin in these processes. Such discrepancies may reflect non-monotonic responses to talin expression-level variation that are not detectable via canonical "binary" comparisons of aggregated control versus knockdown cell populations. Here, we deployed an "analogue" approach to map talin influence across a continuous expression-level spectrum, which we extended with sub-maximal RNAi-mediated talin depletion. Applying correlative imaging to link live cell and fixed immunofluorescence data on a single cell basis, we related per cell talin levels to per cell measures quantitatively defining an array of cellular properties. This revealed both linear and non-linear correspondences between talin expression and cellular properties, including non-monotonic influences over cell shape, adhesion complex-F-actin association and adhesion localization. Furthermore, we demonstrate talin level-dependent changes in networks of correlations among adhesion/migration properties, particularly in relation to cell migration speed. Importantly, these correlation networks were strongly affected by talin expression heterogeneity within the natural range, implying that this endogenous variation has a broad, quantitatively detectable influence. Overall, we present an accessible analogue method that reveals complex dependencies on talin expression-level, thereby establishing a framework for considering non-linear and non-monotonic effects of protein expression-level heterogeneity in cellular systems.

Intracellular signalling during neutrophil recruitment

Recruitment of leucocytes such as neutrophils to the extravascular space is a critical step of the inflammation process and plays a major role in the development of various diseases including several cardiovascular diseases. Neutrophils themselves play a very active role in that process by sensing their environment and responding to the extracellular cues by adhesion and de-adhesion, cellular shape changes, chemotactic migration, and other effector functions of cell activation. Those responses are co-ordinated by a number of cell surface receptors and their complex intracellular signal transduction pathways. Here, we review neutrophil signal transduction processes critical for recruitment to the site of inflammation. The two key requirements for neutrophil recruitment are the establishment of appropriate chemoattractant gradients and the intrinsic ability of the cells to migrate along those gradients. We will first discuss signalling steps required for sensing extracellular chemoattractants such as chemokines and lipid mediators and the processes (e.g. PI3-kinase pathways) leading to the translation of extracellular chemoattractant gradients to polarized cellular responses. We will then discuss signal transduction by leucocyte adhesion receptors (e.g. tyrosine kinase pathways) which are critical for adhesion to, and migration through the vessel wall. Finally, additional neutrophil signalling pathways with an indirect effect on the neutrophil recruitment process, e.g. through modulation of the inflammatory environment, will be discussed. Mechanistic understanding of these pathways provide better understanding of the inflammation process and may point to novel therapeutic strategies for controlling excessive inflammation during infection or tissue damage.

Homer3 regulates the establishment of neutrophil polarity

Most chemoattractants rely on activation of the heterotrimeric G-protein Gαi to regulate directional cell migration, but few links from Gαi to chemotactic effectors are known. Through affinity chromatography using primary neutrophil lysate, we identify Homer3 as a novel Gαi2-binding protein. RNA interference-mediated knockdown of Homer3 in neutrophil-like HL-60 cells impairs chemotaxis and the establishment of polarity of phosphatidylinositol 3,4,5-triphosphate (PIP3) and the actin cytoskeleton, as well as the persistence of the WAVE2 complex. Most previously characterized proteins that are required for cell polarity are needed for actin assembly or activation of core chemotactic effectors such as the Rac GTPase. In contrast, Homer3-knockdown cells show normal magnitude and kinetics of chemoattractant-induced activation of phosphoinositide 3-kinase and Rac effectors. Chemoattractant-stimulated Homer3-knockdown cells also exhibit a normal initial magnitude of actin polymerization but fail to polarize actin assembly and intracellular PIP3 and are defective in the initiation of cell polarity and motility. Our data suggest that Homer3 acts as a scaffold that spatially organizes actin assembly to support neutrophil polarity and motility downstream of GPCR activation.

Network representations of immune system complexity

The mammalian immune system is a dynamic multiscale system composed of a hierarchically organized set of molecular, cellular, and organismal networks that act in concert to promote effective host defense. These networks range from those involving gene regulatory and protein–protein interactions underlying intracellular signaling pathways and single‐cell responses to increasingly complex networks of in vivo cellular interaction, positioning, and migration that determine the overall immune response of an organism. Immunity is thus not the product of simple signaling events but rather nonlinear behaviors arising from dynamic, feedback‐regulated interactions among many components. One of the major goals of systems immunology is to quantitatively measure these complex multiscale spatial and temporal interactions, permitting development of computational models that can be used to predict responses to perturbation. Recent technological advances permit collection of comprehensive datasets at multiple molecular and cellular levels, while advances in network biology support representation of the relationships of components at each level as physical or functional interaction networks. The latter facilitate effective visualization of patterns and recognition of emergent properties arising from the many interactions of genes, molecules, and cells of the immune system. We illustrate the power of integrating ‘omics’ and network modeling approaches for unbiased reconstruction of signaling and transcriptional networks with a focus on applications involving the innate immune system. We further discuss future possibilities for reconstruction of increasingly complex cellular‐ and organism‐level networks and development of sophisticated computational tools for prediction of emergent immune behavior arising from the concerted action of these networks. WIREs Syst Biol Med 2015, 7:13–38. doi: 10.1002/wsbm.1288

This article is categorized under: 1

Analytical and Computational Methods > Computational Methods

2

Laboratory Methods and Technologies > Macromolecular Interactions, Methods

A neutrophil treadmill to decouple spatial and temporal signals during chemotaxis

After more than 50 years of debates, the role of spatial and temporal gradients during cell chemotaxis is still a contentious matter. One major challenge is that when cells move in response to a heterogeneous chemical environment they are exposed to both spatial and temporal concentration changes. Even in the presence of perfectly stable chemical gradients, moving cells experience temporal changes of concentration simply by moving between locations with different chemical concentrations in a heterogeneous environment. Thus, the effects of the spatial and temporal stimuli cannot be dissociated and studied independently, hampering progress towards understanding the mechanisms of cell chemotaxis. Here we employ microfluidic and other engineering tools to build a system that accomplishes a function analogous to a treadmill at the cellular scale, holding a moving cell at a specified, unchanging location in a chemical gradient. Using this system, we decouple the spatial and temporal gradients around moving human neutrophils and find that temporal gradients are necessary for the directional persistence of human neutrophils during chemotaxis. Our results suggest that temporal chemoattractant changes are important during neutrophil migration and should be taken into account when deciphering the signalling pathways of cell chemotaxis.

On comparing heterogeneity across biomarkers

Microscopy reveals complex patterns of cellular heterogeneity that can be biologically informative. However, a limitation of microscopy is that only a small number of biomarkers can typically be monitored simultaneously. Thus, a natural question is whether additional biomarkers provide a deeper characterization of the distribution of cellular states in a population. How much information about a cell's phenotypic state in one biomarker is gained by knowing its state in another biomarker? Here, we describe a framework for comparing phenotypic states across biomarkers. Our approach overcomes the current limitation of microscopy by not requiring costaining biomarkers on the same cells; instead, we require staining of biomarkers (possibly separately) on a common collection of phenotypically diverse cell lines. We evaluate our approach on two image datasets: 33 oncogenically diverse lung cancer cell lines stained with 7 biomarkers, and 49 less diverse subclones of one lung cancer cell line stained with 12 biomarkers. We first validate our method by comparing it to the "gold standard" of costaining. We then apply our approach to all pairs of biomarkers and use it to identify biomarkers that yield similar patterns of heterogeneity. The results presented in this work suggest that many biomarkers provide redundant information about heterogeneity. Thus, our approach provides a practical guide for selecting independently informative biomarkers and, more generally, will yield insights into both the connectivity of biological networks and the complexity of the state space of biological systems.

Mst1 directs Myosin IIa partitioning of low and higher affinity integrins during T cell migration

Chemokines promote T cell migration by transmitting signals that induce T cell polarization and integrin activation and adhesion. Mst1 kinase is a key signal mediator required for both of these processes; however, its molecular mechanism remains unclear. Here, we present a mouse model in which Mst1 function is disrupted by a hypomorphic mutation. Microscopic analysis of Mst1-deficient CD4 T cells revealed a necessary role for Mst1 in controlling the localization and activity of Myosin IIa, a molecular motor that moves along actin filaments. Using affinity specific LFA-1 antibodies, we identified a requirement for Myosin IIa-dependent contraction in the precise spatial distribution of low and higher affinity LFA-1 on the membrane of migrating T cells. Mst1 deficiency or Myosin inhibition resulted in multipolar cells, difficulties in uropod detachment and mis-localization of low affinity LFA-1. Thus, Mst1 regulates Myosin IIa dynamics to organize high and low affinity LFA-1 to the anterior and posterior membrane during T cell migration.

How to understand and outwit adaptation

Adaptation is the ability of a system to respond and reset itself even in the continuing presence of a stimulus. On one hand, adaptation is a physiological necessity that enables proper neuronal signaling and cell movement. On the other hand, adaptation can be a source of annoyance, as it can make biological systems resistant to experimental perturbations. Here we speculate where adaptation might live in eukaryotic chemotaxis and how it can be encoded in the signaling network. We then discuss tools and strategies that can be used to both understand and outwit adaptation in a wide range of cellular contexts.

Moving towards a paradigm: common mechanisms of chemotactic signaling in Dictyostelium and mammalian leukocytes

Chemotaxis, or directed migration of cells along a chemical gradient, is a highly coordinated process that involves gradient sensing, motility, and polarity. Most of our understanding of chemotaxis comes from studies of cells undergoing amoeboid-type migration, in particular the social amoeba Dictyostelium discoideum and leukocytes. In these amoeboid cells the molecular events leading to directed migration can be conceptually divided into four interacting networks: receptor/G protein, signal transduction, cytoskeleton, and polarity. The signal transduction network occupies a central position in this scheme as it receives direct input from the receptor/G protein network, as well as feedback from the cytoskeletal and polarity networks. Multiple overlapping modules within the signal transduction network transmit the signals to the actin cytoskeleton network leading to biased pseudopod protrusion in the direction of the gradient. The overall architecture of the networks, as well as the individual signaling modules, is remarkably conserved between Dictyostelium and mammalian leukocytes, and the similarities and differences between the two systems are the subject of this review.

Plasticity in the macromolecular-scale causal networks of cell migration

Heterogeneous and dynamic single cell migration behaviours arise from a complex multi-scale signalling network comprising both molecular components and macromolecular modules, among which cell-matrix adhesions and F-actin directly mediate migration. To date, the global wiring architecture characterizing this network remains poorly defined. It is also unclear whether such a wiring pattern may be stable and generalizable to different conditions, or plastic and context dependent. Here, synchronous imaging-based quantification of migration system organization, represented by 87 morphological and dynamic macromolecular module features, and migration system behaviour, i.e., migration speed, facilitated Granger causality analysis. We thereby leveraged natural cellular heterogeneity to begin mapping the directionally specific causal wiring between organizational and behavioural features of the cell migration system. This represents an important advance on commonly used correlative analyses that do not resolve causal directionality. We identified organizational features such as adhesion stability and adhesion F-actin content that, as anticipated, causally influenced cell migration speed. Strikingly, we also found that cell speed can exert causal influence over organizational features, including cell shape and adhesion complex location, thus revealing causality in directions contradictory to previous expectations. Importantly, by comparing unperturbed and signalling-modulated cells, we provide proof-of-principle that causal interaction patterns are in fact plastic and context dependent, rather than stable and generalizable.

Disentangling biological signaling networks by dynamic coupling of signaling lipids to modifying enzymes

An unresolved problem in biological signal transduction is how particular branches of highly interconnected signaling networks can be decoupled, allowing activation of specific circuits within complex signaling architectures. Although signaling dynamics and spatiotemporal mechanisms serve critical roles, it remains unclear if these are the only ways cells achieve specificity within networks. The transcription factor Steroidogenic Factor-1 (SF-1) is an excellent model to address this question, as it forms dynamic complexes with several chemically distinct lipid species (phosphatidylinositols, phosphatidylcholines and sphingolipids). This property is important since lipids bound to SF-1 are modified by lipid signaling enzymes (IPMK & PTEN), regulating SF-1 biological activity in gene expression. Thus, a particular SF-1/lipid complex can interface with a lipid signaling enzyme only if SF-1 has been loaded with a chemically compatible lipid substrate. This mechanism permits dynamic downstream responsiveness to constant upstream input, disentangling specific pathways from the full network. The potential of this paradigm to apply generally to nuclear lipid signaling is discussed, with particular attention given to the nuclear receptor superfamily of transcription factors and their phospholipid ligands.

Independence of symmetry breaking on Bem1-mediated autocatalytic activation of Cdc42

Rather than acting directly on Cdc42, Bem1 works in concert with the Cdc42 binding partner Rdi1 to relocalize Cdc42 to the cytosol during symmetry breaking in the absence of an intact actin cytoskeleton.

The ability to break symmetry and polarize through self-organization is a fundamental feature of cellular systems. A prevailing theory in yeast posits that symmetry breaking occurs via a positive feedback loop, wherein the adaptor protein Bem1 promotes local activation and accumulation of Cdc42 by directly tethering Cdc42^GTP with its guanine nucleotide exchange factor (GEF) Cdc24. In this paper, we find that neither Bem1 nor the ability of Bem1 to bind Cdc42^GTP is required for cell polarization. Instead, Bem1 functions primarily by boosting GEF activity, a role critical for polarization without actin filaments. In the absence of actin-based transport, polarization of Cdc42 is accomplished through Rdi1, the Cdc42 guanine nucleotide dissociation inhibitor. A mathematical model is constructed describing cell polarization as a product of distinct pathways controlling Cdc42 activation and protein localization. The model predicts a nonmonotonic dependence of cell polarization on the concentration of Rdi1 relative to that of Cdc42.

p21-Activated kinase (PAK) regulates cytoskeletal reorganization and directional migration in human neutrophils

Neutrophils serve as a first line of defense in innate immunity owing in part to their ability to rapidly migrate towards chemotactic factors derived from invading pathogens. As a migratory function, neutrophil chemotaxis is regulated by the Rho family of small GTPases. However, the mechanisms by which Rho GTPases orchestrate cytoskeletal dynamics in migrating neutrophils remain ill-defined. In this study, we characterized the role of p21-activated kinase (PAK) downstream of Rho GTPases in cytoskeletal remodeling and chemotactic processes of human neutrophils. We found that PAK activation occurred upon stimulation of neutrophils with f-Met-Leu-Phe (fMLP), and PAK accumulated at the actin-rich leading edge of stimulated neutrophils, suggesting a role for PAK in Rac-dependent actin remodeling. Treatment with the pharmacological PAK inhibitor, PF3758309, abrogated the integrity of RhoA-mediated actomyosin contractility and surface adhesion. Moreover, inhibition of PAK activity impaired neutrophil morphological polarization and directional migration under a gradient of fMLP, and was associated with dysregulated Ca(2+) signaling. These results suggest that PAK serves as an important effector of Rho-family GTPases in neutrophil cytoskeletal reorganization, and plays a key role in driving efficient directional migration of human neutrophils.

Mathematical modeling of eukaryotic cell migration: insights beyond experiments

A migrating cell is a molecular machine made of tens of thousands of short-lived and interacting parts. Understanding migration means understanding the self-organization of these parts into a system of functional units. This task is one of tackling complexity: First, the system integrates numerous chemical and mechanical component processes. Second, these processes are connected in feedback interactions and over a large range of spatial and temporal scales. Third, many processes are stochastic, which leads to heterogeneous migration behaviors. Early on in the research of cell migration it became evident that this complexity exceeds human intuition. Thus, the cell migration community has led the charge to build mathematical models that could integrate the diverse experimental observations and measurements in consistent frameworks, first in conceptual and more recently in molecularly explicit models. The main goal of this review is to sift through a series of important conceptual and explicit mathematical models of cell migration and to evaluate their contribution to the field in their ability to integrate critical experimental data.

Pannexin 1 channels link chemoattractant receptor signaling to local excitation and global inhibition responses at the front and back of polarized neutrophils

Neutrophil chemotaxis requires excitatory signals at the front and inhibitory signals at the back of cells, which regulate cell migration in a chemotactic gradient field. We have previously shown that ATP release via pannexin 1 (PANX1) channels and autocrine stimulation of P2Y2 receptors contribute to the excitatory signals at the front. Here we show that PANX1 also contributes to the inhibitory signals at the back, namely by providing the ligand for A2A adenosine receptors. In resting neutrophils, we found that A2A receptors are uniformly distributed across the cell surface. In polarized cells, A2A receptors redistributed to the back where their stimulation triggered intracellular cAMP accumulation and protein kinase A (PKA) activation, which blocked chemoattractant receptor signaling. Inhibition of PANX1 blocked A2A receptor stimulation and cAMP accumulation in response to formyl peptide receptor stimulation. Treatments that blocked endogenous A2A receptor signaling impaired the polarization and migration of neutrophils in a chemotactic gradient field and resulted in enhanced ERK and p38 MAPK signaling in response to formyl peptide receptor stimulation. These findings suggest that chemoattractant receptors require PANX1 to trigger excitatory and inhibitory signals that synergize to fine-tune chemotactic responses at the front and back of neutrophils. PANX1 channels thus link local excitatory signals to the global inhibitory signals that orchestrate chemotaxis of neutrophils in gradient fields.

Biophysical description of multiple events contributing blood leukocyte arrest on endothelium

Blood leukocytes have a remarkable capacity to bind to and stop on specific blood vessel areas. Many studies have disclosed a key role of integrin structural changes following the interaction of rolling leukocytes with surface-bound chemoattractants. However, the functional significance of structural data and mechanisms of cell arrest are incompletely understood. Recent experiments revealed the unexpected complexity of several key steps of cell-surface interaction: (i) ligand-receptor binding requires a minimum amount of time to proceed and this is influenced by forces. (ii) Also, molecular interactions at interfaces are not fully accounted for by the interaction properties of soluble molecules. (iii) Cell arrest depends on nanoscale topography and mechanical properties of the cell membrane, and these properties are highly dynamic. Here, we summarize these results and we discuss their relevance to recent functional studies of integrin-receptor association in cells from a patient with type III leukocyte adhesion deficiency. It is concluded that an accurate understanding of all physical events listed in this review is needed to unravel the precise role of the multiple molecules and biochemical pathway involved in arrest triggering.

Identifying network motifs that buffer front-to-back signaling in polarized neutrophils

Neutrophil polarity relies on local, mutual inhibition to segregate incompatible signaling circuits to the leading and trailing edges. Mutual inhibition alone should lead to cells having strong fronts and weak backs or vice versa. However, analysis of cell-to-cell variation in human neutrophils revealed that back polarity remains consistent despite changes in front strength. How is this buffering achieved? Pharmacological perturbations and mathematical modeling revealed a functional role for microtubules in buffering back polarity by mediating positive, long-range crosstalk from front to back; loss of microtubules inhibits buffering and results in anticorrelation between front and back signaling. Furthermore, a systematic, computational search of network topologies found that a long-range, positive front-to-back link is necessary for back buffering. Our studies suggest a design principle that can be employed by polarity networks: short-range mutual inhibition establishes distinct signaling regions, after which directed long-range activation insulates one region from variations in the other.

Envisioning migration: mathematics in both experimental analysis and modeling of cell behavior

The complex nature of cell migration highlights the power and challenges of applying mathematics to biological studies. Mathematics may be used to create model equations that recapitulate migration, which can predict phenomena not easily uncovered by experiments or intuition alone. Alternatively, mathematics may be applied to interpreting complex data sets with better resolution--potentially empowering scientists to discern subtle patterns amid the noise and heterogeneity typical of migrating cells. Iteration between these two methods is necessary in order to reveal connections within the cell migration signaling network, as well as to understand the behavior that arises from those connections. Here, we review recent quantitative analysis and mathematical modeling approaches to the cell migration problem.