Cell motility: Bioelectrical control of behavior without neurons

Single cells are capable of remarkably sophisticated, sometimes animal-like, behaviors. New work demonstrates bioelectric control of motility through the differential regulation of appendage movements in a unicellular organism that walks across surfaces using leg-like bundles of cilia.

The nuclear transport factor CSE1 drives macronuclear volume increase and macronuclear node coalescence in Stentor coeruleus

Stentor coeruleus provides a unique opportunity to study how cells regulate nuclear shape because its macronucleus undergoes a rapid, dramatic, and developmentally regulated shape change. We found that the volume of the macronucleus increases during coalescence, suggesting an inflation-based mechanism. When the nuclear transport factor, CSE1, is knocked down by RNAi, the shape and volume changes of the macronucleus are attenuated, and nuclear morphology is altered. CSE1 protein undergoes a dynamic relocalization correlated with nuclear shape changes, being mainly cytoplasmic prior to nuclear coalescence, and accumulating inside the macronucleus during coalescence. At the end of regeneration, CSE1 protein levels are reduced as the macronucleus returns to its pre-coalescence volume. We propose a model in which nuclear transport via CSE1 is required to increase the volume of the macronucleus, thereby decreasing the surface-to-volume ratio and driving coalescence of the nodes into a single mass.

Mitochondrial networks through the lens of mathematics

Mitochondria serve a wide range of functions within cells, most notably via their production of ATP. Although their morphology is commonly described as bean-like, mitochondria often form interconnected networks within cells that exhibit dynamic restructuring through a variety of physical changes. Further, though relationships between form and function in biology are well established, the extant toolkit for understanding mitochondrial morphology is limited. Here, we emphasize new and established methods for quantitatively describing mitochondrial networks, ranging from unweighted graph-theoretic representations to multi-scale approaches from applied topology, in particular persistent homology. We also show fundamental relationships between mitochondrial networks, mathematics, and physics, using ideas of graph planarity and statistical mechanics to better understand the full possible morphological space of mitochondrial network structures. Lastly, we provide suggestions for how examination of mitochondrial network form through the language of mathematics can inform biological understanding, and vice versa.

Homologous chromosome recognition via nonspecific interactions

In many organisms, most notably Drosophila, homologous chromosomes in somatic cells associate with each other, a phenomenon known as somatic homolog pairing. Unlike in meiosis, where homology is read out at the level of DNA sequence complementarity, somatic homolog pairing takes place without double strand breaks or strand invasion, thus requiring some other mechanism for homologs to recognize each other. Several studies have suggested a "specific button" model, in which a series of distinct regions in the genome, known as buttons, can associate with each other, presumably mediated by different proteins that bind to these different regions. Here we consider an alternative model, which we term the "button barcode" model, in which there is only one type of recognition site or adhesion button, present in many copies in the genome, each of which can associate with any of the others with equal affinity. An important component of this model is that the buttons are non-uniformly distributed, such that alignment of a chromosome with its correct homolog, compared with a non-homolog, is energetically favored; since to achieve nonhomologous alignment, chromosomes would be required to mechanically deform in order to bring their buttons into mutual register. We investigated several types of barcodes and examined their effect on pairing fidelity. We found that high fidelity homolog recognition can be achieved by arranging chromosome pairing buttons according to an actual industrial barcode used for warehouse sorting. By simulating randomly generated non-uniform button distributions, many highly effective button barcodes can be easily found, some of which achieve virtually perfect pairing fidelity. This model is consistent with existing literature on the effect of translocations of different sizes on homolog pairing. We conclude that a button barcode model can attain highly specific homolog recognition, comparable to that seen in actual cells undergoing somatic homolog pairing, without the need for specific interactions. This model may have implications for how meiotic pairing is achieved.

The flagellar length control system: exploring the physical biology of organelle size

How cells build and maintain dynamic structures of defined size is currently an important unsolved problem in quantitative cell biology. The flagella of the unicellular green algaChlamydomonasprovide a highly tractable model system to investigate this general question, but while the powerful genetics of this organism have revealed numerous genes required for proper flagellar length, in most cases we do not understand their mechanistic role in length control. Flagellar length can be viewed as the steady state solution of a dynamical system involving assembly and disassembly of axonemal microtubules, with assembly depending on an active transport process known as intraflagellar transport (IFT). The inherent length dependence of IFT gives rise to a family of simple models for length regulation that can account for many previously described phenomena such as the ability of flagella to maintain equal lengths. But these models requires that the cell has a way to measure flagellar length in order to adjust IFT rates accordingly. Several models for length sensing have been modeled theoretically and evaluated experimentally, allowing them to be ruled out. Current data support a model in which the diffusive return of the kinesin motor driving IFT provides a length dependence that ultimately is the basis for length regulation. By combining models of length sensing with a more detailed representation of cargo transport and availability, it is now becoming possible to formulate concrete hypotheses to explain length altering mutants.

Single-cell analysis of habituation in Stentor coeruleus

Although learning is often viewed as a unique feature of organisms with complex nervous systems, single-celled organisms also demonstrate basic forms of learning. The giant ciliate Stentor coeruleus responds to mechanical stimuli by contracting into a compact shape, presumably as a defense mechanism. When a Stentor cell is repeatedly stimulated at a constant level of force, it will learn to ignore that stimulus but will still respond to stronger stimuli. Prior studies of habituation in Stentor reported a graded response, suggesting that cells transition through a continuous range of response probabilities. By analyzing single cells using an automated apparatus to deliver calibrated stimuli, we find that habituation occurs via a single step-like switch in contraction probability within each cell, with the graded response in a population arising from the random distribution of switching times in individual cells. This step-like response allows Stentor behavior to be represented by a simple two-state model whose parameters can be estimated from experimental measurements. We find that transition rates depend on stimulus force and also on the time between stimuli. The ability to measure the behavior of the same cell to the same stimulus allowed us to quantify the functional heterogeneity among single cells. Together, our results suggest that the behavior of Stentor is governed by a two-state stochastic machine whose transition rates are sensitive to the time series properties of the input stimuli.

Testing the ion-current model for flagellar length sensing and IFT regulation

