Nanotopography modulates intracellular excitable systems through cytoskeleton actuation.

Cellular sensing of most environmental cues involves receptors that affect a signal-transduction excitable network (STEN), which is coupled to a cytoskeletal excitable network (CEN). We show that the mechanism of sensing of nanoridges is fundamentally different. CEN activity occurs preferentially on nanoridges, whereas STEN activity is constrained between nanoridges. In the absence of STEN, waves disappear, but long-lasting F-actin puncta persist along the ridges. When CEN is suppressed, wave propagation is no longer constrained by nanoridges. A computational model reproduces these experimental observations. Our findings indicate that nanotopography is sensed directly by CEN, whereas STEN is only indirectly affected due to a CEN-STEN feedback loop. These results explain why texture sensing is robust and acts cooperatively with multiple other guidance cues in complex, in vivo microenvironments.

Contractility, focal adhesion orientation, and stress fiber orientation drive cancer cell polarity and migration along wavy ECM substrates.

Contact guidance is a powerful topographical cue that induces persistent directional cell migration. Healthy tissue stroma is characterized by a meshwork of wavy extracellular matrix (ECM) fiber bundles, whereas metastasis-prone stroma exhibit less wavy, more linear fibers. The latter topography correlates with poor prognosis, whereas more wavy bundles correlate with benign tumors. We designed nanotopographic ECM-coated substrates that mimic collagen fibril waveforms seen in tumors and healthy tissues to determine how these nanotopographies may regulate cancer cell polarization and migration machineries. Cell polarization and directional migration were inhibited by fibril-like wave substrates above a threshold amplitude. Although polarity signals and actin nucleation factors were required for polarization and migration on low-amplitude wave substrates, they did not localize to cell leading edges. Instead, these factors localized to wave peaks, creating multiple "cryptic leading edges" within cells. On high-amplitude wave substrates, retrograde flow from large cryptic leading edges depolarized stress fibers and focal adhesions and inhibited cell migration. On low-amplitude wave substrates, actomyosin contractility overrode the small cryptic leading edges and drove stress fiber and focal adhesion orientation along the wave axis to mediate directional migration. Cancer cells of different intrinsic contractility depolarized at different wave amplitudes, and cell polarization response to wavy substrates could be tuned by manipulating contractility. We propose that ECM fibril waveforms with sufficiently high amplitude around tumors may serve as "cell polarization barriers," decreasing directional migration of tumor cells, which could be overcome by up-regulation of tumor cell contractility.

Elongate mineral particles (EMP) characteristics and mesothelioma: Summary and resolution for session I of the Monticello II conference.

The summary contains a consensus opinion regarding the current state of the science about the dimensions of Elongate Mineral Particles (EMPs) as a factor impacting their carcinogenicity.

Excitable systems: A new perspective on the cellular impact of elongate mineral particles.

We investigate how the geometry of elongate mineral particles (EMPs) in contact with cells influences esotaxis, a recently discovered mechanism of texture sensing. Esotaxis is based on cytoskeletal waves and oscillations that are nucleated, shaped, and steered by the texture of the surroundings. We find that all EMPs studied trigger an esotactic response in macrophages, and that this response dominates cytoskeletal activity in these immune cells. In contrast, epithelial cells show little to no esotactic response to the EMPs. These results are consistent with the distinct interactions of both cell types with ridged nanotopographies of dimensions comparable to those of asbestiform EMPs. Our findings raise the question of whether narrow, asbestiform EMPs may also dominate cytoskeletal activity in other types of immune cells that exhibit similar esotactic effects. These findings, together with prior studies of esotaxis, lead us to the hypothesis that asbestiform EMPs suppress the migration of immune cells and activate immune signaling, thereby outcompeting signals that would normally stimulate the immune system in nearby tissue.

Experimentally Measuring Rolling and Sliding in Three-Dimensional Dense Granular Packings.

We experimentally measure a three-dimensional (3D) granular system's reversibility under cyclic compression. We image the grains using a refractive-index-matched fluid, then analyze the images using the artificial intelligence of variational autoencoders. These techniques allow us to track all the grains' translations and 3D rotations with accuracy sufficient to infer sliding and rolling displacements. Our observations reveal unique roles played by 3D rotational motions in granular flows. We find that rotations and contact-point motion dominate the dynamics in the bulk, far from the perturbation's source. Furthermore, we determine that 3D rotations are irreversible under cyclic compression. Consequently, contact-point sliding, which is dissipative, accumulates throughout the cycle. Using numerical simulations whose accuracy our experiment supports, we discover that much of the dissipation occurs in the bulk, where grains rotate more than they translate. Our observations suggest that the analysis of 3D rotations is needed for understanding granular materials' unique and powerful ability to absorb and dissipate energy.

