The brain, like most other organs, is formed by the coordinated growth of a few unspecialized cells in the embryo, which give rise to billions of neurons. For the brain to work properly, it is crucial that, during embryonic development, each neuron ends up in the correct location. This migration to the right spot has to happen while the brain grows and changes shape, which affects how and how far neurons and their precursor cells need to move to reach their final position. If these movements and changes in shape are not coordinated correctly, neurons can end up in the wrong place, form the wrong connections, and ultimately impact how the brain works.

Previous work done in fruit flies and zebrafish resulted in three-dimensional maps of these animals’ healthy brains, which allowed scientists to have a holistic view of how brains are organized. Although these maps are a valuable resource to study the structure of the brain, they do not provide information on how the brain transforms over time, especially during embryonic development. To get a clearer picture of how a few precursor cells give rise to the incredibly complex tissue that is the brain, a three-dimensional map spanning the entire developmental process is needed.

To fill this gap in knowledge, Blanc et al. developed a digital atlas-maker pipeline (DAMAKER) that allows scientists to generate three-dimensional models of the embryonic brain from microscopy images of several individuals. They then used this pipeline to construct a three-dimensional digital atlas of how a part of the brain called the hindbrain develops in the zebrafish embryo.

First, they collected images of the hindbrain showing neurons born at different times and matched these images to the existing static maps. Next, DAMAKER was used to follow neurons from the time of their birth to their final location, allowing Blanc et al. to create a map showing where neurons born at different stages during development end up. This type of map allows users to accurately determine when different populations of mature neurons are born, which allows scientists to estimate when different defects in brain development might originate.

Based on these data, Blanc et al. concluded that in zebrafish most of the cells that will end up forming the hindbrain acquire their specialized neuronal identities very early in development, between 24 and 48 hours post fertilization. These temporal maps of healthy hindbrains were then compared to maps of brains in which the birth of neurons was disrupted, thus changing the final number of neurons in the brain. This experiment showed that changing the number of neurons that are born early in development alters the final positions of neurons and the overall shape of the brain. Therefore, for the brain to grow to its correct size, there must be a balance between the number of unspecialized cells in the developing brain, and the rate at which these cells become neurons.

The DAMAKER pipeline not only provides scientists with a tool to study neurodevelopmental disorders, but also serves as a method that can be adjusted to map growth and shaping of other organs.

Generating a precise temporal cartography of cell populations in the three-dimensional (3D) embryonic brain is essential for understanding how brain morphogenesis impacts the position, and therefore function, of neuronal progenitors and neuronal clusters in normal and pathological conditions. Pioneer work in Drosophila provided 3D-registration atlases and software that can automatically find the corresponding landmarks in a subject brain, and map them to the coordinate system of a target brain (Heckscher et al., 2014; Peng et al., 2011). However, up to now, the incorporation of time —and therefore information of tissue morphogenesis— in 3D-atlases is still one of the main challenges that needs to be solved.

Digital 3D-atlases enable the exploration of brain structures in the 3D-context and can be updated to incorporate new information, which is essential for the exponential increase of knowledge. Existing gene expression 3D-atlases in larval and adult fish provide valuable tools (Chow et al., 2020; Jaggard et al., 2020; Kenney et al., 2021; Kunst et al., 2019; Ronneberger et al., 2012; Tabor et al., 2019). However, to incorporate the temporal component has been challenging, most probably because whole brain images are particularly hard to register across space and time, even though there are recent efforts aiming to make this more accessible (Dsilva et al., 2015; Fernandez and Moisy, 2020). Temporal digital 3D-atlases are essential to understand how the generation of brain cell diversity occurs concomitantly with brain morphogenesis, which results in a dramatic transformation from a simple tubular structure (e.g. the neural tube) to a highly convoluted structure (e.g. the brain). Thus, our goal was to provide an atlas allowing (i) the incorporation of time to enable morphogenesis studies, (ii) the quantification of 3D-patterns with a standardized method, (iii) the comparison of tissue shapes, and (iv) an easy access to the atlas data with user-friendly upgradability, without the need to send the data to the atlas-makers for processing. Such features of a digital atlas are crucial for studies that rely on 3D-topologies to understand how specific regions develop; moreover, to do this in zebrafish would enhance the use of zebrafish as human avatars in neurodevelopmental disorders by enabling to understand how specific territories respond to injury or gene disruptions.

Here, we developed a new ‘Digital Atlas-MAKER’ pipeline (DAMAKER) that fulfills the previous requirements and can be used to create an expandable atlas of virtually any tissue in any organism. Further, it is able to perform data analysis and quantification mostly on a single platform, FIJI, which is already established as the standard for image analysis in the field. The expandable zebrafish digital 3D-atlas was created by combining in vivo imaging data from distinct transgenic embryos with the Fijiyama registration tool (Fernandez and Moisy, 2020), which allows multi-modal rigid registration, and the Trainable Weka Segmentation (TWS) (Arganda-Carreras et al., 2017) that leverages a limited number of manual annotations and permits data to be automatically segmented after a short training and to quantify tissue volumes. We applied DAMAKER to the monitoring of the growth of the neuronal differentiation domain in the context of the whole zebrafish hindbrain —the embryonic brainstem—, to study how cell differentiation and morphogenesis might be intertwined. We made use of the hindbrain, which is the most conserved vertebrate brain vesicle (Murakami et al., 2005) and undergoes tissue segmentation along the anteroposterior (AP) axis leading to the formation of seven rhombomeres (Fraser et al., 1990; Jimenez-Guri et al., 2010; Kiecker and Lumsden, 2005; Krumlauf and Wilkinson, 2021; Pujades, 2020). This morphogenetic process is accompanied with the initial neuronal differentiation, which starts prior 24hours post-fertilization (hpf) with a pronounced increase from 30hpf onwards (Voltes et al., 2019). Neuronal differentiation involves changes in cell proliferative capacity and an extraordinary displacement of neurons from their birth site, and occurs while other dramatic changes take place, such as the generation of the brain ventricle (Belzunce et al., 2020; Gutzman et al., 2008; Hevia et al., 2022; Lyons et al., 2003). Although previous work evoked the importance of the temporal order of neuronal differentiation in ascribing final neuronal position within specific circuits (Kinkhabwala et al., 2011; Pujala and Koyama, 2019; Wan et al., 2019), little is known about its impact in the final organization of differentiated neurons within the hindbrain.

