A classic problem in microbiology is that bacteria display two types of growth behavior when cultured on a mixture of two carbon sources: the two sources are sequentially consumed one after another (diauxie) or they are simultaneously consumed (co-utilization). The search for the molecular mechanism of diauxie led to the discovery of the lac operon. However, questions remain as why microbes would bother to have different strategies of taking up nutrients. Here we show that diauxie versus co-utilization can be understood from the topological features of the metabolic network. A model of optimal allocation of protein resources quantitatively explains why and how the cell makes the choice. In case of co-utilization, the model predicts the percentage of each carbon source in supplying the amino acid pools, which is quantitatively verified by experiments. Our work solves a long-standing puzzle and provides a quantitative framework for the carbon source utilization of microbes.

While most processes in biology are highly deterministic, stochastic mechanisms are sometimes used to increase cellular diversity, such as in the specification of sensory receptors. In the human and Drosophila eye, photoreceptors sensitive to various wavelengths of light are distributed randomly across the retina. Mechanisms that underlie stochastic cell fate specification have been analysed in detail in the Drosophila retina. In contrast, the retinas of another group of dipteran flies exhibit highly ordered patterns. Species in the Dolichopodidae, the "long-legged" flies, have regular alternating columns of two types of ommatidia (unit eyes), each producing corneal lenses of different colours. Individual flies sometimes exhibit perturbations of this orderly pattern, with "mistakes" producing changes in pattern that can propagate across the entire eye, suggesting that the underlying developmental mechanisms follow local, cellular-automaton-like rules. We hypothesize that the regulatory circuitry patterning the eye is largely conserved among flies such that the difference between the Drosophila and Dolichopodidae eyes should be explicable in terms of relative interaction strengths, rather than requiring a rewiring of the regulatory network. We present a simple stochastic model which, among its other predictions, is capable of explaining both the random Drosophila eye and the ordered, striped pattern of Dolichopodidae.

We investigate the geometrical and mechanical properties of adherent cells characterized by a highly anisotropic actin cytoskeleton. Using a combination of theoretical work and experiments on micropillar arrays, we demonstrate that the shape of the cell edge is accurately described by elliptical arcs, whose eccentricity expresses the degree of anisotropy of the internal cell stresses. This results in a spatially varying tension along the cell edge, that significantly affects the traction forces exerted by the cell on the substrate. Our work highlights the strong interplay between cell mechanics and geometry and paves the way towards the reconstruction of cellular forces from geometrical data.

Recent experimental studies have suggested the ratio between T-helper and T-suppressor lymphocytes as an index of immunosuppression in HIV, cancer, immunosenescence and inflammatory and auto-immune diseases. However, a quantitative understanding of the impact of this ratio on the immune response has lagged behind data and its validity as a tool for prognostic monitoring or therapeutic target remains an open question. In this work, we use statistical physics and dynamical systems approaches to analyze the time-dependent response to an antigen, of a simplified model of the adaptive immune system, which comprises B, T-helper and T-suppressor lymphocytes. The model is remarkably robust against changes in the noise level and kinetic parameters, but it is very sensitive to changes in the ratio between T-helper and T-suppressor lymphocytes, exhibiting, in particular, a transition from a responsive to an immuno-suppressed phase, as the ratio is lowered below a critical value, which is in line with experiments. This result supports the validity of the T-helper/T-suppressor ratio as an index of immunosuppression and may provide a useful theoretical benchmark to interpret and compare experiments.

Gene regulatory network (GRN)-based morphogenetic models have recently gained an increasing attention. However, the relationship between microscopic properties of intracellular GRNs and macroscopic properties of morphogenetic systems has not been fully understood yet. Here we propose a theoretical morphogenetic model representing an aggregation of cells, and reveal the relationship between criticality of GRNs and morphogenetic pattern formation. In our model, the positions of the cells are determined by spring-mass-damper kinetics. Each cell has an identical Kauffman's $NK$ random Boolean network (RBN) as its GRN. We varied the properties of GRNs from ordered, through critical, to chaotic by adjusting node in-degree $K$. We randomly assigned four cell fates to the attractors of RBNs for cellular behaviors. By comparing diverse morphologies generated in our morphogenetic systems, we investigated what the role of the criticality of GRNs is in forming morphologies. We found that nontrivial spatial patterns were generated most frequently when GRNs were at criticality. Our finding indicates that the criticality of GRNs facilitates the formation of nontrivial morphologies in GRN-based morphogenetic systems.