Signaling through polymerization and degradation: Analysis and simulations of T cell activation mediated by Bcl10.

The adaptive immune system serves as a potent and highly specific defense mechanism against pathogen infection. One component of this system, the effector T cell, facilitates pathogen clearance upon detection of specific antigens by the T cell receptor (TCR). A critical process in effector T cell activation is transmission of signals from the TCR to a key transcriptional regulator, NF-κB. The transmission of this signal involves a highly dynamic process in which helical filaments of Bcl10, a key protein constituent of the TCR signaling cascade, undergo competing processes of polymeric assembly and macroautophagy-dependent degradation. Through computational analysis of three-dimensional, super-resolution optical micrographs, we quantitatively characterize TCR-stimulated Bcl10 filament assembly and length dynamics, and demonstrate that filaments become shorter over time. Additionally, we develop an image-based, bootstrap-like resampling method that demonstrates the preferred association between autophagosomes and both Bcl10-filament ends and punctate-Bcl10 structures, implying that autophagosome-driven macroautophagy is directly responsible for Bcl10 filament shortening. We probe Bcl10 polymerization-depolymerization dynamics with a stochastic Monte-Carlo simulation of nucleation-limited filament assembly and degradation, and we show that high probabilities of filament nucleation in response to TCR engagement could provide the observed robust, homogeneous, and tunable response dynamic. Furthermore, we demonstrate that the speed of filament disassembly preferentially at filament ends provides effective regulatory control. Taken together, these data suggest that Bcl10 filament growth and degradation act as an excitable system that provides a digital response mechanism and the reliable timing critical for T cell activation and regulatory processes.

Dynamics phenotyping across length and time scales in collective cell migration.

Processes in collective migration span many length and time scales. In this review, we focus on length scales ranging from tens of microns (single cells) to a few millimeters (cell clusters) and the motion of these cells and cell groups on time scales of minutes to hours. We focus on epithelial cell sheets and metrics of motion developed to measure migration phenotypes in this system. Comparisons between cell motion and fluid flows, facilitated by the popular image analysis technique particle image velocimetry, yield metrics that can be used to study migration across a range of length and time scales. Measuring collective cell migration across these scales provides a complex, quantitative phenotype useful for migration models, in particular those that compare and contrast collective cell migration to movement of particles near a transition to jamming. Contrasting the motion of epithelial cells and the jamming transition illustrates aspects of collective motion that can be attributed to the jammed character of cell clusters, and highlights aspects of collective behavior that likely involve active motility and cell-cell guidance. The application of multiple migration metrics, which span multiple scales of the system, thus allows us to link cell-scale signals and mechanics to collective behavior.

The social shape of sperm: using an integrative machine-learning approach to examine sperm ultrastructure and collective motility.

Sperm is one of the most morphologically diverse cell types in nature, yet they also exhibit remarkable behavioural variation, including the formation of collective groups of cells that swim together for motility or transport through the female reproductive tract. Here, we take advantage of natural variation in sperm traits observed across Peromyscus mice to test the hypothesis that the morphology of the sperm head influences their sperm aggregation behaviour. Using both manual and automated morphometric approaches to quantify their complex shapes, and then statistical modelling and machine learning to analyse their features, we show that the aspect ratio of the sperm head is the most distinguishing morphological trait and statistically associates with collective sperm movements obtained from in vitro observations. We then successfully use neural network analysis to predict the size of sperm aggregates from sperm head morphology and show that species with relatively wider sperm heads form larger aggregates, which is consistent with the theoretical prediction that an adhesive region around the equatorial region of the sperm head mediates these unique gametic interactions. Together these findings advance our understanding of how even subtle variation in sperm design can drive differences in sperm function and performance.

Gravity Governs Shear Localization in Confined Dense Granular Flows.

The prediction of flow profiles of slowly sheared granular materials is a major geophysical and industrial challenge. Understanding the role of gravity is particularly important for future planetary exploration in varying gravitational environments. Using the principle of minimization of energy dissipation, and combining experiments and variational analysis, we disentangle the contributions of the gravitational acceleration, confining pressure, and layer thickness on shear strain localization induced by moving fault boundaries at the bottom of a granular layer. The flow profile is independent of the gravity for geometries with a free top surface. However, under a confining pressure or if the sheared layer withstands the weight of the upper layers, increasing gravity promotes the transition from closed shear zones buried in the bulk to open ones that intersect the top surface. We show that the center position and width of the shear zone and the axial angular velocity at the top surface follow universal scaling laws when properly scaled by the gravity, applied pressure, and layer thickness. Our finding that the flow profiles lie on a universal master curve opens the possibility to predict the quasistatic shear flow of granular materials in varying gravitational environments.

