Distribution of neurosensory progenitor pools during inner ear morphogenesis unveiled by cell lineage reconstruction
Abstract
Reconstructing the lineage of cells is central to understanding how the wide diversity of cell types develops. Here, we provide the neurosensory lineage reconstruction of a complex sensory organ, the inner ear, by imaging zebrafish embryos in vivo over an extended timespan, combining cell tracing and cell fate marker expression over time. We deliver the first dynamic map of early neuronal and sensory progenitor pools in the whole otic vesicle. It highlights the remodeling of the neuronal progenitor domain upon neuroblast delamination, and reveals that the order and place of neuroblasts’ delamination from the otic epithelium prefigure their position within the SAG. Sensory and non-sensory domains harbor different proliferative activity contributing distinctly to the overall growth of the structure. Therefore, the otic vesicle case exemplifies a generic morphogenetic process where spatial and temporal cues regulate cell fate and functional organization of the rudiment of the definitive organ.
https://doi.org/10.7554/eLife.22268.001eLife digest
Our ears, eyes and other sensory organs collect information about the world around us. In the inner ear – which is responsible for balance and hearing – specialized cells known as hair cells detect sounds and body position. This information is passed on to other cells called sensory neurons, which relay the information to the brain. Both of these cell types originate from a pool of cells known as “progenitors” located in a structure called the otic vesicle within the embryo, where the progenitor cells follow different sets of instructions to make hair cells or sensory neurons. Although we have identified many of the genes that are important for setting these instructions, it is not known how the progenitor cells behave in the inner ear or how they follow these guidelines.
Dyballa et al. used high-resolution imaging to reconstruct the histories of individual hair cells and sensory neurons in the inner ear of zebrafish embryos. The experiments combined techniques that allowed individual cells to be tracked over time and showed what types of cell they developed into. Dyballa et al. used these data to develop a map of the progenitor cells in the whole otic vesicle.
The findings help to fill the void in our understanding of the link between gene activity and tissue architecture in the inner ear. The next challenge is to use these findings as a basis to further explore models for how sensory neurons and hair cells develop, and to understand how to overcome the inner ear degenerati as we age.
https://doi.org/10.7554/eLife.22268.002Introduction
A major challenge in developmental biology is to explain how spatiotemporally controlled cell specification and differentiation occur alongside morphogenesis in the construction of functional organs. The inner ear is an attractive model to take on this challenge since it is accessible to manipulation, and it undergoes dynamic self-organization. It contains a manageable number of distinct cell types, which develop rapidly on an organized schedule to generate the functional units of the mature organ – the sensory patches. The key cell types of the inner ear, the supporting cells and the hair cells of the sensory patches, and the sensory neurons that innervate them, originate early during embryonic development from progenitors located in the otic vesicle, a 3D-structure arising from the otic placode adjacent to the developing hindbrain (Durruthy-Durruthy et al., 2014; Raft et al., 2007; Sapède et al., 2012; Satoh and Fekete, 2005), and they are easy to score by morphology, position and specific markers (Haddon and Lewis, 1996; Raft and Groves, 2015; Whitfield, 2015; Wu and Kelley, 2012). All these features have established the inner ear as a model widely used for the study of patterning and cell specification (Atkinson et al., 2015; Cai et al., 2015; Fritzsch et al., 2006a, 2006b; Whitfield and Hammond, 2007; Wu and Kelley, 2012).
Despite a good understanding about the molecular hierarchies, the knowledge of how individual progenitors behave throughout patterning, proliferation, and morphogenesis remains elusive. Neuronal vs. sensory specification is achieved through well-defined bHLH transcription factors: atoh1 for hair cell formation (Millimaki et al., 2007; Bermingham et al., 1999), neurog1 for sensory neuron determination (Andermann et al., 2002; Ma et al., 1998), and neuroD for sensory neuron differentiation and survival (Jahan et al., 2010; Kim et al., 2001). Signals arising in the surrounding tissues regionalize the otic vesicle along axes (Maier and Whitfield, 2014; Radosevic et al., 2011; Riccomagno et al., 2002, 2005; Sapède and Pujades, 2010), and this multiple step process implies a gradual restriction of cell fates over time (Whitfield and Hammond, 2007; Wu and Kelley, 2012). However, the phenotypes of targeted mutants for these signaling pathways are not always easy to reconcile (Raft and Groves, 2015), due to the limited comprehension of how developmental gene regulatory networks are integrated. For this, cellular data are needed as it can address how patterns can be achieved while the cells proliferate and the tissue undergoes morphogenesis, which may affect cell positioning and exposure to signals, and therefore cell specification.
Recent developments in 4D-microscopy imaging and cell tracking tools permit now simultaneous measurements at high spatial-temporal coverage and resolution, and therefore the assessment of cell lineages and cell behaviors including displacements and proliferations (Amat et al., 2014; Blanpain and Simons, 2013; Faure et al., 2016; Keller, 2013; Li et al., 2015; Olivier et al., 2010; Truong et al., 2011). Thus, it is time to progress in filling the void between gene regulatory networks and tissue architecture. With this purpose, we reconstructed the otic neurosensory lineage and investigated their single cell behavior by using in vivo imaging technologies paired with image processing tools (Figure 1, Figure 1—figure supplement 1; Faure et al., 2016). Our dynamic analyses revealed some surprising results such as the impact of neuroblast delamination and otic vesicle morphogenesis on the size and shape of this progenitor domain, and further that place and order of neuroblast delamination foreshadow their position within the statoacoustic ganglion (SAG). The comparative map of neuronal and sensory progenitors in the whole otic vesicle allows understanding how their distribution changes over time, being largely segregated with a small region of putative overlap. These findings provide the cellular data helping to understand how gene regulatory networks may work during development, tissue degeneration and regeneration.
Results
Genetic requirements for specification of otic neuroblasts are rather clear (Andermann et al., 2002; Ma et al., 1998); however, little is known about the cellular mechanisms underlying neuroblast development. Specifically we want to understand (i) how the neuroblast progenitor population is altered upon delamination of cells, (ii) how delamination coordinates in space and time, and (iii) how delaminated cells arrange to form the ganglion. We addressed these questions by exploring in depth and detail a selected number of embryos and support our findings with less extensive analyses of additional embryos, which we provide in supplementary form. Data from selected experiments can be downloaded from URL: http://bioemergences.eu/eLife2016.
Expansion of the neuroblast delamination domain
In order to study the spatiotemporal dynamics of neuroblast delamination we in vivo imaged Tg[neuroD:GFP] embryos, which express GFP in neuronal progenitors just prior to and after delamination from the otic epithelium. We observed the first delamination events at 18 hpf in the most anterolateral region of the otic floor (Figure 1—figure supplement 2a–a’). This domain expanded towards middle and posterior regions with cells delaminating very close to the neural tube (Figure 1—figure supplement 2b–c´), as had been described by analysis of serial transverse sections from 22 hpf onwards (Haddon and Lewis, 1996). Segmentation of the epithelial neuroD-domain allowed us to illustrate the quick expansion of this territory from anterolateral to posteromedial (Figure 1b–d), and when comparing three different embryos this emerges as a common trend (Figure 1—figure supplement 2d). After the delamination domain is established many more cells delaminate from this territory, accumulating in the SAG just below the epithelium in close intimacy with the anterior lateral line ganglion (ALLg; Figure 1i, Video 1).
Neuroblast delamination implies that cells exit the otic epithelium basally and therefore undergo cell shape changes. Segmentation of individual cells -delineation of the cell contours by computational tools- when transiting from the epithelium to the SAG showed that the neuroblast cell body moves basally in less than one hour after undergoing apical constriction (see red cell in Video 2; Figure 1e–g’). Neighboring neuroblasts often delaminate consecutively to accumulate in the SAG, which quickly becomes an adjacent compact mass extending beneath the ventral floor of the otic vesicle along the anteroposterior axis by 22 hpf (Figure 1h–i; Video 1). To understand the origin of SAG-neuroblasts and how delamination coordinates in space and time we backtracked these neuroblasts to their progenitor state and followed their dynamics.
Order and place of neuroblasts delamination from the otic epithelium prefigure their position within the SAG
As early as 20 hpf the SAG rudiment already innervates the different maculae, and by 48 hpf the SAG is composed of two segregated neuronal populations that although displaying the same molecular signature innervate different sensory patches (Sapède and Pujades, 2010; Zecca et al., 2015). We were eager to explore the spatiotemporal features controlling these neuronal populations and to gather information on population dynamics and lineage relationships. We first investigated whether neuroblasts were primed to different SAG-populations by backtracking SAG-neuroblasts to their progenitor state, and assessing individual cell position and tissue movements over time systematically. For this, we used the Tg[cldnb:lynGFP]Tg[Isl3:GFP] embryo injected with H2B-mCherry (ID Dataset: 140426aX, Tables 1–2).
Analysis of the dynamics of this process revealed that once the first delamination events are observed, delamination occurs massively. The lineage tree shows that within six hours a big bulk of neuroblasts delaminates (see white, yellow and orange lines in Figure 2a corresponding to single delaminated neuroblasts), and many neuroblasts can divide after exiting (see lines branching in Figure 2a) as previously described (Vemaraju et al., 2012). This lineage analysis allowed ascribing delamination time (see hpf in x-axis in Figure 2a) and delamination position to each neuroblast. Then, we color-coded the reconstructed cell centers according to these two criteria and displayed them (i) early when in the otic vesicle (Figure 2b,d,f), and (ii) late when within the SAG (c-c’,e-e’,g-g’). The first observation was that neuroblasts delaminated randomly in space and time (see intermingled colors of the reconstructed cell centers in Figure 2b), but their organization within the SAG relies on the delamination order as delaminated cells aggregate laterally to the preexisting ganglion (Video 3), thereby generating a mediolateral (ML) growth pattern with early-delaminated cells (white) located more medially than late-delaminated ones (red) (Figure 2c–c´; see Figure 2—figure supplement 1a–b for additional analysis). While the inner ear grows and undergoes morphogenesis, the SAG becomes squeezed in the space between the ventromedial wall of the ear and the neural tube, such that its organization is gradually converted from ML to dorsoventral (DV), with early-delaminated sensory neurons located in close contact with the neural tube (Video 3; Figure 2c). Thus, the time of neuroblast delamination foreshadows the ML gradient of neuronal differentiation in the SAG.
Then, we analyzed the epithelial origin of two functionally distinct SAG populations: anterior, mainly involved in vestibular, and posterior responsible for vestibular and acoustic functions (Haddon and Lewis, 1996). We found that the relative position of progenitors foresees their final location in the SAG: neuroblasts from the anterior portion of the SAG delaminated from the anterolateral floor of the otic vesicle, while the ones of the posterior SAG derived from the posteromedial otic epithelium (Video 4). Next, we addressed whether the epithelial coordinates prefigure the position of neuroblasts in the SAG primordium. We observed that the anteroposterior (AP) coordinates of neuroblasts delamination (see differently colored epithelial cell center populations in Figure 2d) generally defined their relative AP position within the SAG (Figure 2e–e’; Video 5; see Figure 2—figure supplement 1c–d’ for additional analysis). These data show that progenitors delaminating within similar spatial regions along the AP axis maintain similar relative positions later on, and indeed the analysis of their trajectories indicates that cells, while in the otic epithelium, maintain their neighbor relationships. No such correlation was observed regarding the ML axis, because cells that were originally separated in the epithelium were found close in the SAG (Figure 2f–g’; Video 5), suggesting that the AP and time cues prevail. These findings establish a link between the place and order of neuroblast delamination and their neuronal (functional) identity within the SAG, and suggest the existence of a spatial and temporal regulation in the otic epithelium of the fate and exit of SAG neuroblasts.
Remodeling of the neuronal progenitor domain upon delamination
We were interested in understanding whether cell proliferation had any role in conferring neuroblasts categories. To investigate this question we used the previous dataset (ID Dataset: 140426aX, Tables 1–2) to undertake the neuroblast clonal analysis (Figure 3a), which showed that SAG-neuroblasts originate from a neuronal progenitor pool with different division behaviors: one third of the analyzed neuroblasts divided before exiting the epithelium (n = 42 red cells in Figure 3a, Figure 3—figure supplement 1c), one third divided after delamination (n = 40 orange cells in Figure 3a; Figure 3—figure supplement 1c), and a similar fraction of neuroblasts did not divide within this time interval (n = 34 blue cells in Figure 3—figure supplement 1c). Thus, cell position within the otic epithelium is not relevant for these cell division behaviors because no specific spatial distribution can be observed (Figure 3—figure supplement 1c). In addition, cell division and delamination are independent events, in contrast to what was proposed for epithelial invagination of the Drosophila tracheal placode (Kondo and Hayashi, 2013).
Neuroblasts dividing before delamination give rise to two daughter cells that will delaminate (Figure 3a–c). In most cases, sister cell delaminations happened within a two hours interval, however in a few cases the time of delamination between two sister cells can be as long as seven hours (Figure 3b). Neuroblasts giving rise to two daughter cells that delaminated within a delay interval of more than two hours were evenly distributed in the otic epithelium (n = 11/116 bicolored cell centers, Figure 3—figure supplement 1d). Overall, these analyses suggest that (i) epithelial neuroblasts’ divisions give rise to cells that will delaminate (Figure 3c), and (ii) cell division is not necessary as a trigger for delamination.
Finally, we used the same dataset (ID Dataset: 140426aX; Tables 1–2) to understand how the neuronal progenitor domain changes upon cell delamination and morphogenesis of the inner ear over time (orange circles in Figure 3d; n = 131 at 18 hpf). We followed as well the behavior of the surrounding otic epithelial cells as a repair for the edge of the progenitor domain (see grey circles in Figure 3d, n = 64 at 26 hpf). Due to the massive cell delamination in a relatively short time (more than 150 delamination events in the period 18-26 hpf), the neuronal progenitor domain undergoes dramatic size and shape changes (Figure 3d, compare the territory containing orange cell circles at 18 hpf and 26 hpf), as described by Haddon and Lewis (1996). Upon delamination a large amount of cells exits from the otic epithelium: about 25% of total otic cells at 18 hpf are lost by delamination between 18 hpf and 22 hpf. The exit of neuroblasts and the growth of the vesicle are part of the morphogenetic processes that progressively restricted the progenitor field to a more lateral region of the ventral floor of the vesicle (orange cell centers in Figure 3d, 26 hpf).
Spatiotemporal development of the sensory patches
Hair cells are key cell types for building a functional sensory patch, and they are specified in the otic vesicle in the same way in all vertebrates. However, one difference is that in the fish the first hair cells differentiate while neuronal precursors are still delaminating, and that hair cells continue to be produced throughout life (Haddon and Lewis, 1996). This raises the question of how neuronal and hair cell production is coordinated from progenitor domains located very close within the otic epithelium, and this led us to compare the progenitor maps for hair cells and neurons.
Up to date, the formation of the anterior (AM) and posterior (PM) maculae (the first sensory patches to arise) has been followed up by gene or transgene expression mainly in 2D, and therefore their proper allocation within the 3D-otocyst was difficult to assess. In order to understand how these sensory patches arose we followed in vivo the incorporation of newly differentiated hair cells. The first hair cells to arise -the so called tether cells (Riley et al., 1997)- constituted the posterior pole of each of the maculae (white-dotted cells in Figure 4a–c; Figure 4—figure supplement 1b,c) and prefigured the position of these sensory patches within the otic vesicle (Video 6). Next, these two sensory patches gradually increased their size incorporating differentiated hair cells at their anterior pole with a specific pattern (see color-coded cell centers in Figure 4b–c; Video 6; ID Dataset: 140507aX; and see Figure 4—figure supplement 1b,c for additional samples exhibiting the same pattern). At the same time, growth and morphogenesis of the vesicle led to a structure in which the AM remained anterior and ventral, and the PM positioned in the posterior ventral edge of the medial wall of the vesicle. Once hair cells differentiate they become postmitotic (Video 6). Thus, the formation of the sensory patches depends on a pool of progenitors providing postmitotic differentiated hair cells. AM and PM develop asynchronously, with the incorporation of new hair cells being delayed in the PM (Figure 4—figure supplement 1; Sapède and Pujades, 2010).
Previous findings indicated that a specific pool of neuronal progenitors switches its fate to hair cells of the PM upon abrogation of neurog1 function (Sapède et al., 2012). However, they were unable to clearly assess the position of this progenitor pool and how it behaved. To do so, we first undertook the dynamical analysis of maculae generation upon neurog1-downregulation by injection of translation-blocking morpholino (ID Dataset: 140519aX, Tables 1–2), which fully recapitulates the phenotype of neurog1hi1059 mutants (Sapède et al., 2012). In the absence of neurog1, sensory neurons do not form and supernumerary hair cells are produced in the PM (Figure 4f–h; Video 6). Tracking the differentiated hair cells allowed us to compare the dynamics of sensory development at the single-cell level: the overall spatial pattern of hair cell generation in the AM/PM was very similar between control and morphant embryos (compare Figure 4b–c with g–h); however the kinetics of hair cell production in the PM differed due to a boost of differentiation in MO-neurog1 from 34 hpf onwards (Figure 4—figure supplement 1a). All supernumerary hair cells in the MO-neurog1 originate from newly differentiated cells, and not from differentiated hair cells re-entering the cell cycle (Video 6).
To unveil the organization of the PM hair cell progenitor pool we generated the map of hair cell progenitors in the whole otic vesicle over time (ID Datasets: 140507aX, 140326aX). For this, differentiated hair cells at 42 hpf were backtracked to their progenitor state, and information about differentiation time (color-code) and progenitor distribution was combined (Figure 4d–e; Figure 4—figure supplement 1c). 3D-reconstructions revealed that progenitors for hair cells of the PM were distributed over the ventromedial domain of the otic vesicle at early stages of embryonic development (see color-coded cell centers in Figure 4d and Figure 4—figure supplement 1c), and that upon neurog1 inhibition the hair cell progenitor pool expanded along the anteroposterior dimension allocating more medially (ID Dataset: 140519aX; Figure 4i, Video 7). Furthermore, the position of hair cell progenitors foreshadows the organization of differentiated cells within the sensory patch, demonstrating that otic epithelial cells did not rearrange during early stages of hair cell differentiation. However, the relative positions of the maculae changed upon morphogenesis, resulting in the growth of the PM towards the anterior region while growth of the AM was mainly anterior and towards dorsal (Figure 4c,h). Finally, our analysis allowed us to compare the progenitor maps for hair cells and neurons at the onset of hair cell differentiation (24 hpf): they were largely segregated with a small region of putative overlap (Video 8), consistent with the existence of a pool of dual progenitors (Sapède et al., 2012).
Sensory and non-sensory domains of the otic vesicle display different proliferative activity
As overall shape and size of the otic vesicle remain relatively constant during neuroblast delamination, we assumed that there was some homeostatic process involving cell proliferation elsewhere and we further validated this hypothesis.
To understand the proliferative behavior of the distinct otic territories, we reconstructed during 16 hr the lineages of 51–64 neighboring cells located either in the non-sensory (dorsolateral region of the otic vesicle; Figure 5a–c, Video 9) or in the sensory domain (ventromedial epithelial region; Figure 5d–f, Video 9). For this we used the Tg[Brn3c:GFP] H2B-mCherry injected embryo (ID Dataset: 140507aX). Cell behavior was assessed according to cell division (dividing/non-dividing) or cell differentiation (progenitor/differentiated) status. Cell position was monitored over time to identify spatial organization related to these features (Video 9). Cells within the non-sensory domain actively proliferated (Figure 5a–c; Video 9); in contrast, half of the analyzed cells in the ventral region did not divide and some differentiated into hair cells (Figure 5d–f; Video 9), suggesting that once hair cell progenitors are committed they divide less. Interestingly, even though the non-sensory domain displayed higher proliferative activity, it was less compacted (note neighboring reconstructed cell centers are more spaced in Figure 5a–b than in d–e), supporting the idea that it contributes to the overall growth of the vesicle during this time window. To quantify this, we calculated the nearest neighbor distance (NN-distance) for all cells in these domains over time and plotted the reconstructed cell centers color-coded for their NN-distance in the whole vesicle (see that green cell centers are more compacted than blue ones in Figure 5—figure supplement 1a; see Figure 5—figure supplement 1b for analyses of additional specimens).
Similar cell behaviors and cell organization were observed in the MO-neurog1 with no increase in proliferative activity in the sensory domain (Figure 5g–l), revealing that upon neurog1 abrogation hair cell progenitors do not divide more actively, and therefore supernumerary hair cells derive either from expansion of the progenitor pool, or from precocious differentiation of progenitors.
Since in MO-neurog1 cells do not undergo delamination, we determined how these non-delaminated cells located within the epithelium. For this, we estimated the distribution of local cell densities over the entire otic epithelium (Figure 5m–n; Video 10). While upon neurog1 abrogation the number of cells was greatly increased (Figure 5o; control = 384 ± 17 cells n = 3 vs. neurog1hi5109/hi5109 = 479 ± 14 cells n = 3 or MO-neurog1 = 506 ± 14 cells n = 2), the global pattern of cell densities did not change: cells in either condition are most densely packed in the anterolateral and ventral region and most spaced in the dorsal domain (see red and blue cells in Figure 5m–n; Video 10). This suggests that the impact of non-delaminated cells in the compaction of the otic epithelium is low. Thus, if in the MO-neurog1 cell density is not altered but cell numbers are higher in the vesicle, the volume of the otic structure should be larger than the control one. Indeed, when assessing the volumes by 3D-point-cloud segmentation the otic vesicles of MO-neurog1 were larger than the control ones (Figure 5—figure supplement 1c–e). Altogether these results show that the otic epithelium is a heterogeneous tissue, where proliferative activity and cell compaction differ between sensory and non-sensory domains during the generation of the neurosensory cellular elements.
Discussion
We provide information about cellular/population dynamics and lineage relationships of neurosensory elements in the inner ear from the reconstruction of their lineage trees from 4D in vivo data. These results enable us to: (i) understand the proportions of the system, (ii) reveal the impact of morphogenesis in the spatiotemporal distribution of neurosensory cell progenitor pools, and (iii) provide the cellular data to complement the well described gene regulatory networks involved in neurosensory specification. Our strategy has been to analyze in depth and detail few selected embryos, and support our findings with lower power analyses of other specimens. This analysis of progenitor populations in their native environment over extended timespans provides the opportunity to understand in vivo cell behaviors, both at the single cell level and at the cell population level. However, due to the technically demanding nature of the experiments and analyses the sample size could not be very high. All biological systems have intrinsic noise, and therefore variability: absolute cell numbers may differ between embryos, the onset of developmental processes might occur at slightly different times, or even distinct transgenic backgrounds may behave different. Therefore, to overcome these possible biases we analyzed additional specimens with lower power analysis and demonstrate that our approach can yield reproducible results in terms of cell type proportions and cellular behaviors.
To assure the robust outcome of a functional inner ear, with sensory patches containing the precise number of mechanosensory hair cells properly innervated by sensory neurons arranged in a precise topology (Video 11), the developmental strategy used by distinct progenitor populations differs: neuronal specification is concomitant with proliferation (before/after delamination), while hair cell specification and differentiation lead to postmitotic cells indicating that the final number of sensory cells relies on the control of the progenitor pool. Interestingly, the position of hair cell progenitors foreshadows the organization of differentiated cells within the sensory patch, demonstrating that although the otic vesicle undergoes morphogenesis and cells do extensively proliferate, progenitors do not rearrange within the epithelium during early stages of hair cell differentiation.
The inner ear comprises two distinct functional modalities –vestibular and acoustic- carried out by different SAG populations. However, their epithelial origin and how SAG neuroblasts acquire their specific identity are still open questions. Here we show the importance of place and time of neuroblast delamination in their allocation within the SAG, and shed light on how distinct information may converge in the progenitor cells. Delamination place confers position along the AP axis of the SAG, and therefore function, most probably as the result of the integration of patterning signals involved in the emergence of the different domains (Maier and Whitfield, 2014; Radosevic et al., 2011). Additionally, delamination time prefigures the relative position of early-delaminated neuroblasts vs. late-delaminated ones and therefore, the gradient of neuronal differentiation within the SAG. Consequently, an epithelial neuroblast will need to integrate spatial (its position in the epithelium) and temporal (its time of delamination) information for its final allocation within the SAG. However, to achieve the fully functional organ, complexity needs to increase later on possibly by some new delamination events, division of SAG neuroblasts and changes in morphology. Given that auditory neurons, which innervate the PM, accumulate later than the vestibular neurons (Sapède and Pujades, 2010), it is possible that the majority of neuroblasts becoming auditory neurons may delaminate from posterior domains, since no specific increase in cell proliferation has been observed at later stages in the posterior part of the SAG (Vemaraju et al., 2012). However, if this is the mechanism and how it is controlled remain to be explored.
The knowledge of how the gradient of neuronal differentiation within the SAG is built up helps the comprehension of its differential gene expression pattern according to the neuronal differentiation state (Zecca et al., 2015), and the organized distribution of differentiated neurons (Vemaraju et al., 2012). Thus, our results provide (i) a model to discuss how the selection of different progenitors and the determination of their relative population size could be regulated, which would not be possible to derive by gene expression only, and (ii) a framework to challenge the system –both in silico and in vivo- upon modification of key factors for cell fate decisions and patterning.
Interesting to note is that during a relatively short period of time, 18-30 hpf, a big proportion of epithelial neuroblasts delaminate. This drives a dramatic remodeling of the neurosensory progenitor domain that may lead to changes in the exposure of progenitor domains to sources of signals, such as Shh from the floor plate and the ventral neural tube. This loss of epithelial cells is probably compensated for by the active proliferation in the dorsal non-sensory territory in order to keep the homeostasis of the organ.
Finally, we provide the dynamic map of neurosensory progenitors based on in vivo cell lineage studies supplying a global and temporal perspective to previous otic neurosensory lineage analyses, which were mainly focused on the spatial dynamics of gene expression (Durruthy-Durruthy et al., 2014). This permits targeting and challenging progenitor pools specifically. Therefore, our findings establish a ground to further explore intrinsic vs. extrinsic models for cell fate determination, and will contribute to the mechanistic understanding of the developmental gene regulatory networks.
Materials and methods
Zebrafish strains
Request a detailed protocolZebrafish embryos were obtained by mating of adult fish using standard methods. All fish strains were maintained individually as inbred lines. All protocols used have been approved by the Institutional Animal Care and Use Ethic Committee (PRBB–IACUEC), and implemented according to national and European regulations. All experiments were carried out in accordance with the principles of the 3Rs. Wild type zebrafish strain was AB/Tu (RRID:ZIRC_ZL1/RRID:ZIRC_ZL57). Tg[neuroD:GFP] expresses GFP in neuronal progenitors (Obholzer et al., 2008) and differentiating neuroblasts (Zecca et al., 2015), and Tg[Isl3:GFP] (also called Isl2b) expresses GFP in afferent sensory neurons of cranial ganglia (Pittman et al., 2008). Tg[cldnb:lynGFP] labels the plasma membranes of the otic and anterior lateral line structures (Haas and Gilmour, 2006) and allows to visualize the kinocilium of differentiated hair cells. Tg[Brn3c:GFP] expresses GFP in differentiated hair cells of the ear and lateral line system, coinciding with the onset of differentiation (Xiao et al., 2005). Embryos homozygous for the neurog1hi1059 mutation (Golling et al., 2002) in the Tg[Isl3:GFP] background were obtained by incross of heterozygous carriers (Sapède et al., 2012); the presence of the neurog1hi1059 allele was identified by PCR genotyping fin-clips or embryo tails genomic DNA.
Phalloidin and DAPI staining for assessing the total number of cells in otic vesicles
Request a detailed protocolTg[neurog1hi1059/+ Isl3:GFP] fish were crossed and embryos fixed overnight at 4°C in 4% Paraformaldehyde (PFA) and washed with 0.1% Tween-20 in PBS (PBST) at room temperature (RT). To label the plasma membranes, embryos were incubated for 12 hr at RT with Phalloidin-Alexa658 (Invitrogen, Carlsbad, CA) diluted 1:20 in PBST containing 1.5% Triton X-100. They were washed with PBST and incubated with DAPI to stain nuclei. Embryos were then used for the analysis of the total number of cells within the otic vesicle (Figure 5o).
Antisense morpholinos and expression constructs
Request a detailed protocolEmbryos were injected at 1cell-stage with translation-blocking morpholino oligomers (MOs) obtained from GeneTools, Inc.: neurog1-MO 5’ -ACG ATC TCC ATT GTT GAT AAC CTG G-3’ (Cornell and Eisen, 2002). 5 ng of MO-neurog1 was injected as previously described (Sapède et al., 2012). Morphants display the same phenotype as neurog1hi1059/hi1059 mutants (Sapède et al., 2012). For in vivo imaging morphant embryos were used for practical reasons: it is not possible to phenotype the mutants as early as the stage we start the imaging (24 hpf), and since neurog1hi5109/hi5109 need to be obtained by incross of neurog1hi5109/+ fish because the line cannot be maintained in homozygosis, only 25% of embryos are neurog1hi5109/hi5109 with the inconveniences that this poses for life imaging.
For mRNA expression, capped H2B-mcherry, H2B-cerulean (Olivier et al., 2010) or lyn-TdTomato mRNAs were synthetized with mMessage mMachine kit (Ambion). Embryos were injected at 1cell-stage and let to develop until the desired stages.
Confocal imaging of stained whole mount samples
Request a detailed protocolStained fixed samples were mounted in 1% Low Melting Point (LMP)-agarose on glass-bottom Petri dishes (Mattek) and imaged on a Leica SP5 inverted confocal microscope using a 20x objective (NA 0.7).
Embryo mounting and SPIM imaging
Request a detailed protocolEmbryos were anesthetized in 0.17 mg/ml tricaine in system water and mounted in 0.75% LMP-agarose in glass capillaries size 2 (volume 20 µl, BRAND GMBH). Time-lapse imaging was performed at 26.5°C (to avoid melting of LMP-agarose) on a Zeiss Lightsheet Z.1 microscope using a 20x or a 40x objective and the developmental stage was corrected accordingly (at 26.5°C development is delayed about 0.7 fold). Nuclei, plasma membranes and cell fate were recorded simultaneously for the entire system (Figure 1a, Figure 1—figure supplement 1). The cohort of embryos and datasets used in this study are depicted in Tables 1 and 2. Each dataset corresponds to the imaging of a distinct embryo inner ear, except for 140402aX in which both ears were imaged. Note that we have performed deep and detailed analyses in few datasets, and supported the observations and conclusions with partial analyses of additional datasets mainly included in supplementary figures.
3D+time image analysis pipeline
Image pre-processing (Figure 1—figure supplement 1)
Request a detailed protocolImage pre-processing was done using the Zeiss ZEN software and involved dual (illumination) side fusion and deconvolution (Regularized Inverse Method). To compensate for morphogenetic movements during image acquisition a semi-automated rigid registration was carried out using developed FIJI-scripts: the user is guided through the time steps of the 16-bit. czi data set and chooses by clicking a fixpoint of the structure (eg. a tether cell or a nucleus on the dorsolateral wall of the otic vesicle) for which the x, y, z coordinates are recorded. After navigating through all time steps of the video and having specified the fixpoints for all time steps the dataset is processed according to the recorded coordinates. This registration allows as well decreasing the size of region of interest, which substantially reduces data size. A separate 8-bit. vtk file with the corresponding transformation is generated for each time step and channel. Changes on the levels of the transgene fluorescent protein expression were compensated for by modifying dynamic ranges upon 16-bit to 8-bit mapping. Different datasets were generated for this work as displayed in Tables 1–2.
Center detection and automatic tracking (Figure 1—figure supplement 1)
Request a detailed protocolThe .vtk files were then uploaded to the Bioemergences platform (Faure et al., 2016) and center detection was launched. Detected centers were validated using the CenterSelect application available in the platform. Subsequently, automatic tracking of the validated centers was launched and once completed the tracking data was validated and curated using MovIT, a custom-made graphical interface that offers the tools for segmentation and tracking of cells to accurately reconstruct their positions, movements and divisions (Faure et al., 2016). Videos of the developmental processes displayed in the manuscript were generated using FIJI (RRID:SCR_002285), MovIT, or a combination of both. Only validated/curated cell tracks were used for further analyses. MovIT has several display modes of time-lapse data such as: (i) orthoslices through the volume showing one xy-plane, one xz-plane and one yz-plane, used to validate cell tracks, (ii) oblique volume slices, which show a slice of a density map that is not parallel to the faces of the volume box and allows orienting the image along the axes, and (iii) volume rendering allowing to display a 2D projection of a 3D discretely sampled dataset; this permits to display cell centers in the context of the raw images in the whole volume (Video 1).
Datasets analyses (Figure 1—figure supplement 1, Tables 1–2)
Request a detailed protocolAll used datasets are clearly stated in the Figure legends, and mostly throughout the text. For the comparison of the posterior expansion of the delamination domain in several specimens, the relative position of delamination events along the AP axis of the otic vesicle was analyzed: the position of posterior-most neuroD:GFP expressing cells in the otic epithelium, and the position of the anterior and the posterior edges of the otic vesicle (dorsal view) were recorded over time (Figure 1—figure supplement 2d). Then the posterior edges of the delamination domains were plotted in percent of the AP otic vesicle length as a function of time.
For manual segmentation of the delamination domain, ITK-Snap software was used on three time steps of an imaging sequence of Tg[neuroD:GFP] embryos injected with lyn-TdTomato mRNA (Figure 1b–d). The delamination domain (identified by otic epithelial neuroD:GFP expression) and the otic vesicle volume (determined from the lyn-TdTomato signal) were segmented, and the resulting. vtk files were used to display neuroD-epithelial neuroblasts and membrane signal of the otic epithelium only (FIJI-3D viewer).
For cell lineage analysis, validated centers were grouped in selections and displayed with different appearances either alone or in context of the original data (Figure 1h–i; Figure 2b–g’; Figure 2—figure supplement 1b–d’; Figure 3—figure supplement 1; Figure 4a–c,f–h; Figure 4—figure supplement 1b–c; Figure 5a–b,d–e,g–h,j–k; Figure 5—figure supplement 1a). Cell lineages were also displayed as lineage trees (Figure 2a, Figure 3a, Figure 2—figure supplement 1a).
For analysis of cell proliferative activities over time, cell centers were grouped in selections according to cell states (eg. dividing or non dividing) and these data were further explored within SciPy, a Python-based ecosystem for scientific computing (Figure 5). Namely the data was imported as Pandas DataFrames; different time points and samples were registered using the iterative closest points transformation algorithm as implemented in. vtk and the results were plotted using matplotlib and seaborn.
To estimate the local cell density the other cell centers within 20 µm of each center were counted, and this number was then divided by the volume of a sphere of that radius (Figure 5m–n, Video 10). The nearest neighbor distance (NN-distance) was used instead when comparing subsets of cell centers to limit boundary effects (Figure 5—figure supplement 1a,b). Both estimates were calculated using k-d-trees as implemented in the SciPy library.
For studying cell shape changes upon delamination, single neuroblasts delaminating from the ventromedial aspect of the otic vesicle were chosen for automatic segmentation (n = 5; Figure 1e–g’, Video 2). Additionally to the reconstructed cell center that defines the lineage of the cell, two more cell centers were added to each cell at each time step to serve as seed points for segmentation. Automatic segmentation of these cells was then generated by the GSUBSURF method (Mikula et al., 2011), which is an image segmentation method based on solving a level set partial differential equation of the form
where is an edge detector function depending on the intensity function of the segmented grayscale image and and are positive real constants steering strength of advective velocity field and curvature regularization. An evolution process of this form converges to a steady state, and it is stopped once we do not observe a significant change in the function . The steady state of evolving level set function provides a segmentation result which can be directly displayed in FIJI together with the raw data or a selected isosurface of function can represent the surface of the cell. The equation is solved numerically on image voxel grid by the semi-implicit finite volume scheme, and the initial condition is constructed so that we have a reasonable initial approximation of the segmented cell (e.g. a union of spheres or ellipsoids centered in the provided cell centers).
To calculate the volume of the otic vesicle a semi-automatic strategy was applied (Figure 5—figure supplement 1c–e), because a fully automatic segmentation of the otic vesicle structure is a rather complicated task – it is not trivial to construct an appropriate initial condition and to distinguish the borders of the vesicle from the borders of the surrounding cells. A set of points (usually containing 100–300 points) lying on the surface of the vesicle was marked manually using FIJI; the inner and the outer surfaces were represented by separate point sets. Then, we used a method for the reconstruction of 3D objects from point clouds based on Lagrangian evolution of surfaces in 3D and the corresponding evolution model reads
Here, we directly evolve the surface represented by the map (parametrization) . The initial condition is an ellipsoid containing the given point cloud in its inside. The function is the distance function to the point cloud and it is the driving force of the evolution. The vector is a unit normal to the surface and is its mean curvature. The parameters and are positive real constants. The vector field is a tangential vector field specifically designed to control the quality of the surface discretization mesh during the evolution. The evolving surface is approximated by a triangulated surface. That means that at the end of the evolution, we directly obtain a triangular representation of the otic vesicle surfaces that can be easily used to compute the volume of the vesicle. As in the case of cell segmentation, the evolution model is discretized semi-implicitly in time and by a finite volume method in space.
References
-
Neurogenin1 defines zebrafish cranial sensory ganglia precursorsDevelopmental Biology 251:45–58.https://doi.org/10.1006/dbio.2002.0820
-
Unravelling stem cell dynamics by lineage tracingNature Reviews Molecular Cell Biology 14:489–502.https://doi.org/10.1038/nrm3625
-
The Role of Atonal Factors in Mechanosensory Cell Specification and FunctionMolecular Neurobiology 52:1315–1329.https://doi.org/10.1007/s12035-014-8925-0
-
Delta/Notch signaling promotes formation of zebrafish neural crest by repressing Neurogenin 1 functionDevelopment 129:2639–2648.
-
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
-
Neurod1 regulates survival and formation of connections in mouse ear and brainCell and Tissue Research 341:95–110.https://doi.org/10.1007/s00441-010-0984-6
-
NeuroD-null mice are deaf due to a severe loss of the inner ear sensory neurons during developmentDevelopment 128:417–426.
-
Segmentation of 3D cell membrane images by PDE methods and its applicationsComputers in Biology and Medicine 41:326–339.https://doi.org/10.1016/j.compbiomed.2011.03.010
-
Vesicular glutamate transporter 3 is required for synaptic transmission in zebrafish hair cellsJournal of Neuroscience 28:2110–2118.https://doi.org/10.1523/JNEUROSCI.5230-07.2008
-
Specification of the mammalian cochlea is dependent on Sonic hedgehogGenes & Development 16:2365–2378.https://doi.org/10.1101/gad.1013302
-
Development of the inner earCurrent Opinion in Genetics & Development 32:112–118.https://doi.org/10.1016/j.gde.2015.02.006
-
Axial patterning in the developing vertebrate inner earThe International Journal of Developmental Biology 51:507–520.https://doi.org/10.1387/ijdb.072380tw
-
Molecular mechanisms of inner ear developmentCold Spring Harbor Perspectives in Biology 4:a008409.https://doi.org/10.1101/cshperspect.a008409
Article and author information
Author details
Funding
Unidad de Excelencia María de Maetzu (2015-19 MDM-2014-0370 to DCEXS-UPF)
- Sylvia Dyballa
- Andrea Zecca
- Cristina Pujades
Becas de la Generalitat de Catalunya (Predoctoral FI-fellowship)
- Sylvia Dyballa
- Andrea Zecca
Agence Nationale de la Recherche (ANR-11-EQPX-0029)
- Thierry Savy
- Nadine Peyriéras
Agence Nationale de la Recherche (ANR-10-INBS-04)
- Thierry Savy
- Nadine Peyriéras
Schweizerischer Nationalfonds zur Förderung der Wissenschaftlichen Forschung (SINERGIA CRSII3 141918)
- Philipp Germann
Centro de Excelencia Severo Ochoa (2013-17 SEV-2012-0208 to CRG)
- Philipp Germann
Ministerio de Economía y Competitividad (BFU2012-31994)
- Cristina Pujades
Institució Catalana de Recerca i Estudis Avançats (ICREA Academia Award)
- Cristina Pujades
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We are grateful to T Pujol for his valuable help in SPIM imaging and Zeiss for use of the Lightsheet Z.1 microscope. We thank members of Peyriéras and Pujades laboratories for insightful discussions and technical help, especially D Fabreges and G Recher. We acknowledge P Daniel for help in experiments related to otic vesicle point cloud segmentation. We thank V Lecaudey, A Nechiporuk and B Riley for kindly providing transgenic fish lines, and S Fraser for his critical reading of the manuscript and valuable comments. This work was supported from Spanish Ministry of Economy and Competitiveness (MINECO) by Grant BFU2012-31994 to CP, Unidad de Excelencia María de Maetzu 2015–19 MDM-2014–0370 to DCEXS-UPF, Centro de Excelencia Severo Ochoa 2013–17, SEV-2012–0208 and Swiss National Science Foundation (SINERGIA CRSII3 141918) to CRG. BioEmergences services were funded by ANR-10-INBS-04 and ANR-11-EQPX-0029. SD and AZ were recipients of predoctoral FI-fellowships from AGAUR (Generalitat de Catalunya). CP is recipient of an ICREA Academia award (Generalitat de Catalunya).
Ethics
Animal experimentation: This study was performed in strict accordance with the European Regulations. Zebrafish embryos were obtained by mating of adult fish using standard methods. All fish strains were maintained individually as inbred lines. All protocols used have been approved by the Institutional Animal Care and Use Ethic Committee (PRBB-IACUEC), and implemented according to national and European regulations. All experiments were carried out in accordance with the principles of the 3Rs. All our experiments were carried out using the CPC16-008/9125 protocol approved by the Generalitat of Catalonia.
Copyright
© 2017, Dyballa 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
-
- 1,768
- views
-
- 399
- downloads
-
- 25
- 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
-
- Cell Biology
- Physics of Living Systems
The endoplasmic reticulum (ER), the largest cellular compartment, harbours the machinery for the biogenesis of secretory proteins and lipids, calcium storage/mobilisation, and detoxification. It is shaped as layered membranous sheets interconnected with a network of tubules extending throughout the cell. Understanding the influence of the ER morphology dynamics on molecular transport may offer clues to rationalising neuro-pathologies caused by ER morphogen mutations. It remains unclear, however, how the ER facilitates its intra-luminal mobility and homogenises its content. It has been recently proposed that intra-luminal transport may be enabled by active contractions of ER tubules. To surmount the barriers to empirical studies of the minuscule spatial and temporal scales relevant to ER nanofluidics, here we exploit the principles of viscous fluid dynamics to generate a theoretical physical model emulating in silico the content motion in actively contracting nanoscopic tubular networks. The computational model reveals the luminal particle speeds, and their impact in facilitating active transport, of the active contractile behaviour of the different ER components along various time–space parameters. The results of the model indicate that reproducing transport with velocities similar to those reported experimentally in single-particle tracking would require unrealistically high values of tubule contraction site length and rate. Considering further nanofluidic scenarios, we show that width contractions of the ER’s flat domains (perinuclear sheets) generate local flows with only a short-range effect on luminal transport. Only contractions of peripheral sheets can reproduce experimental measurements, provided they are able to contract fast enough.
-
- Cell Biology
Cell migration towards stiff substrates has been coined as durotaxis and implicated in development, wound healing, and cancer, where complex interplays between immune and non-immune cells are present. Compared to the emerging mechanisms underlying the strongly adhesive mesenchymal durotaxis, little is known about whether immune cells - migrating in amoeboid mode - could follow mechanical cues. Here, we develop an imaging-based confined migration device with a stiffness gradient. By tracking live cell trajectory and analyzing the directionality of T cells and neutrophils, we observe that amoeboid cells can durotax. We further delineate the underlying mechanism to involve non-muscle myosin IIA (NMIIA) polarization towards the soft-matrix-side but may not require differential actin flow up- or down-stiffness gradient. Using the protista Dictyostelium, we demonstrate the evolutionary conservation of amoeboid durotaxis. Finally, these experimental phenomena are theoretically captured by an active gel model capable of mechanosensing. Collectively, these results may shed new lights on immune surveillance and recently identified confined migration of cancer cells, within the mechanically inhomogeneous tumor microenvironment or the inflamed fibrotic tissues.