Topological floppy modes in models of epithelial tissues.

Recent advances in topological mechanics have revealed unusual phenomena such as topologically protected floppy modes and states of self-stress that are exponentially localized at boundaries and interfaces of mechanical networks. In this paper, we explore the topological mechanics of epithelial tissues, where the appearance of these boundary and interface modes could lead to localized soft or stressed spots and play a role in morphogenesis. We consider both a simple vertex model (VM) governed by an effective elastic energy and its generalization to an active tension network (ATN) which incorporates active adaptation of the cytoskeleton. By analyzing spatially periodic lattices at the Maxwell point of mechanical instability, we find topologically polarized phases with exponential localization of floppy modes and states of self-stress in the ATN when cells are allowed to become concave, but not in the VM.

Defect patterns on the curved surface of fish retinae suggest a mechanism of cone mosaic formation.

The outer epithelial layer of zebrafish retinae contains a crystalline array of cone photoreceptors, called the cone mosaic. As this mosaic grows by mitotic addition of new photoreceptors at the rim of the hemispheric retina, topological defects, called "Y-Junctions", form to maintain approximately constant cell spacing. The generation of topological defects due to growth on a curved surface is a distinct feature of the cone mosaic not seen in other well-studied biological patterns like the R8 photoreceptor array in the Drosophila compound eye. Since defects can provide insight into cell-cell interactions responsible for pattern formation, here we characterize the arrangement of cones in individual Y-Junction cores as well as the spatial distribution of Y-junctions across entire retinae. We find that for individual Y-junctions, the distribution of cones near the core corresponds closely to structures observed in physical crystals. In addition, Y-Junctions are organized into lines, called grain boundaries, from the retinal center to the periphery. In physical crystals, regardless of the initial distribution of defects, defects can coalesce into grain boundaries via the mobility of individual particles. By imaging in live fish, we demonstrate that grain boundaries in the cone mosaic instead appear during initial mosaic formation, without requiring defect motion. Motivated by this observation, we show that a computational model of repulsive cell-cell interactions generates a mosaic with grain boundaries. In contrast to paradigmatic models of fate specification in mostly motionless cell packings, this finding emphasizes the role of cell motion, guided by cell-cell interactions during differentiation, in forming biological crystals. Such a route to the formation of regular patterns may be especially valuable in situations, like growth on a curved surface, where the resulting long-ranged, elastic, effective interactions between defects can help to group them into grain boundaries.

Physical limits to sensing material properties.

All materials respond heterogeneously at small scales, which limits what a sensor can learn. Although previous studies have characterized measurement noise arising from thermal fluctuations, the limits imposed by structural heterogeneity have remained unclear. In this paper, we find that the least fractional uncertainty with which a sensor can determine a material constant λ0 of an elastic medium is approximately [Formula: see text] for a â‰« d â‰« ξ, [Formula: see text], and D > 1, where a is the size of the sensor, d is its spatial resolution, ξ is the correlation length of fluctuations in λ0, Δλ is the local variability of λ0, and D is the dimension of the medium. Our results reveal how one can construct devices capable of sensing near these limits, e.g. for medical diagnostics. We use our theoretical framework to estimate the limits of mechanosensing in a biopolymer network, a sensory process involved in cellular behavior, medical diagnostics, and material fabrication.

Apical stress fibers enable a scaling between cell mechanical response and area in epithelial tissue.

Biological systems tailor their properties and behavior to their size throughout development and in numerous aspects of physiology. However, such size scaling remains poorly understood as it applies to cell mechanics and mechanosensing. By examining how the Drosophila pupal dorsal thorax epithelium responds to morphogenetic forces, we found that the number of apical stress fibers (aSFs) anchored to adherens junctions scales with cell apical area to limit larger cell elongation under mechanical stress. aSFs cluster Hippo pathway components, thereby scaling Hippo signaling and proliferation with area. This scaling is promoted by tricellular junctions mediating an increase in aSF nucleation rate and lifetime in larger cells. Development, homeostasis, and repair entail epithelial cell size changes driven by mechanical forces; our work highlights how, in turn, mechanosensitivity scales with cell size.

Statistics of noisy growth with mechanical feedback in elastic tissues.