Receptor-mediated endocytosis is an ubiquitous process through which cells internalize biological or synthetic nanoscale objects, including viruses, unicellular parasites, and nanomedical vectors for drug or gene delivery. In passive endocytosis the cell plasma membrane wraps around the "invader" particle driven by ligand-receptor complexation. By means of theory and numerical simulations, here we demonstrate how particles decorated by freely diffusing and non-mutually-interacting (ideal) ligands are significantly more difficult to wrap than those where ligands are either immobile or interact sterically with each other. Our model rationalizes the relationship between uptake mechanism and structural details of the invader, such as ligand size, mobility and ligand/receptor affinity, providing a comprehensive picture of pathogen endocytosis and helping the rational design of efficient drug delivery vectors.

In this paper we propose and study a hybrid discrete in continuous mathematical model of collective motion under alignment and chemotaxis effect. Starting from the paper by Di Costanzo et al (2015a), in which the Cucker-Smale model (Cucker and Smale, 2007) was coupled with other cell mechanisms, to describe the cell migration and self-organization in the zebrafish lateral line primordium, we introduce a simplified model in which the coupling between an alignment and chemotaxis mechanism acts on a system of interacting particles. In particular we rely on a hybrid description in which the agents are discrete entities, while the chemoattractant is considered as a continuous signal. The proposed model is then studied both from an analytical and a numerical point of view. From the analytic point of view we prove, globally in time, existence and uniqueness of the solution. Then, the asymptotic behaviour of a linearised version of the system is investigated. Through a suitable Lyapunov functional we show that for $t\rightarrow +\infty$, the migrating aggregate exponentially converges to a state in which all the particles have a same position with zero velocity. Finally, we present a comparison between the analytical findings and some numerical results, concerning the behaviour of the full nonlinear system.

Multiple myeloma (MM), a plasma cell cancer, is associated with many health challenges, including damage to the kidney by tubulointerstitial fibrosis. We develop a mathematical model which captures the qualitative behavior of the cell and protein populations involved. Specifically, we model the interaction between cells in the proximal tubule of the kidney, free light chains, renal fibroblasts, and myeloma cells. We analyze the model for steady-state solutions to find a mathematically and biologically relevant stable steady-state solution. This foundational model provides a representation of dynamics between key populations in tubulointerstitial fibrosis that demonstrates how these populations interact to affect patient prognosis in patients with MM and renal impairment.

Laboratory models are often used to understand the interaction of related pathogens via host immunity. For example, recent experiments where ferrets were exposed to two influenza strains within a short period of time have shown how the effects of cross-immunity vary with the time between exposures and the specific strains used. On the other hand, studies of the workings of different arms of the immune response, and their relative importance, typically use experiments involving a single infection. However, inferring the relative importance of different immune components from this type of data is challenging. Using simulations and mathematical modelling, here we investigate whether the sequential infection experiment design can be used not only to determine immune components contributing to cross-protection, but also to gain insight into the immune response during a single infection. We show that virological data from sequential infection experiments can be used to accurately extract the timing and extent of cross-protection. Moreover, the broad immune components responsible for such cross-protection can be determined. Such data can also be used to infer the timing and strength of some immune components in controlling a primary infection, even in the absence of serological data. By contrast, single infection data cannot be used to reliably recover this information. Hence, sequential infection data enhances our understanding of the mechanisms underlying the control and resolution of infection, and generates new insight into how previous exposure influences the time course of a subsequent infection.

Recent research has shown that motile cells can adapt their mode of propulsion to the mechanical properties of the environment in which they find themselves--crawling in some environments while swimming in others. The latter can involve movement by blebbing or other cyclic shape changes, and both highlysimplified and more realistic models of these modes have been studied previously. Herein we study swimming that is driven by membrane tension gradients that arise from flows in the actin cortex underlying the membrane, and does not involve imposed cyclic shape changes. Such gradients can lead to a number of different characteristic cell shapes, and our first objective is to understand how different distributions of membrane tension influence the shape of cells in an inviscid quiescent fluid. We then analyze the effects of spatial variation in other membrane properties, and how they interact with tension gradients to determine the shape. We also study the effect of fluid--cell interactions and show how tension leads to cell movement, how the balance between tension gradients and a variable bending modulus determine the shape and direction of movement, and how the efficiency of movement depends on the properties of the fluid and the distribution of tension and bending modulus in the membrane.