Eukaryotic cilia and flagella are microtubule-based organelles whose relatively simple shape makes them ideal for investigating the fundamental question of organelle size regulation. Most of the flagellar materials are transported from the cell body via an active transport process called intraflagellar transport (IFT). The rate of IFT entry into flagella, known as IFT injection, has been shown to negatively correlate with flagellar length. However, it remains unknown how the cell measures the length of its flagella and controls IFT injection. One of the most-discussed theoretical models for length sensing to control IFT is the ion-current model, which posits that there is a uniform distribution of Ca²⁺ channels along the flagellum and that the Ca²⁺ current from the flagellum into the cell body increases linearly with flagellar length. In this model, the cell uses the Ca²⁺ current to negatively regulate IFT injection. The recent discovery that IFT entry into flagella is regulated by the phosphorylation of kinesin through a calcium-dependent protein kinase has provided further impetus for the ion-current model. To test this model, we measured and manipulated the levels of Ca²⁺ inside of Chlamydomonas flagella and quantified IFT injection. Although the concentration of Ca²⁺ inside of flagella was weakly correlated with the length of flagella, we found that IFT injection was reduced in calcium-deficient flagella, rather than increased as the model predicted, and that variation in IFT injection was uncorrelated with the occurrence of flagellar Ca²⁺ spikes. Thus, Ca²⁺ does not appear to function as a negative regulator of IFT injection, hence it cannot form the basis of a stable length control system.

Role of intraflagellar transport in transcriptional control during flagellar regeneration in Chlamydomonas

Biosynthesis of organelle precursors is a central part of the organelle size control problem, but what systems are required to control precursor production? Genes encoding flagellar proteins are upregulated during flagellar regeneration in Chlamydomonas, and this upregulation is critical for flagella to reach their final length, but it not known how the cell triggers these genes during regeneration. We present two models based on transcriptional repressor that is either produced in the flagellum, or else is produced in the cell body and sequestered in the growing flagellum. Both models lead to stable flagellar length control and can reproduce the observed dynamics of gene expression. The two models make opposite predictions regarding the effect of mutations that block intraflagellar transport (IFT). Using quantitative measurements of gene expression, we show that gene expression during flagellar regeneration is greatly reduced in mutations of the heterotrimeric kinesin-2 that drives IFT. This result is consistent with the predictions of the model in which a repressor is sequestered in the flagellum by IFT. Inhibiting axonemal assembly has much less effect on gene expression. The repressor sequestration model allows precursor production to occur when flagella are growing rapidly, representing a form of derivative control.

A unicellular walker controlled by a microtubule-based finite-state machine

Cells are complex biochemical systems whose behaviors emerge from interactions among myriad molecular components. Computation is often invoked as a general framework for navigating this cellular complexity. However, it is unclear how cells might embody computational processes such that the theories of computation, including finite-state machine models, could be productively applied. Here, we demonstrate finite-state-machine-like processing embodied in cells using the walking behavior of Euplotes eurystomus, a ciliate that walks across surfaces using fourteen motile appendages (cirri). We found that cellular walking entails regulated transitions among a discrete set of gait states. The set of observed transitions decomposes into a small group of high-probability, temporally irreversible transitions and a large group of low-probability, time-symmetric transitions, thus revealing stereotypy in the sequential patterns of state transitions. Simulations and experiments suggest that the sequential logic of the gait is functionally important. Taken together, these findings implicate a finite-state-machine-like process. Cirri are connected by microtubule bundles (fibers), and we found that the dynamics of cirri involved in different state transitions are associated with the structure of the fiber system. Perturbative experiments revealed that the fibers mediate gait coordination, suggesting a mechanical basis of gait control.

Modular, cascade-like transcriptional program of regeneration in Stentor

The giant ciliate Stentor coeruleus is a classical model system for studying regeneration and morphogenesis in a single cell. The anterior of the cell is marked by an array of cilia, known as the oral apparatus, which can be induced to shed and regenerate in a series of reproducible morphological steps, previously shown to require transcription. If a cell is cut in half, each half regenerates an intact cell. We used RNA sequencing (RNAseq) to assay the dynamic changes in Stentor’s transcriptome during regeneration, after both oral apparatus shedding and bisection, allowing us to identify distinct temporal waves of gene expression including kinases, RNA -binding proteins, centriole biogenesis factors, and orthologs of human ciliopathy genes. By comparing transcriptional profiles of different regeneration events, we identified distinct modules of gene expression corresponding to oral apparatus regeneration, posterior holdfast regeneration, and recovery after wounding. By measuring gene expression after blocking translation, we show that the sequential waves of gene expression involve a cascade mechanism in which later waves of expression are triggered by translation products of early-expressed genes. Among the early-expressed genes, we identified an E2F transcription factor and the RNA-binding protein Pumilio as potential regulators of regeneration based on the expression pattern of their predicted target genes. RNAi-mediated knockdown experiments indicate that Pumilio is required for regenerating oral structures of the correct size. E2F is involved in the completion of regeneration but is dispensable for earlier steps. This work allows us to classify regeneration genes into groups based on their potential role for regeneration in distinct cell regeneration paradigms, and provides insight into how a single cell can coordinate complex morphogenetic pathways to regenerate missing structures.

Modeling cell biological features of meiotic chromosome pairing to study interlock resolution

During meiosis, homologous chromosomes become associated side by side in a process known as homologous chromosome pairing. Pairing requires long range chromosome motion through a nucleus that is full of other chromosomes. It remains unclear how the cell manages to align each pair of chromosomes quickly while mitigating and resolving interlocks. Here, we use a coarse-grained molecular dynamics model to investigate how specific features of meiosis, including motor-driven telomere motion, nuclear envelope interactions, and increased nuclear size, affect the rate of pairing and the mitigation/resolution of interlocks. By creating in silico versions of three yeast strains and comparing the results of our model to experimental data, we find that a more distributed placement of pairing sites along the chromosome is necessary to replicate experimental findings. Active motion of the telomeric ends speeds up pairing only if binding sites are spread along the chromosome length. Adding a meiotic bouquet significantly speeds up pairing but does not significantly change the number of interlocks. An increase in nuclear size slows down pairing while greatly reducing the number of interlocks. Interestingly, active forces increase the number of interlocks, which raises the question: How do these interlocks resolve? Our model gives us detailed movies of interlock resolution events which we then analyze to build a step-by-step recipe for interlock resolution. In our model, interlocks must first translocate to the ends, where they are held in a quasi-stable state by a large number of paired sites on one side. To completely resolve an interlock, the telomeres of the involved chromosomes must come in close proximity so that the cooperativity of pairing coupled with random motion causes the telomeres to unwind. Together our results indicate that computational modeling of homolog pairing provides insight into the specific cell biological changes that occur during meiosis.

