Size control of the inner ear via hydraulic feedback
Abstract
Animals make organs of precise size, shape, and symmetry but how developing embryos do this is largely unknown. Here, we combine quantitative imaging, physical theory, and physiological measurement of hydrostatic pressure and fluid transport in zebrafish to study size control of the developing inner ear. We find that fluid accumulation creates hydrostatic pressure in the lumen leading to stress in the epithelium and expansion of the otic vesicle. Pressure, in turn, inhibits fluid transport into the lumen. This negative feedback loop between pressure and transport allows the otic vesicle to change growth rate to control natural or experimentally-induced size variation. Spatiotemporal patterning of contractility modulates pressure-driven strain for regional tissue thinning. Our work connects molecular-driven mechanisms, such as osmotic pressure driven strain and actomyosin tension, to the regulation of tissue morphogenesis via hydraulic feedback to ensure robust control of organ size.
Editorial note: This article has been through an editorial process in which the authors decide how to respond to the issues raised during peer review. The Reviewing Editor's assessment is that all the issues have been addressed (see decision letter).
https://doi.org/10.7554/eLife.39596.001Introduction
A fundamental question in developmental biology is how different organs acquire their proper sizes, which are necessary for their healthy function. The existence of control mechanisms is evident in the consistency of organ size in the face of intrinsic noise in biological reactions such as gene expression, and in the observed recovery from size perturbations during development (Waddington, 1959; Debat and Peronnet, 2013; Rao et al., 2002; Lestas et al., 2010). However, unlike in engineered systems, where there is often a clear distinction and hierarchy between the controller and the system, in organ growth one may not have a clear hierarchy—instead there may be control mechanisms distributed across tissues and across scales. Furthermore, in developmental biology, we observe an evolved system that is not necessarily robust to all experimental perturbations that we apply when trying to understand their control networks. Consequently, it can be difficult to distinguish what is necessary for growth from what controls size.
Identifying specific mechanisms that coordinate growth—to ultimately control organ size—has been difficult because the phenomenon of growth encompasses regulatory networks that can span the molecular to organismic. Classical organ transplantation and regeneration studies in the fly (Bryant and Levinson, 1985; Hariharan, 2015), mouse (Metcalf, 1963; Metcalf, 1964), and salamander (Twitty and Schwind, 1931) have indicated that both organ-autonomous and non-autonomous mechanisms control size. In his ‘chalone’ model, Bullough proposed growth duration to be regulated by an inhibitor of proliferation that is secreted by the growing organ and upon crossing a concentration threshold stops organ growth at the target size (Bullough and Laurence, 1964). Modern evidence for organ intrinsic chalones exists in myostatin for skeletal muscle, GDF11 for the nervous system, BMP3 for bone, and BMP2/4 for hair (McPherron et al., 1997; Wu et al., 2003; Plikus et al., 2008; Gamer et al., 2009). Several existing models for size control are based on global positional information regulating cell proliferation based on a morphogen gradient until final organ size is achieved (Day and Lawrence, 2000; Rogulja and Irvine, 2005; Wartlick et al., 2011). Other models emphasize the role of local cell-cell interactions in regulating cell proliferation or cell lineages to make tissues of the correct proportions (García-Bellido, 2009; Kunche et al., 2016). Given that cells are coupled to each other through cell-cell and cell-substrate contacts, physical constraints and tissue geometry provide tissue-level feedback. More recent models emphasize the role of tissue mechanics in regulating cell proliferation via anisotropic stresses and strain rates (Shraiman, 2005; Ingber, 2005; Savin et al., 2011; Hufnagel et al., 2007; Behrndt et al., 2012; Irvine and Shraiman, 2017; Nelson et al., 2017; Pan et al., 2016). From a molecular perspective, the insulin, Hippo (Dupont et al., 2011; Legoff et al., 2013; Pan et al., 2016) and TOR signaling pathways (Colombani et al., 2003; Zhang et al., 2000) have been well-established as regulators of organ size. Several studies have demonstrated that genetic mutation in these pathways is sufficient to alter organ or body size through increases in cell number, cell size, or both (Tumaneng et al., 2012), but the mechanisms that control size in the engineering sense (e.g. feedback of size on growth rate) are generally not known.
Most size control theories have focused on regulation of cell proliferation. Control may also arise from regulation of other parameters such as cell shape, material properties, transepithelial transport, adhesion, and the extracellular environment. In particular, fluid accumulation is a feature of developmental growth for several luminized organs including the embryonic brain (Desmond and Jacobson, 1977; Lowery and Sive, 2005), eye (Coulombre, 1956), gut (Bagnat et al., 2007), Kupffer’s vesicle (Navis et al., 2013; Dasgupta et al., 2018), the inner ear (Abbas and Whitfield, 2009; Hoijman et al., 2015), and the whole mammalian embryo (Chan et al., 2019). Water transport across an epithelium underlies these phenomenon (Frömter and Diamond, 1972; Günzel and Yu, 2013; Rubashkin et al., 2006; Fischbarg, 2010), and for the developing brain and eye it was shown that fluid accumulation coincides with increased hydrostatic pressure (Desmond and Jacobson, 1977; Coulombre, 1956). Just as water is fundamental to the size and function of a cell’s cytoplasm, the fluids filling the lumens of these organs, which are central to their development and physiological function, are fundamental components of these organs. Although specific ion transporters necessary for fluid accumulation have been identified (Lowery and Sive, 2005; Bagnat et al., 2007; Navis et al., 2013; Abbas and Whitfield, 2009), it is only very recently that we are beginning to get a glimpse of how how ion transport and transepithelial fluid flow are regulated, and their role in growth control, and much still remains to be explored .
Catch-up growth during development is the phenomenon where, after growth delay or perturbation, an organ transiently elevates its growth rate relative to other organs to get back on course. During fly development, if the growth of one imaginal disc is perturbed then a hormone, ecdysone, signals to the other imaginal discs to slow their growth such that the perturbed organ can catch-up and the animal’s coordinated growth can resume (Parker and Shingleton, 2011). The phenomenon of catch-up growth clarified ecdysone’s activity as being important for size control. Catch-up growth also occurs in vertebrates: if an infant heart or kidney is transplanted into an adult, it grows faster than the surrounding tissue to catch-up to a target size (Dittmer et al., 1974; Silber, 1976). Recently, the related phenomenon of organ symmetry has been addressed in the context of tails and the inner ear; but, the control mechanism underlying catch-growth was not clearly identified (Das et al., 2017; Green et al., 2017). Catch-up growth also occurs during bone growth and its study has clarified insulin signaling activity as being important for bone size control (Roselló-Díez and Joyner, 2015; Roselló-Díez et al., 2017; Roselló-Díez et al., 2018). Nonetheless, catch-up growth has been underused in the study of vertebrate-specific mechanisms of organ size control (Roselló-Díez et al., 2018).
Here we use a newly revealed instance of catch-up growth combined with physical theory to uncover how size control is achieved in the zebrafish otic vesicle, a fluid-filled closed epithelium that develops into the inner ear. We postulate that fluid pressure is a fundamental regulator of developmental growth in lumenized organs and hydraulic feedback can give rise to robust control of size.
Results
In toto imaging of otic vesicle development shows lumenal inflation dominates growth, not cell proliferation
We sought to determine how size control is achieved in the zebrafish otic vesicle, a 3D lumenized epithelial cyst that becomes the inner ear. Prior studies used qualitative observations and 4D imaging to examine the formation of the otic vesicle (Haddon and Lewis, 1996; Hoijman et al., 2015; Dyballa et al., 2017). To systematically investigate inner ear morphogenesis at longer timescales between 12–45 hours post-fertilization (hpf), we used high-resolution 3D+t confocal imaging combined with automated algorithms for quantifying cell and tissue morphology (Figure 1—figure supplement 1A–F) (Megason, 2009). Beginning at 12 hpf, bilateral regions of ectoderm adjacent to the hindbrain proliferate and subcutaneously accumulate to form the otic placodes (Figure 1—video 1, Figure 1A). The complex morphology of the inner ear arises from progressive changes in cell number, size, shape, and arrangement along with tissue-level patterns of polarization (12–14 hpf, Figure 1A), mesenchymal-to-epithelial transition (14–16 hpf, Figure 1B) and cavitation (16–24 hpf, Figure 1C). These steps build a closed ovoid epithelial structure, the otic vesicle, filled with a fluid called endolymph (Figure 1—figure supplement 1G–H). After assembly, the otic vesicle undergoes a period of rapid growth (16–45 hpf, Figure 1D) prior to the development of more complex substructures such as the semicircular canals and endolymphatic sac.
To evaluate growth kinetics, we used 3D image analysis (Figure 1—figure supplement 1I–M) to quantify a number of morphodynamic parameters between 16 and 45 hpf. During this period, cell number increased nearly three-fold from 415 ± 26 to 1106 ± 52 cells (blue curve, Figure 1E, for all data-points in Figure 1 n = 10 otic vesicles, data spread is the standard deviation). However, cell proliferation was offset by a decrease in average cell size from 0.55 ± 0.02 pl at 16 hpf to 0.34 ± 0.03 pl at 28 hpf and stayed constant thereafter (red curve). Tissue volume, the product of cell number and average cell size, remained effectively constant (230.6 ± 7.4 pl) until 28 hpf and subsequently increased linearly by 132 pl to 45 hpf (green curve, Figure 1F). The volume of the otic vesicle increased dramatically, by 572 pl from 235 ± 16 pl to 807 ± 23 pl (orange curve). The majority of the increase in size of the otic vesicle is due to an increase in lumen volume (blue curve) from 0 to 440 ± 18 pl (77% of the total increase) while tissue growth contributed only 23% to the increase in size.
Pressure inflates the otic vesicle and stretches tissue viscoelastically
A mismatch between the volumetric growth of the lumen and the tissue enclosing the lumen indicates a potential role for otic tissue remodeling. Since the size of the luminal vesicle, which scales as the cube root of lumen volume, increases more rapidly than the surface area of the vesicle, which scales as the square root of the cell number enclosing that volume, we investigated how epithelial cell shape changes to accommodate growth. We observed a monotonic increase in average cell apical surface area ( is lumenal surface area, is the cell number, Figure 1G). Since the otic epithelium is not uniform in thickness, we examined regions of the epithelium that contribute to the stretch. Except for the future sensory patches at the anterior and posterior ends (poles), epithelial thickness of the remaining otic vesicle significantly decreased from 20 µm to 4 µm during otic vesicle growth (lateral and medial regions, Figure 1H). The increase in lumenal volume and large cell stretching rates suggested that the vesicle is pressurized and the epithelium is under tension.
In the absence of extrinsic forces, cells round-up to a spherical morphology during mitosis by balancing internal osmotic pressure with tension provided by cortical actomyosin (Stewart et al., 2011). To investigate the development of pressure derived stress in the epithelium, we used mitotic cells as 'strain gauges’ by measuring their deviation in aspect ratios from spheres. We observed that mitotic cells fail to round up fully in regions where the otic epithelium pushes against the hindbrain and ectoderm, in contrast to non-contact regions at the anterioposterior poles (Figure 2A,C–G). Furthermore, cell division planes are closely aligned with the surface-normal to the epithelium (red markers, Figure 2H) in comparison to the broader distribution exhibited by the non-contact cell populations (blue). The overall alignment progressively increases in developmental time as cells become more stretched (Figure 2I), consistent with mechanical stress driven spindle alignment previously observed in various systems including the zebrafish gastrula (Campinho et al., 2013), fly imaginal disc (Legoff et al., 2013), and zebrafish pre-enveloping layer (Xiong et al., 2014). Given that the otic vesicle is wedged between the hindbrain and skin (Figure 1—figure supplement 1G–H), we examined the impact of its volumetric growth on these tissues. We reasoned that if pressure is present, the vesicle would exhibit higher rigidity and consequently deform neighboring structures as it increased in size. To test this idea, we quantified the indentation of the hindbrain and otic vesicle interface. We observed that as the vesicle grows, the initially planar hindbrain surface indents in, and the skin bulges out (Figure 2B–F,J).
To directly determine the presence of pressure within the otic vesicle, we developed a novel pressure probe able to accurately measure small pressures in small volumes of liquid. This probe consists of a solid-state piezo-resistive sensor coupled to a glass capillary needle filled with water (Figure 3A–B, see Materials and methods). This device is capable of measuring pressure differences of 5 Pascals (≈0.5 mm of water depth) across the range of 50–400 Pascals (Figure 3—figure supplement 1A). Prior to 30 hpf, lumenal pressure is too low for the needle to penetrate the epithelium. From 30 hpf onwards, we observed that the needle can penetrate into the otic vesicle with no observable volume change due to leakage around the needle (Figure 3B). Pressure is transmitted from the otic vesicle lumen through the needle tip to the sensor. Readings after puncture increased gradually before reaching a stable pressure level (Figure 3D). The positive pressure remained until the glass capillary was withdrawn from the otic vesicle, after which the pressure reading dropped to the baseline value (hydrostatic pressure of the buffer due to its depth in the petri dish), further indicating that a pressure difference exists across the epithelium (Figure 3—figure supplement 1F). We measured the pressure level at 30, 36 and 48 hpf and found that the pressure level gradually increases from 100 Pa to upwards of 300 Pa (Figure 3C and D). We are uncertain whether there is a drop in pressure upon insertion of the pressure probe into the otic vesicle because there is no alternate measuring device. These values are similar to prior measurements of the much larger inner ear of adult guinea pigs (Feldman et al., 1979).
To directly test if lumenal pressure ‘inflates’ the otic vesicle to drive inner ear growth, we punctured otic vesicles at different stages between 25–45 hpf (right vesicle in Figure 3E). Immediately following puncture, we observed a significant decrease in vesicle diameter (white arrows, Figure 3E) and loss of lumenal volume (≈30–40%) (Figure 3F). Examination of punctured vesicles showed that as the vesicle shrunk, the epithelium became thicker (Figure 3E and G). Indeed, the excess surface area of the lumenal cavity was absorbed by a significant change in epithelial cell shape to become more columnar while preserving cell volume (in-plane:normal diameter change from 6.7 ± 0.2 µm:13.2 ± 2.9 µm to 5.9 ± 0.2 µm:18.4 ± 4.1 µm at 30hpf) (Figure 3H). A similar transition in cell shape is seen when puncturing was conducted at later stages in development (Figure 3—figure supplement 2E), but importantly to a less columnar resting state suggesting a viscous component. Together, the puncturing experiments provided three insights into the mechanics of the otic vesicle: (i) the lumenal fluid is under hydrostatic pressure that is released when the vesicle is punctured, (ii) lumenal pressure generates stress in the epithelium that alters the shape of epithelial cells, causing them to stretch and become flatter, and (iii) the epithelial tissue response is viscoelastic, being elastic on short time scales, consistent with the epithelium becoming thicker immediately after puncturing, and viscous at longer time scales, consistent with long-term irreversible deformations.
Theoretical framework linking tissue geometry, fluid flux, and osmotic pressure
Given the complex interplay of lumenal pressure, geometry, and viscoelastic mechanics associated with growth, we sought to develop a mathematical model that accounts for these features (Figure 3I, See Materials and methods for mathematical model) (Ruiz-Herrero et al., 2017). In a spherically symmetric setting, the relationship between average vesicle radius (), wall thickness (), and tissue growth rate () can be specified as . Similarly, the relationship between growth in lumenal volume () and transport across the lumenal surface (of area ) is related to fluid influx per unit surface area . Defining as the homeostatic pressure required to balance the chemi-osmotic potential driving fluid flux and as the flux in the absence of a pressure differential, we may write the wild-type fluid flux as = - where is the permeability coefficient. Intuitively, is the fluid flux that maintains the homeostatic pressure , which in turn, remodels tissue to accommodate the incoming fluid.
Changes in luminal volume can be used to directly determine fluid flux because water is incompressible at low pressures. By using population-averaged measurements of lumenal volume and surface area, we calculate that after an initial rapid expansion (16–20 hpf), flux was approximately constant () throughout the period 21–45 hpf (Figure 4A). The flux is initially high when there is no pressure but then quickly goes down as pressure builds. Thereafter, fluid accumulation can only occur through viscous expansion of the vesicle. Interestingly, analysis of the system of equations in our model shows that the vesicle will adjust endolymph flux to account for perturbations to vesicle size, via a mechanical feedback loop that links pressure to flux (Figure 3J). Such a control system could be useful to correct natural as well as experimentally induced asymmetry across the left-right axis as we have found in early zebrafish inner ear development (Green et al., 2017), and in the whole mammalian embryo (Chan et al., 2019).
Model prediction and validation: pressure negatively regulates fluid flux
The model predicts that loss of fluid from the lumen (such as by puncture) should lead to a pressure drop and an increased rate of fluid flux back into the lumen and thus a higher than normal growth rate until size is restored, a phenomenon called catch-up growth. To test these model predictions, we experimentally examined whether pressure and fluid flux couple to each other to result in force-based feedback control of development. We first examined the response of the otic epithelium to puncturing perturbations between 25–45 hpf. We injected fluorescent dye (Alexa Fluor 594, 759 MW) into the fluid outside the inner ear (the perilymph) and tracked its movement into the lumen (Figure 3—figure supplement 2A). Puncturing the otic vesicle and withdrawing the needle created a loss in lumenal volume and allowed the dye from the perilymph to move into the lumen immediately (Figure 3—figure supplement 2B,C). However, when the dye was injected into the perilymph 5 min after the puncture, there was no rapid movement of dye into the lumen (Figure 3—figure supplement 2D). This showed that the otic epithelium rapidly seals after puncture and restores the epithelial barrier.
Next, we punctured the vesicle at 30 hpf, withdrew the needle, and evaluated its growth relative to the unpunctured contralateral vesicle (control) from 30 to 45 hpf by simultaneously imaging both otic vesicles. Interestingly, we observed the complete regeneration of lumenal volume in punctured vesicles by an increased growth rate relative to wild type to restore bilateral symmetry (Figure 4B–E, Figure 4—video 1). The cell shape changes were also reversed (Figure 4F) suggesting that the vesicle was re-pressurized. The rate of regeneration from fluid flux () was high immediately after puncture with a slow, gradual decay as bilateral symmetry is restored (Figure 4G). During the rapid recovery phase, in the punctured vesicle (blue curve) was 2-5X higher than that in the unpunctured vesicle (red curve). Our model predicts that upon loss of pressure from puncturing, the vesicle dynamically adjusts the fluid flux in linear proportion to volume lost (Equation 13 Materials and methods). To test this prediction, we pooled data from multiple punctured embryos regenerating from varying levels of pressure loss, to measure how fluid flux related to the volume loss. Consistently, fluid flux in the punctured vesicle correlates with the difference in lumenal volumes between the left and right vesicles (blue markers in Figure 4H and Figure 4—figure supplement 1 for other developmental stages). These data together show that vesicle pressure negatively regulates fluid flux and suggest that this feedback could buffer variations in size and drive catch-up growth in the otic vesicle.
Ion pumps are required for lumenal expansion
The transport of salts and fluid across an epithelium can occur through a variety of mechanisms involving transcytosis, electrogenic pumps/transporters and aquaporins for transcellular or paracellular flow (Preston et al., 1992; Hill and Shachar-Hill, 2006; Fischbarg, 2010). Paracellular transport refers to the transfer of fluid across an epithelium by passing through the intercellular space between the cells. This is in contrast to transcellular transport, where fluid travels through the cell, passing through both the apical membrane and basolateral membrane. Previous work in the chick otic vesicle identified the activity of Na+-K+-ATPase in setting up a transmural potential (Represa et al., 1986) to drive the selective movement of water and ions. In the zebrafish, a role for ion pumps in ear growth is supported by the previous identification of the Na+-K+-Cl− transporter Slc12a2 as defective in little ear mutants (Abbas and Whitfield, 2009). We administered ouabain, an inhibitor of Na+-K+-ATPase pump activity, to embryos at the 20 hpf stage and quantified vesicle morphology at 30 hpf. We observed a dose-dependent decrease in otic vesicle volume (blue) and wall thickness deformation (red) compared to the wild-type values (Figure 5A) consistent with previous work (Hoijman et al., 2015). In punctured embryos at 25 hpf, adding 500 µM ouabain to the buffer completely inhibited further growth (left vesicle) and post-puncture regeneration (right vesicle in Figure 5B–C, Figure 5—video 1). Knockdown of Na+-K+-ATPase expression in morpholino-injected embryos inhibited lumenal fluid transport in a dose-dependent manner (Figure 5D). We additionally find that otic vesicle growth is sensitive to variations in extracellular pH and blockers of chloride channel activity (Figure 5E–F). Together, these data argue that a network of ion transporters for and is required for fluid flux into the lumen.
Patterning of tissue material properties causes local differences in epithelial thinning
Our minimal mathematical model assumes that the vesicle is spherical allowing us to understand and predict the qualitative trends of our experiments. For pressure acting inside a thin spherical shell, the tensional tissue-stress—the force pulling cells apart that arises from the radially-outward pushing force of hydrostatic pressure—is . Since the tissue is elastic on short timescales and viscous on long timescales, the radial strain-rate—the change in radius of the otic vesicle, ()—may be related to via the constitutive relation for a Maxwell fluid given by , where is the tissue shear modulus—the material property that relates force experienced to deformation—and is the ratio of the tissue viscosity μ and elasticity .
For long timescales, we can use Stokes’ law and force balance (Equations 14, 16, Materials and methods) to derive an effective tissue viscosity where
Using this relationship, our morphodynamic measurements, and pressure measurements we estimate the effective viscosity of the otic vesicle tissue to be about 6.3 ± 0.30 × 106 Pa*s from 24 to 36 hpf and then 2.2 ± 0.13 × 107 Pa*s from 36 to 48 hpf (see Materials and methods for error propagation calculations). These values are within the range of tissue viscosities that had been experimentally measured, (Gordon et al., 1972), and indicate that the otic vesicle’s tissue becomes more viscous through development.
Since the vesicle is not actually spherical (Figure 1—figure supplement 1N,O) and the epithelium is not uniform in thickness (Figure 1H; Hoijman et al., 2015), we examined whether (i) the non-spherical vesicle shape creates non-uniform stress distribution as in Laplace’s law and/or (ii) non-uniform patterning of the material properties produce differential strain among cells. To test the first of these possibilities, we tracked anteroposterior pole cells (future sensory) and examined their shapes as they moved from high-curvature to low-curvature regions of the lumenal surface due to epithelial tread milling caused by regional differences in proliferation and emigration (Figure 1—video 1). We found that cells retained their columnar shapes independent of the underlying tissue curvature, suggesting that material property patterning may instead contribute to differential cell strain-rates (Figure 6—figure supplement 1A,B). To test the second possibility, we quantified spatial differences in elasticity () and viscosity () of the otic vesicle using puncturing experiments (Figure 6A). We observed that by eliminating pressure by puncture, medial and lateral cells deformed significantly more compared to the pole cells, indicating that they are softer (smaller ) (Figure 6B). Likewise, the resting shapes (post-puncture) of medial and lateral cells were more stretched out as the otic vesicle progressed in time, indicative of their lower viscosity (, Figure 6C). The observation that the medial and lateral cells become more viscous during development (Figure 6C) agrees with the increased effective viscosity we independently derived from our model using measurements of unperturbed otic vesicle growth. Thus, puncturing perturbations to eliminate pressure forces allowed us to measure spatial differences in cell viscoelasticity.
Cell stretching and steady proliferation contribute to tissue viscosity
Our model, calculations, and experimental measurements of cell material properties show that cells progressively become more rigid and more viscous. With diminished ability to remodel cell shape, we examined how otic vesicle growth can be sustained with lumenal pressure. To sustain the same growth rates, our model predicts that overall tissue viscoelasticity should be invariant to changes in cell material properties. While otic tissue elasticity arises from reversible cell stretching (), tissue viscosity is the net result of irreversible cellular stretching () as well as proliferation-driven increase in tissue surface area. Thus, we speculated that cell stretching has a more significant role in early growth while proliferation plays a more important role in later stages of growth.
To test these predictions, we evaluated the growth in lumenal surface area () in terms of individual contributions from division () and cell stretching () (Figure 6D). Our analysis shows that lumenal surface area growth is linear through time (blue markers). To support this growth, the contribution from cell-stretching is high initially but monotonically decreasing (green markers) and buffered by division (red markers). A break-even point occurs at around 33 hpf when the contribution to tissue viscosity from cell proliferation exceeds that from cell stretching (dashed black line). Interestingly, our data also show that cell shape stabilizes by this time (Figure 1H). Thus, cell stretching and proliferation play complementary roles through time to sustain a uniform increase in lumenal surface area.
Tissue material properties are patterned through actomyosin regulation
To identify how cell material properties are patterned, we examined localization patterns of F-actin and Myosin II using transgenic zebrafish (Tg(actb2:myl12.1-eGFP)e2212 for visualizing myosin II distribution, and Tg(actb2:GFP-Hsa.UTRN)e116 for visualizing F-actin distribution (Behrndt et al., 2012). Both, F-actin (Figure 6E and Figure 6—video 1) and Myosin II (Figure 6F and Figure 6—video 2) were apically localized prior to lumenization through 12–16 hpf to form a band around the cavity. Through early growth between 16–22 hpf, gradual and non-uniform changes in the apical density of these molecules are observed. By 30 hpf, we find that these proteins are localized to apicolateral junctions inside cells (Figure 6G–H). Expression levels are retained at pole cells but reduced in medial and lateral cells. Using the movies, we tracked individual cells to understand the relationship between cell shape change and apical marker intensity. We find that wild-type cell deformation during normal growth is positively correlated to localized accumulation of F-actin and Myosin II (Figure 6I). In the transgenic embryos, we used a mosaic labeling strategy for tracking cells to measure the relationship between apical localization and deformation of individual cells immediately following puncture (Figure 6—figure supplement 1C). We find that upon puncture, the instantaneous deformation observed in individual cells is linearly correlated with the levels of apical localization of F-actin and Myosin II suggesting that actomyosin tension sets effective tissue elasticity (Figure 6J). As it is unclear what contribution the neighboring tissue has to the effective material properties of the growing otic vesicle, we are unable to distinguish whether the correlation between actomyosin patterns and tissue thinning is organ autonomous or whether elastic forces from neighboring tissue are influencing these behaviors.
To further link spatial patterning of actomyosin localization with epithelial thickness, we conducted loss-of-function experiments and used our model to interpret experimental results. Upon reducing cell elasticity (), our model predicts: (i) an increase in strain-rate () to equilibrate with pressure forces, and (ii) an increase in lumenal dimensions to accomodate increased strain and secretion rates. To decrease cell elasticity, we inhibited actin dynamics by treating embryos at different stages between 16–35 hpf with 100 μM cytochalasin D and used high frame-rate imaging (one frame/s) to measure vesicle deformations. As predicted by theory, we observed an increase in by a factor of 2-5X over the wild-type values (Figure 6K and Figure 6—video 3). The decrease in apical myosin fluorescence positively correlated with the increase in secretion rates (Figure 6L and Figure 6—video 3). In these embryos, the DV diameter was found to increase most, compared to LR and AP diameters that experience reactionary forces from the hindbrain and skin (Figure 6M and Figure 6—video 3). We also observed that embryos between 16–25 hpf lost volume, presumably, due to a loss in epithelial connectivity and lack of pressure needed for deformation (Figure 6—video 4). Together, these data show how spatial patterning of the actomyosin cytokskeleton can lead to spatially varied strain in responses to spatially uniform pressure, and thus contribute to regional differences in the otic vesicle epithelium during growth.
Discussion
Here, we show that hydrostatic lumenal pressure develops in the zebrafish otic vesicle in response to fluid transport across the otic epithelium to drive growth. We used in toto imaging and newly developed quantitative image analysis tools to track changes in cell number, tissue volume, and vesicle lumen volume—which is fluid flux because water is nearly incompressible. We developed a pressure probe device that is amenable to low pressures and small volumes, which enabled us to quantify a developmental increase in hydrostatic pressure. Furthermore, we identified and characterized a new instance of catch-up growth that we leveraged to develop a theoretical framework for otic vesicle size control. With the aid of a multiscale mathematical model, we hypothesized and experimentally confirmed the presence of a hydraulic negative feedback loop between pressure and fluid transport for achieving size control. Modeling helped us systematically integrate the individual contributions of cell physioligical mechanisms underlying pressure and fluid flux, cell proliferation and shape, vesicle geometry, tissue strain-rate and viscoelasticity, to show how growth of the early otic vesicle is controlled. The negative feedback architecture that we found is similar to the chalone model of size control in that the act of growth feeds back to inhibit the rate of growth. However, compared to chalones or morphogen based growth control strategies which are limited in speed by diffusion, the pressure based strategy allows nearly instant communication between different parts of a tissue via hydraulic coupling to allow for 'course corrections’ to developmental trajectories. Indeed, our study is in line with the fact that hydraulic interactions are relevant to the developmental growth of many internal organs with vesicular and tubular origins (Ruiz-Herrero et al., 2017), and most recently the entire mammalian embryo (Chan et al., 2019).
Cause of cell thinning and origin of endolymph
Prior work found that when cells enter into mitosis and round up, their neighbors are stretched and become thinner (Hoijman et al., 2015). This was interpreted as being the mechanism by which cells thin to increase the surface area of the otic vesicle. Here we show that the otic epithelium is fairly elastic and cells can re-thicken following loss of lumenal pressure when the vesicle is punctured (Figure 3G,H), so stretching by neighboring mitotic cells may be short lived. Rather, in our model we postulate that sustained cell thinning during otic vesicle development is caused by in-plane epithelial tension in response to lumenal pressure. It was also noted that cells decrease in volume during early stages of otic vesicle growth and suggested that volume lost from cells is used to inflate the lumen (Hoijman et al., 2015). We attribute the reduction in cell size during early otic development to be due to cell division, and note that the net tissue volume (number of cells multiplied by average cell size) is constant at this stage (Figure 1F). Further, our measurements show that the lumen volume continues to grow until it exceeds the tissue volume (Figure 1F). We thus infer that transepithelial fluid flow is the primary if not sole source of fluid accumulation in the lumen.
Potential application of the lumen growth model to other systems
Our model for vesicle growth integrates a range of cellular behaviors including division, transport, force generation, material property patterning, and tissue thinning. By tailoring our model equations to different geometries and growth parameters, a unified mathematical framework can be realized to understand size control in hollow organs including the eyes, brain, kidneys, vasculature, and heart. The advantages of such a mesoscale model are several. First, a mesoscale model can be more easily applied to other contexts since the level of abstraction is higher making it less dependent on the specifics of the original context (fewer parameters). Second, growth kinetics and geometry parameters can be experimentally measured using in toto imaging approaches. New optical technologies for measuring tissue stresses in vivo using oil droplets (Campàs et al., 2014) and laser-ablation (Campinho et al., 2013; Hoijman et al., 2015), ionic gradients using fluorecent ionophores (Adams and Levin, 2012), and pressure (Link et al., 2004) via probes—like the one developed here—hold promise in providing reliable biophysical measurements necessary for understanding morphogenesis. Indeed, some of these very approaches have recently been taken in helping understand how the size of a mammalian embryo is controlled (Chan et al., 2019). And finally, the contribution of different molecular pathways in regulating model parameters can be prioritized for experimental investigation.
Boundary conditions of the otic vesicle and our model
The otic vesicle is not growing in isolation. In the embryo, it is immediately surrounded by extracellular matrix, mesenchymal cells, skin, and the brain. Within our model, these influences are abstracted as the effective material properties of the otic vesicle tissue. In fact, they may set limits to growth where the tension within the tissue begins to increase rapidly. We are likely observing an influence of these boundary conditions when we observe the spatial patterning of actinomyosin localization and regional tissue thinning (Figure 6). This boundary condition may accelerate cellular and molecular feedback mechanisms that were beyond the scope of this work. For instance, the cells within the tissue may respond to elevated tension by modulating proliferation rates, which may effectively alter the material properties of the tissue and alter strain (Halder and Johnson, 2011; Gudipaty et al., 2017; Gnedeva et al., 2017).
Comparison of pressure-regulation mechanisms
Our integrated approach combining quantitative imaging and theory-guided experimentation allowed us to identify a novel hydraulic-based mechanism for regulatory control of 3D vesicle growth. This mechanism enables long-range, fast, and uniform transmission of force and connects effects at multiple scales from global pressure forces to supracellular tension and cell stretching mechanics to molecular-scale actions of ion pumps and actomyosin regulation. Later, in the adult ear, tight control of inner ear fluid pressure and ionic composition is necessary to properly detect sound, balance, and body position. Pressure is also important to maintain the structural integrity of organs. Dysregulation of pressure homeostasis can give rise to diseases including hypertension in the vasculature, Ménière’s disease in the inner ear, glaucoma in the eye, and hydrocephalus in the brain. Pressure homeostasis mechanisms may vary. In the inner ear, pressure is initially regulated by feedback between pressure and lumenal fluid flux. However later in development, we have found that a physical pressure relief valve is necessary for pressure homeostasis in the inner ear (Swinburne et al., 2018). Together, we expect that these new insights on ear development and physiology will be critical to the development of effective clinical therapies for hearing and balance disorders, and for understanding size control in closed epithelial tissues.
Materials and methods
Contact for reagent and resource sharing
Request a detailed protocolFurther information and requests for resources and reagents should be directed to and will be fulfilled by Lead Contact, Sean Megason (megason@hms.harvard.edu).
Experimental model and subject details
Request a detailed protocolEmbryos were collected using natural spawning methods and the time of fertilization was recorded according to the single cell stage of each clutch. Embryos are incubated at 28°C during imaging and all other times except room temperature during injection and dechorionation steps. Staging was recorded using hours post-fertilization (hpf) as a measure and aligned to the normal table (Kimmel et al., 1995).
Zebrafish strains and maintenance
Request a detailed protocolThe following fluorescent transgenic strains were used in this study: (i) nuclear-localized tomato and membrane-localized citrine (Tg(actb2:Hsa.H2B-tdTomato); Tg(actb2:mem-citrine)hm32,33), Tg(actb2:mem-citrine-citrine)hm30 (ii) membrane-localized mCherry2 (Tg(actb2:mem-mcherry2)hm29) (iii) Tg(actb2:myl12.1-eGFP)e2212 for visualizing myosin II distribution, and (iv) Tg(actb2:GFP-Hsa.UTRN)e116 for visualizing F-actin distribution (Aguet et al., 2016; Behrndt et al., 2012; Xiong et al., 2014; Xiong et al., 2013). All fish are housed in fully equipped and regularly maintained and inspected aquarium facilities. All fish related procedures were carried out with the approval of Institutional Animal Care and Use Committee (IACUC) at Harvard University under protocol 04877. Full details of procedures are given in Extended Experimental Procedures.
Method details
Timelapse confocal imaging
Request a detailed protocolA canyon mount was cast in 1% agarose from a Lucite-plexiglass template and filled with 1X Danieau buffer (Figure 1—figure supplement 1A). The composition of the Danieau buffer is 14.4 mM sodium chloride, 0.21 mM potassium chloride, 0.12 mM magnesium sulfate, 0.18 mM calcium nitrate, and 1.5 mM HEPES buffered to pH 7.6. The template created three linear-ridges of width 400 µm, depth of 1.5 mm, and length 5 mm (Figure 1—figure supplement 1B). Canyon-mounted embryos developed normally for at least 3 days with a consistent orientation (lateral or dorsal or dorso-lateral) and can be continuously imaged during this time. Embryos at 15 hpf stage were dechorionated using sharp tweezers (Dumont 55) and mounted dorsally or dorso-laterally (Figure 1—figure supplement 1C,D) into the immersed canyon mount with a stereoscope (Leica MZ12.5). Multiple embryos for concurrent imaging were mounted in arrays (Figure 1—figure supplement 1B). Live imaging was performed using a Zeiss 710 confocal microscope (objectives: Plan-Apochromat 20 × 1.0 NA, C-Apochromat 40 × 1.2 NA) with a home-made heating chamber maintaining 28°C. For experiments requiring the imaging of both left and right ears in an embryo simultaneously, embryos were mounted dorsally and a Plan-Apochromat 20 × 1.0 NA objective was used. The inner ear is situated closest to the embryo surface when viewed along the dorso-lateral axis. Dorso-lateral mounting permitted the imaging of the entire ear structure with the best resolution and minimizes the depth of imaging. High-resolution imaging with a C-Apochromat 40 × 1.2 NA objective facilitated the use of automated image analysis scripts for cell and lumen segmentation, and tracking the movement of fluorescent dyes. Laser wavelengths 488 nm, 514 nm, 561 nm and 594 nm lasers were used for confocal time-courses and other single Z-stacks. Embryos were immobilized by injecting 2.3 nl of 20ng/µl -bungarotoxin mRNA (paralytic) at the 1 cell stage for experiments requiring long-term time-lapse imaging (Swinburne et al., 2015).
Wild-type growth curves
Request a detailed protocolThe process of anaesthetizing an embryo to prevent twitching and preparing an embryo for continuous imaging can alter wild-type growth dynamics in the long-term. Therefore, we collected single Z-stacks for separate sets of embryos (n = 10–15) to establish the wild-type growth curves at hourly intervals between 16–45 hpf. These sets of embryos were immobilized rapidly by soaking in 1% tricaine.
Confocal microscope settings
Request a detailed protocolImage settings vary by brightness of signal from maternal deposit. For example, (please see Figure 1—video 1): labels:membrane-citrine; lasers: 514 nm (20 mW, 3%); objective: C-Apochromat 40 × 1.2 NA at 1.0 zoom; pixel dwell time: 1.58 µs; pinhole size: 89 µm; line averaging: 1; image spacing: 0.2 × 0.2 µm, and 1024 × 1024 pixels per image, with an interval of 1.0 µm through Z for 80 µm and a temporal resolution of 2 min. The starting Z location for the embryo is ≈20 µm above the top of animal pole to allow sufficient space for it to stay in the field of view or sink in the agarose (Figure 1—figure supplement 1E,F). A total of 25 control time-lapse (covering the developmental time-period of 15–45 hpf), 450 control Z-stacks (covering the developmental time-period of 15–45 hpf), 45 perturbation-related time-lapse datasets, and 65 perturbation-related Z-stacks were collected for the current report. Embryos were screened for their health before imaging.
Vesicular fluid pressure probe
Request a detailed protocolOur pressure probe design was inspired by capillary-based pressure sensing techniques (Hüsken et al., 1978; Tomos and Leigh, 1999). A piezo-resistive solid-state pressure sensor (Honeywell, HSC series) with high resolution (≈2 Pa, 2 kHz) and minimal mechanical deformation (detailed below) was chosen as the sensing element. The sensor was coupled via a high pressure fitting to a ≈2 cm long glass capillary (World Precision Instruments) with a conical tip of 6–13 µm inner diameter (Figure 3A). Before they were coupled, both the capillary and sensor were filled with deionized water and the sensor was carefully degassed to ensure the entire interior volume is filled with water. Thus, the fabricated pressure probe had a water-filled dead-end cavity with the only opening at the capillary tip. A detailed fabrication procedure will be available in a separate publication. The digital output from the sensor was sampled by a developer board at 10 Hz (Arduino Uno-R3, and in-house Matlab program). We calibrated our fabricated pressure probe by measuring the hydrostatic pressure at different water depths (Figure 3—figure supplement 1A) that matched with the sensor calibration provided by the manufacturer. Similar tests were conducted with various capillary diameters and ionic concentration within the bath (deionized water and a solution that resembles mature endolymph) to ensure there was no additional effect (Figure 3—figure supplement 1B). The composition of our synthetic endolymph was 5 mM sodium chloride, 150 mM potassium chloride, 0.2 mM calcium chloride, 0.5 mM glucose, 10 mM tris, buffered to pH 7.5.
To measure the lumenal pressure in otic vesicle, zebrafish embryos were immobilized by injecting -bungarotoxin mRNA and dorsolaterally positioned in a canyon mount as before. The pressure probe was mounted on a micro-manipulator typically used for injection. Puncturing was done under a stereo-microscope (Figure 3B, Figure 3—figure supplement 1F): the capillary tip was first placed next to the otic vesicle to measure the reference hydrostatic pressure, and then the tip was punctured into the otic vesicle from the lateral direction. A tight sealing was indicated by the vesicle being intact while the capillary tip was inside the vesicle. The pressure profiles are shown in Figure 3D. We took the mean pressure value at the plateau, that is after the initial pressure rise and before any drop due to leakage (Stage III in Figure 3—figure supplement 1F), as the fluid pressure inside the vesicle. The data is presented in Figure 3C in the main text. As a negative control, we also punctured in bulk tissue regions such as in the neural tube and measured no pressure rise.
Since the sensor, glass capillary and otic vesicle form a closed volume after puncture, any deformation on the sensing element is compensated by out-flux of the luminal fluid from the vesicle. To verify that the volume change of the piezoresistive element is negligible, and therefore does not significantly reduce the luminal pressure, we calculated the elastic deformation of the sensing membrane and compare against the otic vesicle volume. We disassembled the sensor housing and measured the membrane area to be a square with edge length =850 µm. The membrane thickness, , is estimated to be 5-50 µm (Ruiz et al., 2012). The material is simplified as an isotropic silicon plate with Young’s modulus =163 GPa (Chicot et al., 1996; Hess, 1996; Dolbow and Gosz, 1996) and Poisson’s ratio =0.27 (Hess, 1996; Gan et al., 1996). The small transverse displacement, , of a thin plate under an uniform transverse hydrostatic pressure, , can be calculated using the Kirchhoff-Love plate theory (Timoshenko and Woinowsky-Krieger, 1959) , where the bending stiffness . The boundary condition for the built-in edges satisfies , where is the in-plane normal of the edges. The transverse deformation field, , was obtained by solving the above partial differential equation with a finite element solver (Matlab, Mathwork). An analytical solution exists for the deformation at the center as (Timoshenko and Woinowsky-Krieger, 1959 Article 44) , which agrees well with the numerical solution (Figure 3—figure supplement 1C). The corresponding volume change, , was depicted in Figure 3—figure supplement 1D. Comparing to the volume of a 200 µm diameter sphere (dotted line in Figure 3—figure supplement 1D), the volume change is at least 2 orders of magnitude smaller and hence the resultant pressure drop is negligible.
We also estimated the time scale at which the endolymph is diluted by diffusive exchange. The vesicle is a spherical domain of µm diameter with the initial ionic concentration, mM, the estimated value in the wild type endolymph. It is connected with a conical tube ( full cone angle, mm long) with variable inner diameter of 5 - 15 µm. The tube is initially filled with liquid of 0 mM, despite the actual concentration being slightly higher due to exchange with the bath (Danieau buffer, mainly composed with 14.4 mM sodium chloride). The boundary condition of all surfaces is zero flux, except for the back end of the tip being fixed at 0 mM. The temporal evolution of concentration field was obtained by numerically solving the diffusion equation , , where is the concentration field and m2/s is the diffusion coefficient of sodium chloride in water (Guggenheim, 1954). Example solutions and the normalized average concentration in the vesicle, , is shown in Figure 3—figure supplement 1E. At the typical rise time (around 0.5 min) of the probing stage, (Stage II in Figure 3—figure supplement 1F), the concentration remains at about 70%-90% for inner diameter of 15 µm - 5 µm. We expect the impact on the pressure reading was small at this time scale. After about 5 to 12 min, the concentration drops to 10%, which may significant modify the chemical potential. Together with the imperfection in sealing, they could contribute to the fluctuation measured at the longer time scale. However, we have ignored some factors that can maintain the ionic concentration such as active transport of ions or higher viscosity in the lumen of the otic vesicle.
Ouabain and cytochalasin D treatment
Request a detailed protocolIn order to inhibit the activity of Na+,K+-ATPase, embryos at the 20 hpf stage were soaked in 1% DMSO + ouabain (Sigma Aldrich, CAS 11018-89-6) across a range of concentrations from 0 to 1 mM. Ouabain was stored at 10 mM in 1% DMSO and diluted to required concentrations in 1X Danieau buffer before use. Ear size was assessed at 30 hpf as an endpoint. For long-term imaging, ouabain was added to 1% agarose mold used for mounting the embryos. To ensure consistent penetration, 2.3 nl of 0.75 mM ouabain was injected (Nanoject) into the cardiac chamber for circulation throughout the embryo. Assuming an average embryo volume of ∼180 nl, this injected dose guarantees an effective concentration of 10 µM. In order to perturb the actin network in the otic vesicle, 2.3 nl of 2 mM cytochalasin D was injected into the cardiac chamber for an effective concentration of 25 µM in the ear.
Buffer pH and niflumic acid perturbation
Request a detailed protocolTo study the effect of pH on otic vesicle size, embryos at the 12 hpf stage were soaked in 1X Danieau buffer titrated to different pH levels ranging from 6.5 to 8.5 at 12 hpf. We chose to use the 12 hpf to ensure that the embryo pH homeostasis was adequately perturbed before ear growth commenced. We assessed sizes at 25 and 30 hpf. In order to inhibit the activity of chloride channels and pH regulation in the embryos, embryos at the 20 hpf stage were soaked in 1% DMSO + niflumic acid (Sigma Aldrich, CAS 4394-00-7) across a range of concentrations from 0 to 1 mM. Ear size and cell shape was assessed at 30 hpf.
Antisense morpholino injection
Request a detailed protocolA total of four α1-like and two β subunit Na,K-ATPase genes have been identified in the inner ear with distinct spatiotemporal patterns of expression (Blasiole et al., 2003). Antisense morpholino (Gene Tools LLC; Philomath, OR) with sequence (5’-gccttctcctcgtcccattttgctg-3’) targeted against the Na,K-ATPase α-subunit gene atp1a1a.1 (α1a.1) was developed to knockdown expression in the early otic vesicle (Blasiole et al., 2006). The ability of the morpholino to act specifically to knockdown translation of only the relevant isoform of the Na,K-ATPase mRNA was previously demonstrated using an in vitro translation assay (Blasiole et al., 2006). To examine the role of the Na,K-ATPase in controlling ear growth, the morpholino was injected into 1 cell wild-type zebrafish embryos. Here, we report our results from using two different doses consisting of 0.25 ng in Figure 5D. In comparison to wild-type phenotypes, 0.5 ng morphants developed smaller otic vesicles, displayed smaller or absent otoliths, curved tails, and lagged in overall development (data not shown). Higher doses of morpholino injection (>1 ng) made embryos unhealthy prior to otic vesicle lumenization.
Puncturing of otic vesicle
Request a detailed protocolTo study the development of pressure in the otic vesicle, embryos were mounted dorso-laterally in a canyon mount (1% agarose by weight) with 1X Danieau buffer and an unclipped glass needle was slowly inserted into the otic vesicle. The needle pierced the vesicle in a lateral direction. Puncturing locally affected epithelial connectivity, causing on average 1–2 otic cell and 1–2 ectodermal cell deaths. Lumenal fluid (endolymph) leaked out along the circumference of the needle and the epithelium (Figure 3—figure supplement 2B,C). Needles were positioned using a micromanipulator. The needle was later slowly withdrawn to study regeneration dynamics. Thereafter, the embryo was re-mounted in a dorsal orientation for imaging both the ears simultaneously.
Quantifying the viscoelastic material properties
Request a detailed protocolTo identify the viscoelastic material properties of otic cells, we punctured vesicles and noted the deformation in cell shape. The puncturing experiment was carried out at 5 hr intervals (25, 30, 35, 40 hpf) to determine trends in material property patterning (Figure 6A). We used a sample size of at each timepoint. The observed deformation in the shape of a cell (before vs. after puncture) located at position and time is inversely proportional to the spring constant (Figure 6B). The rate of change in the resting shapes (after puncture) is inversely proportional to the viscosity coefficient (Figure 6C).
Mathematical model: linking geometry, growth, mechanics and regulation
To understand otic vesicle growth, we developed a compact mathematical model that links vesicle geometry, tissue mechanics, and cellular behavior. Quantitative imaging identified aspects of cell behavior including cell division, cell size, cell shape, and material properties as being relevant to the size control problem. The process of realizing a multiscale model enabled the identification of two fundamental mechanisms regulating growth: (i) We identified a negative feedback signal linking the development of hydrostatic pressure to the inhibition of lumenal fluid flux for robust control of size; and, (ii) Spatial patterning of actomyosin contractility affected tissue response to pressure forces to shape the ear. In this section, we elaborate on the development of model equations and explain how theory-guided experimentation allowed us to arrive at the correct representation of the process.
Conservation of mass links tissue flux and fluid flux to geometry
Request a detailed protocolQuantitative analysis of vesicle geometry pointed to two critically changing parameters, namely, vesicle radii and tissue wall thickness (Figure 1—figure supplement 1N,O and Figure 1H). A simple model treats the geometry of the otic vesicle as a spherical shell (Figure 3I) of average radius (µm) and wall thickness (µm). As notation, variable subscripts , , and refer to entire otic vesicle, otic lumen, and otic tissue components respectively. To account for changes in geometry parameters from growth, we represent otic tissue growth rate as () and fluid flux (per unit surface-area) as (µm/hr).
Conservation of tissue mass implies the rate of change of the otic tissue volume is related to the change in tissue geometry. Since the average thickness µm and the smallest otic vesicle radius µm at 30 hpf, we can make the assumption that wall-thickness is relatively small compared to the average radius (), so that
We quantified the parameters , , and and found that average radius increased linearly (Figure 1—figure supplement 1N,O), wall-thickness (Figure 1H) reduced asymptotically, and tissue volume (Figure 1F) was constant initially (16–28 hpf, ) but linearly increased thereafter (28–45 hpf, ). Since tissue volume is a product of cell number () and average cell-size (), we further investigated the role of these cellular parameters during growth.
From Figure 1E, cell number monotonically increased non-linearly between 16–45 hpf. In contrast, average cell-size monotonically decreased till 28 hpf and stabilized to a constant value () thereafter. Therefore, between 16–28 hpf, the increase in cell number was offset by a decrease in cell-size (Figure 1E), effectively keeping tissue mass constant. In addition, between 28–45 hpf, increase in tissue mass occurred due to increase in cell number alone with a constant average cell-size.
Similar to volume, we note that lumenal surface area is the product of cell number () and average cell apical surface area (, µm2). Therefore, we investigated their individual contribution in driving the increase in surface area.
From Figure 1G and E, we find that lumenal surface area and cell number increase monotonically from 16 to 45 hpf. Average apical cell surface showed a saturating response instead. By analyzing both the terms in Figure 6D, we find that both terms contribute significantly to lumenal surface area growth in a temporally complementary manner. Cell stretching () harbors a more significant role during early growth (≤32 hpf) with a slow diminishing influence. In contrast, division () is more significant role during later growth (> 32 hpf). Note that this is a second role of division during growth, in addition to the earlier role in regulating tissue volume increase (Equation 3).
Conservation of fluid volume implies the rate of change in lumenal volume () is equal to product of surface area () and fluid flux .
For a spherical geometry, lumenal volume is and surface-area is , so that
In other words, fluid flux is the same as the rate of change of lumenal radius. From Figure 4A and Figure 1—figure supplement 1N, we find that average radius increased linearly and lumenal fluid flux was approximately constant () between 16–45 hpf.
Thus, our analysis of growth kinetics using simple conservation laws points to the role of mechanisms involved in regulating cell division, cell volume, fluid transport, and cell shape in controlling size.
Modeling pressure generation and feedback to fluid transport mechanisms
Request a detailed protocolSince lumenal volume growth contributes the most to vesicle growth, we examined the phenomenon of fluid transport into a closed cavity (Figure 1F). Prevailing theories of fluid transport suggest the movement of charged ions and fluid through intercellular junctions and channels on cell surfaces (Hill and Shachar-Hill, 2006; Fischbarg, 2010). To model the phenomenon of fluid transport into a closed cavity, we denote the rate of transport of solute as and its concentration in the lumen as . Then, fluid flux required to retain this concentration is given by the relationship:
Earlier, our analysis showed that lumenal fluid flux is a constant throughout growth. Under the assumption of isotonic transport ( is a constant), we inferred that lumenal solute flux is also constant throughout growth. Since solute and fluid flux are coupled, we investigated which of these parameters are regulated. Two possible scenarios exist: (i) Otic epithelium ensures a solute flux of exactly leading to a fluid flux of , or (ii) Otic epithelium restricts fluid flux to on account of wall distensibility and pressure gradient, which effectively retains a net solute flux of . Indeed, the existence and presence of lumenal pressure is evident in terms of causing shape deformation in hindbrain tissue and mitotic cells (Figure 2). To concretely show the existence of pressure and how it restricts fluid transport, we performed simple perturbation experiments:
Puncturing experiments to verify the presence of pressure by assaying the loss of lumenal volume (Figure 3).
Regeneration experiments to demonstrate that the otic vesicle is capable of increased rates of solute and fluid transport in the absence of pressure. We measured regeneration dynamics (Figure 4) and found that vesicles were able to increasing flux by a factor the wild-type values (Figure 4—figure supplement 1).
Experimental outcomes suggest that pressure creates an opposing force to the osmotic potential forces to inhibit the movement of fluid, thus setting growth rates appropriate to the developmental stage. Thus, in unperturbed wild-type embryos, the observed fluid flux () is a function of the difference between the osmotically-derived flux (, arising purely due the osmotic potential alone) and the deviation from a set-point pressure, that we term the homeostatic pressure in the vesicle . Making the assumption that the functional form is linear then suggests the equation
In punctured embryos, the loss in pressure () leads to a larger flux following regeneration ().
that is the increase in flux relative to the unperturbed value is proportional to the loss in pressure relative to its homeostatic value, with
For small changes in the radius in the neighborhood of the fixed point, the loss-in-pressure is proportional to the loss-in-volume , where and are unperturbed and perturbed volumes respectively), so that we arrive at the following relationship:
We validate model predictions of a correlation between regeneration flux () and volume-loss between punctured and unpunctured ears (). Our experimental data (regeneration in blue curve in Figure 4H and Figure 4—figure supplement 1E,F) depicts a proportional relationship, thus validating the accuracy of our model. Our data shows that the unperturbed flux and slope is approximately 3–5.
Modeling pressure-driven growth and deformation of the vesicle shape
Request a detailed protocolThe otic vesicle is expected to deform under the action of pressure, thus leading to growth and reshaping of the tissue. We therefore sought to assess how these forces affect individual cells. Given the vesicle geometry, cells in the otic vesicle experience three types of pressure-derived forces: (i) pressure force (, ) normal to the apical membranes in a direction that flatten cells, (ii) tissue stress (, ) distributed normal to lateral membranes, and (iii) reactionary or support forces from hindbrain and skin tissue.
To formalize this, we modify a recent framework introduced to study the growth of cellular cysts (Ruiz-Herrero et al., 2017). In a closed geometry, tissue stress is related to lumenal pressure by a simple force balance equation. For a spherical pressure vessel, the pressure force pushing one hemisphere is counterbalanced by tissue tension.
To link tissue stress to deformation, we investigated the material properties of the tissue (). We initially modeled otic tissue as being similar to a purely viscous fluid flowing under a tangential shear stress with a strain-rate of .
For a purely viscous fluid with viscosity coefficient () deforming in a spherically-restricted geometry, the shear stress () in the fluid layers relates to rate of radius change using Stokes’ law as:
Since a viscous material deforms continuously under the action of a force and retains the deformation when force is removed, we validated our assumption by studying the temporal trends in resting shapes of cells after pressure forces are eliminated. Indeed, we observed that the resting shape of cells (after puncture) had undergone irreversible deformation in developmental time (Figure 3G and D (red bars progressively decrease)). In other words, a cell at 35 hpf has a relatively squamous morphology at rest compared to the same cell at 25 hpf. When compared to the contralateral unpunctured ear, however, we observed an elastic response in punctured vesicles wherein cells had dynamically reshaped to more columnar morphologies (Figure 3E and F,H (difference between red and blue bars)). This suggested that the tissue behaved elastically in short time-scales and viscous in long time-scales. To incorporate elastic behavior on short timescales, we modeled the overall deformation as:
We next quantified cell material properties by quantifying the short (elastic) and long timescale (viscous) responses in response to puncturing perturbations (Figure 6A). The reversible component of a cell’s deformation, when all acting forces are eliminated, is inversely proportional to its elasticity (). The irreversible component of a cell’s deformation that accumulates over time is inversely proportional to the viscosity coefficient (). We also investigated and identified the molecular origins of material property patterning in terms of actomyosin networks that were found to be apically expressed in the otic cells (Figure 6). Thus, our model identifies the role of pressure forces in deforming cells and driving growth. Material property patterning was found to locally modulate the deformation to anisotropically shape the vesicle through growth.
Quantification and statistical analysis
Request a detailed protocolEach data point of morphodynamic quantification was obtained using our automated bioimage informatics pipeline (ITIAT) on images from 10 different embryos that were immobilized with -bungarotoxin mRNA. 30 time-points were analyzed—every hour between 16 and 45 hpf—which included 300 otic vesicle lumen measurements and more than 200,000 segmented cells. For Figure 1E–H, Figure 1—figure supplement 1N–O the translucent spread of the data is the standard deviation.
For analysis of mitotic cell deformation, we identified 54 mitotic cells. The spread of the data presented as translucent overlays in Figure 2G is the calculated standard deviation. For measuring the hindbrain deformation in Figure 2J, n = 10 embryos/timepoint were used.
Pressure measurements were acquired from 5, 8, and nine different embryos at 30, 36, and 48 hpf (Figure 3C).
For puncture experiments, 10 otic vesicles were measured that were punctured and there morphometrics were compared to 10 unpunctured otic vesicles (Figures 3E–H and 4E–H, Figure 4—figure supplement 1). The error bars are standard deviation and the p-values were obtained using a student t-test.
For ion-channel and pH perturbations, each data point is five embryos and error bars are standard deviation (Figure 5).
For calculating the effective tissue viscosity from Equation 1, the time-averaged data of all variables were obtained from the datasets shown in Figure 3C,E–F. The instantaneous growth rate, , was calculated from the difference of radius of each time point before averaging. A formula for error propagation (Ku, 1966) was derived from Equation 1:
where is the standard deviation of variable . The mean and standard deviations are summarized in the table below:
Mean (Standard deviation) | |||||
---|---|---|---|---|---|
hpf | Pa | µm | µm | 10-4 µm/s | 106 Pa⋅s |
24-36 | 116 (23.1) | 36 (0.18) | 13 (0.12) | 2.31 (0.73) | 6.26 (0.30) |
36-48 | 166 (23.4) | 48 (0.22) | 10 (0.21) | 2.15 (0.95) | 22.2 (1.29) |
For the derivation of relative differences in tissue elasticity and viscosity, each data point was obtained from 10 different punctured embryos (Figure 6A–C) and the translucent spread represents the standard deviation.
Figure 6I-J, 22 utrophin expressing cells and 22 myosin expressing cells were analyzed. Figure 6K–L–M, scatterplots with 15, 16, and 12 embryos were used respectively.
Data and software availability
Image analysis
Request a detailed protocolOur entire bioimage informatics pipeline called ‘In toto image analysis toolkit (ITIAT)’ is available online at: https://wiki.med.harvard.edu/SysBio/Megason/GoFigureImageAnalysis. The pipeline can be used for generating 3D surface reconstructions, automatic whole-cell and nuclei segmentations, and cell tracking. The code is open-source and written using the C++ programming language. The code can be downloaded and compiled in any platform using CMake, a tool for generating native Makefiles. The code is built by linking to pre-compiled open-source libraries, namely, ITK (http://www.itk.org) (Yoo et al., 2002) and VTK toolkit (http://www.vtk.org) (Schroeder et al., 2006) for image analysis and visualization respectively. We also used the open-source and cross-platform GoFigure2 application software for the visualization, interaction, and semi-automated analysis of image data (http://www.gofigure2.org) (Gelas, Mosaliganti, and Megason et al., in preparation). Measurement of mitotic cell aspect-ratios was carried out using ZEN (Carl Zeiss) software 3D distance functionality. Measurements were analysed and plotted with Matlab (Mathworks) and Microsoft Excel. To obtain 3D models of the otic vesicles, 2D contours were first placed along regularly-sampled z-planes in GoFigure2 (Figure 1I). 3D reconstructions were obtained using the Powercrust reconstruction algorithm (https://github.com/krm15/Powercrust) (Amenta et al., 2001). Automatic cell and lumen-segmentation were performed using ACME software for whole-cell segmentations (Mosaliganti et al., 2012; Xiong et al., 2013) (Figure 1—figure supplement 1J–M).
Data availability
All data generated or analyzed during this study are included in the manuscript and supporting files.
References
-
General principles for measuring resting membrane potential and ion concentration using fluorescent bioelectricity reportersCold Spring Harbor Protocols 2012:pdb.top067710.https://doi.org/10.1101/pdb.top067710
-
Membrane dynamics of dividing cells imaged by lattice light-sheet microscopyMolecular Biology of the Cell 27:3418–3435.https://doi.org/10.1091/mbc.e16-03-0164
-
The power crust, unions of balls, and the medial Axis transformComputational Geometry 19:127–153.https://doi.org/10.1016/S0925-7721(01)00017-7
-
Genetic control of single lumen formation in the zebrafish gutNature Cell Biology 9:954–960.https://doi.org/10.1038/ncb1621
-
Mitotic control by internal secretion: the role of the chalone-adrenalin complexExperimental Cell Research 33:176–194.https://doi.org/10.1016/S0014-4827(64)81025-9
-
The role of intraocular pressure in the development of the chick eye. I. control of eye sizeJournal of Experimental Zoology 133:211–225.https://doi.org/10.1002/jez.1401330202
-
Patterned disordered cell motion ensures vertebral column symmetryDevelopmental Cell 42:170–180.https://doi.org/10.1016/j.devcel.2017.06.020
-
Measuring dimensions: the regulation of size and shapeDevelopment 127:2977–2987.
-
Embryonic brain enlargement requires cerebrospinal fluid pressureDevelopmental Biology 57:188–198.https://doi.org/10.1016/0012-1606(77)90364-5
-
The growth of infant hearts grafted to young and adult ratsAmerican Journal of Anatomy 141:155–160.https://doi.org/10.1002/aja.1001410112
-
Measurement of the hydrostatic pressures of the cochlear compartmentsNeurological Research 1:11–18.https://doi.org/10.1080/01616412.1979.11739537
-
Route of passive ion permeation in epitheliaNature New Biology 235:9–13.https://doi.org/10.1038/newbio235009a0
-
The cellular and genetic bases of organ size and shape in DrosophilaThe International Journal of Developmental Biology 53:1291–1303.https://doi.org/10.1387/ijdb.072459ag
-
A rheological mechanism sufficient to explain the kinetics of cell sortingJournal of Theoretical Biology 37:43–73.https://doi.org/10.1016/0022-5193(72)90114-2
-
Recovery of shape and size in a developing organ pairDevelopmental Dynamics 246:451–465.https://doi.org/10.1002/dvdy.24498
-
The diffusion coefficient of sodium chlorideTransactions of the Faraday Society 50:1048–1051.https://doi.org/10.1039/tf9545001048
-
Claudins and the modulation of tight junction permeabilityPhysiological Reviews 93:525–569.https://doi.org/10.1152/physrev.00019.2012
-
Early ear development in the embryo of the zebrafish, Danio rerioThe Journal of Comparative Neurology 365:113–128.https://doi.org/10.1002/(SICI)1096-9861(19960129)365:1<113::AID-CNE9>3.0.CO;2-6
-
Organ size control: lessons from DrosophilaDevelopmental Cell 34:255–265.https://doi.org/10.1016/j.devcel.2015.07.012
-
Laser diagnostics of mechanical and elastic properties of silicon and carbon filmsApplied Surface Science 106:429–437.https://doi.org/10.1016/S0169-4332(96)00369-8
-
A new approach to epithelial isotonic fluid transport: an osmosensor feedback modelJournal of Membrane Biology 210:77–90.https://doi.org/10.1007/s00232-005-0847-3
-
Mitotic cell rounding and epithelial thinning regulate lumen growth and shapeNature Communications 6:7355.https://doi.org/10.1038/ncomms8355
-
Mechanical control of growth: ideas, facts and challengesDevelopment 144:4238–4248.https://doi.org/10.1242/dev.151902
-
Stages of embryonic development of the zebrafishDevelopmental Dynamics 203:253–310.https://doi.org/10.1002/aja.1002030302
-
Notes on the use of propagation of error formulasJournal of Research of the National Bureau of Standards 70:263–273.
-
Feedback, lineages and Self-Organizing morphogenesisPLOS Computational Biology 12:e1004814.https://doi.org/10.1371/journal.pcbi.1004814
-
Intraocular pressure in zebrafish: comparison of inbred strains and identification of a reduced melanin mutant with raised IOPInvestigative Opthalmology & Visual Science 45:4415–4422.https://doi.org/10.1167/iovs.04-0557
-
In toto imaging of embryogenesis with confocal time-lapse microscopyMethods in Molecular Biology 546:317–332.https://doi.org/10.1007/978-1-60327-977-2_19
-
The autonomous behaviour of normal Thymus graftsAustralian Journal of Experimental Biology and Medical Science 41:437–447.https://doi.org/10.1038/icb.1963.64
-
Restricted growth capacity of multiple spleen graftsTransplantation 2:387–392.https://doi.org/10.1097/00007890-196405000-00008
-
ACME: automated cell morphology extractor for comprehensive reconstruction of cell membranesPLOS Computational Biology 8:e1002780.https://doi.org/10.1371/journal.pcbi.1002780
-
Electrical properties of the otic vesicle epithelium in the chick embryoJournal of Embryology and Experimental Morphology 97:125–139.
-
Regulation of long bone growth in vertebrates; It is time to catch upEndocrine Reviews 36:646–680.https://doi.org/10.1210/er.2015-1048
-
Organ size control via hydraulically gated oscillationsDevelopment 144:4422–4427.https://doi.org/10.1242/dev.153056
-
Growth of baby kidneys transplanted into adultsArchives of Surgery 111:75–77.https://doi.org/10.1001/archsurg.1976.01360190077014
-
THE PRESSURE PROBE: a versatile tool in plant cell physiologyAnnual Review of Plant Physiology and Plant Molecular Biology 50:447–472.https://doi.org/10.1146/annurev.arplant.50.1.447
-
Organ size control by hippo and TOR pathwaysCurrent Biology 22:R368–R379.https://doi.org/10.1016/j.cub.2012.03.003
-
The growth of eyes and limbs transplanted heteroplastically between two species of AmblystomaJournal of Experimental Zoology 59:61–86.https://doi.org/10.1002/jez.1400590105
-
Engineering and algorithm design for an image processing api: a technical report on ITK--the insight toolkitStudies in Health Technology and Informatics 85:586–678.
-
Regulation of cellular growth by the Drosophila target of rapamycin dTORGenes & Development 14:2712–2724.https://doi.org/10.1101/gad.835000
Article and author information
Author details
Funding
National Institutes of Health (K25 HD071969)
- Kishore R Mosaliganti
Hearing Health Foundation
- Ian A Swinburne
National Institutes of Health (5F32HL097599)
- Ian A Swinburne
National Institutes of Health (DC010791)
- Sean G Megason
National Institutes of Health (DC015478)
- Sean G Megason
John D. and Catherine T. MacArthur Foundation
- L Mahadevan
National Science Foundation (BMMB 15-36616)
- L Mahadevan
Agency for Science, Technology and Research
- Chon U Chan
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We thank members of Megason lab for feedback, Mr. Dante D’India for fish care, and Ms. Suzanne Mosaliganti for help in proof-reading the manuscript. This work was supported by the National Institute of Health grant K25 HD071969 (KRM), Novartis Fellowship for Systems Biology (IAS), National Institute of Health grant 5F32HL097599 (IAS), A*STAR International Fellowship (CUC), Hearing Health Foundation (IAS), the MacArthur Foundation (LM), National Science Foundation grant BMMB 15–36616 (LM), National Institute of Health grant R01 DC010791 (SGM), and National Institute of Health Grant R01 DC015478 (SGM).
Ethics
Animal experimentation: This study was performed in strict accordance with the recommendations in the Guide for the Care and Use of Laboratory Animals of the National Institutes of Health. The Harvard Medical Area Standing Committee on Animals approved zebrafish work under protocol number 04487.
Copyright
© 2019, Mosaliganti et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 4,144
- views
-
- 334
- downloads
-
- 58
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Computational and Systems Biology
- Microbiology and Infectious Disease
Bacterial membranes are complex and dynamic, arising from an array of evolutionary pressures. One enzyme that alters membrane compositions through covalent lipid modification is MprF. We recently identified that Streptococcus agalactiae MprF synthesizes lysyl-phosphatidylglycerol (Lys-PG) from anionic PG, and a novel cationic lipid, lysyl-glucosyl-diacylglycerol (Lys-Glc-DAG), from neutral glycolipid Glc-DAG. This unexpected result prompted us to investigate whether Lys-Glc-DAG occurs in other MprF-containing bacteria, and whether other novel MprF products exist. Here, we studied protein sequence features determining MprF substrate specificity. First, pairwise analyses identified several streptococcal MprFs synthesizing Lys-Glc-DAG. Second, a restricted Boltzmann machine-guided approach led us to discover an entirely new substrate for MprF in Enterococcus, diglucosyl-diacylglycerol (Glc2-DAG), and an expanded set of organisms that modify glycolipid substrates using MprF. Overall, we combined the wealth of available sequence data with machine learning to model evolutionary constraints on MprF sequences across the bacterial domain, thereby identifying a novel cationic lipid.
-
- Computational and Systems Biology
- Neuroscience
Perception can be highly dependent on stimulus context, but whether and how sensory areas encode the context remains uncertain. We used an ambiguous auditory stimulus – a tritone pair – to investigate the neural activity associated with a preceding contextual stimulus that strongly influenced the tritone pair’s perception: either as an ascending or a descending step in pitch. We recorded single-unit responses from a population of auditory cortical cells in awake ferrets listening to the tritone pairs preceded by the contextual stimulus. We find that the responses adapt locally to the contextual stimulus, consistent with human MEG recordings from the auditory cortex under the same conditions. Decoding the population responses demonstrates that cells responding to pitch-changes are able to predict well the context-sensitive percept of the tritone pairs. Conversely, decoding the individual pitch representations and taking their distance in the circular Shepard tone space predicts the opposite of the percept. The various percepts can be readily captured and explained by a neural model of cortical activity based on populations of adapting, pitch and pitch-direction cells, aligned with the neurophysiological responses. Together, these decoding and model results suggest that contextual influences on perception may well be already encoded at the level of the primary sensory cortices, reflecting basic neural response properties commonly found in these areas.