Recent studies show the existing clinical tests to detect Cardio/cerebrovascular diseases (CVD) are ineffectual as they do not consider different stages of platelet activation or the molecular dynamics involved in platelet interactions. Further they are also incapable to consider inter-individual variability. A physical description of platelets deposition was introduced recently in Chopard et. al. [2017], by integrating fundamental understandings of how platelets interact in a numerical model, parameterized by five parameters. These parameters specify the deposition process and are relevant for a biomedical understanding of the phenomena. One of the main intuition is that these parameters are precisely the information needed for a pathological test identifying CVD captured and that they capture the inter-individual variability. Following this intuition, here we devise a Bayesian inferential scheme for estimation of these parameters. As the likelihood function of the numerical model is intractable due to the complex stochastic nature of the model, we use a likelihood-free inference scheme approximate Bayesian computation (ABC) to calibrate the parameters in a data-driven manner. As ABC requires the generation of many pseudo-data by expensive simulation runs, we use a high performance computing (HPC) framework for ABC to make the inference possible for this model. We illustrate that our mean posterior prediction of platelet deposition pattern matches the experimental dataset closely with a tight posterior prediction error margin for a collective dataset of 7 volunteers. The present approach can be used to build a new generation of personalized platelet functionality tests for CVD detection, using numerical modeling of platelet deposition, Bayesian uncertainty quantification and High performance computing.

Cellular phenotype is characterized by different components such as cell size, protein content and cell cycle time. These are global variables that are the outcome of multiple internal microscopic processes. Accordingly, they display some universal statistical properties and scaling relations, such as distribution collapse and relation between moments. Cell size statistics and its relation to growth and division has been mostly studied separately from proteins and other cellular variables. Here we present experimental and theoretical analyses of these phenotype components in a unified framework that reveals their correlations and interactions inside the cell. We measure these components simultaneously in single cells over dozens of generations, quantify their correlations, and compare population to temporal statistics. We find that cell size and highly expressed proteins have very similar dynamics over growth and division cycles, which result in parallel statistical properties, both universal and individual. In particular, while distribution shapes of fluctuations along time are common to all cells and components, other properties are variable and remain distinct in individual cells for a surprisingly large number of generations. These include temporal averages of cell size and protein content, and the structure of their auto-correlation functions. We explore possible roles of the different components in controlling cell growth and division. We find that in order to stabilize exponential accumulation and division of all components across generations, coupled dynamics among them is required. Finally, we incorporate effective coupling within the cell cycle with a phenomenological mapping across consecutive cycles, and show that this model reproduces the entire array of experimental observations.

To estimate the time, many organisms, ranging from cyanobacteria to animals, employ a circadian clock which is based on a limit-cycle oscillator that can tick autonomously with a nearly 24h period. Yet, a limit-cycle oscillator is not essential for knowing the time, as exemplified by bacteria that possess an 'hourglass': a system that when forced by an oscillatory light input exhibits robust oscillations from which the organism can infer the time, but that in the absence of driving relaxes to a stable fixed point. Here, using models of the Kai system of cyanobacteria, we compare a limit- cycle oscillator with two hourglass models, one that without driving relaxes exponentially and one that does so in an oscillatory fashion. In the limit of low input-noise, all three systems are equally informative on time, yet in the regime of high input-noise the limit-cycle oscillator is far superior. The same behavior is found in the Stuart-Landau model, indicating that our result is universal.

Only about 8 percent of patients with pancreatic ductal adenocarcinoma, PDAC, live up to 5 years following diagnosis; by 2020 PDAC will become the second most lethal cancer in the United States. PDAC includes an anisotropic fibrous-like stroma, desmoplasia, encompassing most of the tumor mass. Desmoplasia is produced by cancer-associated fibroblasts, CAFs, and their cell-derived extracellular matrices, CDMs. Since elimination of CAFs is detrimental to patients, CDM reprogramming, as opposed to desmoplasia ablation, is therapeutically desirable. In this study we used a human mimetic three-dimensional CAF producing CDM system and proceeded to study its dynamic architectural modifications, following underlying substrate stiffness alterations, using digital imaging analyses, atomic force microscopy, mathematical modeling, cell biology, biochemistry and human tissue quantitative simultaneous multiplex immunofluorescence. Results suggested that the architecture of CDMs can be manipulated to render a tumor-suppressive microenvironment. We posit that perhaps treatments that could reprogram desmoplasia to become tumor-restrictive or that could target tumoral ERK2 might provide future new means for treating PDAC patients.

