An image-based data-driven analysis of cellular architecture in a developing tissue

Systematic data analysis pipeline for quantitative morphological cell phenotyping

Quantitative morphological phenotyping (QMP) is an image-based method used to capture morphological features at both the cellular and population level. Its interdisciplinary nature, spanning from data collection to result analysis and interpretation, can lead to uncertainties, particularly among those new to this actively growing field. High analytical specificity for a typical QMP is achieved through sophisticated approaches that can leverage subtle cellular morphological changes. Here, we outline a systematic workflow to refine the QMP methodology. For a practical review, we describe the main steps of a typical QMP; in each step, we discuss the available methods, their applications, advantages, and disadvantages, along with the R functions and packages for easy implementation. This review does not cover theoretical backgrounds, but provides several references for interested researchers. It aims to broaden the horizons for future phenome studies and demonstrate how to exploit years of endeavors to achieve more with less.

Spherical harmonics analysis reveals cell shape-fate relationships in zebrafish lateral line neuromasts

Cell shape is a powerful readout of cell state, fate and function. We describe a custom workflow to perform semi-automated, 3D cell and nucleus segmentation, and spherical harmonics and principal components analysis to distill cell and nuclear shape variation into discrete biologically meaningful parameters. We apply these methods to analyze shape in the neuromast cells of the zebrafish lateral line system, finding that shapes vary with cell location and identity. The distinction between hair cells and support cells accounted for much of the variation, which allowed us to train classifiers to predict cell identity from shape features. Using transgenic markers for support cell subpopulations, we found that subtypes had different shapes from each other. To investigate how loss of a neuromast cell type altered cell shape distributions, we examined atoh1a mutants that lack hair cells. We found that mutant neuromasts lacked the cell shape phenotype associated with hair cells, but did not exhibit a mutant-specific cell shape. Our results demonstrate the utility of using 3D cell shape features to characterize, compare and classify cells in a living developing organism.

Comprehensive and quantitative analysis of intracellular structure polarization at the apical-basal axis in elongating Arabidopsis zygotes

A comprehensive and quantitative evaluation of multiple intracellular structures or proteins is a promising approach to provide a deeper understanding of and new insights into cellular polarity. In this study, we developed an image analysis pipeline to obtain intensity profiles of fluorescent probes along the apical-basal axis in elongating Arabidopsis thaliana zygotes based on two-photon live-cell imaging data. This technique showed the intracellular distribution of actin filaments, mitochondria, microtubules, and vacuolar membranes along the apical-basal axis in elongating zygotes from the onset of cell elongation to just before asymmetric cell division. Hierarchical cluster analysis of the quantitative data on intracellular distribution revealed that the zygote may be compartmentalized into two parts, with a boundary located 43.6% from the cell tip, immediately after fertilization. To explore the biological significance of this compartmentalization, we examined the positions of the asymmetric cell divisions from the dataset used in this distribution analysis. We found that the cell division plane was reproducibly inserted 20.5% from the cell tip. This position corresponded well with the midpoint of the compartmentalized apical region, suggesting a potential relationship between the zygote compartmentalization, which begins with cell elongation, and the position of the asymmetric cell division.

Post-injury hydraulic fracturing drives fissure formation in the zebrafish basal epidermal cell layer

The skin epithelium acts as the barrier between an organism's internal and external environments. In zebrafish and other freshwater organisms, this barrier function requires withstanding a large osmotic gradient across the epidermis. Wounds breach this epithelium, causing a large disruption to the tissue microenvironment due to the mixing of isotonic interstitial fluid with the external hypotonic fresh water. Here, we show that, following acute injury, the larval zebrafish epidermis undergoes a dramatic fissuring process that resembles hydraulic fracturing, driven by the influx of external fluid. After the wound has sealed-preventing efflux of this external fluid-fissuring starts in the basal epidermal layer at the location nearest to the wound and then propagates at a constant rate through the tissue, spanning over 100 μm. During this process, the outermost superficial epidermal layer remains intact. Fissuring is completely inhibited when larvae are wounded in isotonic external media, suggesting that osmotic gradients are required for fissure formation. Additionally, fissuring partially depends on myosin II activity, as myosin II inhibition reduces the distance of fissure propagation away from the wound. During and after fissuring, the basal layer forms large macropinosomes (with cross-sectional areas ranging from 1 to 10 μm²). We conclude that excess external fluid entry through the wound and subsequent closure of the wound through actomyosin purse-string contraction in the superficial cell layer causes fluid pressure buildup in the extracellular space of the zebrafish epidermis. This excess fluid pressure causes tissue to fissure, and eventually the fluid is cleared through macropinocytosis.

Rosemary essential oil and its components 1,8-cineole and α-pinene induce ROS-dependent lethality and ROS-independent virulence inhibition in Candida albicans