Determining protein polarization proteome-wide using physical dissection of individual Stentor coeruleus cells

Cellular components are non-randomly arranged with respect to the shape and polarity of the whole cell.^1-4 Patterning within cells can extend down to the level of individual proteins and mRNA.^5,6 But how much of the proteome is actually localized with respect to cell polarity axes? Proteomics combined with cellular fractionation^7-11 has shown that most proteins localize to one or more organelles but does not tell us how many proteins have a polarized localization with respect to the large-scale polarity axes of the intact cell. Genome-wide localization studies in yeast^12-15 found that only a few percent of proteins have a localized position relative to the cell polarity axis defined by sites of polarized cell growth. Here, we describe an approach for analyzing protein distribution within a cell with a visibly obvious global patterning-the giant ciliate Stentor coeruleus.^16,17 Ciliates, including Stentor, have highly polarized cell shapes with visible surface patterning.^1,18 A Stentor cell is roughly 2 mm long, allowing a "proteomic dissection" in which microsurgery is used to separate cellular fragments along the anterior-posterior axis, followed by comparative proteomic analysis. In our analysis, 25% of the proteome, including signaling proteins, centrin/SFI proteins, and GAS2 orthologs, shows a polarized location along the cell's anterior-posterior axis. We conclude that a large proportion of all proteins are polarized with respect to global cell polarity axes and that proteomic dissection provides a simple and effective approach for spatial proteomics.

Biosynthesis of Linear Protein Nanoarrays Using the Flagellar Axoneme

Applications in biotechnology and synthetic biology often make use of soluble proteins, but there are many potential advantages of anchoring enzymes to a stable substrate, including stability and the possibility for substrate channeling. To avoid the necessity of protein purification and chemical immobilization, there has been growing interest in bio-assembly of protein-containing nanoparticles, exploiting the self-assembly of viral capsid proteins or other proteins that form polyhedral structures. However, these nanoparticles are limited in size, which constrains the packaging and the accessibility of the proteins. An axoneme, the insoluble protein core of the eukaryotic flagellum or cilium, is a highly ordered protein structure that can be several microns in length, orders of magnitude larger than other types of nanoparticles. We show that when proteins of interest are fused to specific axonemal proteins and expressed in living Chlamydomonas reinhardtii cells, they become incorporated into linear arrays, which have the advantages of high protein loading capacity and single-step purification with retention of biomass. The arrays can be isolated as membrane-enclosed vesicles or as exposed protein arrays. The approach is demonstrated for both a fluorescent protein and an enzyme (beta-lactamase), showing that incorporation into axonemes retains protein function in a stable, easily isolated array form.

A simple method to generate human airway epithelial organoids with externally orientated apical membranes

Organoids, which are self-organizing three-dimensional cultures, provide models that replicate specific cellular components of native tissues or facets of organ complexity. We describe a simple method to generate organoid cultures using isolated human tracheobronchial epithelial cells grown in mixed matrix components and supplemented at day 14 with the Wnt pathway agonist R-spondin 2 (RSPO2) and the bone morphogenic protein antagonist Noggin. In contrast to previous reports, our method produces differentiated tracheobronchospheres with externally orientated apical membranes without pretreatments, providing an epithelial model to study cilia formation and function, disease pathogenesis, and interaction of pathogens with the respiratory mucosa. Starting from 3 × 10⁵ cells, organoid yield at day 28 was 1,720 ± 302. Immunocytochemistry confirmed the cellular localization of airway epithelial markers, including CFTR, Na⁺/K⁺ ATPase, acetylated-α-tubulin, E-cadherin, and ZO-1. Compared to native tissues, expression of genes related to bronchial differentiation and ion transport were similar in organoid and air-liquid interface (ALI) cultures. In matched primary cultures, mean organoid cilia length was 6.1 ± 0.2 µm, similar to that of 5.7 ± 0.1 µm in ALI cultures, and ciliary beating was vigorous and coordinated with frequencies of 7.7 ± 0.3 Hz in organoid cultures and 5.3 ± 0.8 Hz in ALI cultures. Functional measurement of osmotically induced volume changes in organoids showed low water permeability. The generation of numerous single testable units from minimal starting material complements prior techniques. This culture system may be useful for studying airway biology and pathophysiology, aiding diagnosis of ciliopathies, and potentially for high-throughput drug screening.

Regeneration in Stentor coeruleus

We often think about regeneration in terms of replacing missing structures, such as organs or tissues, with new structures generated via cell proliferation and differentiation. But at a smaller scale, single cells, themselves, are capable of regenerating when part of the cell has been removed. A classic model organism that facilitates the study of cellular regeneration in the giant ciliate Stentor coeruleus. These cells, which can grow to more than a millimeter in size, have the ability to survive after extensive wounding of their surface, and are able to regenerate missing structures. Even a small piece of a cell can regenerate a whole cell with normal geometry, in a matter of hours. Such regeneration requires cells to be able to trigger organelle biogenesis in response to loss of structures. But subcellular regeneration also relies on intracellular mechanisms to create and maintain global patterning within the cell. These mechanisms are not understood, but at a conceptual level they involve processes that resemble those seen in animal development and regeneration. Here we discuss single-celled regeneration in Stentor from the viewpoint of standard regeneration paradigms in animals. For example, there is evidence that regeneration of the oral apparatus in Stentor follows a sender-receiver model similar to crustacean eyestalk regeneration. By drawing these analogies, we find that many of the concepts already known from the study of animal-scale regeneration and development can be applied to the study of regeneration at the cellular level, such as the concepts of determination, induction, mosaic vs. regulative development, and epimorphosis vs. morphallaxis. We propose that the similarities may go beyond analogy, and that some aspects of animal development and regeneration may have evolved by exploiting pre-existing subcellular developmental strategies from unicellular ancestors.