Piscine orthoreovirus Strain PRV-1 is the causative agent of heart and skeletal muscle inflammation (HSMI) in Atlantic salmon (Salmo salar). Given its high prevalence in net pen salmon, debate has arisen on whether PRV poses a risk to migratory salmon, especially in British Columbia (BC) where commercially important wild Pacific salmon are in decline. Various strains of PRV have been associated with diseases in Pacific salmon, including erythrocytic inclusion body syndrome (EIBS), HSMI-like disease, and jaundice/anemia in Japan, Norway, Chile and Canada. We examine the developmental pathway of HSMI and jaundice/anemia associated with PRV-1 in farmed Atlantic and Chinook (Oncorhynchus tshawytscha) salmon in BC, respectively. In situ hybridization localized PRV-1 within developing lesions in both diseases. The two diseases showed dissimilar pathological pathways, with inflammatory lesions in heart and skeletal muscle in Atlantic salmon, and degenerative-necrotic lesions in kidney and liver in Chinook salmon, plausibly explained by differences in PRV load tolerance in red blood cells. Viral genome sequencing revealed no consistent differences in PRV-1 variants intimately involved in the development of both diseases, suggesting that migratory Chinook salmon may be at more than a minimal risk of disease from exposure to the high levels of PRV occurring on salmon farms.

Acute lymphoblastic leukemia (ALL) and chronic lymphocytic leukemia (CLL) are two major forms of leukemia that arise from lymphoid cells (LCs). ALL occurs mostly in children and CLL occurs mainly in old people. However, the Philadelphia-chromosome-positive ALL (Ph+-ALL) and the Ph-like ALL occur in both children and adults. To understand childhood leukemia/lymphoma, we have recently proposed two hypotheses on the causes and the mechanism of cell transformation of a LC. Hypothesis A is: repeated bone-remodeling during bone-growth and bone-repair may be a source of cell injuries of marrow cells including hematopoietic stem cells (HSCs), myeloid cells, and LCs. Hypothesis B is: a LC may have three pathways on transformation: a slow, a rapid, and an accelerated. We discuss in the present paper the developing mechanisms of ALL and CLL by these hypotheses. Having a peak incidence in young children, ALL may develop mainly as a result of rapid cell transformation of a lymphoblast (or pro-lymphocyte). Differently, Ph+-ALL and Ph-like ALL may develop as results of transformation of a lymphoblast via accelerated pathway. Occurring mainly in adults, CLL may be a result of transformation of a memory B-cell via slow pathway. By causing cell injuries of HSCs and LCs, repeated bone-remodeling during bone-growth and bone-repair may be related to the cell transformation of a LC. In conclusion, ALL may develop as a result of cell transformation of a lymphoblast via rapid or accelerated pathway; and repeated bone-remodeling during bone-growth may be a trigger for the cell transformation of a lymphoblast in a child.

Lymphomas are a large group of neoplasms developed from lymphoid cells (LCs) in lymph nodes (LNs) or lymphoid tissues (LTs). Some forms of lymphomas, including Burkitt lymphoma (BL), ALK+ anaplastic large cell lymphoma (ALK+-ALCL), and T-cell lymphoblastic lymphoma/leukemia (T-LBL), occur mainly in children and teenagers. Hodgkin's lymphoma (HL) has a peak incidence at age 20s. To understand pediatric lymphoma, we have recently proposed two hypotheses on the causes and the mechanism of cell transformation of a LC. Hypothesis A is: repeated bone-remodeling during bone-growth and bone-repair may be a source of cell injuries of marrow cells including hematopoietic stem cells (HSCs), myeloid cells, and LCs, and thymic involution may be a source of damage to the developing T-cells in thymus. Hypothesis B is: a LC may have three pathways on transformation: a slow, a rapid, and an accelerated. In this paper, we discuss pediatric lymphomas by this hypothesis. Having a peak incidence at young age, BL, T-LBL, ALK+-ALCL, and HL develop more likely as a result of rapid transformation of a LC. In BL, ALK+-ALCL, and HL, the cell transformations may be triggered by severe viral infections. In T-LBL, the cell transformation may be related to thymic involution. Occurring in both adults and children, diffuse large B-cell lymphoma (DLBCL) may develop via slow or accelerated pathway. In conclusion, pediatric lymphoma may develop as a result of "one-step" cell transformation of a LC, and severe viral infections may be the main trigger for the rapid transformation of a LC in a LN/LT.