Tissue growth is a fundamental aspect of development and is intrinsically noisy. Stochasticity has important implications for morphogenesis, precise control of organ size, and regulation of tissue composition and heterogeneity. However, the basic statistical properties of growing tissues, particularly when growth induces mechanical stresses that can in turn affect growth rates, have received little attention. Here, we study the noisy growth of elastic sheets subject to mechanical feedback. Considering both isotropic and anisotropic growth, we find that the density-density correlation function shows power law scaling. We also consider the dynamics of marked, neutral clones of cells. We find that the areas (but not the shapes) of two clones are always statistically independent, even when they are adjacent. For anisotropic growth, we show that clone size variance scales like the average area squared and that the mode amplitudes characterizing clone shape show a slow [Formula: see text] decay, where n is the mode index. This is in stark contrast to the isotropic case, where relative variations in clone size and shape vanish at long times. The high variability in clone statistics observed in anisotropic growth is due to the presence of two soft modes-growth modes that generate no stress. Our results lay the groundwork for more in-depth explorations of the properties of noisy tissue growth in specific biological contexts.

Robustness of Clocks to Input Noise.

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 24 h 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.

Optimal entrainment of circadian clocks in the presence of noise.

Circadian clocks are biochemical oscillators that allow organisms to estimate the time of the day. These oscillators are inherently noisy due to the discrete nature of the reactants and the stochastic character of their interactions. To keep these oscillators in sync with the daily day-night rhythm in the presence of noise, circadian clocks must be coupled to the dark-light cycle. In this paper, we study the entrainment of phase oscillators as a function of the intrinsic noise in the system. Using stochastic simulations, we compute the optimal coupling strength, intrinsic frequency, and shape of the phase-response curve, that maximize the mutual information between the phase of the clock and time. We show that the optimal coupling strength and intrinsic frequency increase with the noise, but that the shape of the phase-response curve varies nonmonotonically with the noise: in the low-noise regime, it features a dead zone that increases in width as the noise increases, while in the high-noise regime, the width decreases with the noise. These results arise from a tradeoff between maximizing stability-noise suppression-and maximizing linearity of the input-output, i.e., time-phase, relation. We also show that three analytic approximations-the linear-noise approximation, the phase-averaging method, and linear-response theory-accurately describe different regimes of the coupling strength and the noise.

Anisotropic Müller glial scaffolding supports a multiplex lattice mosaic of photoreceptors in zebrafish retina.

BACKGROUND: The multiplex, lattice mosaic of cone photoreceptors in the adult fish retina is a compelling example of a highly ordered epithelial cell pattern, with single cell width rows and columns of cones and precisely defined neighbor relationships among different cone types. Cellular mechanisms patterning this multiplex mosaic are not understood. Physical models can provide new insights into fundamental mechanisms of biological patterning. In earlier work, we developed a mathematical model of photoreceptor cell packing in the zebrafish retina, which predicted that anisotropic mechanical tension in the retinal epithelium orients planar polarized adhesive interfaces to align the columns as cone photoreceptors are generated at the retinal margin during post-embryonic growth. METHODS: With cell-specific fluorescent reporters and in vivo imaging of the growing retinal margin in transparent juvenile zebrafish we provide the first view of how cell packing, spatial arrangement, and cell identity are coordinated to build the lattice mosaic. With targeted laser ablation we probed the tissue mechanics of the retinal epithelium. RESULTS: Within the lattice mosaic, planar polarized Crumbs adhesion proteins pack cones into a single cell width column; between columns, N-cadherin-mediated adherens junctions stabilize Müller glial apical processes. The concentration of activated pMyosin II at these punctate adherens junctions suggests that these glial bands are under tension, forming a physical barrier between cone columns and contributing to mechanical stress anisotropies in the epithelial sheet. Unexpectedly, we discovered that the appearance of such parallel bands of Müller glial apical processes precedes the packing of cones into single cell width columns, hinting at a possible role for glia in the initial organization of the lattice mosaic. Targeted laser ablation of Müller glia directly demonstrates that these glial processes support anisotropic mechanical tension in the planar dimension of the retinal epithelium. CONCLUSIONS: These findings uncovered a novel structural feature of Müller glia associated with alignment of photoreceptors into a lattice mosaic in the zebrafish retina. This is the first demonstration, to our knowledge, of planar, anisotropic mechanical forces mediated by glial cells.

Period Robustness and Entrainability of the Kai System to Changing Nucleotide Concentrations.

