Search

Applications in biotechnology and synthetic biology often make use of soluble proteins, but there are many potential advantages to 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. But these nanoparticle are limited in size which constrains the packaging and the accessibility of the proteins. The 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 cells, they become incorporated into linear arrays which have the advantages of high protein loading capacity, high stability, and single-step purification with retention of biomass. The arrays can be isolated as membrane enclosed vesicle or as exposed protein arrays. The approach is demonstrated for both fluorescent proteins and enzymes, and in the latter case it is found that incorporation into axoneme arrays provides increased stability for the enzyme.

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 × 105 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.

The flagellar axoneme is a complex protein-based machine capable of generating motile forces by coordinating the action of thousands of dynein motors. A key element of the axoneme is the central pair apparatus, consisting of a pair of microtubules surrounded by additional structures. In an effort to better understand the organization of the central pair, we used 2D DIGE to identify proteins that are depleted from flagella isolated from two differentChlamydomonas reinhardtiimutants,pf15andpf18, that lack the central pair. The set of proteins contained almost no components of the central apparatus. We find that three proteins of the radial spoke head RSP1, RSP9, and RSP10, as well as a number of other protein components associated with the outer doublets, are depleted from flagella of mutants lacking the central apparatus. Two of the other proteins depleted frompf15andpf18flagella, the microtubule inner proteins (MIPs) FAP21 and FAP161, are missing from the genome ofThalassiosira, an organism that lacks a central pair and radial spokes, and RNAi of FAP21 in planaria shows that it has a role in ciliary motility. Based on the depletion of radial spoke head proteins, as well as MIPs and other axonemal components, from flagella lacking the central pair apparatus, we hypothesize that the central apparatus may play a role in scaffolding the assembly or retention of radial spokes and other axonemal structures.

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. Here 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. We show that both models lead to stable flagellar length control, can reproduce the observed dynamics of gene expression, and are consistent with the effects of protein synthesis inhibitors on 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 FLA3, FLA8, and FLA10, which encode the three components of the heterotrimer 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. In contrast to the effects of IFT mutants, we find that inhibiting axonemal assembly has much less effect on gene expression, suggesting that transport is more important than axonemal assembly. The repressor sequestration model allows precursor production to occur when flagella are growing rapidly, representing a form of derivative control.

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 Ca2+channels along the flagellum and that the Ca2+current from the flagellum into the cell body increases linearly with flagellar length. In this model, the cell uses the Ca2+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 Ca2+inside ofChlamydomonasflagella and quantified IFT injection. Although the concentration of Ca2+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 Ca2+spikes. Thus, Ca2+does not appear to function as a negative regulator of IFT injection, hence it cannot form the basis of a stable length control system.

Mucus clearance constitutes the primary defence of the respiratory system against viruses, bacteria and environmental insults. This transport across the entire airway emerges from the integrated activity of thousands of multiciliated cells, each containing hundreds of cilia, which together must coordinate their spatial arrangement, alignment and motility. The mechanisms of fluid transport have been studied extensively at the level of an individual cilium, collectively moving metachronal waves and, more generally, the hydrodynamics of active matter. However, the connection between local cilia architecture and the topology of the flows they generate remains largely unexplored. Here, we image the mouse airway from subcellular (nm) to organ (mm) scales, characterizing quantitatively its ciliary arrangement and the generated flows. Locally, we measure heterogeneity in both cilia organization and flow structure, but, across the trachea, fluid transport is coherent. To examine this result, a hydrodynamic model was developed for a systematic exploration of different tissue architectures. Surprisingly, we find that disorder enhances particle clearance, whether it originates from fluctuations, heterogeneity in multiciliated cell arrangement or ciliary misalignment. This resembles elements of ‘stochastic resonance’, in the sense that noise can improve the function of the system. Taken together, our results shed light on how the microstructure of an active carpet determines its emergent dynamics. Furthermore, this work is also directly applicable to human airway pathologies, which are the third leading cause of deaths worldwide. Fluid flow through airways—necessary to keep lungs healthy and free from particles—occurs thanks to moving cilia. Here the authors show that defects in the arrangement of these cilia can facilitate particle clearance through the lungs.