Here, we registered experimental whole zebrafish hindbrain images and mapped neuronal differentiation patterns with a unified spatial representation. We build up a digital 3D-hindbrain atlas in which signals corresponding to specific territories and differentiated neurons were merged. With these combinations, we could associate specific neuronal populations with anatomical landmarks and quantify complex volumetric data. Next, we performed a tissue wide neuronal birthdate analysis by using a method based in photoconverted fluorescent protein tracing in vivo (Caron et al., 2008), which permitted to distinguish early-born vs. late-born neuronal territories in live embryos. This allowed us to develop neuronal differentiation maps that revealed the remodeling of the neuronal progenitor and neuronal cluster domains over time. By generating these maps, we unveiled that early-born neurons were always located in the inner part of the differentiation domain and surrounded by younger neurons, generating an inner–outer gradient of older vs. younger neurons. As a proof-of-concept, we demonstrated that this temporal 3D-atlas could be used as a proxy to infer the birthdate of given neuronal populations, as well as to analyze neurogenic defects. We believe that this tool —the digital 3D-atlas maker, DAMAKER— will allow the integration of the information generated in other laboratories and will help ascribing neuronal birthdate and neuronal differentiation times upon gene disruption in zebrafish avatars for human diseases. This collaborative integration is essential for expanding our overall knowledge about how the brain is functionally organized, and therefore advancing brain research, medicine and brain-inspired information technology.

Our scientific interest was to monitor the dynamic growth of the neuronal differentiation domain in the context of the whole zebrafish hindbrain with the emphasis in how neuronal differentiation order or birthdate may prefigure cell position within the differentiation domain. For this, we developed a standardized and reliable 3D-imaging processing method —the digital 3D-atlas maker or DAMAKER— that allowed us to generate digital 3D-models, to quantify volumes and compare tissues shape, and most importantly for us to incorporate time.

Developing a standardized, reliable 3D-image processing and quantification method

Building upon previous work that reconstructed prototype 3D-atlases of whole embryos in Drosophila (Heckscher et al., 2014; Peng et al., 2011) and of larval or adult zebrafish (Chow et al., 2020; Jaggard et al., 2020; Kenney et al., 2021; Kunst et al., 2019; Ronneberger et al., 2012; Tabor et al., 2019), we aimed to create a digital 3D-atlas of the hindbrain neuronal differentiation domain that could incorporate time. To build it up, we designed a pipeline requiring the following steps: (i) registration of in vivo 3D-images displaying specific signals in order to align raw samples from different embryos (Figure 1A), (ii) segmentation of the 3D-images to generate labeling masks, each corresponding to a given marker or cell population (Figure 1B), and (iii) processing these masks to create average digital 3D-models (Figure 1C). In addition, the pipeline should allow to isolate specific signals from z-stack images in a semi-automatic manner in order to (i) quantify them along any body axes (Figure 1D) or calculate signal volumes upon time (Figure E), (ii) generate a digital 3D-atlas by intersection of several reference 3D-models (Figure 1F), and (iii) perform tissue shape comparisons (Figure 1G).

Figure 1 with 1 supplement see all

Download asset Open asset

Since we were interested in neuronal differentiation patterns, to establish a reference tissue signal we made use of Tg[HuC:GFP] embryos at 72hpf, a stage in which most of the cells have undergone neuronal differentiation and GFP is uniformly displayed in the whole neuronal differentiation domain (Hevia et al., 2022). Acquired 3D-images from different Tg[HuC:GFP] embryos (e1–eY; n=5each) were registered to the reference by Fijiyama registration, in which rigid registration was favored, although non-rigid registration could be considered at the cost of sample shape variability. Then, the obtained transformation matrix could be applied to other signals, such as those labeling specific tissue domains, cell populations, etc (Figure 1A). After alignment of 3D-images, tissue segmentation was applied in order to generate binary masks for each labeling/staining (Figure 1B). Here, segmentation settings were at tissue scale, both in an effort to keep computing power low for pipeline accessibility and for our interest in cell populations, although single cell resolution could be achieved with nuclear reporter lines. In order to discern noise vs. signal in the raw data from each embryo, different users (u1–uX; n=5) were invited to employ the Trainable Weka Segmentation (TWS) tool to train classifiers (one per user) (Figure 1B). In the case of low or sparse signals (or specific structures within a signal) specific training of the algorithm with different settings would be better to leverage the structural recognition capabilities of TWS. To reduce the potential impact of the human bias from the trained classifiers, we averaged the resulting segmented images from multiple users’ trained classifiers (Figure 1B). In the case of having an outlier user, this classifier was not considered (as example, see u4 in Figure 1B as an outlier that was discarded). Accordingly, we generated as many binary 3D-masks as users (n=4; Figure 1B). The users’ averaged embryo 3D-masks displayed a gradient, which represented the spatial similarity between users’ segmentations. We selected the signals agreed by at least three users (Figure 1B), although this could be differently thresholded to display signal consensus among users, with the strength of the threshold setting up the diversity. For image processing, these embryo 3D-masks were averaged to merge information from all embryos (e1–eY; n=5) to obtain the ‘label 3D-model’. As many label 3D-models as desired signals were generated, by reiteration of the whole process (Figure 1C). Each digital label 3D-model displayed a gradient, which represented the spatial similarity between embryos. To establish the most accurate consensus selection, we compared label 3D-models thresholded for different levels of consensus between embryos. We selected the signals displayed by at least 3 embryos, compared the distribution of the labeling between embryos, and built up the digital reference 3D-atlas (Figure 1C). To evaluate the accuracy of the registration, we assessed (i) the interindividual variability of different signals (Figure 1—figure supplement 1A–E,A’–E’), (ii) the size of the neuronal differentiation territory (HuC-domain) in different transgenic backgrounds (Figure 1—figure supplement 1F), and (iii) quantified the interindividual variability between the HuC volume of the individuals and the HuC volume of the generated 3D-model (Figure 1—figure supplement 1G). Thus, this pipeline allowed us to (i) assess signal similarity between segmentations limiting the bias due to the users’ subjectivity, (ii) obtain reliable 3D-tissue quantifications, (iii) produce a representative 3D-model of a label, (iv) estimate interindividual variability, and (v) build up a reference 3D-brain to be further used for spatiotemporal comparisons.

Building a multilayered digital 3D-atlas

Next, we wanted to add multiple layers to the neuronal differentiation 3D-atlas. For this, we made use of double transgenic embryos displaying a fluorescence reporter in the neuronal differentiation domain that was used as a reference for alignment in one channel, and another fluorescence reporter of the specific cell population or territory of interest (Figure 2). In addition, we made use of the rhombomeres (r) 3 and 5 as anatomical landmarks to have a positional coordinate system on which to align the quantitative representation of the given signal distribution across the tissue. Thus, we imaged Tg[Mu4127;HuC:GFP] embryos, which displayed DsRed in r3 and r5 (Distel et al., 2009), and GFP in the neuronal differentiation domain (Figure 2A–A’). Upon rendering the processed 3D-model of the signal, we could observe the 3D-view of the neuronal differentiation r3 and r5 volumes (Figure 2—video 1A). To monitor differences in the specific signal spatial distribution provided by Mu4127 within the neuronal differentiation territory, we obtained the volumes’ value across body axes by taking sample segmentations, created ‘virtual slices’ along axes and quantified the signal ratio at each ‘slice’ (Figure 1D; Figure 2B). Accordingly, the signal distribution was concentrated in r3 and r5 along the anteroposterior (AP) axis, whereas it was similarly distributed along the mediolateral (ML) axis (Figure 2B).