Circadian clocks must be able to entrain to time-varying signals to keep their oscillations in phase with the day-night rhythm. On the other hand, they must also exhibit input compensation: their period must remain approximately one day in different constant environments. The posttranslational oscillator of the Kai system can be entrained by transient or oscillatory changes in the ATP fraction, yet is insensitive to constant changes in this fraction. We study in three different models of this system how these two seemingly conflicting criteria are met. We find that one of these (our recently published Paijmans model) exhibits the best tradeoff between input compensation and entrainability: on the footing of equal phase-response curves, it exhibits the strongest input compensation. Performing stochastic simulations at the level of individual hexamers allows us to identify a new, to our knowledge, mechanism, which is employed by the Paijmans model to achieve input compensation: at lower ATP fraction, the individual hexamers make a shorter cycle in the phosphorylation state space, which compensates for the slower pace at which they traverse the cycle.

Robustness of synthetic oscillators in growing and dividing cells.

Synthetic biology sets out to implement new functions in cells, and to develop a deeper understanding of biological design principles. Elowitz and Leibler [Nature (London) 403, 335 (2000)NATUAS0028-083610.1038/35002125] showed that by rational design of the reaction network, and using existing biological components, they could create a network that exhibits periodic gene expression, dubbed the repressilator. More recently, Stricker et al. [Nature (London) 456, 516 (2008)NATUAS0028-083610.1038/nature07389] presented another synthetic oscillator, called the dual-feedback oscillator, which is more stable. Detailed studies have been carried out to determine how the stability of these oscillators is affected by the intrinsic noise of the interactions between the components and the stochastic expression of their genes. However, as all biological oscillators reside in growing and dividing cells, an important question is how these oscillators are perturbed by the cell cycle. In previous work we showed that the periodic doubling of the gene copy numbers due to DNA replication can couple not only natural, circadian oscillators to the cell cycle [Paijmans et al., Proc. Natl. Acad. Sci. (USA) 113, 4063 (2016)PNASA60027-842410.1073/pnas.1507291113], but also these synthetic oscillators. Here we expand this study. We find that the strength of the locking between oscillators depends not only on the positions of the genes on the chromosome, but also on the noise in the timing of gene replication: noise tends to weaken the coupling. Yet, even in the limit of high levels of noise in the replication times of the genes, both synthetic oscillators show clear signatures of locking to the cell cycle. This work enhances our understanding of the design of robust biological oscillators inside growing and diving cells.

A thermodynamically consistent model of the post-translational Kai circadian clock.

The principal pacemaker of the circadian clock of the cyanobacterium S. elongatus is a protein phosphorylation cycle consisting of three proteins, KaiA, KaiB and KaiC. KaiC forms a homohexamer, with each monomer consisting of two domains, CI and CII. Both domains can bind and hydrolyze ATP, but only the CII domain can be phosphorylated, at two residues, in a well-defined sequence. While this system has been studied extensively, how the clock is driven thermodynamically has remained elusive. Inspired by recent experimental observations and building on ideas from previous mathematical models, we present a new, thermodynamically consistent, statistical-mechanical model of the clock. At its heart are two main ideas: i) ATP hydrolysis in the CI domain provides the thermodynamic driving force for the clock, switching KaiC between an active conformational state in which its phosphorylation level tends to rise and an inactive one in which it tends to fall; ii) phosphorylation of the CII domain provides the timer for the hydrolysis in the CI domain. The model also naturally explains how KaiA, by acting as a nucleotide exchange factor, can stimulate phosphorylation of KaiC, and how the differential affinity of KaiA for the different KaiC phosphoforms generates the characteristic temporal order of KaiC phosphorylation. As the phosphorylation level in the CII domain rises, the release of ADP from CI slows down, making the inactive conformational state of KaiC more stable. In the inactive state, KaiC binds KaiB, which not only stabilizes this state further, but also leads to the sequestration of KaiA, and hence to KaiC dephosphorylation. Using a dedicated kinetic Monte Carlo algorithm, which makes it possible to efficiently simulate this system consisting of more than a billion reactions, we show that the model can describe a wealth of experimental data.

Vertex stability and topological transitions in vertex models of foams and epithelia.

In computer simulations of dry foams and of epithelial tissues, vertex models are often used to describe the shape and motion of individual cells. Although these models have been widely adopted, relatively little is known about their basic theoretical properties. For example, while fourfold vertices in real foams are always unstable, it remains unclear whether a simplified vertex model description has the same behavior. Here, we study vertex stability and the dynamics of T1 topological transitions in vertex models. We show that, when all edges have the same tension, stationary fourfold vertices in these models do indeed always break up. In contrast, when tensions are allowed to depend on edge orientation, fourfold vertices can become stable, as is observed in some biological systems. More generally, our formulation of vertex stability leads to an improved treatment of T1 transitions in simulations and paves the way for studies of more biologically realistic models that couple topological transitions to the dynamics of regulatory proteins.

