Cell movement has essential functions in development, immunity, and cancer. Various cell migration patterns have been reported, but no general rule has emerged so far. Here, we show on the basis of experimental data in vitro and in vivo that cell persistence, which quantifies the straightness of trajectories, is robustly coupled to cell migration speed. We suggest that this universal coupling constitutes a generic law of cell migration, which originates in the advection of polarity cues by an actin cytoskeleton undergoing flows at the cellular scale. Our analysis relies on a theoretical model that we validate by measuring the persistence of cells upon modulation of actin flow speeds and upon optogenetic manipulation of the binding of an actin regulator to actin filaments. Beyond the quantitative prediction of the coupling, the model yields a generic phase diagram of cellular trajectories, which recapitulates the full range of observed migration patterns.

The distribution of exit times is computed for a Brownian particle in spherically symmetric twodimensional domains (disks, angular sectors, annuli) and in rectangles that contain an exit on their boundary. The governing partial differential equation of Helmholtz type with mixed Dirichlet-Neumann boundary conditions is solved analytically. We propose both an exact solution relying on a matrix inversion, and an approximate explicit solution. The approximate solution is shown to be exact for an exit of vanishing size and to be accurate even for large exits. For angular sectors, we also derive exact explicit formulas for the moments of the exit time. For annuli and rectangles, the approximate expression of the mean exit time is shown to be very accurate even for large exits. The analysis is also extended to biased diffusion. Since the Helmholtz equation with mixed boundary conditions is encountered in microfluidics, heat propagation, quantum billiards, and acoustics, the developed method can find numerous applications beyond exit processes.

We present an exact calculation of the mean first-passage time to a target on the surface of a 2D or 3D spherical domain, for a molecule alternating phases of surface diffusion on the domain boundary and phases of bulk diffusion. We generalize the results of Bénichou et al. in (J. Stat. Phys. 142:657, 2011) and consider a biased diffusion in a general annulus with an arbitrary number of regularly spaced targets on a partially reflecting surface. The presented approach is based on an integral equation which can be solved analytically. Numerically validated approximation schemes, which provide more tractable expressions of the mean first-passage time are also proposed. In the framework of this minimal model of surface-mediated reactions, we show analytically that the mean reaction time can be minimized as a function of the desorption rate from the surface.

During development, many mutations cause increased variation in phenotypic outcomes, a phenomenon termed decanalization. Phenotypic discordance is often observed in the absence of genetic and environmental variations, but the mechanisms underlying such inter-individual phenotypic discordance remain elusive. Here, using the anterior-posterior (AP) patterning of the Drosophila embryo, we identified embryonic geometry as a key factor predetermining patterning outcomes under decanalizing mutations. With the wild-type AP patterning network, we found that AP patterning is robust to variations in embryonic geometry; segmentation gene expression remains reproducible even when the embryo aspect ratio is artificially reduced by more than twofold. In contrast, embryonic geometry is highly predictive of individual patterning defects under decanalized conditions of either increased bicoid (bcd) dosage or bcd knockout. We showed that the phenotypic discordance can be traced back to variations in the gap gene expression, which is rendered sensitive to the geometry of the embryo under mutations.

We introduce and investigate the escape problem for random walkers that may eventually die, decay, bleach, or lose activity during their diffusion towards an escape or reactive region on the boundary of a confining domain. In the case of a first-order kinetics (i.e., exponentially distributed lifetimes), we study the effect of the associated death rate onto the survival probability, the exit probability, and the mean first passage time. We derive the upper and lower bounds and some approximations for these quantities. We reveal three asymptotic regimes of small, intermediate and large death rates. General estimates and asymptotics are compared to several explicit solutions for simple domains, and to numerical simulations. These results allow one to account for stochastic photobleaching of fluorescent tracers in bio-imaging, degradation of mRNA molecules in genetic translation mechanisms, or high mortality rates of spermatozoa in the fertilization process. This is also a mathematical ground for optimizing storage containers and materials to reduce the risk of leakage of dangerous chemicals or nuclear wastes.

We present an exact expression for the mean exit time through the cap of a confining sphere for particles alternating phases of surface and of bulk diffusion. The present approach is based on an integral equation which can be solved analytically. In contrast to the statement of Berezhkovskii and Barzykin [J. Chem. Phys. 136, 54115 (2012)], we show that the mean exit time can be optimized with respect to the desorption rate, under analytically determined criteria.

In curved environments, geometric constraints can result in cellular competition for space. Quantitative information about cell geometry was extracted to explore how columnar-like cells fit within the curved anterior of the Drosophila embryo. Cell deformations and cell rearrangements enable packing in this highly curved region.

Organ formation is an inherently biophysical process, requiring large-scale tissue deformations. Yet, understanding how complex organ shape emerges during development remains a major challenge. During zebrafish embryogenesis, large muscle segments, called myotomes, acquire a characteristic chevron morphology, which is believed to aid swimming. Myotome shape can be altered by perturbing muscle cell differentiation or the interaction between myotomes and surrounding tissues during morphogenesis. To disentangle the mechanisms contributing to shape formation of the myotome, we combine single-cell resolution live imaging with quantitative image analysis and theoretical modeling. We find that, soon after segmentation from the presomitic mesoderm, the future myotome spreads across the underlying tissues. The mechanical coupling between the future myotome and the surrounding tissues appears to spatially vary, effectively resulting in spatially heterogeneous friction. Using a vertex model combined with experimental validation, we show that the interplay of tissue spreading and friction is sufficient to drive the initial phase of chevron shape formation. However, local anisotropic stresses, generated during muscle cell differentiation, are necessary to reach the acute angle of the chevron in wild-type embryos. Finally, tissue plasticity is required for formation and maintenance of the chevron shape, which is mediated by orientated cellular rearrangements. Our work sheds light on how a spatiotemporal sequence of local cellular events can have a nonlocal and irreversible mechanical impact at the tissue scale, leading to robust organ shaping.