Acute myeloid leukemia (AML) and chronic myeloid leukemia (CML) are two major formsof leukemia developed from myeloid cells (MCs). To understand why AML and CML occurin children, we analyzed the causes and the mechanism of cell transformation of a MC. I. Forthe MCs in marrow cavity, repeated bone-remodeling during bone-growth may be a source ofcell injuries. II. As a type of blood cell, a MC may have higher survivability from DNAchanges and require obtaining fewer cancerous properties for cell transformation than a tissuecell. III. Point DNA mutations (PDMs) and chromosome changes (CCs) are the two majortypes of DNA changes. CCs have three subtypes by effects on a cell: great effect CCs(GECCs), mild-effect CCs (MECCs), and intermediate-effect CCs (IECCs). A GECC affectsone or more genes and can alone trigger cell transformation. PDMs/MECCs are mostly mildand can accumulate in cells. Some of the PDMs/MECCs contribute to cell transformation. AnIECC affects one or more genes and participates in cell transformation. IV. Based on II andIII, we hypothesize that a MC may have two pathways on transformation: a slow and anaccelerated. Slow pathway is driven by accumulation of PDMs/MECCs. Accelerated pathwayis driven by accumulation of PDMs/MECCs/IECC(s). A transformation via slow pathwayoccurs at old age; whereas that via accelerated pathway occurs at any age. Thus, CML andpediatric AML may develop via accelerated pathway, and adult AML may develop via bothpathways. In conclusion, pediatric AML and CML may develop as a result of transformationof a MC via accelerated pathway; and repeated bone-remodeling for bone-growth may be atrigger for the transformation of a MC in a child.

One of the many ways cells transmit information within their volume is through steady spatial gradients of different proteins. However, the mechanism through which proteins without any sources or sinks form such single-cell gradients is not yet fully understood. One of the models for such gradient formation, based on differential diffusion, is limited to proteins with large ratios of their diffusion constants or to specific protein-large molecule interactions. We introduce a novel mechanism for gradient formation via the coupling of the proteins within a single cell with a molecule, that we call a "pronogen", whose action is similar to that of morphogens in multi-cell assemblies, the pronogen is produced with a fixed flux at one side of the cell. This coupling results in an effectively non-linear diffusion degradation model for the pronogen dynamics within the cell, which leads to a steady-state gradient of the protein concentration. We use a stability analysis to show that these gradients are linearly stable with respect to perturbations.

Lymphoid leukemia (LL) and lymphoma are neoplasms developed from lymphoid cells (LCs). To understand why different forms of LL/lymphoma occur at different ages, we analyzed the effects of different types of DNA changes on a LC and the cellular characteristics of LCs. Point DNA mutations (PDMs) and chromosome changes (CCs) are the two major types of DNA changes. CCs have three subtypes by their effects on a LC: great-effect CCs (GECCs), mild-effect CCs (MECCs), and intermediate-effect CCs (IECCs). PDMs and MECCs are mostly mild thus can accumulate in cells. Some of the PDMs/MECCs contribute to cell transformation. A GECC affects one or more genes and can alone drive cell transformation. An IECC affects one or more genes and participates in cell transformation. Due to cellular characteristics, a LC may have higher survivability from DNA changes and require obtaining fewer cancerous properties for transformation than a tissue cell. Hence, a LC can be more rapidly transformed by a CC. On this basis, we hypothesize that a LC may have three pathways on transformation: a slow, a rapid, and an accelerated. Slow pathway is driven by accumulation of PDMs/MECCs. Rapid pathway is driven by a GECC in "one step". Accelerated pathway is driven by accumulation of PDMs/MECCs/IECC(s). Cell transformations of a LC via different pathways occur at different ages. A transformation via slow pathway occurs mainly in adults. A transformation via rapid pathway occurs at any age and has no increasing incidence with age. A transformation via accelerated pathway occurs also at any age but has increasing incidence with age. In conclusion, a LC may have three pathways on cell transformation, and the occurring age of LL/lymphoma may be determined by the transforming pathway of a LC.