The essential oil from Rosmarinus officinalis L., a composite mixture of plant-derived secondary metabolites, exhibits antifungal activity against virulent candidal species. Here we report the impact of rosemary oil and two of its components, the monoterpene α-pinene and the monoterpenoid 1,8-cineole, against Candida albicans, which induce ROS-dependent cell death at high concentrations and inhibit hyphal morphogenesis and biofilm formation at lower concentrations. The minimum inhibitory concentrations (100% inhibition) for both rosemary oil and 1,8-cineole were 4500 μg/ml and 3125 μg/ml for α-pinene, with the two components exhibiting partial synergy (FICI = 0.55 ± 0.07). At MIC and 1/2 MIC, rosemary oil and its components induced a generalized cell wall stress response, causing damage to cellular and organelle membranes, along with elevated chitin production and increased cell surface adhesion and elasticity, leading to complete vacuolar segregation, mitochondrial depolarization, elevated reactive oxygen species, microtubule dysfunction, and cell cycle arrest mainly at the G1/S phase, consequently triggering cell death. Interestingly, the same oils at lower fractional MIC (1/8-1/4) inhibited virulence traits, including reduction of mycelium (up to 2-fold) and biofilm (up to 4-fold) formation, through a ROS-independent mechanism.

Computational modeling and analysis of the morphogenetic domain signaling networks regulating C. elegans embryogenesis