Figure 2 with 2 supplements see all

Download asset Open asset

Next, we investigated the distribution of the differentiated neurons deriving from specific neuronal committed progenitors, according to their proneural gene expression (Guillemot, 2007), in this case neurog1. In Tg[neurog1:DsRed;HuC:GFP] embryos, we observed the striped-pattern distribution of neurog1-derivatives in the neuronal differentiated domain (Figure 2C–C’; Figure 2—video 1B). The automatic quantification of the neurog1-volume displayed differences in the ratio of labeling along the AP axis (Figure 2D), supporting the observation of neuronal differentiated regions with less neurog1-derivatives when compared with more posterior domains (Figure 2C–C’, see white asterisk in the lateral r5). In addition, the more medial distribution of neurog1-derivatives in the hindbrain was observed (Figure 2D) when compared with the spinal cord as it was previously described (Belzunce et al., 2020). Thus, this quantitative analysis allowed to associate signal distribution to specific territories or cell populations. For mapping distinct neuronal differentiation programs, we first analyzed the spatial distribution of motoneurons by using Tg[isl1:GFP;HuC:Kaede^Red] embryos and observed their ventral and medial position (Figure 2E–E’; Figure 2—video 2A). Upon quantification, higher signal was observed in discrete patterns along the AP axis with a clear increase in the spinal cord (Figure 2F). We explored the specific neurotransmitter-expressing cell populations, considering the mainly expressed markers of excitatory (glutamatergic neurons) and inhibitory populations (GABAergic neurons) in the zebrafish hindbrain (Higashijima et al., 2004). Using Tg[vglut2a:DsRed;HuC:GFP] and Tg[gad1b:DsRed;HuC:GFP] embryos we observed the main contribution of glutamatergic neurons to the most dorsolateral neuronal differentiation domain (Figure 2G–G’), whereas GABAergic neurons allocated more medially (Figure 2I–I’). We built up the digital 3D-maps (Figure 2—video 2B–C) and quantified the distribution along the axes. We observed a smaller contribution of glutamatergic neurons to the r3-r5 region (Figure 2H), whereas GABAergic signal at this stage was enriched in r2 (Figure 2J) when compared to the rest of the hindbrain. This quantification recapitulated the spatially restricted domains of neurotransmitter expression (Figure 2G’ and I’) previously described (Kinkhabwala et al., 2011; Koyama et al., 2011; Pujala and Koyama, 2019). The overlay of the distinct digital 3D-maps provided information about the segregation of populations (Figure 2—video 2D). Overall, the expression analysis of fluorescent reporter transgenic lines allows the virtual colocalization of signals, and therefore the model intersection in order to study specific subsets of neuronal populations and their spatial relationships.

Building the digital 3D-atlas of neuronal birthdate

Previous work demonstrated the importance of neuronal birthdate in ascribing neuronal position and function. The stripe patterning of alx-expressing neurons, which are involved in swimming behavior, is age-related (Kinkhabwala et al., 2011), and the position of the V2a cell body and the order of spinal projections in the hindbrain relies on their neuronal birth (Pujala and Koyama, 2019). Thus, we wanted to explore whether neuronal birthdate was a general rule that dictated the spatial distribution of differentiated neurons within the hindbrain.

To build birthdating 3D-maps of neuronal differentiation, we performed spatiotemporal analyses of neuron birthdating using Kaede photoconversion experiments as a tool to register temporality (Caron et al., 2008). In this case, Kaede^Green in Tg[HuC:Kaede] embryos was in vivo photoconverted to Kaede^Red at different developmental stages (24, 36 or 48hpf), and embryos were imaged either at 48hpf or 72hpf to assess the relative position of neurons born at different times within the differentiation domain (Figure 3; earlier-born neurons Kaede^Red-cells vs. older-born as Kaede^Green-cells). To generate the birthdate 3D-map at 48hpf, we first photoconverted the HuC:Kaede^Green at 24hpf and imaged embryos at 48hpf. We observed that the onset of neuronal differentiation was before 24hpf and that the first differentiated neurons located more ventral and medial than late-born cells (Figure 3A–A’; see magenta vs. green cells in 3 A’). When we labeled the neurons born before 36hpf we appreciated that neuronal differentiation dramatically increased between 24hpf and 36hpf, since the magenta territory in this case was larger, and that later differentiated neurons piled up on the top of the previously differentiated ones (Figure 3B–B’; see magenta vs. green cells in 3B’). We further analyzed the remodeling of the neuronal differentiation territory, by generating the reference 3D-models of the early-born (24hpf) and late-born (36hpf) neurons in the context of the whole differentiated neuronal domain at 48hpf (Figure 3—video 1A–B; see in green neurons born between 24 and 36hpf, and in blue those born between 36 and 48hpf). We compared them to determine whether the temporal order of neuronal differentiation prefigured the spatial distribution of neurons in the tissue. We overlapped and color-coded them according to the neuronal birthdate, and observed that neurons born at different stages were precisely organized at 48hpf (Figure 3C—C’; see youngest neurons in blue, neurons differentiated between 24 and 36hpf in green, and the oldest ones in magenta). Many neurons differentiated between 24 and 36hpf (compare Figure 3A’ and B’), and they surrounded the existing differentiated neurons generating an inner-outer differentiation gradient (Figure 3C’). To demonstrate that indeed newly born cells were orderly incorporating into the neuronal differentiation domain, we combined HuC:Kaede^Green photoconversion with an additional age-labeling method such as EdU-incorporation. We photoconverted HuC:Kaede^Green at 24hpf, pulsed embryos at 30hpf with EdU to label proliferating cells, and let them grow until 48hpf (Figure 3D). We observed cells with EdU-incorporation exclusively in the newly differentiated domain, stacking dorsally, with no signal indicating a migration of newly born neurons in the older neuronal domain, therefore confirming the establishment of segregated domains according to neuronal differentiation timings (Figure 3E—E’).

Figure 3 with 4 supplements see all

Download asset Open asset

To seek the impact of morphogenesis in the position of different birthdated neurons, we expanded our study to 72hpf, in which most of the cells have already differentiated. As before, we observed that early differentiated neurons were located more ventrally and medially, and that late-born cells piled up (Figure 3F–H, 3F’–H’; see magenta vs. green cells in 3F’–H’). To monitor possible differences in the spatiotemporal distribution within the neuronal differentiation domain, we generated the neuronal birthdate 3D-map at 72hpf and observed that the differentiation capacity diminished with time (Figure 3I–I’’; see small yellow domain). Indeed, by 48hpf, almost the whole differentiation domain was already occupied by previously differentiated neurons (Figure 3I–I’’; see blue, green and magenta domains). This continuous birth of differentiated neurons could be clearly observed by assessing the relative contribution of early- vs. late-born neurons to the neuronal differentiation domain (Figure 3—video 2A–C; see youngest neurons in yellow, those differentiated between 36hpf and 48hpf in blue, and between 24hpf and 36hpf in green). The intersection of the different neuronal differentiation intervals allowed to better appreciate the inner-outer neuronal differentiation gradient (Figure 3I’).