Length control of flagella represents a simple and tractable system to investigate the dynamics of organelle size. Models for flagellar length control in the model organism, Chlamydomonas reinhardtii have focused on the length-dependence of the intraflagellar transport (IFT) system which manages the delivery and removal of axonemal subunits at the tip of the flagella. One of these cargoes, tubulin, is the major axonemal subunit, and its frequency of arrival at the tip plays a central role in size control models. However, the mechanisms determining tubulin dynamics at the tip are still poorly understood. We discovered a loss-of-function mutation that leads to shortened flagella, and found that this was an allele of a previously described gene, SHF1, whose molecular identity had not previously been determined. We found that SHF1 encodes a Chlamydomonas ortholog of Crescerin, previously identified as a cilia-specific TOG-domain array protein that can bind tubulin via its TOG domains and increase tubulin polymerization rates. In this mutant, flagellar regeneration occurs with the same initial kinetics as wild-type cells, but plateaus at a shorter length. Using a computational model in which the flagellar microtubules are represented by a differential equation for flagellar length combined with a stochastic model for cytoplasmic microtubule dynamics, we found that our experimental results are best described by a model in which Crescerin/SHF1 binds tubulin dimers in the cytoplasm and transports them into the flagellum. We suggest that this TOG-domain protein is necessary to efficiently and preemptively increase intra-flagella tubulin levels to offset decreasing IFT cargo at the tip as flagellar assembly progresses.

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.Author SummaryEarly in meiosis, homologous chromosomes must find each other within the crowded nuclear space and become aligned along their entire length in a process known as homologous chromosome pairing. It remains unclear how the cell manages to align each pair of chromosomes quickly while mitigating and resolving interlocks. Here, we study this process by using a computational model. Our model attempts to capture the large-scale cell biological picture of meiotic pairing including the random initial 3D search, active motion of the chromosome ends, and meiosis specific constraints such as telomere attachment to the nuclear envelope. We use our model to study how these different features of meiosis affect the rate of pairing and the mitigation/resolution of interlocks. Importantly, our model gives us detailed movies of interlock resolution events, which we then analyze to build a step-by-step recipe for interlock resolution. We believe computational modeling of homolog pairing provides valuable insight into this complex biological process.

SummaryMany individual proteins have been identified as having defined positions relative to cell polarity axes, raising the question of what fraction of all proteins may have polarized localizations. We took advantage of the giant ciliate Stentor coeruleus to quantify the extent of polarized localization proteome-wide. This trumpet-shaped unicellular organism shows a clear morphological anterior-posterior axis defined by a circular array of cilia known as a membranellar band at one end, and a holdfast at the other end. Because individual Stentor cells are over a millimeter in length, we were able to cut the cells into three pieces along the anterior-posterior axis, followed by proteomic analysis of proteins enriched in each piece. We find that approximately 30% of all detected proteins show a polarized location relative to the anterior-posterior cell axis. Proteins with polarized enrichment include centrin-like proteins, calcium-regulated kinases, orthologs of SFI1 and GAS2, and proteases. At the organelle level, nuclear and mitochondrial proteins are enriched in the anterior half of the cell body, but not in the membranellar band itself, while ribosome related proteins are apparently uniformly distributed. RNAi of signaling proteins enriched in the membranellar band, which is the anterior-most structure in the cell, revealed a protein phosphatase 2 subunit b ortholog required for closure of the membranellar band into the ring shape characteristic of Stentor. These results suggest that a large fraction of the Stentor proteome has a polarized localization, and provide a protein-level framework for future analysis of pattern formation and regeneration in Stentor as well as defining a general strategy for subcellular spatial proteomics based on physical dissection of cells.

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 alga Chlamydomonas provide 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.

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.

We use computational modeling to study within-host viral infection and evolution. In our model, viruses exhibit variable binding to cells, with better infection and replication countered by a stronger immune response and a high rate of mutation. By varying host conditions (permissivity to viral entry T and immune clearance intensity A) for large numbers of cells and viruses, we study the dynamics of how viral populations evolve from initial infection to steady state and obtain a phase diagram of the range of cell and viral responses. We find three distinct replicative strategies corresponding to three physiological classes of viral infections: acute, chronic, and opportunistic. We show similarities between our findings and the behavior of real viral infections such as common flu, hepatitis, and SARS-CoV-2019. The phases associated with the three strategies are separated by a phase transition of primarily first order, in addition to a crossover region. Our simulations also reveal a wide range of physical phenomena, including metastable states, periodicity, and glassy dynamics. Lastly, our results suggest that the resolution of acute viral disease in patients whose immunity cannot be boosted can only be achieved by significant inhibition of viral infection and replication.Author summaryVirus, in particular RNA viruses, often produce offspring with slightly altered genetic composition. This process occurs both across host populations and within a single host over time. Here, we study the interactions of viruses with cells inside a host over time. In our model, the viruses encounter host cell defenses characterized by two parameters: permissivity to viral entry T and immune response A). The viruses then mutate upon reproduction, eventually resulting in a distribution of related viral types termed a quasi-species distribution. Across varying host conditions (T, A), three distinct viral quasi-species types emerge over time, corresponding to three classes of viral infections: acute, chronic and opportunistic. We interpret these results in terms of real viral types such as common flu, hepatitis, and also SARS-CoV-2019. Analysis of viral of viral mutant populations over a wide range of permissivity and immunity, for large numbers of cells and viruses, reveals phase transitions that separate the three classes of viruses, both in the infection-cycle dynamics and at steady state. We believe that such a multiscale approach for the study of within-host viral infections, spanning individual proteins to collections of cells, can provide insight into developing more effective therapies for viral disease.