Caenorhabditis elegans, often referred to as the ‘roundworm’, provides a powerful model for studying cell autonomous and cell–cell interactions through the direct observation of embryonic development in vivo. By leveraging the precisely mapped cell lineage at single cell resolution, we are able to study at a systems level how early embryonic cells communicate across morphogenetic domains for the coordinated processes of gene expressions and collective cellular behaviors that regulate tissue morphogenesis. In this study, we developed a computational framework for the exploration of the morphogenetic domain cell signaling networks that may regulate C. elegans gastrulation and embryonic organogenesis. We demonstrated its utility by producing the following results, i) established a virtual reference model of developing C. elegans embryos through the spatiotemporal alignment of individual embryo cell nuclear imaging samples; ii) integrated the single cell spatiotemporal gene expression profile with the established virtual embryo model by data pooling; iii) trained a Machine Learning model (Random Forest Regression), which predicts accurately the spatial positions of the cells given their gene expression profiles for a given developmental time (e.g. total cell number of the embryo); iv) enabled virtual 4-dimensional tomographic graphical modeling of single cell data; v) inferred the biology signaling pathways that act in each of morphogenetic domains by meta-data analysis. It is intriguing that the morphogenetic domain cell signaling network seems to involve some crosstalk of multiple biology signaling pathways during the formation of tissue boundary pattern. Lastly, we developed the Software tool ‘Embryo aligner version 1.0’ and provided it as an Open Source program to the research community for virtual embryo modeling, and phenotype perturbation analyses (https://github.com/csniuben/embryo_aligner/wiki and https://bioinfo89.github.io/C.elegansEmbryonicOrganogenesisweb/).

Dynamics of Drosophila endoderm specification

During early Drosophila embryogenesis, a network of gene regulatory interactions orchestrates terminal patterning, playing a critical role in the subsequent formation of the gut. We utilized CRISPR gene editing at endogenous loci to create live reporters of transcription and light-sheet microscopy to monitor the individual components of the posterior gut patterning network across 90 min prior to gastrulation. We developed a computational approach for fusing imaging datasets of the individual components into a common multivariable trajectory. Data fusion revealed low intrinsic dimensionality of posterior patterning and cell fate specification in wild-type embryos. The simple structure that we uncovered allowed us to construct a model of interactions within the posterior patterning regulatory network and make testable predictions about its dynamics at the protein level. The presented data fusion strategy is a step toward establishing a unified framework that would explore how stochastic spatiotemporal signals give rise to highly reproducible morphogenetic outcomes.

A role for the centrosome in regulating the rate of neuronal efferocytosis by microglia in vivo

During brain development, many newborn neurons undergo apoptosis and are engulfed by microglia, the tissue-resident phagocytes of the brain, in a process known as efferocytosis. A hallmark of microglia is their highly branched morphology characterized by the presence of numerous dynamic extensions that these cells use for scanning the brain parenchyma and engulfing unwanted material. The mechanisms driving branch formation and apoptotic cell engulfment in microglia are unclear. By taking a live-imaging approach in zebrafish, we show that while microglia generate multiple microtubule-based branches, they only successfully engulf one apoptotic neuron at a time. Further investigation into the mechanism underlying this sequential engulfment revealed that targeted migration of the centrosome into one branch is predictive of phagosome formation and polarized vesicular trafficking. Moreover, experimentally doubling centrosomal numbers in microglia increases the rate of engulfment and even allows microglia to remove two neurons simultaneously, providing direct supporting evidence for a model where centrosomal migration is a rate-limiting step in branch-mediated efferocytosis. Conversely, light-mediated depolymerization of microtubules causes microglia to lose their typical branched morphology and switch to an alternative mode of engulfment, characterized by directed migration towards target neurons, revealing unexpected plasticity in their phagocytic ability. Finally, building on work focusing on the establishment of the immunological synapse, we identified a conserved signalling pathway underlying centrosomal movement in engulfing microglia.

Origami: Single-cell 3D shape dynamics oriented along the apico-basal axis of folding epithelia from fluorescence microscopy data

A common feature of morphogenesis is the formation of three-dimensional structures from the folding of two-dimensional epithelial sheets, aided by cell shape changes at the cellular-level. Changes in cell shape must be studied in the context of cell-polarised biomechanical processes within the epithelial sheet. In epithelia with highly curved surfaces, finding single-cell alignment along a biological axis can be difficult to automate in silico. We present 'Origami', a MATLAB-based image analysis pipeline to compute direction-variant cell shape features along the epithelial apico-basal axis. Our automated method accurately computed direction vectors denoting the apico-basal axis in regions with opposing curvature in synthetic epithelia and fluorescence images of zebrafish embryos. As proof of concept, we identified different cell shape signatures in the developing zebrafish inner ear, where the epithelium deforms in opposite orientations to form different structures. Origami is designed to be user-friendly and is generally applicable to fluorescence images of curved epithelia.

Organ geometry channels reproductive cell fate in the Arabidopsis ovule primordium

In multicellular organisms, sexual reproduction requires the separation of the germline from the soma. In flowering plants, the female germline precursor differentiates as a single spore mother cell (SMC) as the ovule primordium forms. Here, we explored how organ growth contributes to SMC differentiation. We generated 92 annotated 3D images at cellular resolution in Arabidopsis. We identified the spatio-temporal pattern of cell division that acts in a domain-specific manner as the primordium forms. Tissue growth models uncovered plausible morphogenetic principles involving a spatially confined growth signal, differential mechanical properties, and cell growth anisotropy. Our analysis revealed that SMC characteristics first arise in more than one cell but SMC fate becomes progressively restricted to a single cell during organ growth. Altered primordium geometry coincided with a delay in the fate restriction process in katanin mutants. Altogether, our study suggests that tissue geometry channels reproductive cell fate in the Arabidopsis ovule primordium.

Data science in cell imaging

Cell imaging has entered the 'Big Data' era. New technologies in light microscopy and molecular biology have led to an explosion in high-content, dynamic and multidimensional imaging data. Similar to the 'omics' fields two decades ago, our current ability to process, visualize, integrate and mine this new generation of cell imaging data is becoming a critical bottleneck in advancing cell biology. Computation, traditionally used to quantitatively test specific hypotheses, must now also enable iterative hypothesis generation and testing by deciphering hidden biologically meaningful patterns in complex, dynamic or high-dimensional cell image data. Data science is uniquely positioned to aid in this process. In this Perspective, we survey the rapidly expanding new field of data science in cell imaging. Specifically, we highlight how data science tools are used within current image analysis pipelines, propose a computation-first approach to derive new hypotheses from cell image data, identify challenges and describe the next frontiers where we believe data science will make an impact. We also outline steps to ensure broad access to these powerful tools - democratizing infrastructure availability, developing sensitive, robust and usable tools, and promoting interdisciplinary training to both familiarize biologists with data science and expose data scientists to cell imaging.

Spatial frequency metrics for analysis of microscopic images of musculoskeletal tissues

Purpose: Imaging-based metrics for analysis of biological tissues are powerful tools that can extract information such as shape, size, periodicity, and many other features to assess the requested qualities of a tissue. Muscular and osseous tissues consist of periodic structures that are directly related to their function, and so analysis of these patterns likely reflects tissue health and regeneration.Methods: A method for assessment of periodic structures is by analyzing them in the spatial frequency domain using the Fourier transform. In this paper, we present two filters which we developed in the spatial frequency domain for the purpose of analyzing musculoskeletal structures. These filters provide information about 1) the angular orientation of the tissues and 2) their periodicity. We explore periodic structural patterns in the mitochondrial network of skeletal muscles that are reflective of muscle metabolism and myogenesis; and patterns of collagen fibers in the bone that are reflective of the organization and health of bone extracellular matrix.Results: We present an analysis of mouse skeletal muscle in healthy and injured muscles. We used a transgenic mouse that ubiquitously expresses fluorescent protein in their mitochondria and performed 2-photon microscopy to image the structures. To acquire the collagen structure of the bone we used non-linear SHG microscopy of mouse flat bone. We analyze and compare juvenile versus adult mice, which have different structural patterns.Conclusions: Our results indicate that these metrics can quantify musculoskeletal tissues during development and regeneration.

Zebrafish Posterior Lateral Line primordium migration requires interactions between a superficial sheath of motile cells and the skin

The Zebrafish Posterior Lateral Line primordium migrates in a channel between the skin and somites. Its migration depends on the coordinated movement of its mesenchymal-like leading cells and trailing cells, which form epithelial rosettes, or protoneuromasts. We describe a superficial population of flat primordium cells that wrap around deeper epithelialized cells and extend polarized lamellipodia to migrate apposed to the overlying skin. Polarization of lamellipodia extended by both superficial and deeper protoneuromast-forming cells depends on Fgf signaling. Removal of the overlying skin has similar effects on superficial and deep cells: lamellipodia are lost, blebs appear instead, and collective migration fails. When skinned embryos are embedded in Matrigel, basal and superficial lamellipodia are recovered; however, only the directionality of basal protrusions is recovered, and migration is not rescued. These observations support a key role played by superficial primordium cells and the skin in directed migration of the Posterior Lateral Line primordium.