To demonstrate that the birthdating 3D-maps displayed low interindividual variability, we verified that photoconversion of the HuC:Kaede^Green at different stages was consistent in all the analyzed embryos. For this, we first checked that the overall HuC-volumes did not dramatically change among embryos that underwent HuC:Kaede^Green photoconversion at different stages (Figure 3—figure supplement 1A, C). Second, we calculated the difference between the HuC volume of the individuals and the generated models and showed that average individuals’ variability from the model did not exceed 15% (Figure 3—figure supplement 1B, D). And third, we overlaid the HuC:Kaede^Red domains in 72hpf embryos that were generated at different times, and displayed them as a gray gradient (Figure 3—figure supplement 1E–G,E’–G’). Next, we verified that the HuC:Kaede^Red cells labeled in an early conversion were included in the HuC:Kaede^Red domain labeled in a late conversion. For this, we overlaid the differentiation domains contributed at each time interval (Figure 3—figure supplement 1H–J,H’–J’). Overall, these observations allowed us to confidently use the generated birthdating 3D-maps to estimate the importance of neuronal birthdate in predicting the neuronal position within the HuC-domain.

Next, we compared the territory occupied by neurons differentiating at distinct embryonic stages by quantifying the volume size of the neuronal differentiation domains. We detected that at 48hpf, of the whole differentiated domain (12.07 [SD ± 1.36]×10⁶ μm³), only a small part comprised neurons differentiated before 24hpf (3.42 [SD ± 1.56]×10⁶ μm³ early-born neurons vs. 8.64 [SD ± 0.86]×10⁶ μm³ late-born neurons), whereas most of the differentiated neurons arose before 36hpf (9.61 [SD ± 0.68]×10⁶ μm³) and the rest differentiated between 36 and 48hpf (3.10 [SD ± 0.35]×10⁶ μm³). This timely growth in the differentiated neuronal domain could be clearly observed when the relative contribution of early- vs- late-born neurons as compared to the whole neuronal differentiation territory at 48hpf was plotted (Figure 3—figure supplement 2A). When quantifying the signal ratio along the axes, we could observe similar dynamics along the AP and ML axes between 24hpf and 36hpf, and a clear difference along the DV demonstrating that neurons born at 24hpf located more ventrally than neurons born at 36hpf (Figure 3—figure supplement 2B). When volumes were quantified at 72hpf, we observed that only a small part of the differentiation domain came from neurons born before 24hpf (2.77 [SD ± 1.23]×10⁶ μm³ for early-born vs. 11.88 [SD ± 1.28]×10⁶ μm³ for late-born). Most neurons underwent differentiation between 24hpf and 48hpf, with half of the differentiated neurons arising before 36hpf (6.75 [SD ± 2.19]×10⁶ μm³ born before 36hpf vs. 7.58 [SD ± 2.37]×10⁶ μm³ after 36hpf) and the rest were differentiated between 36 and 48hpf (11.19 [SD ± 1.09]×10⁶ μm³ born before 48hpf vs. 2.99 [SD ± 0.90]×10⁶ μm³ after 48hpf). When the ratio of neurons born at different stages over the whole neuronal differentiation domain was assessed, a steady increase of neuronal differentiation domain was observed (Figure 3—figure supplement 2C; 24hpf: 0.21 ± 0.04; 36hpf: 0.47 ± 0.16; 48hpf: 0.79 ± 0.07). Upon quantifying the signal ratio along the axes, the domains corresponding to early- vs. late-born neurons followed the same trend with no main differences along AP and ML axes (Figure 3—figure supplement 2D; e.g., the magenta, green and blue lines did not display major changes in the multiple virtual slices). As expected, when the analysis was performed along the DV axis, differences in cell allocation relying on the neuronal birthdate could be appreciated, with early-born neurons (24-36hpf) displaying a more ventral enrichment (Figure 3—figure supplement 1D, see DV-axis graph).

Overall, the overlaying of birthdating 3D-maps showed that most of the neuronal differentiation events occurred between 24 and 48hpf, and early-born neurons were always located in the inner part of the differentiation domain and surrounded by younger neurons, generating an inner-outer gradient of early-born vs. late-born neurons. This suggested that the position of neurons relied on their neuronal birthdate, and revealed the the importance of morphogenesis for remodeling the differentiation domain.

The birthdating 3D-atlas as a proxy to infer neuronal differentiation order

Our next goal was to investigate whether position in the neuronal differentiation domain at late embryonic stages could be used as a proxy for neuronal birthdating. As a proof-of-concept, we used the Tg[vglut:DsRed] and Tg[gad1b:DsRed] lines, which display DsRed in the glutamatergic and GABAergic neurons, respectively (Satou et al., 2013). To dissect when glutamatergic neurons were born, we generated the 3D-models of neuronal differentiation of Tg[vglut2a:DsRed] at 48 and 72hpf, and intersected them with the birthdate 3D-maps created in Figure 3, in such a manner that glutamatergic neurons were color-coded according to their birthdate (Figure 4A—A’ and B—B’). Glutamatergic neurons contributed in a big proportion to the differentiation domain, and were born during all the analyzed temporal windows (Figure 4—video 1A–B). They displayed the inner-outer distribution depending on birth-time, with the oldest cells cells in the inner domain (see magenta cells in Figure 4A’–B’; Figure 4—video 1A–B), and the youngest neurons in the most outer position (see blue and yellow domains in Figure 4A’–B’; Figure 4—video 1A–B). In the case of GABAergic neurons, we took the same approach and intersected the 3D-model of Tg[gad1b:DsRed] embryos on the birthdate 3D-map, and color-coded them according to the differentiation time. Some GABAergic neurons in r3 and r5 were born as early as 24 hpf (Figure 4C–C’; Figure 4—video 1C–D). These observations suggested that an important part of glutamatergic and GABAergic neurons was born between 24 and 36hpf, although some were born earlier and later too. Accordingly, when assessing the relative contribution of neurons born at different intervals to the corresponding neuronal populations, we observed that of the total glutamatergic and GABAergic neurons, approximately 20% and 30%, respectively, were contributed by cells born before 24hpf, 60% by cells born between 26 and 36hpf, and only by a small percentage of neurons born between 36 and 48hpf (Figure 4E–F). Comparison of the growth dynamics of both populations at 72hpf showed that they were slightly different, with glutamatergic neurons produced earlier than GABAergic neurons (Figure 4G). Overall, this analysis indicated that position of differentiated neurons at 72hpf could recapitulate neuronal birthdates.