Analysis of Motility Patterns of Stentor During and After Oral Apparatus Regeneration using Cell Tracking

Stentor coeruleus is a well-known model organism for the study of unicellular regeneration. Transcriptomic analysis of individual cells revealed hundreds of genes-many not associated with the oral apparatus (OA)-that are differentially regulated in phases throughout the regeneration process. It was hypothesized that this systemic reorganization and mobilization of cellular resources towards growth of a new OA will lead to observable changes in movement and behavior corresponding in time to the phases of differential gene expression. However, the morphological complexity of S. coeruleus necessitated the development of an assay to capture the statistics and timescale. A custom script was used to track cells in short videos, and statistics were compiled over a large population (N ~100). Upon loss of the OA, S. coeruleus initially loses the ability for directed motion; then starting at ~4 h, it exhibits a significant drop in speed until ~8 h. This assay provides a useful tool for the screening of motility phenotypes and can be adapted for the investigation of other organisms.

Analysis of biological noise in the flagellar length control system

Any proposed mechanism for organelle size control should be able to account not only for average size but also for the variation in size. We analyzed cell-to-cell variation and within-cell variation of length for the two flagella in Chlamydomonas, finding that cell-to-cell variation is dominated by cell size, whereas within-cell variation results from dynamic fluctuations. Fluctuation analysis suggests tubulin assembly is not directly coupled with intraflagellar transport (IFT) and that the observed length fluctuations reflect tubulin assembly and disassembly events involving large numbers of tubulin dimers. Length variation is increased in long-flagella mutants, an effect consistent with theoretical models for flagellar length regulation. Cells with unequal flagellar lengths show impaired swimming but improved gliding, raising the possibility that cells have evolved mechanisms to tune biological noise in flagellar length. Analysis of noise at the level of organelle size provides a way to probe the mechanisms determining cell geometry.

Microfluidic guillotine reveals multiple timescales and mechanical modes of wound response in Stentor coeruleus

Background

Wound healing is one of the defining features of life and is seen not only in tissues but also within individual cells. Understanding wound response at the single-cell level is critical for determining fundamental cellular functions needed for cell repair and survival. This understanding could also enable the engineering of single-cell wound repair strategies in emerging synthetic cell research. One approach is to examine and adapt self-repair mechanisms from a living system that already demonstrates robust capacity to heal from large wounds. Towards this end, Stentor coeruleus, a single-celled free-living ciliate protozoan, is a unique model because of its robust wound healing capacity. This capacity allows one to perturb the wounding conditions and measure their effect on the repair process without immediately causing cell death, thereby providing a robust platform for probing the self-repair mechanism.

Results

Here we used a microfluidic guillotine and a fluorescence-based assay to probe the timescales of wound repair and of mechanical modes of wound response in Stentor. We found that Stentor requires ~ 100–1000 s to close bisection wounds, depending on the severity of the wound. This corresponds to a healing rate of ~ 8–80 μm²/s, faster than most other single cells reported in the literature. Further, we characterized three distinct mechanical modes of wound response in Stentor: contraction, cytoplasm retrieval, and twisting/pulling. Using chemical perturbations, active cilia were found to be important for only the twisting/pulling mode. Contraction of myonemes, a major contractile fiber in Stentor, was surprisingly not important for the contraction mode and was of low importance for the others.

Conclusions

While events local to the wound site have been the focus of many single-cell wound repair studies, our results suggest that large-scale mechanical behaviors may be of greater importance to single-cell wound repair than previously thought. The work here advances our understanding of the wound response in Stentor and will lay the foundation for further investigations into the underlying components and molecular mechanisms involved.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-021-00970-0.

Deep Convolutional and Recurrent Neural Networks for Cell Motility Discrimination and Prediction

Cells in culture display diverse motility behaviors that may reflect differences in cell state and function, providing motivation to discriminate between different motility behaviors. Current methods to do so rely upon manual feature engineering. However, the types of features necessary to distinguish between motility behaviors can vary greatly depending on the biological context, and it is not always clear which features may be most predictive in each setting for distinguishing particular cell types or disease states. Convolutional neural networks (CNNs) are machine learning models allowing for relevant features to be learned directly from spatial data. Similarly, recurrent neural networks (RNNs) are a class of models capable of learning long term temporal dependencies. Given that cell motility is inherently spacio-temporal data, we present an approach utilizing both convolutional and long- short-term memory (LSTM) recurrent neural network units to analyze cell motility data. These RNN models provide accurate classification of simulated motility and experimentally measured motility from multiple cell types, comparable to results achieved with hand-engineered features. The variety of cell motility differences we can detect suggests that the algorithm is generally applicable to additional cell types not analyzed here. RNN autoencoders based on the same architecture are capable of learning motility features in an unsupervised manner and capturing variation between myogenic cells in the latent space. Adapting these RNN models to motility prediction, RNNs are capable of predicting muscle stem cell motility from past tracking data with performance superior to standard motion prediction models. This advance in cell motility prediction may be of practical utility in cell tracking applications.

Scaling of Subcellular Structures

As cells grow, the size and number of their internal organelles increase in order to keep up with increased metabolic requirements. Abnormal size of organelles is a hallmark of cancer and an important aspect of diagnosis in cytopathology. Most organelles vary in either size or number, or both, as a function of cell size, but the mechanisms that create this variation remain unclear. In some cases, organelle size appears to scale with cell size through processes of relative growth, but in others the size may be set by either active measurement systems or genetic programs that instruct organelle biosynthetic activities to create organelles of a size appropriate to a given cell type.

Simple Rules Determine Distinct Patterns of Branching Morphogenesis

Many metazoan organs are comprised of branching trees of epithelial tubes; how patterning occurs in these trees is a fundamental problem of development. Commonly, branch tips fill the volume of the organ approximately uniformly, e.g., in mammalian lung, airway branch tips are dispersed roughly uniformly throughout the volume of the lung. In contrast, in the developing metanephric kidney, the tips of the ureteric bud tree are located close to the outer surface of the kidney rather than filling the kidney. Here, we describe a simple alteration in the branching rules that accounts for the difference between the kidney pattern that leads to tips near the organ surface versus previously known patterns that lead to the branch tips being dispersed throughout the organ. We further use a simple toy model to deduce from first principles how this rule change accounts for the differences in the two types of trees.