Gaining access to the cell interior is fundamental for many applications, such as electrical recording, drug and biomolecular delivery. A very promising technique consists of culturing cells on nano/micro pillars. The tight adhesion and high local deformation of cells in contact with nanostructures can promote the permeabilization of lipids at the plasma membrane, providing access to the internal compartment. However, there is still much experimental controversy regarding when and how the intracellular environment is targeted and the role of the geometry and interactions with surfaces. Consequently, we investigated, by coarse-grained molecular dynamics simulations of the cell membrane, the mechanical properties of the lipid bilayer under high strain and bending conditions. We found out that a high curvature of the lipid bilayer dramatically lowers the traction force necessary to achieve membrane rupture. Afterwards, we experimentally studied the permeabilization rate of cell membrane by pillars with comparable aspect ratios but different sharpness values at the edges. The experimental data support the simulation results: even pillars with diameters in the micron range may cause local membrane disruption when their edges are sufficiently sharp. Therefore, the permeabilization likelihood is connected to the local geometric features of the pillars rather than diameter or aspect ratio. The present study can also provide significant contributions to the design of 3D biointerfaces for tissue engineering and cellular growth.

Lymphoid leukemia (LL) and lymphoma are blood cancers developed from lymphoid cells (LCs). To understand the cause and the mechanism of cell transformation of a LC, we studied the potential sources of cell injuries of LCs and analyzed how DNA changes are generated and accumulate in LCs. I. The DNA changes that contribute to cell transformation of a LC can be generated in the LCs in marrow, thymus, lymph nodes (LNs), and/or lymphoid tissues (LTs). In LNs/LTs, pathogen-infections may be the main cause for cell injuries of LCs. In marrow cavity, repeated bone-remodeling during bone-growth and bone-repair, by producing toxic substances, may be a source of damage to hematopoietic cells, including hematopoietic stem cells (HSCs) and developing LCs. In thymus, thymic involution and death of stromal cells may be a damaging factor for the developing T-cells. II. Point DNA mutation (PDM) and chromosome change (CC) are the two major types of DNA changes. CCs include numerical CCs (NCCs) and structural CCs (SCCs). Generation of a PDM/SCC may be a result of Misrepair of DNA on DNA breaks. Generation of NCC is rather a consequence of dysfunction of cell division. III. Repeated cell injuries and cell proliferation drive the accumulation of DNA changes in LCs and HSCs. However, long-term accumulation of DNA changes occurs mainly in long-living stem cells including HSCs and memory cells. In conclusion, the DNA changes in LCs are generated and accumulate as a consequence of repeated cell injuries and repeated cell proliferation; and three potential sources of cell injuries of LCs may be: repeated bone-remodeling, long-term thymic involution, and repeated pathogen-infections.

We present a theory of pattern formation in growing domains inspired by biological examples of tissue development. Gradients of signaling molecules regulate growth, while growth changes these graded chemical patterns by dilution and advection. We identify a critical point of this feedback dynamics, which is characterized by spatially homogeneous growth and proportional scaling of patterns with tissue length. We apply this theory to the biological model system of the developing wing of the fruit fly \textit{Drosophila melanogaster} and quantitatively identify signatures of the critical point.

Asymmetric segregation of key proteins at cell division -- be it a beneficial or deleterious protein -- is ubiquitous in unicellular organisms and often considered as an evolved trait to increase fitness in a stressed environment. Here, we provide a general framework to describe the evolutionary origin of this asymmetric segregation. We compute the population fitness as a function of the protein segregation asymmetry $a$, and show that the value of $a$ which optimizes the population growth manifests a phase transition between symmetric and asymmetric partitioning phases. Surprisingly, the nature of phase transition is different for the case of beneficial proteins as opposed to proteins which decrease the single-cell growth rate. Our study elucidates the optimization problem faced by evolution in the context of protein segregation, and motivates further investigation of asymmetric protein segregation in biological systems.

Despite their vast morphological diversity, many invertebrates have similar larval forms characterized by ciliary bands, innervated arrays of beating cilia that facilitate swimming and feeding. Hydrodynamics suggests that these bands should tightly constrain the behavioral strategies available to the larvae; however, their apparent ubiquity suggests that these bands also confer substantial adaptive advantages. Here, we use hydrodynamic techniques to investigate "blinking," an unusual behavioral phenomenon observed in many invertebrate larvae in which ciliary bands across the body rapidly change beating direction and produce transient rearrangement of the local flow field. Using a general theoretical model combined with quantitative experiments on starfish larvae, we find that the natural rhythm of larval blinking is hydrodynamically optimal for inducing strong mixing of the local fluid environment due to transient streamline crossing, thereby maximizing the larvae's overall feeding rate. Our results are consistent with previous hypotheses that filter feeding organisms may use chaotic mixing dynamics to overcome circulation constraints in viscous environments, and it suggests physical underpinnings for complex neurally-driven behaviors in early-divergent animals.