Adenylyl cyclase A mRNA localized at the back of cells is actively translated in live chemotaxing Dictyostelium.

Dictyostelium discoideum cells transport adenylyl cyclase A (ACA)-containing vesicles to the back of polarized cells to relay exogenous cAMP signals during chemotaxis. Fluorescence in situ hybridization (FISH) experiments showed that ACA mRNA is also asymmetrically distributed at the back of polarized cells. By using the MS2 bacteriophage system, we now visualize the distribution of ACA mRNA in live chemotaxing cells. We found that the ACA mRNA localization is not dependent on the translation of the protein product and requires multiple cis-acting elements within the ACA-coding sequence. We show that ACA mRNA is associated with actively translating ribosomes and is transported along microtubules towards the back of cells. By monitoring the recovery of ACA-YFP after photobleaching, we observed that local translation of ACA-YFP occurs at the back of cells. These data represent a novel functional role for localized translation in the relay of chemotactic signals during chemotaxis.

RefCell: multi-dimensional analysis of image-based high-throughput screens based on 'typical cells'.

BACKGROUND: Image-based high-throughput screening (HTS) reveals a high level of heterogeneity in single cells and multiple cellular states may be observed within a single population. Currently available high-dimensional analysis methods are successful in characterizing cellular heterogeneity, but suffer from the "curse of dimensionality" and non-standardized outputs. RESULTS: Here we introduce RefCell, a multi-dimensional analysis pipeline for image-based HTS that reproducibly captures cells with typical combinations of features in reference states and uses these "typical cells" as a reference for classification and weighting of metrics. RefCell quantitatively assesses heterogeneous deviations from typical behavior for each analyzed perturbation or sample. CONCLUSIONS: We apply RefCell to the analysis of data from a high-throughput imaging screen of a library of 320 ubiquitin-targeted siRNAs selected to gain insights into the mechanisms of premature aging (progeria). RefCell yields results comparable to a more complex clustering-based single-cell analysis method; both methods reveal more potential hits than a conventional analysis based on averages.

Leishmania infection inhibits macrophage motility by altering F-actin dynamics and the expression of adhesion complex proteins.

Leishmania is an intracellular protozoan parasite that causes a broad spectrum of clinical manifestations, ranging from self-healing skin lesions to fatal visceralizing disease. As the host cells of choice for all species of Leishmania, macrophages are critical for the establishment of infections. How macrophages contribute to parasite homing to specific tissues and how parasites modulate macrophage function are still poorly understood. In this study, we show that Leishmania amazonensis infection inhibits macrophage roaming motility. The reduction in macrophage speed is not dependent on particle load or on factors released by infected macrophages. L. amazonensis-infected macrophages also show reduced directional migration in response to the chemokine MCP-1. We found that infected macrophages have lower levels of total paxillin, phosphorylated paxillin, and phosphorylated focal adhesion kinase when compared to noninfected macrophages, indicating abnormalities in the formation of signaling adhesion complexes that regulate motility. Analysis of the dynamics of actin polymerization at peripheral sites also revealed a markedly enhanced F-actin turnover frequency in L. amazonensis-infected macrophages. Thus, Leishmania infection inhibits macrophage motility by altering actin dynamics and impairing the expression of proteins that function in plasma membrane-extracellular matrix interactions.

Collective cell migration has distinct directionality and speed dynamics.

When a constraint is removed, confluent cells migrate directionally into the available space. How the migration directionality and speed increase are initiated at the leading edge and propagate into neighboring cells are not well understood. Using a quantitative visualization technique-Particle Image Velocimetry (PIV)-we revealed that migration directionality and speed had strikingly different dynamics. Migration directionality increases as a wave propagating from the leading edge into the cell sheet, while the increase in cell migration speed is maintained only at the leading edge. The overall directionality steadily increases with time as cells migrate into the cell-free space, but migration speed remains largely the same. A particle-based compass (PBC) model suggests cellular interplay (which depends on cell-cell distance) and migration speed are sufficient to capture the dynamics of migration directionality revealed experimentally. Extracellular Ca2+ regulated both migration speed and directionality, but in a significantly different way, suggested by the correlation between directionality and speed only in some dynamic ranges. Our experimental and modeling results reveal distinct directionality and speed dynamics in collective migration, and these factors can be regulated by extracellular Ca2+ through cellular interplay. Quantitative visualization using PIV and our PBC model thus provide a powerful approach to dissect the mechanisms of collective cell migration.