Discrete gene replication events drive coupling between the cell cycle and circadian clocks.

Many organisms possess both a cell cycle to control DNA replication and a circadian clock to anticipate changes between day and night. In some cases, these two rhythmic systems are known to be coupled by specific, cross-regulatory interactions. Here, we use mathematical modeling to show that, additionally, the cell cycle generically influences circadian clocks in a nonspecific fashion: The regular, discrete jumps in gene-copy number arising from DNA replication during the cell cycle cause a periodic driving of the circadian clock, which can dramatically alter its behavior and impair its function. A clock built on negative transcriptional feedback either phase-locks to the cell cycle, so that the clock period tracks the cell division time, or exhibits erratic behavior. We argue that the cyanobacterium Synechococcus elongatus has evolved two features that protect its clock from such disturbances, both of which are needed to fully insulate it from the cell cycle and give it its observed robustness: a phosphorylation-based protein modification oscillator, together with its accompanying push-pull read-out circuit that responds primarily to the ratios of different phosphoform concentrations, makes the clock less susceptible to perturbations in protein synthesis; the presence of multiple, asynchronously replicating copies of the same chromosome diminishes the effect of replicating any single copy of a gene.

Epithelial tricellular junctions act as interphase cell shape sensors to orient mitosis.

The orientation of cell division along the long axis of the interphase cell--the century-old Hertwig's rule--has profound roles in tissue proliferation, morphogenesis, architecture and mechanics. In epithelial tissues, the shape of the interphase cell is influenced by cell adhesion, mechanical stress, neighbour topology, and planar polarity pathways. At mitosis, epithelial cells usually adopt a rounded shape to ensure faithful chromosome segregation and to promote morphogenesis. The mechanisms underlying interphase cell shape sensing in tissues are therefore unknown. Here we show that in Drosophila epithelia, tricellular junctions (TCJs) localize force generators, pulling on astral microtubules and orienting cell division via the Dynein-associated protein Mud independently of the classical Pins/Gαi pathway. Moreover, as cells round up during mitosis, TCJs serve as spatial landmarks, encoding information about interphase cell shape anisotropy to orient division in the rounded mitotic cell. Finally, experimental and simulation data show that shape and mechanical strain sensing by the TCJs emerge from a general geometric property of TCJ distributions in epithelial tissues. Thus, in addition to their function as epithelial barrier structures, TCJs serve as polarity cues promoting geometry and mechanical sensing in epithelial tissues.

Patterning the cone mosaic array in zebrafish retina requires specification of ultraviolet-sensitive cones.

Cone photoreceptors in teleost fish are organized in precise, crystalline arrays in the epithelial plane of the retina. In zebrafish, four distinct morphological/spectral cone types occupy specific, invariant positions within a regular lattice. The cone lattice is aligned orthogonal and parallel to circumference of the retinal hemisphere: it emerges as cones generated in a germinal zone at the retinal periphery are incorporated as single-cell columns into the cone lattice. Genetic disruption of the transcription factor Tbx2b eliminates most of the cone subtype maximally sensitive to ultraviolet (UV) wavelengths and also perturbs the long-range organization of the cone lattice. In the tbx2b mutant, the other three cone types (red, green, and blue cones) are specified in the correct proportion, differentiate normally, and acquire normal, planar polarized adhesive interactions mediated by Crumbs 2a and Crumbs 2b. Quantitative image analysis of cell adjacency revealed that the cones in the tbx2b mutant primarily have two nearest neighbors and align in single-cell-wide column fragments that are separated by rod photoreceptors. Some UV cones differentiate at the dorsal retinal margin in the tbx2b mutant, although they are severely dysmorphic and are eventually eliminated. Incorporating loss of UV cones during formation of cone columns at the margin into our previously published mathematical model of zebrafish cone mosaic formation (which uses bidirectional interactions between planar cell polarity proteins and anisotropic mechanical stresses in the plane of the retinal epithelium to generate regular columns of cones parallel to the margin) reproduces many features of the pattern disruptions seen in the tbx2b mutant.

Exchange of stability as a function of system size in a nonequilibrium system.