Cellular Cognition: Sequential Logic in a Giant Protist

Quantitative analysis of the giant ciliate Stentor roeselii shows that a single cell can make decisions, based on the ability to switch between several different behaviors in a non-random order.

Aging induces aberrant state transition kinetics in murine muscle stem cells

Murine muscle stem cells (MuSCs) experience a transition from quiescence to activation that is required for regeneration, but it remains unknown if the trajectory and dynamics of activation change with age. Here, we use time-lapse imaging and single cell RNA-seq to measure activation trajectories and rates in young and aged MuSCs. We find that the activation trajectory is conserved in aged cells, and we develop effective machine-learning classifiers for cell age. Using cell-behavior analysis and RNA velocity, we find that activation kinetics are delayed in aged MuSCs, suggesting that changes in stem cell dynamics may contribute to impaired stem cell function with age. Intriguingly, we also find that stem cell activation appears to be a random walk-like process, with frequent reversals, rather than a continuous linear progression. These results support a view of the aged stem cell phenotype as a combination of differences in the location of stable cell states and differences in transition rates between them.

Speed and Diffusion of Kinesin-2 Are Competing Limiting Factors in Flagellar Length-Control Model

Flagellar length control in Chlamydomonas is a tractable model system for studying the general question of organelle size regulation. We have previously proposed that the diffusive return of the kinesin motor that powers intraflagellar transport can play a key role in length regulation. Here, we explore how the motor speed and diffusion coefficient for the return of kinesin-2 affect flagellar growth kinetics. We find that the system can exist in two distinct regimes, one dominated by motor speed and one by diffusion coefficient. Depending on length, a flagellum can switch between these regimes. Our results indicate that mutations can affect the length in distinct ways. We discuss our theory's implication for flagellar growth influenced by beating and provide possible explanations for the experimental observation that a beating flagellum is usually longer than its immotile mutant. These results demonstrate how our simple model can suggest explanations for mutant phenotypes.

Aurora kinase inhibitors delay regeneration in Stentor coeruleus at an intermediate step

The giant unicellular ciliate Stentor coeruleus can be cut into pieces and each piece will regenerate into a healthy, full-sized individual. The molecular mechanism for how Stentor regenerates is a complete mystery, however, the process of regeneration shows striking similarities to the process of cell division. On a morphological level, the process of creating a second mouth in division or a new oral apparatus in regeneration have the same steps and occur in the same order. On the transcriptional level, genes encoding elements of the cell division and cell cycle regulatory machinery, including Aurora kinases, are differentially expressed during regeneration. This suggests that there may be some common regulatory mechanisms involved in both regeneration and cell division. If the cell cycle machinery really does play a role in regeneration, then inhibition of proteins that regulate the timing of cell division may also affect the timing of regeneration in Stentor. Here we show that two well-characterized Aurora kinase A+B inhibitors that affect the timing of regeneration. ZM447439 slows down regeneration by at least one hour. PF03814735 completely suppresses regeneration until the drug is removed. Here we provide the first direct experimental evidence that Stentor may harness the cell division machinery to regulate the sequential process of regeneration.

Towards computer-aided design of cellular structure

Cells are complex machines with tremendous potential for applications in medicine and biotechnology. Although much effort has been devoted to engineering the metabolic, genetic, and signaling pathways of cells, methods for systematically engineering the physical structure of cells are less developed. Here we consider how coarse-grained models for cellular geometry at the organelle level can be used to build computer-aided design (CAD) tools for cellular structure.

Reorganization of complex ciliary flows around regenerating Stentor coeruleus

The phenomenon of ciliary coordination has garnered increasing attention in recent decades and multiple theories have been proposed to explain its occurrence in different biological systems. While hydrodynamic interactions are thought to dictate the large-scale coordinated activity of epithelial cilia for fluid transport, it is rather basal coupling that accounts for synchronous swimming gaits in model microeukaryotes such as Chlamydomonas. Unicellular ciliates present a fascinating yet understudied context in which coordination is found to persist in ciliary arrays positioned across millimetre scales on the same cell. Here, we focus on the ciliate Stentor coeruleus, chosen for its large size, complex ciliary organization, and capacity for cellular regeneration. These large protists exhibit ciliary differentiation between cortical rows of short body cilia used for swimming, and an anterior ring of longer, fused cilia called the membranellar band (MB). The oral cilia in the MB beat metachronously to produce strong feeding currents. Remarkably, upon injury, the MB can be shed and regenerated de novo. Here, we follow and track this developmental sequence in its entirety to elucidate the emergence of coordinated ciliary beating: from band formation, elongation, curling and final migration towards the cell anterior. We reveal a complex interplay between hydrodynamics and ciliary restructuring in Stentor, and highlight for the first time the importance of a ring-like topology for achieving long-range metachronism in ciliated structures.

This article is part of the Theo Murphy meeting issue ‘Unity and diversity of cilia in locomotion and transport’.

Testing the role of intraflagellar transport in flagellar length control using length-altering mutants of Chlamydomonas

Cilia and flagella are ideal model organelles in which to study the general question of organelle size control. Flagellar microtubules are steady-state structures whose size is set by the balance of assembly and disassembly. Assembly requires intraflagellar transport (IFT), and measurements of IFT have shown that the rate of entry of IFT particles into the flagellum is a decreasing function of length. It has been proposed that this length dependence of IFT may be the basis for flagellar length control. Here, we test this idea by showing that three different long-flagella mutations in Chlamydomonas all cause increased IFT injection, thus confirming that IFT can influence length control. However, quantitative comparisons with mathematical models suggest that the increase in injection is not sufficient to explain the full increase in length seen in these mutants; hence, some other mechanism may be at work. One alternative mechanism that has been proposed is length-regulated binding of tubulin to the IFT particles. However, we find that the apparent length dependence of tubulin loading that has previously been reported may actually reflect length-dependent organization of IFT trains.

This article is part of the Theo Murphy meeting issue ‘Unity and diversity of cilia in locomotion and transport’.

Microtubules are necessary for proper Reticulon localization during mitosis