SummaryThe giant ciliate,Stentor coeruleus, provides a unique opportunity to study nuclear shape because its macronucleus undergoes a rapid, dramatic, and developmentally regulated shape change. During a 2 hour time period within cell division and regeneration, the 400 um long moniliform macronucleus coalesces into a single mass, elongates into a vermiform shape, and then renodulates, returning to its original beads-on-a-string morphology.1Previous work from the 1960’s - 1980’s demonstrated that the macronuclear shape change is a highly regulated part of cell division and regeneration2,3, but the molecular pathways driving these changes are unknown. With the recent availability of a sequencedStentorgenome, a transcriptome during regeneration, and molecular tools like RNAi4–6, it is now possible to investigate the molecular mechanisms that drive macronuclear 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 is degraded during the period when the macronucleus returns to its pre-coalescence volume. We propose a model in which nuclear transport via CSE1 increases the volume of the macronucleus, thereby decreasing the surface to volume ratio and driving coalescence of the many nodes into a single mass.

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.

The giant ciliate Stentor coeruleus is a classical model system for studying regeneration and morphogenesis at the level of 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 will regenerate an intact cell, including a new oral apparatus in the posterior half. We used 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 implicated in Meckel and Joubert syndromes. By comparing transcriptional profiles of different regeneration events in the same species, we were able to identify distinct modules of gene expression corresponding to oral apparatus regeneration, posterior holdfast regeneration, and recovery after wounding. By measuring gene expression in cells in which translation is blocked, 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 conserved 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. We show that 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.

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. Contrastingly, 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.

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.

Stentor coeruleus has long been employed as a model organism to understand unicellular regeneration. Transcriptomic analysis of individual cells revealed hundreds of genes -- many not associated with the oral apparatus -- differentially regulated in phases throughout the regeneration process. We hypothesized that this systemic reorganization and mobilization of cellular resources towards growth of a new oral apparatus will lead to observable changes in movement and behavior corresponding in time to the phases of differential gene expression, thus allowing an assay of functional regeneration based on motility analysis. However, the morphological complexity of Stentor coeruleus meant the development of an assay to capture the statistics and timescale was necessary. Using a custom script to track cells in short videos, we compile statistics over a large population (N(~)100) to show that, upon loss of the OA, Stentor coeruleus initially loses the ability for directed motion, then starting at approximately four hours exhibits a significant drop in speed until approximately eight hours. This assay provides a useful tool for the screening of motility phenotypes, and can be adapted for the investigation of other organisms.

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.

BackgroundWound 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.ResultsHere we used a microfluidic guillotine and a fluorescence-based assay to probe the timescales of wound repair and of mechanical modes of wound response inStentor. We found thatStentorrequires ~ 100–1000 s to close bisection wounds, depending on the severity of the wound. This corresponds to a healing rate of ~ 8–80 μm2/s, faster than most other single cells reported in the literature. Further, we characterized three distinct mechanical modes of wound response inStentor: 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 inStentor, was surprisingly not important for the contraction mode and was of low importance for the others.ConclusionsWhile 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 inStentorand will lay the foundation for further investigations into the underlying components and molecular mechanisms involved.

Summary 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., Graphical abstract, Highlights • Chlamydomonas flagella show dynamic length fluctuations • Length variation involves processes inside the cell body and in the flagella • Long-flagella mutants have increased length variability • Flagellar length variation leads to variation in cell motility, Biological Sciences; Molecular Biology; Cell Biology; Biophysics; Mathematical Biosciences