Asymmetric nanotopography biases cytoskeletal dynamics and promotes unidirectional cell guidance.

Many biological and physiological processes depend upon directed migration of cells, which is typically mediated by chemical or physical gradients or by signal relay. Here we show that cells can be guided in a single preferred direction based solely on local asymmetries in nano/microtopography on subcellular scales. These asymmetries can be repeated, and thereby provide directional guidance, over arbitrarily large areas. The direction and strength of the guidance is sensitive to the details of the nano/microtopography, suggesting that this phenomenon plays a context-dependent role in vivo. We demonstrate that appropriate asymmetric nano/microtopography can unidirectionally bias internal actin polymerization waves and that cells move with the same preferred direction as these waves. This phenomenon is observed both for the pseudopod-dominated migration of the amoeboid Dictyostelium discoideum and for the lamellipod-driven migration of human neutrophils. The conservation of this mechanism across cell types and the asymmetric shape of many natural scaffolds suggest that actin-wave-based guidance is important in biology and physiology.

Adenylyl cyclase mRNA localizes to the posterior of polarized DICTYOSTELIUM cells during chemotaxis.

BACKGROUND: In Dictyostelium discoideum, vesicular transport of the adenylyl cyclase A (ACA) to the posterior of polarized cells is essential to relay exogenous 3',5'-cyclic adenosine monophosphate (cAMP) signals during chemotaxis and for the collective migration of cells in head-to-tail arrangements called streams. RESULTS: Using fluorescence in situ hybridization (FISH), we discovered that the ACA mRNA is asymmetrically distributed at the posterior of polarized cells. Using both standard estimators and Monte Carlo simulation methods, we found that the ACA mRNA enrichment depends on the position of the cell within a stream, with the posterior localization of ACA mRNA being strongest for cells at the end of a stream. By monitoring the recovery of ACA-YFP after cycloheximide (CHX) treatment, we observed that ACA mRNA and newly synthesized ACA-YFP first emerge as fluorescent punctae that later accumulate to the posterior of cells. We also found that the ACA mRNA localization requires 3' ACA cis-acting elements. CONCLUSIONS: Together, our findings suggest that the asymmetric distribution of ACA mRNA allows the local translation and accumulation of ACA protein at the posterior of cells. These data represent a novel functional role for localized translation in the relay of chemotactic signal during chemotaxis.

Analysis of microtubule growth dynamics arising from altered actin network structure and contractility in breast tumor cells.

The periphery of epithelial cells is shaped by opposing cytoskeletal physical forces generated predominately by two dynamic force generating systems-growing microtubule ends push against the boundary from the cell center, and the actin cortex contracts the attached plasma membrane. Here we investigate how changes to the structure and dynamics of the actin cortex alter the dynamics of microtubules. Current drugs target actin polymerization and contraction to reduce cell division and invasiveness; however, the impacts on microtubule dynamics remain incompletely understood. Using human MCF-7 breast tumor cells expressing GFP-tagged microtubule end-binding-protein-1 (EB1) and coexpression of cytoplasmic fluorescent protein mCherry, we map the trajectories of growing microtubule ends and cytoplasmic boundary respectively. Based on EB1 tracks and cytoplasmic boundary outlines, we calculate the speed, distance from cytoplasmic boundary, and straightness of microtubule growth. Actin depolymerization with Latrunculin-A reduces EB1 growth speed as well as allows the trajectories to extend beyond the cytoplasmic boundary. Blebbistatin, a direct myosin-II inhibitor, reduced EB1 speed and yielded less straight EB1 trajectories. Inhibiting signaling upstream of myosin-II contractility via the Rho-kinase inhibitor, Y-27632, altered EB1 dynamics differently from Blebbistatin. These results indicate that reduced actin cortex integrity can induce distinct alterations in microtubule dynamics. Given recent findings that tumor stem cell characteristics are increased by drugs which reduce actin contractility or stabilize microtubules, it remains important to clearly define how cytoskeletal drugs alter the interactions between these two filament systems in tumor cells.

svclassify: a method to establish benchmark structural variant calls.