Figure 4 with 2 supplements see all

Download asset Open asset

To demonstrate that indeed this could be used as a proxy to infer birthdate, and that there were glutamatergic and GABAergic neurons at 24hpf, we performed in situ hybridizations with vglut2a and gad1b probes at early stages of embryonic development while monitoring the neuronal differentiation domain. Indeed, glutamatergic and GABAergic neurons were already specified at 24hpf and they were located within the HuC territory (Figure 4H–H’’), although in our hands the GABAergic transgenic line did not display the reporter expression at this stage. Thus, these results demonstrate that the neuronal birthdating 3D-map can be used as a proxy to infer neuronal differentiation order, and widens the use of the digital 3D-atlas. To go deeper in the differentiation dynamics of the glutamatergic and GABAergic populations, we quantified the ratio of neurons differentiated in distinct time windows along the axes (Figure 4—figure supplement 1). Analysis at 48hpf showed that the glutamatergic populations born between 24 and 36hpf and between 36 and 48hpf underwent a drastic displacement along the dorsoventral (DV) axis, when compared to the same populations observed at 72hpf. The 24-36hpf neuronal subset, which displayed a relatively homogenous distribution along the DV axis when observed at 48hpf, was medially located at 72hpf (engulfing the 0-24hpf subset). The 36-48hpf subset, mostly dorsal when observed at 48hpf, was distributed both ventrally and dorsally when observed at 72hpf (Figure 4—figure supplement 1A–B). These observations support the idea that newly produced neurons initially located dorsally were displaced towards ventral, resulting in the formation of the inner to outer gradient of neuronal differentiation. In addition, the medio-lateral distribution of the neurons produced at different time intervals recapitulated the different allocations of late born glutamatergic and GABAergic neurons, with the 48-72hpf subset of glutamatergic neurons more laterally located than the GABAergic 48-72hpf subset, which was enriched medially (Figure 4—figure supplement 1B, D).

Tissue shape analysis reveals defects in neurogenesis and the remodeling of the neuronal differentiation domains

To take this newly developed pipeline one step further we tested whether it could serve as a proof-of-concept to analyze aberrant neurogenesis. For this, we disturbed embryonic neurogenesis by conditionally inhibiting the Notch-pathway, which has been demonstrated to play an important role during hindbrain neurogenesis (Hevia et al., 2022; Nikolaou et al., 2009). We made use of photoconversion experiments in Tg[HuC:Kaede] embryos to discern early-born vs. late-born neurons, conditionally inhibited Notch-signaling for a short time (only affecting to newly generated neurons), and embryos were then imaged at 72hpf (Figure 5A). When comparing control embryos with those in which Notch-signaling was downregulated, we observed an increase in newly differentiated neurons that dramatically impacted in the early-born neuronal domain, since newer neurons occupied the most dorsal part of the hindbrain displacing the older neurons (Figure 5B–B’ and C–C’) as previously described (Belzunce et al., 2020; Hevia et al., 2022). Accordingly, upon quantification of the neuronal differentiation domain we observed a significant increase of the HuC volume in treated embryos (Figure 5D), and newly produced neurons were mainly localized in the medial domain (Figure 5E). To deepen our knowledge on the impact of the neurogenesis increase and the remodeling of the hindbrain, we developed a protocol for tissue shape comparison. We generated the 3D-models of the control embryos and of those with higher neurogenesis (in which Notch-signaling was inhibited), extracted the mesh from them and measured the signed distance between the two mesh surfaces in more than 400,000 points across the models. Then, the measurements were displayed as a color-coded map on the control model. When we compared the rendered hindbrains of both groups of embryos we observed a consistent difference in the overall tissue shape, mainly in the dorsomedial differentiation domain (Figure 5F–F’; see color-coded legend showing in blue the highest differences; Figure 5—video 1). This result was supported by the previous observation in which the brain ventricle was distorted due to the increase of neurons. The increased differentiated domain volume was visible mainly in the dorsomedial part, which was occupied by new neurons and remodeled the neuronal differentiation domain.

Figure 5 with 1 supplement see all

Download asset Open asset

One of the main current challenges in biology is the 3D-anatomical reconstruction of the tissues’ growth over time together with the cell lineage histories, which requires to extract meaningful biological information from large imaging datasets. Some of the bottlenecks are the required computing power for the visualization and analysis of the data volumes, and the need of user-friendly tools to be paired with imaging techniques. Here, we developed a digital 3D-atlas maker (DAMAKER), which is a dynamic and expandable 3D-tool for monitoring the temporal changes in growth and shape during tissue morphogenesis. We built a digital representation of the developing vertebrate hindbrain from 3D-imaging data collected at different times, which can be used for further quantitative analysis and multilevel modeling. Our goal is to share an established protocol using open-access tools and algorithms that allow standardized and accessible tissue-wide cell population atlas construction and can be used to create an expandable atlas of any tissue in any organism.

In this work, our interest was to explore how neuronal differentiation and morphogenesis might be intertwined during hindbrain development. For this, we registered experimental whole hindbrain images and mapped neuronal differentiation patterns with a unified spatial representation, using confocal imaging setup and leveraging cross-platform segmentation and registration tools. To improve comparison of various signals across embryos that requires an alignment/registration, which is often done manually and is very time-consuming, we used the Fijiyama plugin in FIJI (Schindelin et al., 2012). This allowed a multi-modal automatic rigid registration process and a large range of customization, and it permitted the application of a previously computed transform; namely, an alignment can be computed using one channel of an embryo matching the chosen reference and this can be applied to another channel that shares no features with the reference signal. This enables the use of virtually any signal previously added to the atlas, as reference to align others.

The quantitative evaluation of datasets often involves manual annotation, which is time-consuming, introduces biases and often constitutes the major constraint in an evaluation pipeline. To overcome this problem, we used the Trainable Weka Segmentation 3D (TWS), (Arganda-Carreras et al., 2017), which combines the image processing toolkit FIJI with the state-of-the-art machine learning algorithms provided in the latest version of the data mining and machine learning toolkit Waikato Environment for Knowledge Analysis (WEKA) (Hall et al., 2009). This tool is modular, allows sophisticated data mining processes to be built from the wide collection of base learning algorithms and tools provided, is open-sourced, and can be easily used as a plugin of FIJI. Users’ trained classifiers, which are a set of filters and thresholds determined by the algorithm to best match input from the user, can be used repeatedly to segment each dataset. By keeping our pipeline in the FIJI ecosystem, we aim at providing a user-friendly tool to share data, upgrade methodology and adapt the pipeline to a specific experimental set up (model organism, organ etc), thereby making it accessible to anyone.