In equilibrium systems with short-ranged interactions, the relative stability of different thermodynamic states generally does not depend on system size (as long as this size is larger than the interaction range). Here, we use a large deviations approach to show that, in contrast, different states can exchange stability as system size is varied in a driven, bistable reaction-diffusion system. This striking effect is related to a shift from a spatially uniform to a nonuniform transition state and should generically be possible in a wide range of nonequilibrium physical and biological systems.

Coupling mechanical deformations and planar cell polarity to create regular patterns in the zebrafish retina.

The orderly packing and precise arrangement of epithelial cells is essential to the functioning of many tissues, and refinement of this packing during development is a central theme in animal morphogenesis. The mechanisms that determine epithelial cell shape and position, however, remain incompletely understood. Here, we investigate these mechanisms in a striking example of planar order in a vertebrate epithelium: The periodic, almost crystalline distribution of cone photoreceptors in the adult teleost fish retina. Based on observations of the emergence of photoreceptor packing near the retinal margin, we propose a mathematical model in which ordered columns of cells form as a result of coupling between planar cell polarity (PCP) and anisotropic tissue-scale mechanical stresses. This model recapitulates many observed features of cone photoreceptor organization during retinal growth and regeneration. Consistent with the model's predictions, we report a planar-polarized distribution of Crumbs2a protein in cone photoreceptors in both unperturbed and regenerated tissue. We further show that the pattern perturbations predicted by the model to occur if the imposed stresses become isotropic closely resemble defects in the cone pattern in zebrafish lrp2 mutants, in which intraocular pressure is increased, resulting in altered mechanical stress and ocular enlargement. Evidence of interactions linking PCP, cell shape, and mechanical stresses has recently emerged in a number of systems, several of which show signs of columnar cell packing akin to that described here. Our results may hence have broader relevance for the organization of cells in epithelia. Whereas earlier models have allowed only for unidirectional influences between PCP and cell mechanics, the simple, phenomenological framework that we introduce here can encompass a broad range of bidirectional feedback interactions among planar polarity, shape, and stresses; our model thus represents a conceptual framework that can address many questions of importance to morphogenesis.

A dynamical model of ommatidial crystal formation.

The crystalline photoreceptor lattice in the Drosophila eye is a paradigm for pattern formation during development. During eye development, activation of proneural genes at a moving front adds new columns to a regular lattice of R8 photoreceptors. We present a mathematical model of the governing activator-inhibitor system, which indicates that the dynamics of positive induction play a central role in the selection of certain cells as R8s. The "switch and template" patterning mechanism we observe is mathematically very different from the well-known Turing instability. Unlike a standard lateral inhibition model, our picture implies that R8s are defined before the appearance of the complete group of proneural cells. The model reproduces the full time course of proneural gene expression and accounts for specific features of the refinement of proneural groups that had resisted explanation. It moreover predicts that perturbing the normal template can lead to eyes containing stripes of R8 cells. We observed these stripes experimentally after manipulation of the Notch and scabrous genes. Our results suggest an alternative to the generally assumed mode of operation for lateral inhibition during development; more generally, they hint at a broader role for bistable switches in the initial establishment of patterns as well as in their maintenance.

Robust circadian clocks from coupled protein-modification and transcription-translation cycles.

The cyanobacterium Synechococcus elongatus uses both a protein phosphorylation cycle and a transcription-translation cycle to generate circadian rhythms that are highly robust against biochemical noise. We use stochastic simulations to analyze how these cycles interact to generate stable rhythms in growing, dividing cells. We find that a protein phosphorylation cycle by itself is robust when protein turnover is low. For high decay or dilution rates (and compensating synthesis rates), however, the phosphorylation-based oscillator loses its integrity. Circadian rhythms thus cannot be generated with a phosphorylation cycle alone when the growth rate, and consequently the rate of protein dilution, is high enough; in practice, a purely posttranslational clock ceases to function well when the cell doubling time drops below the 24-h clock period. At higher growth rates, a transcription-translation cycle becomes essential for generating robust circadian rhythms. Interestingly, although a transcription-translation cycle is necessary to sustain a phosphorylation cycle at high growth rates, a phosphorylation cycle can dramatically enhance the robustness of a transcription-translation cycle at lower protein decay or dilution rates. In fact, the full oscillator built from these two tightly intertwined cycles far outperforms not just each of its two components individually, but also a hypothetical system in which the two parts are coupled as in textbook models of coupled phase oscillators. Our analysis thus predicts that both cycles are required to generate robust circadian rhythms over the full range of growth conditions.