During mitosis, the structure of the Endoplasmic Reticulum (ER) displays a dramatic reorganization and remodeling, however, the mechanism driving these changes is poorly understood. Hairpin-containing ER transmembrane proteins that stabilize ER tubules have been identified as possible factors to promote these drastic changes in ER morphology. Recently, the Reticulon and REEP family of ER shaping proteins have been shown to heavily influence ER morphology by driving the formation of ER tubules, which are known for their close proximity with microtubules. Here, we examine the role of microtubules and other cytoskeletal factors in the dynamics of a Drosophila Reticulon, Reticulon-like 1 (Rtnl1), localization to spindle poles during mitosis in the early embryo. At prometaphase, Rtnl1 is enriched to spindle poles just prior to the ER retention motif KDEL, suggesting a possible recruitment role for Rtnl1 in the bulk localization of ER to spindle poles. Using image analysis-based methods and precise temporal injections of cytoskeletal inhibitors in the early syncytial Drosophila embryo, we show that microtubules are necessary for proper Rtnl1 localization to spindles during mitosis. Lastly, we show that astral microtubules, not microfilaments, are necessary for proper Rtnl1 localization to spindle poles, and is largely independent of the minus-end directed motor protein dynein. This work highlights the role of the microtubule cytoskeleton in Rtnl1 localization to spindles during mitosis and sheds light on a pathway towards inheritance of this major organelle.

Cell learning

Although we often think of cells as small, simple building blocks of life, in fact they are highly complex and can perform a startling variety of functions. In our bodies, cells are programmed by complex differentiation pathways and are capable of responding to a bewildering range of chemical and physical signals. Free-living single-celled organisms, such as bacteria or protists, have to cope with varying environments, locate prey and potential mates, and escape from predators - all of the same tasks that a free-living animal is faced with. When animals face complex behavioral challenges, they rely on their cognitive abilities - the ability to learn from experience, to analyse a situation and choose an appropriate course of action. This ability is essential for survival and should, in principle, be a ubiquitous feature of all living things regardless of the complexity of the organism.

Modeling meiotic chromosome pairing: a tug of war between telomere forces and a pairing-based Brownian ratchet leads to increased pairing fidelity

Meiotic homolog pairing involves associations between homologous DNA regions scattered along the length of a chromosome. When homologs associate, they tend to do so by a processive zippering process, which apparently results from avidity effects. Using a computational model, we show that this avidity-driven processive zippering reduces the selectivity of pairing. When active random forces are applied to telomeres, this drop in selectivity is eliminated in a force-dependent manner. Further simulations suggest that active telomere forces are engaged in a tug-of-war against zippering, which can be interpreted as a Brownian ratchet with a stall force that depends on the dissociation constant of pairing. When perfectly homologous regions of high affinity compete with homeologous regions of lower affinity, the affinity difference can be amplified through this tug of war effect provided the telomere force acts in a range that is strong enough to oppose zippering of homeologs while still permitting zippering of correct homologs. The degree of unzippering depends on the radius of the nucleus, such that complete unzippering of homeologous regions can only take place if the nucleus is large enough to pull the two chromosomes completely apart. A picture of meiotic pairing thus emerges that is fundamentally mechanical in nature, possibly explaining the purpose of active telomere forces, increased nuclear diameter, and the presence of 'Maverick' chromosomes in meiosis.

Diffusion as a Ruler: Modeling Kinesin Diffusion as a Length Sensor for Intraflagellar Transport

An important question in cell biology is whether cells are able to measure size, either whole cell size or organelle size. Perhaps cells have an internal chemical representation of size that can be used to precisely regulate growth, or perhaps size is just an accident that emerges due to constraint of nutrients. The eukaryotic flagellum is an ideal model for studying size sensing and control because its linear geometry makes it essentially one-dimensional, greatly simplifying mathematical modeling. The assembly of flagella is regulated by intraflagellar transport (IFT), in which kinesin motors carry cargo adaptors for flagellar proteins along the flagellum and then deposit them at the tip, lengthening the flagellum. The rate at which IFT motors are recruited to begin transport into the flagellum is anticorrelated with the flagellar length, implying some kind of communication between the base and the tip and possibly indicating that cells contain some mechanism for measuring flagellar length. Although it is possible to imagine many complex scenarios in which additional signaling molecules sense length and carry feedback signals to the cell body to control IFT, might the already-known components of the IFT system be sufficient to allow length dependence of IFT? Here we investigate a model in which the anterograde kinesin motors unbind after cargo delivery, diffuse back to the base, and are subsequently reused to power entry of new IFT trains into the flagellum. By mathematically modeling and simulating such a system, we are able to show that the diffusion time of the motors can in principle be sufficient to serve as a proxy for length measurement. We found that the diffusion model can not only achieve a stable steady-state length without the addition of any other signaling molecules or pathways, but also is able to produce the anticorrelation between length and IFT recruitment rate that has been observed in quantitative imaging studies.

Electron cryo-tomography provides insight into procentriole architecture and assembly mechanism

Centriole is an essential structure with multiple functions in cellular processes. Centriole biogenesis and homeostasis is tightly regulated. Using electron cryo-tomography (cryoET) we present the structure of procentrioles from Chlamydomonas reinhardtii. We identified a set of non-tubulin components attached to the triplet microtubule (MT), many are at the junctions of tubules likely to reinforce the triplet. We describe structure of the A-C linker that bridges neighboring triplets. The structure infers that POC1 is likely an integral component of A-C linker. Its conserved WD40 β-propeller domain provides attachment sites for other A-C linker components. The twist of A-C linker results in an iris diaphragm-like motion of the triplets in the longitudinal direction of procentriole. Finally, we identified two assembly intermediates at the growing ends of procentriole allowing us to propose a model for the procentriole assembly. Our results provide a comprehensive structural framework for understanding the molecular mechanisms underpinning procentriole biogenesis and assembly.

A Dilution Model for Embryonic Scaling

A recent study in Nature Cell Biology (Almuedo-Castillo et al., 2018) describes a mechanism for tissue scaling in zebrafish embryos. The authors show that a fixed relative amount of the diffusible Nodal inhibitor Lefty produces an extended gradient in larger embryos, ensuring proportionally scaled germ layers, irrespective of embryo size.