The digital 3D-brain atlases are essential for neurobiology because they allow to address new questions that can not be answered by the use of the available 2D-atlases. The addition of time as a key parameter is an asset, since it will permit to follow cells upon time, and therefore cell displacements and tissue growth. For this, we made use of an in vivo cell birthdating system that allowed us to monitor intervals of neuronal differentiation (Caron et al., 2008). Thus, this temporal registration using data acquired in vivo at different developmental stages provides important advantages, as it allowed us to (i) dissect the final location of distinct birthdated cell populations, and (ii) obtain a whole-cell population coverage by analyzing the whole differentiated territory, while maintaining the in vivo positioning of cells within the structure. Until now, most of the currently used approaches either lack whole population coverage or are limited by expression timing of the gene or transgene; more importantly, although they permit coverage of the whole cell population, they do not allow the user to foreshadow the birthdate of a given cell population at the time of interest. Here, with the combination of photoconversion experiments as a tool to register temporality of neuronal differentiation and DAMAKER, we provide neuronal birthdating 3D-maps and revealed that the order of neuronal differentiation prefigured the spatial distribution of neurons in the hindbrain, with an inner-outer gradient of temporal differentiation.

Usually, expression of the reporter genes in transgenic lines is delayed as compared to the gene expression onset, unless the reporter has been inserted within the given gene locus. Now, we can specifically knock-in genes within a targeted genome site thanks to the CRISPR-Cas technology, but most of the available —and very valuable—transgenic lines were (and still are) generated with random insertion systems that do not always recapitulate the right onset of gene expression. Thus, our tool can be very useful for predicting neuronal birthdates just by analyzing the position of neurons in the differentiation domain as we showed with glutamatergic and GABAergic neurons. Besides that, this developed digital platform will be useful to interrogate phenotypes in neurological diseases, and zebrafish avatars from human neurological disorders could be analyzed in a very simple manner. As a proof-of-concept, we show tissue shape comparisons after modulating neurogenesis using pharmacological treatments. Thus, the combination of several features can provide computational models allowing to merge information from several high-throughput experiments to enlarge our reductionist view, and implement these digital tools for neuronal disorder studies.

In addition, brain morphogenesis studies could use our digital atlas-maker tool, which maps neuronal differentiation domains in the embryonic brain over time and can be easily shared. In the future, DAMAKER can be further improved by using batch processing algorithms (e.g., to assess the effects of pharmacological treatments or gene disruptions) or by adjusting for more geometrically complex tissues or for tissues with less alignment features (such as heart). Last but not least, this pipeline could be applied to time-lapse videos by repeating the process for each frame. This approach will help us to fill the knowledge gap between gene regulatory networks, cell populations and tissue architecture, as it will allow to follow the dynamic events that pattern the early embryonic organs as they happen within the intact system.

All data generated or analysed during this study are included in the manuscript or in the supplementary files. All code is available in https://github.com/cristinapujades/Blanc-et-al.-2022, (copy archived at swh:1:rev:82c8199b774af1b4ff68d7c0ccd2cb8d4197df56).

1. 1. Arganda-Carreras I
   2. Kaynig V
   3. Rueden C
   4. Eliceiri KW
   5. Schindelin J
   6. Cardona A
   7. Sebastian Seung H

   (2017) Trainable weka segmentation: a machine learning tool for microscopy pixel classification

   Bioinformatics 33:2424–2426.

   https://doi.org/10.1093/bioinformatics/btx180

   • PubMed
   • Google Scholar
2. 1. Belzunce I
   2. Belmonte-Mateos C
   3. Pujades C

   (2020) The interplay of Atoh1 genes in the lower rhombic lip during hindbrain morphogenesis

   PLOS ONE 15:e0228225.

   https://doi.org/10.1371/journal.pone.0228225

   • PubMed
   • Google Scholar
3. Software

   1. Blanc M

   (2022) Cristinapujades, version swh:1:rev:82c8199b774af1b4ff68d7c0ccd2cb8d4197df56

   Software Heritage.

   https://archive.softwareheritage.org/swh:1:dir:133b092a8bfd4e61983aec841b3ff9b91c6d9763;origin=https://github.com/cristinapujades/Blanc-et-al.-2022;visit=swh:1:snp:f7d2a08f163b0c87c925def76acfd880868b8e6c;anchor=swh:1:rev:82c8199b774af1b4ff68d7c0ccd2cb8d4197df56
4. 1. Bourge M
   2. Fort C
   3. Soler M-N
   4. Satiat-Jeunemaître B
   5. Brown SC

   (2015) A pulse-chase strategy combining click-edu and photoconvertible fluorescent reporter: tracking Golgi protein dynamics during the cell cycle

   The New Phytologist 205:938–950.

   https://doi.org/10.1111/nph.13069

   • PubMed
   • Google Scholar
5. 1. Caron SJC
   2. Prober D
   3. Choy M
   4. Schier AF

   (2008) In vivo birthdating by BAPTISM reveals that trigeminal sensory neuron diversity depends on early neurogenesis

   Development 135:3259–3269.

   https://doi.org/10.1242/dev.023200

   • PubMed
   • Google Scholar
6. 1. Chow DM
   2. Sinefeld D
   3. Kolkman KE
   4. Ouzounov DG
   5. Akbari N
   6. Tatarsky R
   7. Bass A
   8. Xu C
   9. Fetcho JR

   (2020) Deep three-photon imaging of the brain in intact adult zebrafish

   Nature Methods 17:605–608.

   https://doi.org/10.1038/s41592-020-0819-7

   • PubMed
   • Google Scholar
7. 1. Distel M
   2. Wullimann MF
   3. Köster RW

   (2009) Optimized GAL4 genetics for permanent gene expression mapping in zebrafish

   PNAS 106:13365–13370.

   https://doi.org/10.1073/pnas.0903060106

   • PubMed
   • Google Scholar
8. 1. Drerup CM
   2. Nechiporuk AV

   (2013) Jnk-Interacting protein 3 mediates the retrograde transport of activated c-Jun N-terminal kinase and lysosomes

   PLOS Genetics 9:e1003303.

   https://doi.org/10.1371/journal.pgen.1003303

   • PubMed
   • Google Scholar
9. 1. Dsilva CJ
   2. Lim B
   3. Lu H
   4. Singer A
   5. Kevrekidis IG
   6. Shvartsman SY

   (2015) Temporal ordering and registration of images in studies of developmental dynamics

   Development 142:1717–1724.

   https://doi.org/10.1242/dev.119396

   • PubMed
   • Google Scholar
10. 1. Fernandez R
    2. Moisy C

    (2020) Fijiyama: a registration tool for 3D multimodal time-lapse imaging

    Bioinformatics 37:1482–1484.

    https://doi.org/10.1093/bioinformatics/btaa846

    • PubMed
    • Google Scholar
