skip to main content

Title: Topological data analysis of zebrafish patterns
Self-organized pattern behavior is ubiquitous throughout nature, from fish schooling to collective cell dynamics during organism development. Qualitatively these patterns display impressive consistency, yet variability inevitably exists within pattern-forming systems on both microscopic and macroscopic scales. Quantifying variability and measuring pattern features can inform the underlying agent interactions and allow for predictive analyses. Nevertheless, current methods for analyzing patterns that arise from collective behavior capture only macroscopic features or rely on either manual inspection or smoothing algorithms that lose the underlying agent-based nature of the data. Here we introduce methods based on topological data analysis and interpretable machine learning for quantifying both agent-level features and global pattern attributes on a large scale. Because the zebrafish is a model organism for skin pattern formation, we focus specifically on analyzing its skin patterns as a means of illustrating our approach. Using a recent agent-based model, we simulate thousands of wild-type and mutant zebrafish patterns and apply our methodology to better understand pattern variability in zebrafish. Our methodology is able to quantify the differential impact of stochasticity in cell interactions on wild-type and mutant patterns, and we use our methods to predict stripe and spot statistics as a function of varying cellular communication. more » Our work provides an approach to automatically quantifying biological patterns and analyzing agent-based dynamics so that we can now answer critical questions in pattern formation at a much larger scale. « less
; ;
Award ID(s):
1714429 1740741 1764421
Publication Date:
Journal Name:
Proceedings of the National Academy of Sciences
Page Range or eLocation-ID:
5113 to 5124
Sponsoring Org:
National Science Foundation
More Like this
  1. ABSTRACT Single mutations frequently alter several aspects of cell behavior but rarely reveal whether a particular statistically significant change is biologically significant. To determine which behavioral changes are most important for multicellular self-organization, we devised a new methodology using Myxococcus xanthus as a model system. During development, myxobacteria coordinate their movement to aggregate into spore-filled fruiting bodies. We investigate how aggregation is restored in two mutants, csgA and pilC , that cannot aggregate unless mixed with wild-type (WT) cells. To this end, we use cell tracking to follow the movement of fluorescently labeled cells in combination with data-driven agent-based modeling. The results indicate that just like WT cells, both mutants bias their movement toward aggregates and reduce motility inside aggregates. However, several aspects of mutant behavior remain uncorrected by WT, demonstrating that perfect recreation of WT behavior is unnecessary. In fact, synergies between errant behaviors can make aggregation robust. IMPORTANCE Self-organization into spatial patterns is evident in many multicellular phenomena. Even for the best-studied systems, our ability to dissect the mechanisms driving coordinated cell movement is limited. While genetic approaches can identify mutations perturbing multicellular patterns, the diverse nature of the signaling cues coupled to significant heterogeneity of individual cellmore »behavior impedes our ability to mechanistically connect genes with phenotype. Small differences in the behaviors of mutant strains could be irrelevant or could sometimes lead to large differences in the emergent patterns. Here, we investigate rescue of multicellular aggregation in two mutant strains of Myxococcus xanthus mixed with wild-type cells. The results demonstrate how careful quantification of cell behavior coupled to data-driven modeling can identify specific motility features responsible for cell aggregation and thereby reveal important synergies and compensatory mechanisms. Notably, mutant cells do not need to precisely recreate wild-type behaviors to achieve complete aggregation.« less
  2. Patil, Kiran Raosaheb (Ed.)
    Chemotactic bacteria form emergent spatial patterns of variable cell density within cultures that are initially spatially uniform. These patterns are the result of chemical gradients that are created from the directed movement and metabolic activity of billions of cells. A recent study on pattern formation in wild bacterial isolates has revealed unique collective behaviors of the bacteria Enterobacter cloacae . As in other bacterial species, Enterobacter cloacae form macroscopic aggregates. Once formed, these bacterial clusters can migrate several millimeters, sometimes resulting in the merging of two or more clusters. To better understand these phenomena, we examine the formation and dynamics of thousands of bacterial clusters that form within a 22 cm square culture dish filled with soft agar over two days. At the macroscale, the aggregates display spatial order at short length scales, and the migration of cell clusters is superdiffusive, with a merging acceleration that is correlated with aggregate size. At the microscale, aggregates are composed of immotile cells surrounded by low density regions of motile cells. The collective movement of the aggregates is the result of an asymmetric flux of bacteria at the boundary. An agent-based model is developed to examine how these phenomena are the result ofmore »both chemotactic movement and a change in motility at high cell density. These results identify and characterize a new mechanism for collective bacterial motility driven by a transient, density-dependent change in motility.« less
  3. Abstract

    The coloring of zebrafish skin is often used as a model system to study biological pattern formation. However, the small number and lack of movement of chromatophores defies traditional Turing-type pattern generating mechanisms. Recent models invoke discrete short-range competition and long-range promotion between different pigment cells as an alternative to a reaction-diffusion scheme. In this work, we propose a lattice-based “Survival model,” which is inspired by recent experimental findings on the nature of long-range chromatophore interactions. The Survival model produces stationary patterns with diffuse stripes and undergoes a Turing instability. We also examine the effect that domain growth, ubiquitous in biological systems, has on the patterns in both the Survival model and an earlier “Promotion” model. In both cases, domain growth alone is capable of orienting Turing patterns above a threshold wavelength and can reorient the stripes in ablated cells, though the wavelength for which the patterns orient is much larger for the Survival model. While the Survival model is a simplified representation of the multifaceted interactions between pigment cells, it reveals complex organizational behavior and may help to guide future studies.

  4. Drying of bacterial suspensions is frequently encountered in a plethora of natural and engineering processes. However, the evaporation-driven mechanical instabilities of dense consolidating bacterial suspensions have not been explored heretofore. Here, we report the formation of two different crack patterns of drying suspensions of Escherichia coli ( E. coli ) with distinct motile behaviors. Circular cracks are observed for wild-type E. coli with active swimming, whereas spiral-like cracks form for immotile bacteria. Using the elastic fracture mechanics and the poroelastic theory, we show that the formation of the circular cracks is determined by the tensile nature of the radial drying stress once the cracks are initiated by the local order structure of bacteria due to their collective swimming. Our study demonstrates the link between the microscopic swimming behaviors of individual bacteria and the mechanical instabilities and macroscopic pattern formation of drying bacterial films. The results shed light on the dynamics of active matter in a drying process and provide useful information for understanding various biological processes associated with drying bacterial suspensions.
  5. Abstract Background

    Repetitive action, resistance to environmental change and fine motor disruptions are hallmarks of autism spectrum disorder (ASD) and other neurodevelopmental disorders, and vary considerably from individual to individual. In animal models, conventional behavioral phenotyping captures such fine-scale variations incompletely. Here we observed male and female C57BL/6J mice to methodically catalog adaptive movement over multiple days and examined two rodent models of developmental disorders against this dynamic baseline. We then investigated the behavioral consequences of a cerebellum-specific deletion in Tsc1 protein and a whole-brain knockout in Cntnap2 protein in mice. Both of these mutations are found in clinical conditions and have been associated with ASD.


    We used advances in computer vision and deep learning, namely a generalized form of high-dimensional statistical analysis, to develop a framework for characterizing mouse movement on multiple timescales using a single popular behavioral assay, the open-field test. The pipeline takes virtual markers from pose estimation to find behavior clusters and generate wavelet signatures of behavior classes. We measured spatial and temporal habituation to a new environment across minutes and days, different types of self-grooming, locomotion and gait.


    Both Cntnap2 knockouts and L7-Tsc1 mutants showed forelimb lag during gait. L7-Tsc1 mutants and Cntnap2 knockouts showed complexmore »defects in multi-day adaptation, lacking the tendency of wild-type mice to spend progressively more time in corners of the arena. In L7-Tsc1 mutant mice, failure to adapt took the form of maintained ambling, turning and locomotion, and an overall decrease in grooming. However, adaptation in these traits was similar between wild-type mice and Cntnap2 knockouts. L7-Tsc1 mutant and Cntnap2 knockout mouse models showed different patterns of behavioral state occupancy.


    Genetic risk factors for autism are numerous, and we tested only two. Our pipeline was only done under conditions of free behavior. Testing under task or social conditions would reveal more information about behavioral dynamics and variability.


    Our automated pipeline for deep phenotyping successfully captures model-specific deviations in adaptation and movement as well as differences in the detailed structure of behavioral dynamics. The reported deficits indicate that deep phenotyping constitutes a robust set of ASD symptoms that may be considered for implementation in clinical settings as quantitative diagnosis criteria.

    « less