Methods for the Study of Regeneration in Stentor

Cells need to be able to regenerate their parts to recover from external perturbations. The unicellular ciliate Stentor coeruleus is an excellent model organism to study wound healing and subsequent cell regeneration. The Stentor genome became available recently, along with modern molecular biology methods, such as RNAi. These tools make it possible to study single-cell regeneration at the molecular level. The first section of the protocol covers establishing Stentor cell cultures from single cells or cell fragments, along with general guidelines for maintaining Stentor cultures. Culturing Stentor in large quantities allows for the use of valuable tools like biochemistry, sequencing, and mass spectrometry. Subsequent sections of the protocol cover different approaches to inducing regeneration in Stentor. Manually cutting cells with a glass needle allows studying the regeneration of large cell parts, while treating cells with either sucrose or urea allows studying the regeneration of specific structures located at the anterior end of the cell. A method for imaging individual regenerating cells is provided, along with a rubric for staging and analyzing the dynamics of regeneration. The entire process of regeneration is divided in three stages. By visualizing the dynamics of the progression of a population of cells through the stages, the heterogeneity in regeneration timing is demonstrated.

Self-repairing cells: How single cells heal membrane ruptures and restore lost structures

Many organisms and tissues display the ability to heal and regenerate as needed for normal physiology and as a result of pathogenesis. However, these repair activities can also be observed at the single-cell level. The physical and molecular mechanisms by which a cell can heal membrane ruptures and rebuild damaged or missing cellular structures remain poorly understood. This Review presents current understanding in wound healing and regeneration as two distinct aspects of cellular self-repair by examining a few model organisms that have displayed robust repair capacity, including Xenopus oocytes, Chlamydomonas, and Stentor coeruleus Although many open questions remain, elucidating how cells repair themselves is important for our mechanistic understanding of cell biology. It also holds the potential for new applications and therapeutic approaches for treating human disease.

Mechanical Forces Program the Orientation of Cell Division during Airway Tube Morphogenesis

Oriented cell division plays a key role in controlling organogenesis. The mechanisms for regulating division orientation at the whole-organ level are only starting to become understood. By combining 3D time-lapse imaging, mouse genetics, and mathematical modeling, we find that global orientation of cell division is the result of a combination of two types of spindles with distinct spindle dynamic behaviors in the developing airway epithelium. Fixed spindles follow the classic long-axis rule and establish their division orientation before metaphase. In contrast, rotating spindles do not strictly follow the long-axis rule and determine their division orientation during metaphase. By using both a cell-based mechanical model and stretching-lung-explant experiments, we showed that mechanical force can function as a regulatory signal in maintaining the stable ratio between fixed spindles and rotating spindles. Our findings demonstrate that mechanical forces, cell geometry, and oriented cell division function together in a highly coordinated manner to ensure normal airway tube morphogenesis.

Inferring cell state by quantitative motility analysis reveals a dynamic state system and broken detailed balance

Cell populations display heterogeneous and dynamic phenotypic states at multiple scales. Similar to molecular features commonly used to explore cell heterogeneity, cell behavior is a rich phenotypic space that may allow for identification of relevant cell states. Inference of cell state from cell behavior across a time course may enable the investigation of dynamics of transitions between heterogeneous cell states, a task difficult to perform with destructive molecular observations. Cell motility is one such easily observed cell behavior with known biomedical relevance. To investigate heterogenous cell states and their dynamics through the lens of cell behavior, we developed Heteromotility, a software tool to extract quantitative motility features from timelapse cell images. In mouse embryonic fibroblasts (MEFs), myoblasts, and muscle stem cells (MuSCs), Heteromotility analysis identifies multiple motility phenotypes within the population. In all three systems, the motility state identity of individual cells is dynamic. Quantification of state transitions reveals that MuSCs undergoing activation transition through progressive motility states toward the myoblast phenotype. Transition rates during MuSC activation suggest non-linear kinetics. By probability flux analysis, we find that this MuSC motility state system breaks detailed balance, while the MEF and myoblast systems do not. Balanced behavior state transitions can be captured by equilibrium formalisms, while unbalanced switching between states violates equilibrium conditions and would require an external driving force. Our data indicate that the system regulating cell behavior can be decomposed into a set of attractor states which depend on the identity of the cell, together with a set of transitions between states. These results support a conceptual view of cell populations as dynamical systems, responding to inputs from signaling pathways and generating outputs in the form of state transitions and observable motile behaviors.

Testing the time-of-flight model for flagellar length sensing

A combination of quantitative imaging, modeling, and genetics has been used to test a proposed mechanism for measuring the size of an organelle. One way to measure distance is to send a clock out on a train and measure the elapsed time when the train returns. We tested a molecular version of this model as a possible regulator of intraflagellar transport by altering the return speed of the transport machinery and probing the effect on a known length-dependent process.

Cilia and flagella are microtubule-based organelles that protrude from the surface of most cells, are important to the sensing of extracellular signals, and make a driving force for fluid flow. Maintenance of flagellar length requires an active transport process known as intraflagellar transport (IFT). Recent studies reveal that the amount of IFT injection negatively correlates with the length of flagella. These observations suggest that a length-dependent feedback regulates IFT. However, it is unknown how cells recognize the length of flagella and control IFT. Several theoretical models try to explain this feedback system. We focused on one of the models, the “time-of-flight” model, which measures the length of flagella on the basis of the travel time of IFT protein in the flagellar compartment. We tested the time-of-flight model using Chlamydomonas dynein mutant cells, which show slower retrograde transport speed. The amount of IFT injection in dynein mutant cells was higher than that in control cells. This observation does not support the prediction of the time-of-flight model and suggests that Chlamydomonas uses another length-control feedback system rather than that described by the time-of-flight model.

Microfluidic guillotine for single-cell wound repair studies