11. 1. Fraser S
    2. Keynes R
    3. Lumsden A

    (1990) Segmentation in the chick embryo hindbrain is defined by cell lineage restrictions

    Nature 344:431–435.

    https://doi.org/10.1038/344431a0

    • PubMed
    • Google Scholar
12. 1. Guillemot F

    (2007) Spatial and temporal specification of neural fates by transcription factor codes

    Development 134:3771–3780.

    https://doi.org/10.1242/dev.006379

    • PubMed
    • Google Scholar
13. 1. Gutzman JH
    2. Graeden EG
    3. Lowery LA
    4. Holley HS
    5. Sive H

    (2008) Formation of the zebrafish midbrain-hindbrain boundary constriction requires laminin-dependent basal constriction

    MECHANISMS OF DEVELOPMENT 125:974–983.

    https://doi.org/10.1016/j.mod.2008.07.004

    • PubMed
    • Google Scholar
14. 1. Hall M
    2. Frank E
    3. Holmes G
    4. Pfahringer B
    5. Reutemann P
    6. Witten IH

    (2009) The WEKA data mining software

    ACM SIGKDD Explorations Newsletter 11:10–18.

    https://doi.org/10.1145/1656274.1656278

    • Google Scholar
15. 1. Harrison C
    2. Wabbersen T
    3. Shepherd IT

    (2014) In vivo visualization of the development of the enteric nervous system using a TG (-8.3bphox2b: kaede) transgenic zebrafish

    Genesis 52:985–990.

    https://doi.org/10.1002/dvg.22826

    • PubMed
    • Google Scholar
16. 1. Heckscher ES
    2. Long F
    3. Layden MJ
    4. Chuang CH
    5. Manning L
    6. Richart J
    7. Pearson JC
    8. Crews ST
    9. Peng H
    10. Myers E
    11. Doe CQ

    (2014) Atlas-builder software and the eneuro atlas: resources for developmental biology and neuroscience

    Development 141:2524–2532.

    https://doi.org/10.1242/dev.108720

    • PubMed
    • Google Scholar
17. 1. Hevia CF
    2. Engel-Pizcueta C
    3. Udina F
    4. Pujades C

    (2022) The neurogenic fate of the hindbrain boundaries relies on notch3-dependent asymmetric cell divisions

    Cell Reports 39:110915.

    https://doi.org/10.1016/j.celrep.2022.110915

    • PubMed
    • Google Scholar
18. 1. Higashijima S
    2. Hotta Y
    3. Okamoto H

    (2000)

    Visualization of cranial motor neurons in live transgenic zebrafish expressing green fluorescent protein under the control of the islet-1 promoter/enhancer

    The Journal of Neuroscience 20:206–218.

    • PubMed
    • Google Scholar
19. 1. Higashijima S-I
    2. Schaefer M
    3. Fetcho JR

    (2004) Neurotransmitter properties of spinal interneurons in embryonic and larval zebrafish

    The Journal of Comparative Neurology 480:19–37.

    https://doi.org/10.1002/cne.20279

    • PubMed
    • Google Scholar
20. 1. Jaggard JB
    2. Lloyd E
    3. Yuiska A
    4. Patch A
    5. Fily Y
    6. Kowalko JE
    7. Appelbaum L
    8. Duboue ER
    9. Keene AC

    (2020) Cavefish brain atlases reveal functional and anatomical convergence across independently evolved populations

    Science Advances 6:eaba3126.

    https://doi.org/10.1126/sciadv.aba3126

    • PubMed
    • Google Scholar
21. 1. Jimenez-Guri E
    2. Udina F
    3. Colas JF
    4. Sharpe J
    5. Padrón-Barthe L
    6. Torres M
    7. Pujades C

    (2010) Clonal analysis in mice underlines the importance of rhombomeric boundaries in cell movement restriction during hindbrain segmentation

    PLOS ONE 5:e10112.

    https://doi.org/10.1371/journal.pone.0010112

    • PubMed
    • Google Scholar
22. 1. Kenney JW
    2. Steadman PE
    3. Young O
    4. Shi MT
    5. Polanco M
    6. Dubaishi S
    7. Covert K
    8. Mueller T
    9. Frankland PW

    (2021) A 3D adult zebrafish brain atlas (AZBA) for the digital age

    eLife 10:e69988.

    https://doi.org/10.7554/eLife.69988

    • PubMed
    • Google Scholar
23. 1. Kiecker C
    2. Lumsden A

    (2005) Compartments and their boundaries in vertebrate brain development

    Nature Reviews. Neuroscience 6:553–564.

    https://doi.org/10.1038/nrn1702

    • PubMed
    • Google Scholar
24. 1. Kinkhabwala A
    2. Riley M
    3. Koyama M
    4. Monen J
    5. Satou C
    6. Kimura Y
    7. Higashijima SI
    8. Fetcho J

    (2011) A structural and functional ground plan for neurons in the hindbrain of zebrafish

    PNAS 108:1164–1169.

    https://doi.org/10.1073/pnas.1012185108

    • PubMed
    • Google Scholar
25. 1. Koyama M
    2. Kinkhabwala A
    3. Satou C
    4. Higashijima S
    5. Fetcho J

    (2011) Mapping a sensory-motor network onto a structural and functional ground plan in the hindbrain

    PNAS 108:1170–1175.

    https://doi.org/10.1073/pnas.1012189108

    • PubMed
    • Google Scholar
26. 1. Krumlauf R
    2. Wilkinson DG

    (2021) Segmentation and patterning of the vertebrate hindbrain

    Development 148:dev186460.

    https://doi.org/10.1242/dev.186460

    • PubMed
    • Google Scholar
27. 1. Kunst M
    2. Laurell E
    3. Mokayes N
    4. Kramer A
    5. Kubo F
    6. Fernandes AM
    7. Förster D
    8. Dal Maschio M
    9. Baier H

    (2019) A cellular-resolution atlas of the larval zebrafish brain

    Neuron 103:21–38.

    https://doi.org/10.1016/j.neuron.2019.04.034

    • PubMed
    • Google Scholar
28. 1. Lyons DA
    2. Guy AT
    3. Clarke JDW

    (2003) Monitoring neural progenitor fate through multiple rounds of division in an intact vertebrate brain

    Development 130:3427–3436.

    https://doi.org/10.1242/dev.00569

    • PubMed
    • Google Scholar
29. 1. Murakami Y
    2. Uchida K
    3. Rijli FM
    4. Kuratani S

    (2005) Evolution of the brain developmental plan: insights from agnathans

    Developmental Biology 280:249–259.

    https://doi.org/10.1016/j.ydbio.2005.02.008

    • PubMed
    • Google Scholar
30. 1. Nikolaou N
    2. Watanabe-Asaka T
    3. Gerety S
    4. Distel M
    5. Köster RW
    6. Wilkinson DG

    (2009) Lunatic fringe promotes the lateral inhibition of neurogenesis

    Development 136:2523–2533.

    https://doi.org/10.1242/dev.034736

    • PubMed
    • Google Scholar
