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. 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.  more » « less
Award ID(s):
1714429 1740741 1764421
NSF-PAR ID:
10160190
Author(s) / Creator(s):
; ;
Date Published:
Journal Name:
Proceedings of the National Academy of Sciences
Volume:
117
Issue:
10
ISSN:
0027-8424
Page Range / eLocation ID:
5113 to 5124
Format(s):
Medium: X
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 cell 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. 
    more » « 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 of 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. 
    more » « 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.

     
    more » « less
  4. Abstract Background

    Sexually dimorphic mating behaviors differ between sexes and involve gonadal hormones and possibly sexually dimorphic gene expression in the brain. However, the associations among the brain, gonad, and sexual behavior in teleosts are still unclear. Here, we utilized germ cells-freetdrd12knockout (KO) zebrafish, and steroid synthesis enzymecyp17a1-deficient zebrafish to investigate the differences and interplays in the brain–gonad–behavior axis, and the molecular control of brain dimorphism and male mating behaviors.

    Methods

    Tdrd12+/−;cyp17a1+/−double heterozygous parents were crossed to obtaintdrd12−/−;cyp17a1+/+(tdrd12 KO),tdrd12+/+;cyp17a1−/−(cyp17a1 KO), andtdrd12−/−;cyp17a1−/−(double KO) homozygous progenies. Comparative analysis of mating behaviors were evaluated using Viewpoint zebrafish tracking software and sexual traits were thoroughly characterized based on anatomical and histological experiments in these KOs and wild types. The steroid hormone levels (testosterone, 11-ketotestosterone and 17β-estradiol) in the brains, gonads, and serum were measured using ELISA kits. To achieve a higher resolution view of the differences in region-specific expression patterns of the brain, the brains of these KOs, and control male and female fish were dissected into three regions: the forebrain, midbrain, and hindbrain for transcriptomic analysis.

    Results

    Qualitative analysis of mating behaviors demonstrated thattdrd12−/−fish behaved in the same manner as wild-type males to trigger oviposition behavior, whilecyp17a1−/−and double knockout (KO) fish did not exhibit these behaviors. Based on the observation of sex characteristics, mating behaviors and hormone levels in these mutants, we found that the maintenance of secondary sex characteristics and male mating behavior did not depend on the presence of germ cells; rather, they depended mainly on the 11-ketotestosterone and testosterone levels secreted into the brain–gonad regulatory axis. RNA-seq analysis of different brain regions revealed that the brain transcript profile oftdrd12−/−fish was similar to that of wild-type males, especially in the forebrain and midbrain. However, the brain transcript profiles ofcyp17a1−/−and double KO fish were distinct from those of wild-type males and were partially biased towards the expression pattern of the female brain. Our results revealed important candidate genes and signaling pathways, such as synaptic signaling/neurotransmission, MAPK signaling, and steroid hormone pathways, that shape brain dimorphism and modulate male mating behavior in zebrafish.

    Conclusions

    Our results provide comprehensive analyses and new insights regarding the endogenous interactions in the brain–gonad–behavior axis. Moreover, this study revealed the crucial candidate genes and neural signaling pathways of different brain regions that are involved in modulating brain dimorphism and male mating behavior in zebrafish, which would significantly light up the understanding the neuroendocrine and molecular mechanisms modulating brain dimorphism and male mating behavior in zebrafish and other teleost fish.

    Graphical Abstract 
    more » « less
  5. 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. 
    more » « less