Modern approaches in quantitative live cell imaging have become an essential tool for exploring cell biology, by enabling the use of statistics and computational modeling to classify and compare biological processes. Although cell culture model systems are great for high content imaging, high throughput studies of cell morphology suggest that ex vivo cultures are limited in recapitulating the morphological complexity found in cells within living organisms. As such, there is a need for a scalable high throughput model system to image living cells within an intact organism. Described here is a protocol for using a high content image analyzer to simultaneously acquire multiple time-lapse videos of embryonic Drosophila melanogaster development during the syncytial blastoderm stage. The syncytial blastoderm has traditionally served as a great in vivo model for imaging biological events; however, obtaining a significant number of experimental replicates for quantitative and high-throughput approaches has been labor intensive and limited by the imaging of a single embryo per experimental repeat. Presented here is a method to adapt imaging and microinjection approaches to suit a high content imaging system, or any inverted microscope capable of automated multipoint acquisition. This approach enables the simultaneous acquisition of 6-12 embryos, depending on desired acquisition factors, within a single imaging session.

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.

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 fractionation7-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 yeast12-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 2mm 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.

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.

Analysis of fluctuation in organelle size provides a new way to probe the mechanisms of organelle size control systems. By analyzing cell-to-cell variation and within-cell fluctuations of flagellar length inChlamydomonas, we show that the flagellar length control system exhibits both types of variation. Cell to cell variation is dominated by cell size, while within-cell variation results from dynamic fluctuations that are subject to a constraint, providing evidence for a homeostatic size control system. We analyzed a series of candidate genes affecting flagella and found that flagellar length variation is increased in mutations which increase the average flagellar length, an effect that we show is consistent with a theoretical model for flagellar length regulation based on length-dependent intraflagellar transport balanced by length-independent disassembly. Comparing the magnitude and time-scale of length fluctuations with simple models suggests that tubulin assembly is not directly coupled with IFT-mediated arrival and that observed fluctuations involve tubulin assembly and disassembly events involving large numbers of tubulin dimers. Cells with greater differences in their flagellar lengths show impaired swimming but improved gliding motility, raising the possibility that cells have evolved mechanisms to tune intrinsic noise in length. Taken together our results show that biological noise exists at the level of subcellular structures, with a corresponding effect on cell function, and can provide new insights into the mechanisms of organelle size control.

In the context of animal or plant development, we tend to think of cells as small, simple, building blocks, such that complex patterns or shapes can only be constructed from large numbers of cells, with cells in different parts of the organism taking on different fates. However, cells themselves are far from simple, and often take on complex shapes with a remarkable degree of intracellular patterning. How do these patterns arise? As in embryogenesis, the development of structure inside a cell can be broken down into a number of basic processes. For each part of the cell, morphogenetic processes create internal structures such as organelles, which might correspond to organs at the level of a whole organism. Given that mechanisms exist to generate parts, patterning processes are required to ensure that the parts are distributed in the correct arrangement relative to the rest of the cell. Such patterning processes make reference to global polarity axes, requiring mechanisms for axiation which, in turn, require processes to break symmetry. These fundamental processes of symmetry breaking, axiation, patterning, and morphogenesis, have been extensively studied in developmental biology but less so at the sub-cellular level. This review will focus on developmental processes that give eukaryotic cells their complex structures, with a focus on cytoskeletal organization in free-living cells, ciliates in particular, in which these processes are most readily apparent.

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 asChlamydomonas.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 ciliateStentor 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 inStentor, 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’.

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.

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., Highlighted Article: Aged muscle stem cells display delayed activation dynamics, but retain a youthful activation trajectory, suggesting that changes to cell state dynamics might contribute to aging pathology.

Cells are complex biochemical systems whose behavior emerges from interactions among myriad molecular components. Computation is often invoked as a general framework for navigating this cellular complexity. However, the manner in which cells might embody computational processes such that theories of computation, including finite state machine models, could be productively applied, remains to be seen. 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 highly regulated transitions between a discrete set of gait states. The set of observed transitions decomposes into a small group of high-probability, temporally irreversible transitions forming a cycle and a large group of low-probability time-symmetric transitions, thus revealing stereotypy in sequential patterns of state transitions. Taken together, these findings implicate a finite state machine-like process. Cirri are connected by microtubule bundles, and we found that the dynamics of cirri involved in different state transitions are associated with the structure of the microtubule bundle system. Perturbative experiments revealed that the fibers mediate gait coordination, suggesting a mechanical basis of gait control. One sentence summary The ciliate Euplotes walks across surfaces with an unusual gait, which involves modular subcellular structure and combines elements of stereotypy and variability according to a computational process actively coordinated by a microtubule fiber system.

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.

Summary 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.

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.