31. 1. Park HC
    2. Kim CH
    3. Bae YK
    4. Yeo SY
    5. Kim SH
    6. Hong SK
    7. Shin J
    8. Yoo KW
    9. Hibi M
    10. Hirano T
    11. Miki N
    12. Chitnis AB
    13. Huh TL

    (2000) Analysis of upstream elements in the huC promoter leads to the establishment of transgenic zebrafish with fluorescent neurons

    Developmental Biology 227:279–293.

    https://doi.org/10.1006/dbio.2000.9898

    • PubMed
    • Google Scholar
32. 1. Peng H
    2. Chung P
    3. Long F
    4. Qu L
    5. Jenett A
    6. Seeds AM
    7. Myers EW
    8. Simpson JH

    (2011) BrainAligner: 3D registration atlases of Drosophila brains

    Nature Methods 8:493–500.

    https://doi.org/10.1038/nmeth.1602

    • PubMed
    • Google Scholar
33. 1. Pujades C

    (2020) The multiple functions of hindbrain boundary cells: tinkering boundaries?

    Seminars in Cell & Developmental Biology 107:179–189.

    https://doi.org/10.1016/j.semcdb.2020.05.002

    • PubMed
    • Google Scholar
34. 1. Pujala A
    2. Koyama M

    (2019) Chronology-based architecture of descending circuits that underlie the development of locomotor repertoire after birth

    eLife 8:e42135.

    https://doi.org/10.7554/eLife.42135

    • PubMed
    • Google Scholar
35. 1. Ronneberger O
    2. Liu K
    3. Rath M
    4. Rueβ D
    5. Mueller T
    6. Skibbe H
    7. Drayer B
    8. Schmidt T
    9. Filippi A
    10. Nitschke R
    11. Brox T
    12. Burkhardt H
    13. Driever W

    (2012) ViBE-Z: a framework for 3D virtual colocalization analysis in zebrafish larval brains

    Nature Methods 9:735–742.

    https://doi.org/10.1038/nmeth.2076

    • PubMed
    • Google Scholar
36. 1. Satou C
    2. Kimura Y
    3. Hirata H
    4. Suster ML
    5. Kawakami K
    6. Higashijima S

    (2013) Transgenic tools to characterize neuronal properties of discrete populations of zebrafish neurons

    Development 140:3927–3931.

    https://doi.org/10.1242/dev.099531

    • PubMed
    • Google Scholar
37. 1. Schindelin J
    2. Arganda-Carreras I
    3. Frise E
    4. Kaynig V
    5. Longair M
    6. Pietzsch T
    7. Preibisch S
    8. Rueden C
    9. Saalfeld S
    10. Schmid B
    11. Tinevez JY
    12. White DJ
    13. Hartenstein V
    14. Eliceiri K
    15. Tomancak P
    16. Cardona A

    (2012) Fiji: an open-source platform for biological-image analysis

    Nature Methods 9:676–682.

    https://doi.org/10.1038/nmeth.2019

    • PubMed
    • Google Scholar
38. 1. Sestelo M
    2. Villanueva NM
    3. Softw L

    (2017) Npregfast: an R package for nonparametric estimation and inference in life sciences

    Journal of Statistical Software 82:1–27.

    https://doi.org/10.18637/jss.v082.i12

    • Google Scholar
39. 1. Tabor KM
    2. Marquart GD
    3. Hurt C
    4. Smith TS
    5. Geoca AK
    6. Bhandiwad AA
    7. Subedi A
    8. Sinclair JL
    9. Rose HM
    10. Polys NF
    11. Burgess HA

    (2019) Brain-Wide cellular resolution imaging of Cre transgenic zebrafish lines for functional circuit-mapping

    eLife 8:e42687.

    https://doi.org/10.7554/eLife.42687

    • PubMed
    • Google Scholar
40. 1. Thisse C
    2. Thisse B

    (2008) High-Resolution in situ hybridization to whole-mount zebrafish embryos

    Nature Protocols 3:59–69.

    https://doi.org/10.1038/nprot.2007.514

    • PubMed
    • Google Scholar
41. 1. Voltes A
    2. Hevia CF
    3. Engel-Pizcueta C
    4. Dingare C
    5. Calzolari S
    6. Terriente J
    7. Norden C
    8. Lecaudey V
    9. Pujades C

    (2019) Yap/taz-TEAD activity links mechanical cues to progenitor cell behavior during zebrafish hindbrain segmentation

    Development 146:dev176735.

    https://doi.org/10.1242/dev.176735

    • PubMed
    • Google Scholar
42. 1. Wan Y
    2. Wei Z
    3. Looger LL
    4. Koyama M
    5. Druckmann S
    6. Keller PJ

    (2019) Single-Cell reconstruction of emerging population activity in an entire developing circuit

    Cell 179:355–372.

    https://doi.org/10.1016/j.cell.2019.08.039

    • PubMed
    • Google Scholar
43. 1. Zhang Q
    2. Xu H
    3. Zhao W
    4. Zheng J
    5. Sun L
    6. Luo C

    (2020) Zygotic VSX1 plays a key role in defining v2a interneuron sub-lineage by directly repressing TAL1 transcription in zebrafish

    International Journal of Molecular Sciences 21:E3600.

    https://doi.org/10.3390/ijms21103600

    • PubMed
    • Google Scholar

Author details

1. Matthias Blanc

   Department of Medicine and Life Sciences, Universitat Pompeu Fabra, Barcelona, Spain

   Contribution

   Conceptualization, Resources, Data curation, Software, Formal analysis, Validation, Investigation, Visualization, Methodology, Writing - original draft, Writing - review and editing

   Competing interests

   No competing interests declared

   "This ORCID iD identifies the author of this article:" 0000-0002-4892-8161

2. Giovanni Dalmasso

   European Molecular Biology Laboratory, Barcelona, Spain

   Contribution

   Software

   Competing interests

   No competing interests declared

3. Frederic Udina

   Department of Economics and Business, Universitat Pompeu Fabra, Barcelona, Spain

   Contribution

   Software

   Competing interests

   No competing interests declared

   "This ORCID iD identifies the author of this article:" 0000-0002-7155-4153

4. Cristina Pujades

   Department of Medicine and Life Sciences, Universitat Pompeu Fabra, Barcelona, Spain

   Contribution

   Conceptualization, Formal analysis, Supervision, Funding acquisition, Investigation, Writing - original draft, Writing - review and editing

   For correspondence

   cristina.pujades@upf.edu

   Competing interests

   No competing interests declared

   "This ORCID iD identifies the author of this article:" 0000-0001-6423-7451

• 941

  views

• 186

  downloads

• 0

  citations

Views, downloads and citations are aggregated across all versions of this paper published by eLife.