BACKGROUND: The human genome contains variants ranging in size from small single nucleotide polymorphisms (SNPs) to large structural variants (SVs). High-quality benchmark small variant calls for the pilot National Institute of Standards and Technology (NIST) Reference Material (NA12878) have been developed by the Genome in a Bottle Consortium, but no similar high-quality benchmark SV calls exist for this genome. Since SV callers output highly discordant results, we developed methods to combine multiple forms of evidence from multiple sequencing technologies to classify candidate SVs into likely true or false positives. Our method (svclassify) calculates annotations from one or more aligned bam files from many high-throughput sequencing technologies, and then builds a one-class model using these annotations to classify candidate SVs as likely true or false positives. RESULTS: We first used pedigree analysis to develop a set of high-confidence breakpoint-resolved large deletions. We then used svclassify to cluster and classify these deletions as well as a set of high-confidence deletions from the 1000 Genomes Project and a set of breakpoint-resolved complex insertions from Spiral Genetics. We find that likely SVs cluster separately from likely non-SVs based on our annotations, and that the SVs cluster into different types of deletions. We then developed a supervised one-class classification method that uses a training set of random non-SV regions to determine whether candidate SVs have abnormal annotations different from most of the genome. To test this classification method, we use our pedigree-based breakpoint-resolved SVs, SVs validated by the 1000 Genomes Project, and assembly-based breakpoint-resolved insertions, along with semi-automated visualization using svviz. CONCLUSIONS: We find that candidate SVs with high scores from multiple technologies have high concordance with PCR validation and an orthogonal consensus method MetaSV (99.7 % concordant), and candidate SVs with low scores are questionable. We distribute a set of 2676 high-confidence deletions and 68 high-confidence insertions with high svclassify scores from these call sets for benchmarking SV callers. We expect these methods to be particularly useful for establishing high-confidence SV calls for benchmark samples that have been characterized by multiple technologies.

Consequences of anomalous diffusion in disordered systems under cyclic forcing.

We study the particle scale response of a 2D frictionless disk system to bulk forcing via cyclic shear with reversal amplitude γ_{r}. We find a subdiffusive γ_{r}-dependent regime, which is consistent with models of anomalous diffusion with scale-invariant cage dynamics, and a crossover to diffusive grain motion at high γ_{r}. Analysis of local displacements of a particle relative to its cage of neighbors reveals a key distinction from thermal systems. Particles are moved by fluctuations of their cage of neighbors rather than rattling in their cage, indicating a distinct cage-breaking mechanism.

Modeling and measuring signal relay in noisy directed migration of cell groups.

We develop a coarse-grained stochastic model for the influence of signal relay on the collective behavior of migrating Dictyostelium discoideum cells. In the experiment, cells display a range of collective migration patterns, including uncorrelated motion, formation of partially localized streams, and clumping, depending on the type of cell and the strength of the external, linear concentration gradient of the signaling molecule cyclic adenosine monophosphate (cAMP). From our model, we find that the pattern of migration can be quantitatively described by the competition of two processes, the secretion rate of cAMP by the cells and the degradation rate of cAMP in the gradient chamber. Model simulations are compared to experiments for a wide range of strengths of an external linear-gradient signal. With degradation, the model secreting cells form streams and efficiently transverse the gradient, but without degradation, we find that model secreting cells form clumps without streaming. This indicates that the observed effective collective migration in streams requires not only signal relay but also degradation of the signal. In addition, our model allows us to detect and quantify precursors of correlated motion, even when cells do not exhibit obvious streaming.

From cellular characteristics to disease diagnosis: uncovering phenotypes with supercells.

Cell heterogeneity and the inherent complexity due to the interplay of multiple molecular processes within the cell pose difficult challenges for current single-cell biology. We introduce an approach that identifies a disease phenotype from multiparameter single-cell measurements, which is based on the concept of "supercell statistics", a single-cell-based averaging procedure followed by a machine learning classification scheme. We are able to assess the optimal tradeoff between the number of single cells averaged and the number of measurements needed to capture phenotypic differences between healthy and diseased patients, as well as between different diseases that are difficult to diagnose otherwise. We apply our approach to two kinds of single-cell datasets, addressing the diagnosis of a premature aging disorder using images of cell nuclei, as well as the phenotypes of two non-infectious uveitides (the ocular manifestations of Behçet's disease and sarcoidosis) based on multicolor flow cytometry. In the former case, one nuclear shape measurement taken over a group of 30 cells is sufficient to classify samples as healthy or diseased, in agreement with usual laboratory practice. In the latter, our method is able to identify a minimal set of 5 markers that accurately predict Behçet's disease and sarcoidosis. This is the first time that a quantitative phenotypic distinction between these two diseases has been achieved. To obtain this clear phenotypic signature, about one hundred CD8(+) T cells need to be measured. Although the molecular markers identified have been reported to be important players in autoimmune disorders, this is the first report pointing out that CD8(+) T cells can be used to distinguish two systemic inflammatory diseases. Beyond these specific cases, the approach proposed here is applicable to datasets generated by other kinds of state-of-the-art and forthcoming single-cell technologies, such as multidimensional mass cytometry, single-cell gene expression, and single-cell full genome sequencing techniques.