Wound repair is a key feature distinguishing living from nonliving matter. Single cells are increasingly recognized to be capable of healing wounds. The lack of reproducible, high-throughput wounding methods has hindered single-cell wound repair studies. This work describes a microfluidic guillotine for bisecting single Stentor coeruleus cells in a continuous-flow manner. Stentor is used as a model due to its robust repair capacity and the ability to perform gene knockdown in a high-throughput manner. Local cutting dynamics reveals two regimes under which cells are bisected, one at low viscous stress where cells are cut with small membrane ruptures and high viability and one at high viscous stress where cells are cut with extended membrane ruptures and decreased viability. A cutting throughput up to 64 cells per minute-more than 200 times faster than current methods-is achieved. The method allows the generation of more than 100 cells in a synchronized stage of their repair process. This capacity, combined with high-throughput gene knockdown in Stentor, enables time-course mechanistic studies impossible with current wounding methods.

Non-model model organisms

Model organisms are widely used in research as accessible and convenient systems to study a particular area or question in biology. Traditionally only a handful of organisms have been widely studied, but modern research tools are enabling researchers to extend the set of model organisms to include less-studied and more unusual systems. This Forum highlights a range of 'non-model model organisms' as emerging systems for tackling questions across the whole spectrum of biology (and beyond), the opportunities and challenges, and the outlook for the future.

Intraflagellar Transport and Ciliary Dynamics

Cilia and flagella are microtubule-based organelles whose assembly requires a motile process, known as intraflagellar transport (IFT), to bring tubulin and other components to the distal tip of the growing structure. The IFT system uses a multiprotein complex with components that appear to be specialized for the transport of different sets of cargo proteins. The mechanisms by which cargo is selected for ciliary import and transport by IFT remain an area of active research. The complex dynamics of cilia and flagella are under constant regulation to ensure proper length control, and this regulation appears to involve regulation at the stage of IFT injection into the flagellum, as well as regulation of flagellar disassembly and, possibly, of cargo binding. Cilia and flagella thus represent a convenient model system to study how multiple motile and signaling pathways cooperate to control the assembly and dynamics of a complex cellular structure.

The Macronuclear Genome of Stentor coeruleus Reveals Tiny Introns in a Giant Cell

The giant, single-celled organism Stentor coeruleus has a long history as a model system for studying pattern formation and regeneration in single cells. Stentor [1, 2] is a heterotrichous ciliate distantly related to familiar ciliate models, such as Tetrahymena or Paramecium. The primary distinguishing feature of Stentor is its incredible size: a single cell is 1 mm long. Early developmental biologists, including T.H. Morgan [3], were attracted to the system because of its regenerative abilities-if large portions of a cell are surgically removed, the remnant reorganizes into a normal-looking but smaller cell with correct proportionality [2, 3]. These biologists were also drawn to Stentor because it exhibits a rich repertoire of behaviors, including light avoidance, mechanosensitive contraction, food selection, and even the ability to habituate to touch, a simple form of learning usually seen in higher organisms [4]. While early microsurgical approaches demonstrated a startling array of regenerative and morphogenetic processes in this single-celled organism, Stentor was never developed as a molecular model system. We report the sequencing of the Stentor coeruleus macronuclear genome and reveal key features of the genome. First, we find that Stentor uses the standard genetic code, suggesting that ciliate-specific genetic codes arose after Stentor branched from other ciliates. We also discover that ploidy correlates with Stentor's cell size. Finally, in the Stentor genome, we discover the smallest spliceosomal introns reported for any species. The sequenced genome opens the door to molecular analysis of single-cell regeneration in Stentor.

Organelles - understanding noise and heterogeneity in cell biology at an intermediate scale

Many studies over the years have shown that non-genetic mechanisms for producing cell-to-cell variation can lead to highly variable behaviors across genetically identical populations of cells. Most work to date has focused on gene expression noise as the primary source of phenotypic heterogeneity, yet other sources may also contribute. In this Commentary, we explore organelle-level heterogeneity as a potential secondary source of cellular 'noise' that contributes to phenotypic heterogeneity. We explore mechanisms for generating organelle heterogeneity and present evidence of functional links between organelle morphology and cellular behavior. Given the many instances in which molecular-level heterogeneity has been linked to phenotypic heterogeneity, we posit that organelle heterogeneity may similarly contribute to overall phenotypic heterogeneity and underline the importance of studying organelle heterogeneity to develop a more comprehensive understanding of phenotypic heterogeneity. Finally, we conclude with a discussion of the medical challenges associated with phenotypic heterogeneity and outline how improved methods for characterizing and controlling this heterogeneity may lead to improved therapeutic strategies and outcomes for patients.

A systematic comparison of mathematical models for inherent measurement of ciliary length: how a cell can measure length and volume

Cells control organelle size with great precision and accuracy to maintain optimal physiology, but the mechanisms by which they do so are largely unknown. Cilia and flagella are simple organelles in which a single measurement, length, can represent size. Maintenance of flagellar length requires an active transport process known as intraflagellar transport, and previous measurements suggest that a length-dependent feedback regulates intraflagellar transport. But the question remains: how is a length-dependent signal produced to regulate intraflagellar transport appropriately? Several conceptual models have been suggested, but testing these models quantitatively requires that they be cast in mathematical form. Here, we derive a set of mathematical models that represent the main broad classes of hypothetical size-control mechanisms currently under consideration. We use these models to predict the relation between length and intraflagellar transport, and then compare the predicted relations for each model with experimental data. We find that three models-an initial bolus formation model, an ion current model, and a diffusion-based model-show particularly good agreement with available experimental data. The initial bolus and ion current models give mathematically equivalent predictions for length control, but fluorescence recovery after photobleaching experiments rule out the initial bolus model, suggesting that either the ion current model or a diffusion-based model is more likely correct. The general biophysical principles of the ion current and diffusion-based models presented here to measure cilia and flagellar length can be generalized to measure any membrane-bound organelle volume, such as the nucleus and endoplasmic reticulum.

Cell Geometry: How Cells Count and Measure Size

The cell represents a highly organized state of living matter in which numerous geometrical parameters are under dynamic regulation in order to match the form of a cell with its function. Cells appear capable of regulating not only the total quantity of their internal organelles, but also the size and number of those organelles. The regulation of three parameters, size, number, and total quantity, can in principle be accomplished by regulating the production or growth of organelles, their degradation or disassembly, and their partitioning among daughter cells during division. Any or all of these steps could in principle be under regulation. But if organelle assembly or disassembly is regulated by number or size, how would the cell know how many copies of an organelle it has, or how big they are?