1. Developmental Biology
  2. Neuroscience

Dermal appendage-dependent patterning of zebrafish atoh1a+ Merkel cells

  1. Tanya L Brown
  2. Emma C Horton
  3. Evan W Craig
  4. Camille EA Goo
  5. Erik C Black
  6. Madeleine N Hewitt
  7. Nathaniel G Yee
  8. Everett T Fan
  9. David W Raible
  10. Jeffrey P Rasmussen  Is a corresponding author
  1. Department of Biology, University of Washington, United States
  2. Molecular and Cellular Biology Program, University of Washington, United States
  3. Department of Biological Structure, University of Washington, United States
  4. Department of Otolaryngology - Head and Neck Surgery, University of Washington, United States
  5. Institute for Stem Cell and Regenerative Medicine, University of Washington, United States
https://doi.org/10.7554/eLife.85800
Share this article
Cite this article
  1. Tanya L Brown
  2. Emma C Horton
  3. Evan W Craig
  4. Camille EA Goo
  5. Erik C Black
  6. Madeleine N Hewitt
  7. Nathaniel G Yee
  8. Everett T Fan
  9. David W Raible
  10. Jeffrey P Rasmussen
(2023)
Dermal appendage-dependent patterning of zebrafish atoh1a+ Merkel cells
eLife 12:e85800.
https://doi.org/10.7554/eLife.85800
  2. Download BibTeX
  3. Download .RIS

Abstract

Touch system function requires precise interactions between specialized skin cells and somatosensory axons, as exemplified by the vertebrate mechanosensory Merkel cell-neurite complex. Development and patterning of Merkel cells and associated neurites during skin organogenesis remain poorly understood, partly due to the in utero development of mammalian embryos. Here, we discover Merkel cells in the zebrafish epidermis and identify Atonal homolog 1a (Atoh1a) as a marker of zebrafish Merkel cells. We show that zebrafish Merkel cells derive from basal keratinocytes, express neurosecretory and mechanosensory machinery, extend actin-rich microvilli, and complex with somatosensory axons, all hallmarks of mammalian Merkel cells. Merkel cells populate all major adult skin compartments, with region-specific densities and distribution patterns. In vivo photoconversion reveals that Merkel cells undergo steady loss and replenishment during skin homeostasis. Merkel cells develop concomitant with dermal appendages along the trunk and loss of Ectodysplasin signaling, which prevents dermal appendage formation, reduces Merkel cell density by affecting cell differentiation. By contrast, altering dermal appendage morphology changes the distribution, but not density, of Merkel cells. Overall, our studies provide insights into touch system maturation during skin organogenesis and establish zebrafish as an experimentally accessible in vivo model for the study of Merkel cell biology.

Editor's evaluation

The authors describe and characterize the touch system in zebrafish as a new model to study Merkel cell development and maintenance. The study demonstrates that the zebrafish touch system shares many characteristics with its mammalian counterpart, including developmental origin, innervation, and molecular characteristics while allowing in vivo analysis of specification, development, and maintenance. This study is the foundation for future detailed cellular and molecular analyses of the touch sensory system and will be of interest to developmental biologists and neuroscientists studying stem cells, regeneration, and aging.

https://doi.org/10.7554/eLife.85800.sa0

Introduction

Skin functions as our primary interface with the physical environment and can distinguish a range of tactile inputs with exquisite acuity. As the skin undergoes organogenesis, the epidermis transforms from a simple, uniform epithelium into a complex, diverse tissue. During these dramatic changes, the skin develops regionally specialized sensory structures and becomes innervated by specific types of somatosensory neurites (reviewed by Jenkins and Lumpkin, 2017). Interactions between somatosensory neurites and cutaneous cell types regulate diverse tactile responses (reviewed by Handler and Ginty, 2021). Altered tactile sensitivity during early mammalian development has been associated with neurodevelopmental disorders (reviewed by Orefice, 2020), underscoring the importance of understanding the cellular and molecular basis of touch system development and function.

Merkel cells (MCs), a specialized mechanosensory cell type found in the vertebrate epidermis (reviewed by Hartschuh et al., 1986), densely populate many highly sensitive regions of skin (Lacour et al., 1991). MCs have several defining cellular characteristics that distinguish them from other epidermal cell types: they are relatively small, extend actin-rich microvilli, contain cytoplasmic granules reminiscent of synaptic vesicles, and form contacts with somatosensory axons (Hartschuh and Weihe, 1980; Mihara et al., 1979; Smith, 1977; Toyoshima et al., 1998). In mammals, a subset of cutaneous somatosensory axons known as Aß slowly adapting type I low-threshold mechanoreceptors (SAI-LTMRs) innervate MCs, forming the MC-neurite complex. MCs detect mechanical inputs via the cation channel Piezo2 (Ikeda et al., 2014; Maksimovic et al., 2014; Woo et al., 2014) and play an active role in touch sensation by releasing neurotransmitters to activate neighboring neurites (Chang et al., 2016; Chang and Gu, 2020; Hoffman et al., 2018). Genetic ablation of rodent MCs indicates they are required for specific aspects of touch system function, including promoting the static phase of the slowly adapting response of Aß SAI-LTMRs and sensory tasks such as texture discrimination (Maricich et al., 2012; Maricich et al., 2009).

Molecular control of MC development has primarily been studied in rodent hairy skin (reviewed by Oss-Ronen and Cohen, 2021). While this system has been useful for understanding many aspects of MC development and function, the rodent system also has several significant limitations that warrant additional models to improve the understanding of MCs. First, vertebrates have diverse types of skin, and MCs are found in both hairy and glabrous (non-hairy) skin, as well as mucocutaneous regions such as the gingiva and palate (Hashimoto, 1972; Lacour et al., 1991; Moayedi et al., 2021). Importantly, MC populations within different skin compartments share similar transcriptional profiles (Nguyen et al., 2019). Thus, the establishment of complementary genetic systems in different types of skin could help reveal both shared and divergent principles of MC development. Second, because in utero development of mammalian skin limits access to the developing touch system—combined with technical limitations of imaging intact mammalian skin—the dynamics of MC development and innervation remain essentially unknown. Third, unbiased screens for regulators of MC development would be difficult or impractical in rodents due to the prohibitive cost of animal housing and difficulty of visualizing MCs in situ.

Anamniote model systems, such as the genetically tractable zebrafish, provide the potential to overcome these limitations. Interestingly, despite the different tactile environments encountered by terrestrial and aquatic vertebrates, MCs have been described by transmission electron microscopy (TEM) in a wide variety of anamniotes, including teleost (ray-finned) fish, lungfish, and lamprey (Fox et al., 1980; Lane and Whitear, 1977; Whitear, 1989; Whitear and Lane, 1981). Here, we identify and characterize a population of zebrafish epidermal cells that we propose are bona fide MCs. Our studies establish the zebrafish as a promising new model to investigate the developmental and cellular biology of MCs.

Results

Ultrastructural identification of presumptive MCs in the adult epidermis

Given the presence of cells with the ultrastructural characteristics of MCs in several teleosts (Lane and Whitear, 1977; Whitear, 1989), we reasoned that the zebrafish epidermis may contain similar cells. Whitear, 1989 defined five ultrastructural criteria for the identification of vertebrate MCs: (1) a relatively small volume of cytoplasm; (2) an association with a nerve fiber; (3) the presence of cytoplasmic granules; (4) desmosomal attachments to neighboring cells; and (5) peripheral microvilli.

We previously demonstrated that somatosensory axons densely innervate the epidermis above scales (Rasmussen et al., 2018), dermal appendages that cover the adult zebrafish trunk (Figure 1A). By TEM, we found that many of the axon endings in the scale epidermis arborize between keratinocyte membranes (Rasmussen et al., 2018). Interestingly, however, we identified additional axon-associated epidermal cells that were distinct from the large, cuboidal keratinocytes that comprise most of the epidermis based on several characteristics. The cells appeared relatively small and spherical with a low cytoplasmic-to-nuclear ratio compared to neighboring keratinocytes (Figure 1B and C; Figure 1—figure supplement 1), contained cytoplasmic vesicles that in some instances localized adjacent to axon contacts (Figure 1B’), and formed desmosomal-like attachments with neighboring keratinocytes (Figure 1B’’ and B’’’). Furthermore, the cells extended spike-like microvillar processes that contacted adjacent cells (Figure 1C and C’). Thus, based on established TEM criteria, we identified presumptive MCs in the adult scale epidermis.

Figure 1 with 3 supplements see all
Download asset Open asset
The adult scale epidermis contains atoh1a+ Merkel cells (MCs).

(A) Illustration of the adult zebrafish trunk anatomy showing the organization of epidermis, scales, and dermis. Scales are flat bony discs arranged in an overlapping, imbricated pattern and coated on their external surface by epidermis. (B) Transmission electron microscopy (TEM) of a presumptive MC from the scale epidermis. Dotted boxes indicate regions of magnification in B’–B’’’. (B’) Magnification of B showing cytoplasmic granules (g, brackets) juxtaposed to a putative axon (a) contact containing a mitochondrion (m). (B’’ and B’’’) Magnifications of B showing desmosomal-like (d, arrows) attachments between keratinocytes (B’’) and between a presumptive MC and keratinocyte (B’’’). (C and C’) TEM of a presumptive MC from the scale epidermis showing a microvillar process (p, arrowhead). (D) Illustration of a cross section of the scale epidermis based on TEM observations. Periderm cells (superficial epidermis; dark blue) are located in the uppermost epidermal stratum, and basal keratinocytes (light blue) are located in the lowermost epidermal stratum. MCs containing cytoplasmic granules, extending microvillar processes, and contacting axons localize between keratinocytes. (E) Lateral confocal micrograph of the trunk epidermis in an adult expressing reporters for keratinocytes (Tg(actb2:LOXP-BFP-LOXP-DsRed)) and atoh1a-expressing cells (Tg(atoh1a:nls-Eos)). Dotted boxes indicate areas of magnification in E’ and E’’. (E’) Magnification of E showing atoh1a+ hair cells (HCs) and progenitors within neuromasts (nm) of the posterior lateral line. (E’’) Magnification of E showing atoh1a+ MCs scattered throughout the scale epidermis. (F) Lateral and reconstructed cross sectional confocal micrographs of the trunk in an adult expressing reporters for keratinocytes (Tg(actb2:LOXP-BFP-LOXP-DsRed)) and atoh1a-expressing cells (Tg(atoh1a:nls-Eos)) and stained with Alizarin Red S (ARS) to label the mineralized scale matrix. Note that atoh1a+ MCs localize to the epidermis above scales (arrowhead). (G) Lateral confocal micrograph of the scale epidermis in an adult expressing reporters for keratinocytes (Tg(krt4:DsRed)) and F-actin within atoh1a+ MCs (Tg(atoh1a:Lifeact-EGFP)). Note that all atoh1a+ MCs extend multiple microvilli. (G’) Magnification of G with arrowheads indicating individual microvillar processes on the surface of MCs. (G’’) Reconstructed cross section along the yellow line in G. MCs localize to the upper epidermal strata as diagrammed in D. Note that Tg(krt4:DsRed) (blue) preferentially labels keratinocytes in the upper epidermal strata, but not in the basal cell layer. Scale bars: 1 µm (B and C), 0.1 µm (B’–B’’’), 0.5 µm (C’), 50 µm (E–E’’ and F), 10 µm (G), and 5 µm (G’ and G’’).

atoh1a reporters label MCs in the adult epidermis

To date, a lack of genetically encoded reagents has hindered in-depth study of anamniote MCs. Since the TEM studies of Whitear and colleagues decades ago, molecular markers have been identified that distinguish mammalian MCs from other epidermal cells. For example, keratins, most notably keratin 8 and keratin 20, have been used extensively as markers of mammalian MCs (Moll et al., 1995; Moll et al., 1984). However, teleost keratins have undergone extensive gene loss and duplication and are not orthologous to mammalian keratin genes (Ho et al., 2022). Thus, we considered alternative molecular markers to label zebrafish MCs. Expression of Atoh1 uniquely identifies MCs in rodent skin and is necessary and sufficient for MC development (Morrison et al., 2009; Ostrowski et al., 2015; Van Keymeulen et al., 2009). The zebrafish genome contains three genes (atoh1a, atoh1b, and atoh1c) encoding Atoh1 homologs (Chaplin et al., 2010; Kani et al., 2010). To determine if the adult epidermis contained cells expressing an Atoh1 homolog, we focused on characterizing the expression pattern of atoh1a due to the availability of an enhancer trap line that expresses a nuclear localized version of the photoconvertible fluorescent protein Eos (nls-Eos) from the endogenous atoh1a locus (Tg(atoh1a:nls-Eos); Pickett et al., 2018). Confocal imaging of the adult trunk revealed that Tg(atoh1a:nls-Eos) labeled hair cells of the posterior lateral line, which formed tight clusters within neuromasts in interscale regions (Figure 1E and E’). In addition to atoh1a+ cells of the lateral line, we identified a second, spatially distinct population of atoh1a+ cells dispersed across the scale surface (Figure 1E, E’’ and F). Reconstructed cross sections showed that this population of atoh1a+ cells resided within the epidermis above scales (Figure 1F), in a similar axial position to the cells we identified by TEM.

The numerous, actin-rich microvilli that emanate from the MC surface morphologically distinguish them from other epidermal cells (Lane and Whitear, 1977; Toyoshima et al., 1998; Yamashita et al., 1993). To determine whether the dispersed epidermal atoh1a+ cells extended microvilli, we created an atoh1a enhancer trap line that expresses Lifeact-EGFP, a reporter for filamentous actin (Riedl et al., 2008). Similar to Tg(atoh1a:nls-Eos), Tg(atoh1a:Lifeact-EGFP) labeled hair cells of the lateral line and inner ear in larvae (Figure 1—figure supplement 2). atoh1a+ cells were notably absent from regions above the larval eye, yolk sac, or caudal fin (Figure 1—figure supplement 2), where neuroepithelial cells (NECs), a morphologically distinct population of sensory cells, have been described in larval skin (Coccimiglio and Jonz, 2012). Confocal microscopy of the scale epidermis in Tg(atoh1a:Lifeact-EGFP) adults revealed actin-rich microvilli densely decorating atoh1a+ cells in close proximity to neighboring keratinocytes (Figure 1G), further suggesting that the epidermal atoh1a+ cell population in the trunk shared key characteristics with the candidate MCs identified by TEM. Immunostaining for Sox2, a transcription factor required for murine MC maturation (Bardot et al., 2013; Perdigoto et al., 2014), demonstrated that the epidermal atoh1a+ cells expressed Sox2 (Figure 1—figure supplement 3). Together, these results define molecular and cellular properties of a previously uncharacterized epidermal cell population in zebrafish and identify genetic reagents for the study of this cell type. Anticipating the conclusion of our analysis below, we shall hereafter refer to the epidermal atoh1a+ cells as MCs, with the majority of the analyses completed on trunk MCs unless stated otherwise.

Somatosensory axons innervate zebrafish MCs, which display neurosecretory and mechanosensory characteristics

We next sought to determine whether zebrafish MCs displayed other key characteristics of MCs defined in mammals, including innervation by somatosensory axons and expression of neurosecretory and mechanosensory machinery.

Our ultrastructural observations suggested that cutaneous axons innervate MCs (Figure 1B). Staining scales with zn-12, a monoclonal antibody that labels several types of peripheral axons (Metcalfe et al., 1990), revealed that >90% of MCs were tightly associated with axons (Figure 2A and C). To determine the type of axon(s) innervating MCs, we examined expression of genetically encoded somatosensory axon reporters that we previously characterized in adult scales (Rasmussen et al., 2018). Analysis of reporters for three somatosensory neuron-expressed genes (p2rx3a, p2rx3b, and trpa1b; Kucenas et al., 2006; Palanca et al., 2013; Pan et al., 2012) demonstrated that somatosensory axons innervated up to 99% of MCs (Figure 2B and C). Consistent with ultrastructural analyses of MCs in the skin of other teleosts (Whitear, 1989), some axons formed ring-like structures that wrapped around MCs with MC-axon contacts containing varicosities or swellings (Figure 2B, inset and Figure 2D-F; Video 1). Additionally, we observed examples of axons forming both bouton- and en passant-like contacts with MCs (Figure 2G and H; Figure 1D).

Figure 2
Download asset Open asset
Somatosensory axons innervate Merkel cells (MCs) in the adult epidermis.

(A) Lateral confocal micrograph of the scale epidermis from an adult expressing an MC reporter immunostained for peripheral axons (zn-12). (B) Lateral confocal micrograph of the scale epidermis showing that somatosensory peripheral axons (Tg(p2rx3b:EGFP)) innervate MCs. Inset of dotted region shows axonal varicosities adjacent to an MC (arrowheads). (C) Quantification of MC innervation in the scale epidermis (17–30 mm standard length [SL]). Each dot represents measurements from an individual scale. Innervation frequencies: zn-12, 91% (284/311 cells; N=3 adults); Tg(p2rx3a>mCherry), 86% (196/228 cells; N=4 adults); Tg(p2rx3b:EGFP), 99% (225/228 cells; N=4 adults); Tg(trpa1b:EGFP), 96% (217/225 cells; N=9 adults). Error bars represent 95% CIs. (D–F) High-magnification confocal micrographs showing examples of somatosensory axons forming extended, ring-like contacts with MCs within the scale epidermis. (G) Three-dimensional (3D) reconstruction of an axon (zn-12 immunostaining, arrowheads) forming a bouton-like ending (asterisk) that terminates in close proximity to an MC. DAPI staining labels epidermal nuclei. (H) 3D reconstruction of a single somatosensory axon (Tg(p2rx3a>mCherry)) that forms en passant-like contacts (asterisks) with multiple MCs. Scale bars: 10 μm (A and B), 5 μm (D–H).

Figure 2—source code 1

ImageJ macro used for Figure 2C.

https://cdn.elifesciences.org/articles/85800/elife-85800-fig2-code1-v2.zip
Figure 2—source data 1

Datasheet for Figure 2C.

https://cdn.elifesciences.org/articles/85800/elife-85800-fig2-data1-v2.xlsx
Video 1
Download asset
Three-dimensional (3D) reconstruction of somatosensory axon and MC interactions.

3D rotation of somatosensory axons (green) and photoconverted MCs (green and magenta) in the adult scale epidermis. Arrows indicate axonal varicosities in close proximity to MCs. Scale bar: 4 µm.

Based on our observation that MCs contained cytoplasmic granules (Figure 1B), we postulated that they would display neurosecretory characteristics. We began by staining scales with an antibody against synaptic vesicle glycoprotein 2 (SV2), which labels secretory vesicle membranes in zebrafish and mammalian neuronal and endocrine cells (Buckley and Kelly, 1985; Jonz and Nurse, 2003). Essentially all MCs contained SV2-positive structures (Figure 3A), suggesting they express neurosecretory machinery that may contain neurotransmitter(s). Indeed, immunostaining revealed that MCs expressed serotonin (5-hydroxytryptamine; 5-HT) (Figure 3B), similar to mammalian MCs (Chang et al., 2016; English et al., 1992; García-Caballero et al., 1989) and zebrafish NECs (Coccimiglio and Jonz, 2012). Both 5-HT and SV2 appeared in a speckled pattern within MCs (Figure 3A and B), consistent with a vesicular localization.

Figure 3 with 1 supplement see all
Download asset Open asset
Merkel cells (MCs) in the adult epidermis express neurosecretory and mechanosensory machinery.

(A and B) Anti-SV2 (A–A’’) or anti-5-hydroxytryptamine (5-HT) (B–B’’) immunostaining of the scale epidermis from an adult expressing an MC reporter. Insets of dotted regions show the punctate localization of SV2 and 5-HT staining in MCs (arrowheads), consistent with a vesicular localization. 96% of MCs (178/185) were SV2+. 98% of MCs (326/332) were 5-HT+. Cells analyzed from n=3 scales from N=2 adults (25–27 mm standard length [SL]). DAPI labels epidermal nuclei. (C) Scale epidermis from an adult expressing an MC reporter stained with AM1-43. 98% of MCs (90/92) were AM1-43+. Cells analyzed from n=6 scales from N=2 adults. Inset of dotted region shows puncta within an MC labeled by AM1-43 (arrowheads). AM1-43 has been reported to stain neurites innervating MCs in murine whisker vibrissae (Meyers et al., 2003). However, our AM1-43 staining regimen did not strongly label cutaneous axons, although we cannot exclude low levels of staining. (D) Scale epidermis from an adult expressing an MC reporter stained with hybridization chain reaction (HCR) probes against piezo2 and an anti-GFP antibody. 99% of MCs (310/314) were piezo2+. Cells analyzed from n=7 scales from N=2 adults. Arrowheads indicate examples of positive staining within an MC. Scale bars: 5µm.

Do zebrafish MCs exhibit properties consistent with mechanosensory function? To address this question, we began by staining scales with AM1-43, an activity-dependent fluorescent styryl dye that labels a variety of sensory cells, including mammalian MCs (Meyers et al., 2003) and hair cells of the zebrafish lateral line (Corey et al., 2004). Following a short preincubation, AM1-43 robustly stained MC membranes and punctate structures reminiscent of vesicular compartments (Figure 3C), suggestive of ion channel expression in MCs (Meyers et al., 2003). Mammalian MCs express the mechanically activated cation channel Piezo2, which is required for MC mechanosensory responses (Ikeda et al., 2014; Maksimovic et al., 2014; Woo et al., 2014). Hybridization chain reaction (HCR) with an antisense probeset against piezo2 labeled MCs in adult scales (Figure 3D). We confirmed this staining pattern by performing fluorescent in situ hybridization (FISH) with a previously described piezo2 probe (Faucherre et al., 2013; Figure 3—figure supplement 1). Together, these data suggest that somatosensory peripheral axons innervate adult MCs, which possess neurosecretory and mechanosensory properties.

MCs arise from basal keratinocyte precursors

What are the precursors of MCs in zebrafish? Analysis of MC progenitors have come to conflicting results in avians and rodents: quail-chick chimeras suggest a neural crest origin for avian MC (Grim and Halata, 2000), whereas Cre-based lineage tracing studies in mouse demonstrate an epidermal origin (Morrison et al., 2009; Van Keymeulen et al., 2009).

To investigate a possible neural crest origin, we crossed a Cre driver expressed in neural crest progenitors (Tg(sox10:Cre); Kague et al., 2012) to a reporter transgene that stably expresses DsRed upon Cre-mediated recombination from a quasi-ubiquitous promoter (Tg(actb2:LOXP-BFP-LOXP-DsRed); Kobayashi et al., 2014; Figure 4A). DsRed+ neural crest-derived cell types, such as Schwann cells, appeared along scales, indicative of successful recombination (Figure 4B). However, we observed <0.5% co-labeling between the neural crest lineage trace and an MC reporter (Figure 4B’ and F). Based on these results, we concluded that zebrafish MCs likely derive from a non-neural crest lineage.

Figure 4 with 1 supplement see all
Download asset Open asset
Merkel cells (MCs) derive from the basal keratinocyte lineage.

(A) Schematic of Cre-based neural crest lineage tracing strategy. (B) Confocal micrograph of the scale epidermis in an adult expressing neural crest lineage (Tg(sox10:Cre); Tg(actb2:LOXP-BFP-LOXP-DsRed)) and MC (Tg(atoh1a:Lifeact-EGFP)) reporters. Brackets denote Schwann cells associated with a nerve along a scale radius. (C) Schematic of Cre-based basal keratinocyte lineage tracing strategy. (D) Confocal micrograph of the scale epidermis in an adult expressing basal keratinocyte lineage (TgBAC(ΔNp63:Cre-ERT2); Tg(actb2:LOXP-BFP-LOXP-DsRed)) and MC (Tg(atoh1a:Lifeact-EGFP)) reporters, which was treated with 4-hydroxytamoxifen (4-OHT) at 1 day post-fertilization (dpf). Arrowheads indicate MCs labeled by the basal keratinocyte lineage reporter. Note that recombination is not complete, possibly explaining why not all MCs express the lineage reporter. (E) Workflow to calculate percentage of MCs expressing lineage reporter and percentage of total cells expressing lineage reporter. (F) Boxplots of the percentage of MCs expressing the lineage tracing reporters diagrammed in panels A and C. Each dot represents an individual scale. Overall percentage of MCs expressing lineage trace reporters: sox10/Lifeact, 0.3% (1/323 cells; N=6 adults, 27.5–31 mm standard length [SL]); ΔNp63/Lifeact, 29.7% (299/1005 cells; N=6 adults, 21–26 mm SL); ΔNp63/nls-Eos, 32.3% (386/1195 cells; N=4 adults, 20–30 mm SL). A one-way ANOVA (F=12.06; p<0.001) with Tukey’s post-hoc honestly significant difference (HSD) test was used to compare groups. **, p<0.01; ***, p<0.001. (G and H) Paired dot plots of the percentage of MCs expressing the indicated atoh1a reporter and the basal keratinocyte lineage reporter compared to the percentage of all cells in the field of view expressing DsRed. Statistical analyses were performed using the Wilcoxon test. Scale bars: 20 µm.

Figure 4—source data 1

Datasheet for Figure 4F.

https://cdn.elifesciences.org/articles/85800/elife-85800-fig4-data1-v2.xlsx
Figure 4—source data 2

Datasheet for Figure 4G and H.

https://cdn.elifesciences.org/articles/85800/elife-85800-fig4-data2-v2.xlsx

To investigate a possible epidermal origin, we considered basal keratinocytes, an epidermal-resident stem cell population, the most likely candidate progenitors. To follow this lineage, we engineered a transgene to express a tamoxifen-inducible Cre recombinase from regulatory sequences of ΔNp63 (TgBAC(ΔNp63:Cre-ERT2)), a basal keratinocyte marker (Bakkers et al., 2002; Lee and Kimelman, 2002). We crossed this transgene to the Cre reporter transgene and treated embryos with 4-hydroxytamoxifen (4-OHT) at 1 day post-fertilization (dpf) to induce Cre-ERT2 activity, which resulted in permanent DsRed expression in basal keratinocytes and their derivatives (Figure 4C and D; Figure 4—figure supplement 1). After raising 4-OHT-treated animals to adulthood, we observed variable (2–81%) co-labeling between the basal keratinocyte lineage trace and MC reporters (Figure 4D’ and F). We note that our lineage tracing strategy did not label all basal keratinocytes (Figure 4D; Figure 4—figure supplement 1), suggestive of incomplete Cre-ERT2 induction and/or transgene recombination. Consistent with the latter possibility, a recent analysis demonstrated Tg(actb2:LOXP-BFP-LOXP-DsRed) has a low recombination efficiency compared to other Cre reporter transgenes (Lalonde et al., 2022). To estimate the local recombination efficiency in imaged regions, we thresholded the DsRed channel and calculated the fraction of skin cells labeled (Figure 4E). Importantly, the proportion of MCs labeled by the basal keratinocyte lineage trace was not significantly different from the local recombination efficiency (Figure 4G–H). These observations support a basal keratinocyte origin of most or all zebrafish MCs.

MCs continuously turn over in adult skin

The longevity and turnover of murine MCs are controversial. Several studies concluded that MC numbers fluctuate with hair cycle stages (Marshall et al., 2016; Moll et al., 1996; Nakafusa et al., 2006), while another found no correlation (Wright et al., 2017). To determine the turnover rate of zebrafish MCs, we photoconverted small regions of the scale epidermis in Tg(atoh1a:nls-Eos) adults and tracked individual cells over time. Exposure to UV light irreversibly photoconverts nls-Eos, allowing us to distinguish pre-existing cells (containing photoconverted nls-Eos) from newly added cells (without photoconverted nls-Eos; Figure 5A–C). By longitudinally tracking individual fish over the course of 28 days, we found a decrease of ~15% of the photoconverted MCs every 7 days (Figure 5D). In addition to the gradual loss of MCs over time, we noted a steady addition of new MCs (Figure 5C’), resulting in a nearly constant total cell number (Figure 5E). Thus, MCs undergo constant cell loss and renewal in adult skin, albeit at a slower rate than atoh1a-expressing hair cells of the lateral line (Cruz et al., 2015).

Figure 5
Download asset Open asset
Homeostatic replacement of Merkel cells (MCs) in the adult epidermis.

(A) Illustration of the photoconversion experiment showing the epidermis (blue), non-converted MCs (green), and converted MCs (magenta) after exposure of a region of the scale epidermis to UV light. (B and C) Representative images of MCs labeled by Tg(atoh1a:nls-Eos) at 0 (B) or 7 (C) days post-conversion (dpc) from a single adult. Cyan dotted box indicates the photoconverted region. White dotted box indicates the area magnified in B’ and C’. (B’ and C) Numbers label examples of individual cells present at 0 and 7 dpc. Arrows indicate examples of newly added cells, which appear green due to the presence of non-converted nls-Eos (green) and absence of converted nls-Eos (magenta). (D and E) Boxplots of the percentage of photoconverted MCs remaining compared to 0 dpc (D) and the total number of MCs (converted+non-converted) present at each day compared to 0 dpc (E). Each dot represents an individual fish. N=5–8 fish (24–32 mm standard length [SL]). Scale bars: 50 μm.

Figure 5—source data 1

Datasheet for Figure 5D and E.

https://cdn.elifesciences.org/articles/85800/elife-85800-fig5-data1-v2.xlsx

MCs are widely distributed across the body, in compartment-specific patterns

MCs localize to specific regions of mammalian skin, such as in crescent-shaped touch domes adjacent to hair follicles in hairy skin and at the bottom of rete ridges in glabrous skin (Boot et al., 1992; Fradette et al., 1995; Iggo and Muir, 1969; Lacour et al., 1991). To determine the distribution pattern of zebrafish MCs, we used confocal microscopy to survey multiple regions of the adult skin. In addition to the MCs found on the trunk, MCs appeared in the epidermis above the eyes, gill covers (opercula), and fins (Figure 6A–F). While MC morphology was similar across the skin compartments (Figure 6G–J), MC densities and spatial distributions varied across skin compartments (Figure 6B). For example, MCs were distributed relatively uniformly across the cornea, although this spatial pattern was not specifically quantified (Figure 6C). By contrast, in the caudal fin, MCs localized to the epidermis above bony rays and in the medial regions of the interray epidermis between bony rays (Figure 6F). Along the trunk, MCs appeared in patches, similar to the pattern of dermal scales beneath the epidermis (Figure 6E). Altogether, our results demonstrate that MCs are widely distributed across the adult zebrafish skin and localize in specific patterns in each skin compartment.

Figure 6
Download asset Open asset
Merkel cells (MCs) are widely distributed across the skin, in compartment-specific patterns.

(A) Illustration indicating the epidermal regions imaged in adult zebrafish. (B) Quantification of MC densities in the specified regions. Each dot represents an individual fish (N=8–18, 20–29.5 mm standard length [SL]). *** indicates p<0.001 using a one-way ANOVA (F=83.94; p<0.001) with post-hoc Tukey’s HSD test. (C–J) Lateral confocal micrographs of MCs in the different skin regions from animals expressing the indicated reporters. The regions imaged are indicated in A. Note that MCs expressing Tg(atoh1a:Lifeact-EGFP) have a similar morphology across skin compartments (G–J). nm, neuromasts of the posterior lateral line. Scale bars: 50 μm (C–F) and 10 μm (G–J).

Figure 6—source data 1

Datasheet for Figure 6B.

https://cdn.elifesciences.org/articles/85800/elife-85800-fig6-data1-v2.xlsx

Trunk MCs develop concomitant with dermal appendage morphogenesis

To examine the mechanisms that generate a compartment-specific MC pattern, we focused on the trunk skin because of its molecular and cellular similarities to murine hairy skin (Aman et al., 2018; Harris et al., 2008). Both during ontogeny and at post-embryonic stages, murine MCs associate with primary (guard) hairs, a subclass of dermal appendages (Jenkins et al., 2019; Nguyen et al., 2018; Perdigoto et al., 2016). Based on these studies in mice, and our previous work showing that epidermal diversification and somatosensory remodeling coincides with scale development in zebrafish (Rasmussen et al., 2018), we postulated that MCs would appear during squamation (scale formation).

Zebrafish post-larval development is staged by standard length (SL) in millimeters (mm; Parichy et al., 2009). Squamation begins at ~9 mm SL (Figure 7A; Aman et al., 2018; Harris et al., 2008; Sire et al., 1997a). Using reporters that label MCs and scale-forming osteoblasts, we rarely observed MCs in the epidermis prior to 8 mm SL (Figure 7B and F). Between 8 and 10 mm SL, MCs appeared at a low density along the trunk (Figure 7C and F). MC density rapidly increased from 10 to 15 mm SL, a period of active scale growth (Figure 7D, E and F). The density and number of MCs positively correlated with scale area (Figure 7G and H), although this trend was less pronounced at stages less than 10 mm SL (Figure 7—figure supplement 1). These data indicate that MC development coincides with dermal appendage growth along the trunk.

Figure 7 with 1 supplement see all
Download asset Open asset
Merkel cells (MCs) develop concomitant with dermal appendage morphogenesis.

(A) Abbreviated zebrafish developmental timeline relative to standard length (SL) in millimeters. Developing scales are drawn in magenta below the approximate corresponding stage. (B–E) Representative lateral confocal micrographs of MCs and osteoblasts along the trunk at the indicated stages. Note that MCs increase in number and density as scale-forming osteoblasts develop below the epidermis. nm, neuromast of the posterior lateral line. (F) Quantification of MC density according to SL. Each dot represents an individual fish. Data represent n=81 scales from N=52 fish. Line indicates segmented linear regression (breakpoint = 8.27 mm SL). (G and H) Quantification of the number (G) or density (H) of MCs relative to scale area. Each dot represents an individual scale. Data represent n=62 scales from N=18 fish. Dot colors represent animal SL as indicated in the legend. Shading indicates a 95% CI around the linear regression lines in G and H. Correlation coefficients (R2): 0.08 (F, slope 1), 0.68 (F, slope 2), 0.73 (G), and 0.31 (H). F-statistics: 3.5 (F, slope 1), 83.9 (F, slope 2),164.6 (G), and 28.31 (H). p-values: 0.07 (F, slope 1), <0.05 (F, slope 2), and <0.05 (G and H). Scale bars: 50 μm.

Figure 7—source data 1

Datasheet for Figure 7F.

https://cdn.elifesciences.org/articles/85800/elife-85800-fig7-data1-v2.xlsx
Figure 7—source data 2

Datasheet for Figure 7G and H.

https://cdn.elifesciences.org/articles/85800/elife-85800-fig7-data2-v2.xlsx

Ectodysplasin signaling promotes trunk MC development

Since the appearance of MCs in the trunk epidermis tightly correlated with scale growth, we examined the consequences of blocking signals required for dermal appendage morphogenesis on MC development. Ectodysplasin (Eda) signaling regulates the formation of many types of skin appendages, including mammalian hair follicles and zebrafish scales (Biggs and Mikkola, 2014; Harris et al., 2008). To determine whether MC development required Eda-dependent signals, we measured MC density in animals homozygous for a presumptive null allele of eda that does not develop scales (edadt1261/dt1261; hereafter eda–/–; Harris et al., 2008). Immediately prior to squamation, we found that there was no difference in MC density between eda mutants and sibling controls (Figure 8A, B and G). However, after the onset of squamation, eda mutants had significantly fewer MCs, a difference that persisted into adulthood (Figure 8C–G). In addition to the decrease in cell density, we observed a dramatic change in the spatial distribution of MCs across the epidermis in eda mutants compared to controls (Figure 8H). Specifically, in siblings, MCs appeared in patches corresponding to the location of the underlying scales (Figure 8C and H). By contrast, the MCs that developed in eda mutants were distributed relatively uniformly across the trunk (Figure 8D and H). Although we found a decrease in MC density in the trunk skin of the mutants, we observed no change in MC density in the corneal or opercular epidermis (Figure 8I), suggesting that the reduced density was specific to the trunk skin. eda mutants lack fins at the stages analyzed (Harris et al., 2008), precluding analysis of these regions in the homozygous mutants.

Figure 8 with 1 supplement see all
Download asset Open asset
Loss of Eda signaling decreases Merkel cell (MC) density in trunk, but not facial skin.

(A–F) Representative confocal images of MCs in the trunk of animals of the indicated genotypes at the indicated stages. Dotted yellow lines indicate posterior scale boundaries. nm, neuromasts of the posterior lateral line. (G) Quantification of MC density in the trunk skin relative to standard length (SL). Gray shading indicates a 95% CI around the linear regression lines. The difference between genotypes was significant above 12.5 mm SL (p<0.05, Johnson-Neyman Technique). Each dot represents an individual fish (N=16–18 fish/genotype). (H) Histograms of the distribution of trunk MCs along a rectangular segment encompassing three scales in a sibling and an identically sized region in an eda mutant (18–19 mm SL). (I) Boxplots of MC densities in the epidermis above the cornea or operculum in animals of the indicated genotypes. ns, not significant (cornea, p=0.21; operculum, p=0.14; Mann-Whitney test). Scale bars: 50 μm (A–F).

Figure 8—source data 1

Datasheet for Figure 8G.

https://cdn.elifesciences.org/articles/85800/elife-85800-fig8-data1-v2.xlsx
Figure 8—source data 2

Datasheet for Figure 8I.

https://cdn.elifesciences.org/articles/85800/elife-85800-fig8-data2-v2.xlsx

The decreased MC density in eda mutant trunk skin could be due to decreased cell addition, increased cell loss, or a combination of the two. Using in vivo photoconversion, we found that the rate of MC addition was significantly reduced in eda mutants compared to siblings (Figure 8—figure supplement 1A–D). Additionally, the rate of cell loss was higher in mutants compared to siblings, although this change was not statistically significant (Figure 8—figure supplement 1E). Thus, our observations indicate that the decrease in MC cell density in eda mutants is mainly due to reduced MC production. Together, these data suggest that Eda signaling, either directly or indirectly, is required for MC development, maintenance, and distribution along the trunk.

MC patterning is not predetermined along the trunk

Since blocking dermal appendage formation through loss of Eda signaling inhibited MC development, we next examined the consequences of altering dermal appendage size and shape on MC patterning. Zebrafish scale morphogenesis is regulated by Fibroblast growth factor (FGF) signaling (Aman et al., 2018; Daane et al., 2016; De Simone et al., 2021; Rohner et al., 2009). To determine whether alterations to scale patterning impacted MC development, we examined animals heterozygous for an allele of hagoromo (hag; fgf8adhiD1Tg/+), which results in fgf8a overexpression in the post-embryonic skin due to a viral insertion near the fgf8a locus (Amsterdam et al., 2009). An independent allele of hag (fgf8adhi4000Tg/+) was previously shown to result in large, disorganized sheets of scale-forming osteoblasts during squamation (Aman et al., 2018). fgf8adhiD1Tg/+ juveniles showed dramatic variability in scale size and shape, with both smaller and larger scales compared to the relatively uniformly patterned scales observed in sibling controls (Figure 9A–D; Figure 9—figure supplement 1A–C). We found no significant differences in MC density between the genotypes (Figure 9—figure supplement 1D,E). Nevertheless, the distribution of MCs tracked with the altered scale size and shape in the mutants (Figure 9E–H), suggesting the MC pattern is not predetermined within the trunk skin compartment.

Figure 9 with 1 supplement see all
Download asset Open asset
Merkel cell (MC) patterning is not predetermined along the trunk.

(A–D) Representative images of juvenile animals of the indicated genotypes expressing an MC reporter and stained with Alizarin Red S (ARS) to visualize scales (A and B) or co-expressing MC and keratinocyte (Tg(krt4:DsRed)) reporters (C and D). Dotted lines indicate scale boundaries. nm, neuromasts of the posterior lateral line. (E–H) Tracings of scale outlines (E and F) and density plots of MC position (G and H) from juvenile animals (n=43–49 scales/genotype from N=10–13 fish/genotype; 11.6–14.7 mm standard length [SL]) of the indicated genotypes. Scale tracings were aligned at the dorsal-ventral midpoint of the posterior scale margin. Note the variability in scale shape and size and corresponding increased spread of MC position in fgf8adhiD1Tg/+ juveniles compared to sibling controls. Scale bars: 100 μm (A–D) and 200 μm (E–H).

Figure 9—source data 1

Datasheet for Figure 9G, H, Figure 9—figure supplement 1E.

https://cdn.elifesciences.org/articles/85800/elife-85800-fig9-data1-v2.xlsx

Discussion

Here, we discover a zebrafish epidermal cell type that we classify as an MC based on ultrastructural criteria (Whitear, 1989). We further present several lines of evidence that suggest zebrafish MCs share molecular, cellular, and lineage properties with mammalian MCs. First, we show that zebrafish MCs express the transcription factors Atoh1a and Sox2, the orthologs of which uniquely mark MCs in mammalian skin (Maricich et al., 2009; Bardot et al., 2013; Perdigoto et al., 2014; Ostrowski et al., 2015; Van Keymeulen et al., 2009). Second, zebrafish MCs extend numerous short, actin-rich microvilli and complex with somatosensory axons, classic morphological hallmarks of MCs (Mihara et al., 1979; Smith, 1977; Toyoshima et al., 1998). Our morphological observations support the interpretation that these cells are MCs rather than Merkel-like cells, which lack axon association and microvillar processes (reviewed by Halata et al., 2003). Third, Cre-based lineage tracing revealed that basal keratinocytes give rise to zebrafish MCs, akin to studies in mouse (Morrison et al., 2009; Van Keymeulen et al., 2009). Fourth, we demonstrate that zebrafish MCs contain neurosecretory machinery and express the neurotransmitter serotonin, the release of which has been proposed to regulate somatosensory responses to touch (Chang et al., 2016; Chang and Gu, 2020; English et al., 1992). Finally, we show that zebrafish MCs express the cation channel Piezo2, which is cell-autonomously required for MC mechanosensory function (Ikeda et al., 2014; Maksimovic et al., 2014; Woo et al., 2014). Importantly, our results extend previous histological studies of MCs in various teleost fish (Lane and Whitear, 1977; Whitear, 1989; Zachar and Jonz, 2012) by identifying the first genetically encoded reagents for the study of this cell type in zebrafish.

While our characterization revealed substantial similarities between mammalian and zebrafish MCs, we did observe anatomical differences in line with previous ultrastructural characterizations of teleost MCs (Lane and Whitear, 1977; Whitear, 1989). For example, the nuclei of mammalian MC are commonly lobulated (Boulais et al., 2009; Cheng Chew and Leung, 1994; Moll et al., 2005; Tachibana and Nawa, 2002). While we did not observe lobulation of zebrafish MC nuclei by TEM, we cannot rule out that serial sectioning or high-resolution reconstruction of nuclear shape would reveal lobulation. Mammalian MCs typically localize adjacent to basal keratinocytes (Boot et al., 1992; Cheng Chew and Leung, 1994; Fradette et al., 1995; Mihara et al., 1979; Moll et al., 1996; Smith, 1977), whereas zebrafish MCs appear in upper strata, typically beneath the periderm (Figure 1D and G’’). As the majority of the analyses completed here focused on MCs found in the trunk epidermis, it will be intriguing to determine whether all MCs in different skin compartments in the juvenile and adult zebrafish share similar properties.

Teleost MCs and somatosensory physiology

In mammalian skin, the MC-neurite complex regulates slowly adapting type I responses to light touch (Iggo and Muir, 1969; Ikeda et al., 2014; Maksimovic et al., 2014; Maricich et al., 2009; Woo et al., 2014). Although physiological studies of somatosensory responses in adult zebrafish have not been reported, extracellular recordings in adult rainbow trout demonstrate that a subset of somatosensory neurons exhibit slowly adapting responses to mechanical skin stimulation (Ashley et al., 2007; Ashley et al., 2006; Sneddon, 2003). We postulate that the slowly adapting responses to mechanical skin stimulation in adults require MCs. Nevertheless, the exact physiological roles of teleost MCs in regulating somatosensory responses and resulting behaviors remain unknown and will require the development of tools to selectively ablate and activate MCs. Interestingly, recordings from zebrafish Rohon-Beard neurons, a transient larval somatosensory population, suggest they have rapidly, but not slowly, adapting mechanosensory responses (Katz et al., 2021). Together these studies correlate with our finding that MCs develop at post-larval stages and suggest that the teleost somatosensory system undergoes significant functional maturation during the juvenile period.

What are the subtypes of somatosensory neurons in fish, and how do they correspond to MC innervation? Several studies have identified molecularly distinct subsets of somatosensory neurons through mRNA, protein, and transgene expression analysis in zebrafish larvae (Faucherre et al., 2013; Gau et al., 2017; Gau et al., 2013; Kucenas et al., 2006; Palanca et al., 2013; Pan et al., 2012; Patten et al., 2007; Slatter et al., 2005). Adult trout somatosensory neurons have been classified based on their responses to mechanical, chemical, and thermal stimuli (Ashley et al., 2007; Ashley et al., 2006; Sneddon, 2003). However, to date, a detailed molecular characterization of the diversity of somatosensory subtypes present in adult fish has not been performed. Our data suggest that somatosensory neurons expressing reporters for p2rx3a, p2rx3b, or trpa1b innervate MCs. Whether these neurons represent a dedicated class of MC-innervating neurons remains unknown. The development of Cre drivers for specific somatosensory subtypes (Bai et al., 2015; Li et al., 2011; Luo et al., 2009; Rutlin et al., 2014; Zylka et al., 2005) and single-cell transcriptional profiling (Sharma et al., 2020; Usoskin et al., 2015; Zeisel et al., 2018) has been fruitful in characterizing the diversity of somatosensory neurons in mammals. The application of these technologies to the teleost somatosensory system is an interesting avenue for further investigation.

MC lineage and homeostasis

The developmental lineage of MCs has been a long-standing question with both epidermal and neural crest origins posited (Hartschuh et al., 1986). Our Cre-based lineage tracing identified basal keratinocytes as MC progenitors. These results extend previous studies in the zebrafish epidermis showing that basal keratinocytes serve as precursors for diverse post-larval cell types, including periderm and immune cells (Lee et al., 2014; Lin et al., 2019). Although previous work in mouse unambiguously identified keratin 14-expressing basal keratinocytes as MC precursors (Morrison et al., 2009; Van Keymeulen et al., 2009), the precise nature of murine MC progenitors varies across skin compartments (Nguyen et al., 2019). Future studies characterizing the molecular properties and cellular behaviors of zebrafish MC precursors will be informative for identifying conserved properties of skin stem cells.

The turnover of MCs in mammalian skin has been a source of controversy. Several studies reported that MC numbers fluctuate with the natural hair cycle in mouse (Marshall et al., 2016; Moll et al., 1996; Nakafusa et al., 2006). By contrast, Wright et al., 2017 found no evidence for changes in MC density based on stages of the hair cycle and demonstrated that MCs could live for months. These types of analyses in murine skin have relied either on histology, which limits tissue sampling, or required the use of advanced (2-photon) microscopy in combination with hair shaving, a mild form of skin injury. Using photoconversion and confocal imaging, we non-invasively tracked individual MCs during normal skin homeostasis in vivo for weeks. We found that trunk MCs have a steady turnover in adult animals, with a half-life of approximately 1 month. Additionally, this further distinguishes MCs from hair cells in the adult lateral line, which have a shorter half-life (Cruz et al., 2015). Whether MC turnover varies at different stages of development, across skin compartments, or following skin insults will require further study.

MC distribution and patterning

Regionally specific sensory structures allow our skin to distinguish tactile inputs with remarkable acuity (Corniani and Saal, 2020). For example, MC densities vary greatly across human skin compartments, with the highest numbers found in particularly sensitive regions such as fingertips and lips (Boot et al., 1992; Lacour et al., 1991). We observed that MCs populate several major skin compartments and have regional-specific densities in adult zebrafish, with the highest densities found in the face (corneal and opercular epidermis). We speculate this may bestow the juvenile and adult skin with the ability to detect innocuous tactile inputs across almost the entire body surface, with perhaps the greatest sensitivity along facial structures.

Although most studies of MC development have centered on the formation of MC aggregates in the touch domes of murine hairy skin, MCs are found in a range of distribution patterns in other types of skin. For example, MCs are found as dispersed, single cells arrayed across the skin of human toe pads (Boot et al., 1992). Similarly, we found that MCs have a dispersed, rather than clustered, pattern in all skin compartments examined. Few studies have addressed how MCs adopt specific distributions, and zebrafish presents a promising model to understand mechanisms of MC pattern formation in vivo.

Dermal appendages and MC development

Our developmental analysis showed that trunk MC density rapidly increases during dermal appendage morphogenesis. Previous genetic analysis in mouse hairy skin revealed that MC development requires Eda signaling (Vielkind et al., 1995; Xiao et al., 2016). We show that zebrafish eda mutants have decreased MC density in the trunk, but not the facial, skin. These observations suggest that MC development in mouse and zebrafish likely share similar genetic pathways, akin to the shared molecular and cellular mechanisms that regulate dermal appendage formation (Aman et al., 2018; Biggs and Mikkola, 2014; Daane et al., 2016; Harris et al., 2008; Rohner et al., 2009). They further support a model whereby MC development requires compartment-specific signals, akin to recent observations on MC development in mouse hairy and glabrous skin (Nguyen et al., 2019). By taking advantage of the ability to image large skin areas in intact zebrafish, we show that eda mutants have altered MC distribution compared to controls. Furthermore, we use in vivo photoconversion to demonstrate the reduction in MC density is largely due to decreased production but also reflects slightly increased cell loss. Further investigations are required to determine whether Eda signaling directly regulates the differentiation of MC progenitors. Alternatively, since eda mutants lack scales (Harris et al., 2008) and have decreased epidermal innervation (Rasmussen et al., 2018), MC development may require scale- and/or somatosensory neuron-derived signals. Finally, we note that trunk MCs are not completely absent in eda mutants, suggesting that a subset of MCs develop independent of Eda signaling.

We found that a gain-of-function allele of fgf8a leads to a change in the overall size and shape of scales. Intriguingly, the MC distribution modifies to accommodate the altered scale size and shape in the fgf8a mutants but retains the same MC density as sibling controls. This result suggests that the number of MCs per scale is not predetermined but rather is titrated relative to appendage size. How are MCs able to populate the much larger scales in fgf8a mutants? Does the size of the MC progenitor domain expand with increases in scale size? Are MCs, or their progenitors, able to migrate to their final destination? Distinguishing between these possibilities will require tracking the behaviors of MCs and their progenitors in vivo.

Summary

Our results establish a promising new system to investigate MC biology. This model will allow for the identification of deeply conserved mechanisms used to regulate vertebrate MC biology. Furthermore, the advantages of zebrafish—such as non-invasive in vivo imaging, genetic and chemical screens, and high regenerative capacity—will complement the strengths of existing rodent models. Specifically, the ability to track individual cells over time has the potential to answer key and long-standing questions surrounding MC biology, including how the MC-neurite relationship is established, how MCs interact with neighboring cell types, and their progenitor dynamics. Addressing these questions, as well as potential novel insights provided by the zebrafish system, represent exciting directions for future research.

Materials and methods

Key resources table
Reagent type (species) or resourceDesignationSource or referenceIdentifiersAdditional information
AntibodyAnti-Serotonin (rabbit polyclonal)MilliporeSigmaCat #: S5545, RRID:AB_477522(1:1000)
AntibodyAnti-Sv2 (mouse monoclonal)DSHB; (Buckley and Kelly, 1985)Cat #: SV2, RRID:AB_2315387(1:50)
AntibodyAnti-Sox2 (rabbit polyclonal)GeneTexCat #: GTX124477, RRID:AB_11178063(1:500)
AntibodyAnti-Fluorescein Polyclonal Antibody, POD Conjugated (sheep polyclonal)RocheCat #: 11426346910, RRID:AB_840257(1:2000)
AntibodyAnti-GFP (rabbit polyclonal)Thermo Fisher ScientificCat #: A11122, RRID:AB_221569(1:1000)
Antibodyzn-12 (mouse monoclonal)Zebrafish International Resource CenterCat #: zn-12, RRID:AB_10013761(1:200)
AntibodyAnti-Rabbit Alexa Fluor 647 (goat polyclonal)Thermo Fisher ScientificCat #: A32733, RRID:AB_2633282(1:500)
AntibodyAnti-Mouse Alexa Fluor 647 (goat polyclonal)Thermo Fisher ScientificCat #: A32728, RRID:AB_2633277(1:500)
AntibodyAnti-Rabbit Alexa Fluor 568 (goat polyclonal)Thermo Fisher ScientificCat #: A-11036, RRID:AB_10563566(1:500)
AntibodyAnti-Mouse Alexa Fluor 568 (goat polyclonal)Thermo Fisher ScientificCat #: A-11031, RRID:AB_144696(1:500)
AntibodyAnti-Rabbit Alexa Fluor 488 (goat polyclonal)Thermo Fisher ScientificCat #: A32731, RRID:AB_2633280(1:500)
Commercial assay and kitTSA Plus Cyanine 5Akoya BiosciencesCat #: NEL705A001KT(1:50)
Sequence-based reagentpiezo2 in situ probe (originally referred to as piezo2b)Faucherre et al., 2013N/A
Sequence-based reagentatoh1a gRNA, 5’-GGA GAC TGA ATA AAG TTA TG-3’Pickett et al., 2018N/A
Sequence-based reagentMbait gRNA, 5’-GGC TGC TGC GGT TCC AGA GGT GG-3’Kimura et al., 2014N/A
Sequence-based reagentZebrafish piezo2 HCR v3.0 probeMolecular InstrumentsN/AUsed at 2 pmol
Chemical compound and drugMS-222MilliporeSigmaCat #: E10521
Chemical compound and drug(Z)–4-Hydroxytamoxifen (4-OHT)MilliporeSigmaCat #: H7904Used at 10 μM
Strain and strain background (Danio rerio)AB (Wild-Type)Zebrafish International Resource CenterZIRC Cat# ZL1, RRID:ZIRC_ZL1
Genetic reagent (Danio rerio)Tg(actb2:LOXP-BFP-LOXP-DsRed)Kobayashi et al., 2014Tg(actb2:LOXP-BFP-LOXP-DsRed)sd27Tg, ZFIN: ZDB-TGCONSTRCT-141111–5
Genetic reagent (Danio rerio)Tg(atoh1a:nls-Eos)Pickett et al., 2018Tg(atoh1a:nls-Eos)w214Tg, ZFIN: ZDB-TGCONSTRCT-190701–2
Genetic reagent (Danio rerio)Tg(atoh1a:lifeact-EGFP)This studyTg(atoh1a:lifeact-EGFP)w259Tg
Genetic reagent (Danio rerio)TgBAC(ΔNp63:Cre-ERT2)This studyTgBAC(ΔNp63:Cre-ERT2)w267Tg
Genetic reagent (Danio rerio)Tg(sox10:Cre)Kague et al., 2012Tg(Mmu.Sox10-Mmu.Fos:Cre)zf384, ZFIN: ZDB-TGCONSTRCT-130614–2
Genetic reagent (Danio rerio)Gt(ctnna-citrine)Trinh et al., 2011Gt(ctnna-citrine)ct3aGt, ZFIN: ZDB-ALT-111010–23
Genetic reagent (Danio rerio)Tg(sp7:mCherry)Singh et al., 2012Tg(Ola.Sp7:mCherry-Eco.NfsB)pd46Tg, ZFIN: ZDB-TGCONSTRCT-120503–4
Genetic reagent (Danio rerio)Tg(p2rx3a>mCherry)Palanca et al., 2013Tg(Tru.P2rx3a:LEXA-VP16,4xLEXOP-mCherry)la207Tg, ZFIN: ZDB-TGCONSTRCT-130307–1
Genetic reagent (Danio rerio)Tg(trpa1b:EGFP)Pan et al., 2012TgBAC(trpa1b:EGFP)a129Tg, ZFIN: ZDB-TGCONSTRCT-120208–2
Genetic reagent (Danio rerio)Tg(p2rx3b:EGFP)Kucenas et al., 2006Tg(p2rx3b:EGFP)sl1Tg, ZFIN: ZDB-TGCONSTRCT-070117–110
Genetic reagent (Danio rerio)Tg(krt4:DsRed)Rieger and Sagasti, 2011Tg(krt4:DsRed)la203Tg, ZFIN: ZDB-TGCONSTRCT-120127–5
Genetic reagent (Danio rerio)edadt1261Harris et al., 2008edadt1261, ZFIN: ZDB-ALT-090324–1
Genetic reagent (Danio rerio)fgf8adhiD1Tg/+Amsterdam et al., 2009fgf8adhiD1Tg/+, ZFIN: ZDB-ALT-010427–4
Software and algorithmFIJIhttp://fiji.scRRID:SCR_002285
Software and algorithmImarisBitplaneRRID:SCR_007370
OtherFetal bovine serumGibcoCat #: 10082–139
OtherNormal goat serumAbcamCat #: ab7481, RRID:AB_2716553
OtherDAPIMilliporeSigmaCat #: 508741Used at 5 ng/μl
OtherAM1-43BiotiniumCat #: 70024Used at 15 μm
OtherAlizarin Red SACROS OrganicsCat #: 400480250Used at 0.01%
OtherProteinase KThermo Fisher ScientificCat #: 100005393Used at 0.1 mg/ml
OtherHoechst 3342Thermo Fisher ScientificCat #: H3570Used at 5 ng/μl

Materials availability

Request a detailed protocol

Strains and DNA constructs generated in this study are available from the corresponding author upon request.

Animals

Zebrafish

Zebrafish were housed at 26–27°C on a 14/10 hr light cycle. See Key resources table for strains used in this study. Animals of either sex were used. Zebrafish were staged according to SL (Parichy et al., 2009). SL of fish was measured using the IC Measure software (The Imaging Source) on images captured on a Stemi 508 stereoscope (Zeiss) equipped with a DFK 33UX264 camera (The Imaging Source). All zebrafish experiments were approved by the Institutional Animal Care and Use Committee at the University of Washington (Protocol: #4439–01).

Creation of Tg(atoh1a:Lifeact-EGFP)

Request a detailed protocol

Tg(atoh1a:Lifeact-EGFP)w259Tg was generated by CRISPR-mediated knock-in as previously described (Kimura et al., 2014). A donor plasmid containing the Mbait, minimal hsp70l promoter, Lifeact-EGFP, and bgh poly(A) sequences was created using Gibson assembly. The insertion was targeted 372 bp upstream of the endogenous atoh1a coding sequence using a previously published guide RNA (gRNA; Pickett et al., 2018). The Mbait-hsp70l-Lifeact-EGFP plasmid, Mbait and atoh1a gRNAs, and Cas9 protein were prepared and injected into single-cell embryos of the AB strain as previously described (Thomas and Raible, 2019). Larvae were screened for Lifeact-EGFP expression at 3 dpf and raised to adulthood. A founder adult was identified and outcrossed to generate a stable transgenic line.

Creation of TgBAC(ΔNp63:Cre-ERT2) and induction with 4-OHT

Request a detailed protocol

The ΔNp63:Cre-ERT2 bacterial artificial chromosome (BAC) was created by modifying the previously generated BAC DKEY-263P13-iTol2-amp (Rasmussen et al., 2015). The predicted ΔNp63 start codon was replaced by a Cre-ERT2-pA-KanR cassette that contained a zebrafish codon-optimized Cre-ERT2 (Kesavan et al., 2018) using a previously described protocol (Suster et al., 2011). TgBAC(ΔNp63:Cre-ERT2)w267Tg was created by injecting tol2 mRNA, which was transcribed from pCS2-zT2TP (Suster et al., 2011), and BAC DNA into one-cell stage embryos and screening adults for germline transmission. To activate Cre-ERT2, 1 dpf embryos were treated with 10 μM 4-OHT for 24 hr. 4-OHT was prepared as described (Felker et al., 2016).

Mutant identification and analysis

Request a detailed protocol

eda mutants and siblings were sorted by visible phenotype starting at 7 mm SL. Mutants were grown separately from siblings. fgf8adhiD1Tg/+ fish were identified based on altered scale patterning and/or pigmentation (Kawakami et al., 2000).

Imaging and photoconversion

Electron microscopy

Request a detailed protocol

Isolated scales were prepared for TEM as described (Sire et al., 1997b), with the following modifications: after dehydration, scales were treated with propylene oxide (PO), infiltrated with PO:Eponate 12, and embedded in Eponate 12. Semithin sections (0.2 μm) stained with toluidine blue were used for orientation. Thin sections (50 nm) were placed on Formvar coated copper slot grids, stained with saturated uranyl acetate and Reynolds’ lead citrate, and examined on a JEOL 100 CX at 60 kV or a Philips CM100 at 80 kV.

Confocal image acquisition

Request a detailed protocol

Confocal z-stacks were collected using a A1R MP+ confocal scanhead mounted on an Ni-E upright microscope (Nikon) using a 16× water dipping objective (NA 0.8) for live imaging or 40× oil immersion objective (NA 1.3) for fixed image acquisition. Images acquired in resonant scanning mode were post-processed using the denoise.ai function in NIS-Elements (Nikon). For live imaging, zebrafish were anesthetized in a solution of 0.006–0.012%buffered MS-222 in system water for 5 min. Anesthetized fish were mounted in a custom imaging chamber, partially embedded in 1% agarose, and covered with MS-222-containing solution. For Video 1, a FLUOVIEW FV3000 scanning confocal microscope (Olympus) equipped with a 100× objective (NA 1.49) was used to collect a z-stack which was 3D rendered with Imaris (Bitplane).

Whole animal photoconversion

Request a detailed protocol

Prior to imaging, Tg(atoh1a:nls-Eos) zebrafish were exposed to light from a UV LED flashlight (McDoer) for 15 min in a reflective chamber constructed from a styrofoam box lined with aluminum foil. A similar lateral region of the trunk was imaged over subsequent days identified by approximate body position below the dorsal fin and relative to underlying pigment stripes.

Regional photoconversion

Request a detailed protocol

After anesthetization and mounting as described above, the Tg(atoh1a:nls-Eos) reporter was photoconverted using the stimulation program of NIS-Elements with the 405 nm laser at 14–18% power for 30–45 s within a 500×500 pixel region of interest with an area of 67,055 µm2. The same lateral region of the trunk was imaged over subsequent days identified by body position under the dorsal fin, position relative to underlying pigment stripes, and presence of photoconverted cells. Animals that died over the course of the experiment were excluded from further analysis.

Staining

Alizarin Red S staining

Request a detailed protocol

To visualize mineralized bone, live animals were stained for 15 min in a solution of 0.01% (wt/vol) Alizarin Red S dissolved in system water and subsequently rinsed 3×5 min in system water prior to imaging as described (Bensimon-Brito et al., 2016).

Antibody staining

Request a detailed protocol

Zebrafish were anesthetized in a solution of 0.012% MS-222 in system water for 5 min. Using metal forceps, up to 10 scales were removed from the lateral side of the trunk in the region below the dorsal fin. Scales were fixed in 4% paraformaldehyde (PFA)/PBS at 4°C overnight. Scales were washed 4×5 min in 1× PBS +0.3% triton-X (PBST) at room temperature and then blocked for 1.5 hr with PBST containing 5% normal goat serum. Incubation with primary antibodies occurred at 4°C overnight, followed by 4×15 min washes in PBST. Scales were incubated in appropriate secondary antibodies for 2 hr at room temperature and washed 4×15 min in PBST. To label nuclei, scales were incubated with DAPI for 5 min at 4°C and washed in PBST 4×5 min at room temperature. Scales were mounted between a microscope slide and coverslip in Prolong gold. All steps were performed on a rotating platform.

AM1-43 staining

Request a detailed protocol

Scales were removed from adult Tg(atoh1a:nls-Eos) zebrafish as described above and placed into the center of a petri dish. 1 ml of L-15 media was added to the dish containing newly plucked scales no longer than 2 min after the scales had been removed. 1.5 µl of 10 mM AM1-43 was added to the dish for a final concentration of 15 µM AM1-43. Scales were incubated for 5 min in this solution to allow for incorporation. Prior to confocal imaging, regional photoconversion of nls-Eos was carried out as described above.

Fluorescent in situ hybridization

Request a detailed protocol

piezo2 antisense RNA was transcribed in vitro from a previously generated plasmid (Faucherre et al., 2013) using SP6 and fluorescein-dUTP. The FISH protocol for adult zebrafish scales was previously described (Lin et al., 2019). Briefly, scales from Tg(atoh1a:Lifeact-EGFP) adults were plucked and fixed in 4% PFA overnight at 4°C then washed three times with 1× PBS +0.1% Tween-20 (PBSTw). Scales were dehydrated in sequential washes of 75% PBSTw:25% methanol (MeOH), 50% PBSTw:50% MeOH, 25%PBSTw:75%MeOH, then placed in 100% MeOH at –20°C overnight. Scales were rehydrated in sequential washes of 25% PBSTw:75% MeOH, 50% PBSTw:50% MeOH, 75%PBSTw:25% MeOH, then washed 3× in PBSTw. Scales were treated with 0.1 mg/ml proteinase K for 5 min, then re-fixed in 4% PFA for 20 min. Scales were washed once in PBSTw, washed once in 50% PBSTw:50% hybridization buffer, then incubated in hybridization buffer for 2 hr at 65°C. Scales were incubated in hybridization buffer with probe (~1 ng/µl) overnight at 65°C. Scales were sequentially washed at 65°C in 75% hybridization buffer:25% 2× SSC +0.1% Tween20 (SSCT), 50% hybridization buffer:50% 2× SSCT, 25% hybridization buffer:75% 2xSSCT, followed by three washes at room temperature in 2× SSCT, followed by three washes in 0.2× SSCT. Scales were then washed 3× in 1× PBS +0.2% Triton X-100 (PBSTr), then blocked for 2 hr in PBSTr +5% fetal bovine serum. Scales were incubated in blocking buffer with an anti-fluorescein peroxidase (POD) conjugated antibody (1:2000) overnight at 4°C. Scales were washed 6× in PBSTr, followed by staining with Tyramide Signal Amplification (TSA) Plus Cyanine 5 (1:50 dilution) for 10 min.

Following FISH, scales were incubated in PBSTr +10% NGS for 2 hr at room temperature. Scales were stained with an anti-GFP antibody (1:1000) in PBSTr +10% NGS overnight at 4°C. Scales were washed in PBSTr, then incubated in secondary antibodies (1:200) for 2 hr at room temperature. Scales were washed in PBSTr, stained with Hoechst (3.24 nM) for 10 min at room temperature, washed in PBSTr, mounted under coverslips in ProLong Gold, and imaged.

Hybridization chain reaction

Request a detailed protocol

A custom piezo2 probe set (set size: 20; amplifier: B3) was ordered using accession number XM_021468270.1. For HCR on adult zebrafish scales, minor alterations were made to a previously described protocol (Ibarra-García-Padilla et al., 2021). Briefly, scales from Tg(atoh1a:Lifeact-EGFP) adults were plucked and fixed in 4% PFA overnight at 4°C, with 20 scales per 1.5 ml tube. Scales were washed 3× in 1× PBS and then dehydrated and permeabilized with 2×10 min washes in 100% methanol. The samples were stored at –20°C overnight. To rehydrate the samples, a series of graded MeOH/PBSTw washes were used for 5 min each: 75% MeOH:25% 1× PBSTw, 50% MeOH:50% 1× PBSTw, 25% MeOH:75% 1× PBSTw, and finally 2× washes in 100% 1× PBSTw. To further permeabilize the scales, samples were incubated in 10 ug/ml Proteinase K diluted in 1× PBSTw for 10 min. Samples were quickly washed 3× in 1× PBSTw, and then post-fixed with 4% PFA for 20 min. After post-fixation, samples underwent 5×5 min washes with 1× PBSTw. Samples were then pre-hybridized with Molecular Instruments HCR hybridization buffer at 37°C for 30 min. After pre-hybridization, samples were incubated with 2 pmol of the probe set diluted in hybridization buffer for 16 hr at 37°C. To remove the probe mixture solution, samples were washed 4× for 15 min each with probe wash buffer at 37°C. Samples were washed 2× for 5 min with 5× SSC +0.1% Tween-20 and then treated with probe amplification buffer for 30 min at room temperature. Samples were washed into hairpin amplification buffer containing snap cooled amplifier hairpins and were incubated at room temperature, protected from light, overnight. Samples were then washed with successive 5× SSC +0.1% Tween-20 washes: 2× washes for 5 min, 2× washes for 30 min, and 1× wash for 5 min. Finally, samples underwent 3×5 min washes with 1× PBSTw. Anti-GFP staining was performed as described above.

Image analysis

Axon contact quantification

Request a detailed protocol

Innervation of Tg(atoh1a:nls-Eos)-expressing cells was scored using a custom ImageJ macro. The macro identified the centroid of each atoh1a+ cell using the ‘3D project’ function and then created a 3D sphere that was 10% larger than the maximum nuclear diameter. A cell was scored as innervated if an axon passed within the sphere. In some cases, the Tg(atoh1a:nls-Eos) reporter was photoconverted prior to image acquisition as described above.

Cell density analysis

Request a detailed protocol

Maximum intensity projections of confocal z-stacks were converted to 8-bit images and thresholded in ImageJ. Cell density was quantified using the ‘Analyze particles’ function of ImageJ. For low-magnification quantification of MC cell density across the trunk of fgf8adhiD1Tg/+ and siblings, tiled images were collected that included multiple scales per region. Cell density was quantified as described above using ImageJ. For high-magnification cell density quantification in the epidermis directly above scales, a small region centered in the epidermis of each full scale in view and positioned based on scale lobe was quantified. For scales with multiple lobes, a density measurement was collected from the center of each lobe and averaged.

Recombination efficiency analysis

Request a detailed protocol

Maximum intensity projections of the DsRed channel of confocal z-stacks for the basal keratinocyte Cre lineage tracing experiments were converted to 8-bit images and thresholded using the ‘Huang’ preset in ImageJ. The ‘%Area’ of the thresholded image was determined using the ‘Measure’ function and taken to represent recombination efficiency.

Statistical analysis

Request a detailed protocol

Statistical tests used are listed in individual figure legends. Plots were created using R or Python. N refers to either the number of individual fish or the number of individual scales where appropriate and is specified in figure legends. For boxplots, each graph possesses a black line within the box that represents the median, two hinges for the first and third quartiles, two whiskers that extend no further than 1.5× (interquartile range) from the adjacent hinge, and outlying points plotted individually beyond the whiskers. Segmented linear regression in Figure 7F was performed using the R function ‘segmented’ with parameter npsi = 1.

Data availability

All data generated or analyses during this study are included in the manuscript and supporting files. Source data files have been provided as indicated.

References

    1. Cheng Chew SB
    2. Leung PY
    (1994)
    Ultrastructural study of the merkel cell and its expression of met-enkephalin immunoreactivity during fetal and postnatal development in mice
    Journal of Anatomy 185 ( Pt 3):511–520.
  26. Book
    1. Fox H
    2. Lane E
    3. Whitear M
    (1980) Sensory nerve endings and receptors in fish and amphibians
    In: Spearman R, Riley P, editors. The Skin of Vertebrates, Linnean Society Symposium Series. Academic Press. pp. 271–281.
    https://doi.org/10.1016/j.acthis.2008.11.002
  36. Book
    1. Hartschuh W
    2. Weihe E
    3. Reinecke M
    (1986) The merkel cell in
    In: Bereiter-Hahn J, Matoltsy AG, Richards KS, editors. Biology of the Integument: 2 Vertebrates. Berlin, Heidelberg: Springer. pp. 605–620.
    https://doi.org/10.1007/978-3-662-00989-5_30
    1. Smith KR
    (1977) The haarscheibe
    The Journal of Investigative Dermatology 69:68–74.
    https://doi.org/10.1111/1523-1747.ep12497883

Decision letter

  1. Marianne E Bronner
    Senior and Reviewing Editor; California Institute of Technology, United States

In the interests of transparency, eLife publishes the most substantive revision requests and the accompanying author responses.

[Editors' note: this paper was reviewed by Review Commons.]

https://doi.org/10.7554/eLife.85800.sa1

Author response

We thank all the reviewers for their positive and constructive comments on the manuscript. We address each of the reviewers’ comments point-by-point below. The major revisions include:

  • Improved statistical details and attention to subjective language throughout.

  • New TEM data included in the new Figure 1—figure supplement 1 to illustrate the drastic ultrastructural differences between MCs and neighboring epidermal cells.

  • Inclusion of an estimate of the “recombination efficiency” of our keratinocyte lineage trace in Figure 4.

  • Additional quantification of MC density in the different body regions (Figure 6) and prior to squamation (Figure 7F). Imaging of the zebrafish oral mucosa (Author response image 1).

  • More nuanced interpretations of the eda and fgf8a mutant phenotypes.

Reviewer #1 (Evidence, reproducibility and clarity (Required)):

The authors describe and characterize the touch system in zebrafish as a new model to study MC development and maintenance. The manuscript is very well written, and the experiments are carefully executed and beautifully illustrated. This study addresses the origin of zebrafish MCs, shows that they are innervated by somatosensory neurons and that they share molecular properties with mammalian MCs. In addition, the authors developed transgenic lines that allow them to study MCs in vivo.

Genetic lineage tracing shows that zf MCs are derived from the epidermis as in mouse and not from neural crest cells, as described for avian MCs. In addition, longevity and turnover of murine MCs was controversial. Here, the authors show that zf MCs constantly turnover and that the distribution and turnover rate in the trunk depends on underlying scales. They show that the loss of scales in eda mutants leads to a decrease in MC production and an increase in MC death showing that scales are required for MC production and maintenance. Using a specific fgf8 mutant allele that causes an increase in Fgf signaling and an increase in scale size they demonstrate that scales are sufficient to induce MCs.

In summary, this manuscript is a rigorous and beautiful characterization of MCs development and maintenance. The authors demonstrate that zebrafish MCs share many characteristics with mammalian MCs. The generation of MCs specific transgenic lines, coupled with existing transgenic lines that label somatosensory cells and cells in the scales sets the stage for detailed further analyses. For example, using these tools one can now study how the size of the MC progenitor domain is controlled, if progenitors migrate and what the identities of the molecular signals from the nerves and scales to progenitors and differentiated MCs are.

Minor comments:

Line 71: Why is the heterogeneity a limitation? Couldn't it also exist in zebrafish?

Thank you for raising this question. The limitations are meant to refer to current limitations of the rodent system and demonstrate an opportunity for a new model system to complement the rodent system. We have rephrased this section to better articulate this point.

Introduction:

“While this system has been useful for understanding many aspects of MC development and function, the rodent system also has several significant limitations that warrant additional models to improve the understanding of MCs.”

We also added the following to the discussion: “As the majority of the analyses completed here focus on MCs found in the trunk epidermis, it will be intriguing to determine whether all MCs in different skin compartments in the juvenile and adult zebrafish share similar molecular, cellular, and functional properties.”

Line 295: The authors write: 'Thus, our observations indicate that the decrease in MC cell density in eda mutants is likely due to both reduced MC production and increased MC turnover'.

It should say: '.. increased MC loss'. In mutants the MCs show poor turnover. I believe the term 'turnover' implies that the cells are being replaced, which is only partially happening here.

Thank you for the clarification. We agree with the reviewer and have changed the wording from “turnover” to “loss” in lines 295 and 301.

Line 301: 'The authors state: 'these data suggest that Eda signaling is required for MC development, maintenance, and distribution along the trunk.'

The authors do not show any data that Eda signaling is involved in MC development but only that scales are needed. The MC inducing signals from scales to the epidermis could be independent from Eda signaling. Please rephrase.

Please discuss that not all MC specification/development depends on scales. Even in the scale-less eda mutants some MCs form (as in the inter scale regions in wt?) and even turnover. Do scales secrete a signal that increases proliferation of existing MC progenitors but scales do not affect specification?

We respectfully disagree with the reviewer on the interpretation of these results. Our experimental manipulation (examination of eda-/- vs. sibling controls) only allows us to conclude that Eda signaling - either directly or indirectly - is required for these processes along the trunk. To determine whether signaling from scales is required would require identification of the signal(s) and/or loss/ablation of scales independent of Eda. We have rephrased the results to more clearly state our interpretation. The corresponding portion of the discussion now reads “Further investigations are required to determine whether Eda signaling directly regulates the differentiation of MC progenitors. Alternatively, since eda mutants lack scales (Harris et al., 2008) and have decreased epidermal innervation (Rasmussen et al., 2018), MC development may require scale-derived and/or somatosensory neuron-derived signals.”

Line 320: The authors describe that the fgf8 allele leads to a redistribution of MCs. Is it really a redistribution, or is it ectopic induction or expansion of existing progenitors? Redistribution implies that the expansion is due to a loss of MCs in another region, which I do not see in the data.

Thank you for raising this point about the potentially poor wording choice relating to “redistribution”. We do not yet know whether the distribution of MCs in fgf8a mutants reflects a redistribution, ectopic induction, or expansion of existing progenitors (these are excellent ideas for future studies). Thus, in response to the reviewer’s comment, we have changed the heading for this results section to “The MC pattern is not predetermined along the trunk” and concluded the section as follows: “... the distribution of MCs tracked with the altered scale size and shape in the mutants, suggesting the MC pattern is not predetermined within the trunk skin compartment (Figure 9E-H).”

Figures:

- Figure 1, panels B-C': EM images are very dark and difficult to see. Letter 'a' is on top of the axon, maybe move to the side and pseudo-color different structures.

In response to these suggestions, we have adjusted the brightness and contrast to lighten the TEM images in Figure 1B-C’ as much as possible. We also moved the ‘a’ off to the side in Figure 1B’ to make the axon more visible. In response to Reviewer #3’s comments (see below), we also added an additional TEM image in the new Figure 1—figure supplement 1 that has presumptive keratinocytes and an MC differentially pseudo-colored. We hesitate to pseudo-color the cells/structures in Figure 1B-C’ for fear of obscuring the underlying TEM images.

- Figure 1, panel D: very difficult to see the magenta axons in the cartoon. Please enlarge and make brighter.

We agree that this needed improvement. In the revised Figure 1D, we made the axons clearer and illustrated the different types of MC-axon associations we observe in Figure 2. We also refer the reader back to this figure in the corresponding axon innervation results section.

- Figure 2, panels A and D: keeping the same antibody stainings in the same color would help with visualization. Matching the bar plots in panel C would be even nicer.

Thank you for the suggestion. The revised Figure 2 now has a consistent color scheme.

- Figure 2, panel C: please identify in the legend if the error bars are SD, SEM or other.

These error bars represent 95% confidence intervals. This information has been added to the figure legend.

- Figure 2, panels G and H: MCs are in cyan in the image, but green in the legend.

This has been corrected.

- Figure 3: include percentages and total number in the image instead of the legend.

The numbers and percentages have now been added to the Figure 3 panels. We have left them in the figure legend for clarity on what was scored.

- Figure 6, panel B: which part of the eye is being depicted?

Thank you for the question. We imaged the corneal epithelium above the lens. This has been clarified in the appropriate parts of Figures 6 and 8 and the corresponding figure legends.

- Figure 6, panel F: please provide error bars and statistics to show that the operculum has a higher density of MCs.

Thank you for the suggestion. In response to the comment, we revised Figure 6F by: (1) increasing the sample size; (2) replotting the data as boxplots rather than bar graphs; and (3) including the results of a one-way ANOVA.

- Figure 7, panels F-H: for simple linear regression, please also provide F and p values.

Thank you for the suggestion. This information has been added to the figure legend.

- Figure 8, panel D: colors for SL do not follow a scale, very hard to understand which is which.

In response to the reviewer’s suggestion, we tried numerous different color palettes. However, we were unable to find a color palette that allowed us to distinguish individual points as well as the rainbow palette used in Figure 8D. Thus, after careful consideration, we have elected to keep the original palette here. For consistency, we have used the same palette in the revised Figure 8–figure supplement 1D and Figure 9–figure supplement 1D.

Methods:

- Line 472: the word "sex" should be used instead of "gender".

Thank you for the correction. This is fixed in the revision.

- Image analysis, line 593. Please provide a more detailed explanation or describe the ImageJ macro used for the analysis.

Our ImageJ macro has been fully annotated and is provided as Figure 2—source code 1 in the revision. The corresponding methods section has also been updated to clarify the methodology.

Reviewer #1 (Significance (Required)):

Soft touch is perceived by Merkel cells (MCs). How MCs develop and are maintained is not well understood because MC development is difficult to study in mammals due to their in utero development. The authors describe and characterize the touch system in zebrafish as a new model to study MC development and maintenance. The study demonstrates that the zebrafish touch system shares many characteristics with its mammalian counterpart, namely its developmental origin, innervation and molecular characteristics. In contrast to mammals, zebrafish transgenic lines that the authors generated, allow the in vivo analysis of Merkel Cell specification, development and maintenance. Therefore this study is the foundation for future detailed cellular and molecular analyses of the touch sensory system and will be of interest to developmental biologists studying stem cells, regeneration and aging, as well as neuroscientists.

We thank the reviewer for their positive assessment of the manuscript.

Reviewer #2 (Evidence, reproducibility and clarity (Required)):

This is a very nice and straightforward paper characterizing mechanosensory Merkel cells in the zebrafish skin. The paper uses a number of criteria, based on our knowledge of Merkel cells in mammals, to identify a population atoh1a expressing cells, with neurosecretory granules and actin rich microvilli as Merkel cells in the zebrafish skin. The authors have used existing transgenic lines and developed some of their own, described in this paper, to follow the development of Merkel cell in zebrafish. They show that Merkel cells are derived from basal keratinocytes not neural crest cells. They have region specific densities that influenced by underlying structures like scales and fin rays. They go to show that Ectodysplasin signaling promotes Merkel cell development in the trunk skin but not above the eye or operculum. Reduction of Merkel cells in eda mutants suggest that Eda signaling is required for their development and maintenance. Finally they show that alteration of zebfrafish scale pattern using a mutant with exaggerated fgf8a expression also alters merkel cell distribution.

The data presented is clear and the conclusions are supported by their observations.

I have no significant issue with the paper as is.

Reviewer #2 (Significance (Required)):

This study will serve as an excellent basis for future work looking at studies of Merkel cell development and function in fish. Though Merkel cells have been studied in mammals, establishing a zebrafish model for their study will help overcome many barriers that make their analysis difficult in mammals.

We thank the reviewer for their positive assessment of the manuscript.

Reviewer #3 (Evidence, reproducibility and clarity (Required)):

In this manuscript, Brown et al. (2022) seek to characterize and address fundamental questions regarding the development and dynamics of Merkel cells (MCs) in zebrafish (Danio rerio). The authors utilize a diverse and complementary suite of methods to characterize presumptive MCs in the epidermis of adult zebrafish, including electron microscopy, novel transgenic lines, confocal imaging, and various immuno- and non-immunohistological staining techniques. These studies demonstrate that zebrafish MCs share many features with vertebrate (including mammalian) MCs, particularly regarding morphology/structure, putative functions, genetic markers, and bodily distribution.

After establishing the identity of zebrafish MCs, the authors employ lineage tracing and cell tracking analyses to determine that trunk MCs derive from basal keratinocytes and exhibit regular cell turnover. Finally, the authors examine how trunk scales may affect MC development by using established scale mutants. These results show that the presence/absence of scales influences trunk MC development, while scale characteristics (e.g. shape, size) change the distribution of MCs.

Major comments:

The key conclusions of the manuscript are convincing, however, several points should be addressed by the authors.

• Throughout the manuscript, the authors make general claims about zebrafish MCs (zMCs) based on the evidence collected. Yet, most of this evidence (particularly claims about MC turnover, development, structure) comes from examination and experimentation of a specific MC population: trunk MCs located in the scale epidermis. The authors remark upon mammalian MC diversity (lines 73-74) and go on to highlight the diversity of MCs throughout the adult zebrafish (Figure 6), which have differing densities and distribution patterns. Any statements that suggest all zebrafish MCs share certain qualities/features should be carefully considered given the evidence presented.

Thank you for raising this important point. We have added wording in the results and discussion to clearly articulate that the majority of our analyses and conclusions are based on trunk MCs:

Results:

“Anticipating the conclusion of our analysis below, we shall hereafter refer to the epidermal atoh1a+ cells as MCs, with the majority of the analyses completed on trunk MCs unless stated otherwise.”

Discussion:

“As the majority of the analyses completed here focus on MCs found in the trunk epidermis, it will be intriguing to determine whether all MCs in different skin compartments in the juvenile and adult zebrafish share similar molecular, cellular, and functional properties.”

• In the manuscript, the authors validate several markers for the identification of zMCs based upon known mammalian markers (e.g. atoh1a, sox2, piezo2, SV2, 5-HT, and AM1-43; Figures 1-3). Yet, another well-known marker for MCs (CK8) is not addressed (Moll, 1995; Moll, 2005). One zebrafish ortholog for CK8 is krt4, a transgene successfully employed in this study to label keratinocytes. Do zMCs express krt4 or other mammalian MC keratins? Answering this question or addressing this discrepancy would further strengthen the authors claims that these cells are bona fide zMCs.

We agree with the reviewer that (1) identification of a keratin(s) that distinguishes MCs from other epidermal cell types in zebrafish would be an excellent reagent; and (2) readers familiar with the mammalian MC literature may similarly wonder why this was not addressed in the manuscript. Indeed, we had considered whether we could identify homologs of CK8, CK20 or other mammalian MC keratins that would label zebrafish MCs. However, despite the confusing nomenclature that would indicate otherwise, the zebrafish keratins share more homology with each other than the corresponding mammalian proteins (Ho et al., 2022; PMID: 34991727). Our revised results section now includes the following to clarify this point: “For example, keratins, most notably keratin 8 and keratin 20, have been used extensively as markers of mammalian MCs (Moll et al., 1995, 1984). However, zebrafish keratins have undergone extensive gene loss and duplication and are not orthologous to mammalian keratin genes (Ho et al., 2022). Thus, we considered alternative molecular markers to label zebrafish MCs.”

• The authors utilize a previously validated eda mutant line to see if ectodysplasin signaling affects zMC development. While the results of these experiments are convincing, the authors need to make clear whether they are claiming that scales, scale-derived Eda signaling, or Eda signaling alone dictate trunk MC development. It appears that there is some conflation of these ideas, particularly with line 306 ("blocking dermal appendage formation inhibited MC development" is a different claim from 'blocking Eda signaling inhibited MC development'). One way to make this differentiation would be to perform a similar experiment as detailed in Xiao et al. 2016: using a Shh agonist in eda mutants. If scale-specific signals are required in addition to Eda, we would expect to see similar MC densities and patterns in both Shh agonist-treated and non-treated eda mutants.

We agree that our interpretation of these results could have been more clearly articulated in our initial submission. As discussed above in response to Reviewer #1, we do not yet know whether Eda signaling directly or indirectly influences MC development. We have revised the results section to clarify our interpretation of the results as follows: “Together, these data suggest that either Eda signaling, or a scale-derived signal, is required for MC development, maintenance, and distribution along the trunk. Further studies are required to determine the specific scale-derived signal that regulates MC development in the trunk.”

The suggestion of using a Shh pathway agonist in eda mutants to attempt to rescue MC differentiation similar to Xiao et al. 2016 is an interesting one. To our knowledge, experiments validating the Smo agonist used by Xiao and colleagues (Hh-Ag1.5) in zebrafish have not been published. We also note that activation of Shh signaling by heat-shock induction of shha expression during squamation led to kyphosis and epidermal migration off of the trunk (Aman et al., 2018; PMID: 30014845). Thus, we respectfully suggest that distinguishing between the various possibilities downstream of Eda is beyond the scope of the current manuscript. We have added a discussion point along these lines: “Further investigations are required to determine whether Eda signaling directly regulates the differentiation of MC progenitors. Alternatively, since eda mutants lack scales (Harris et al., 2008) and have decreased epidermal innervation (Rasmussen et al., 2018), MC development may require scale- and/or somatosensory neuron-derived signals. Finally, we note that trunk MCs are not completely absent in eda mutants, suggesting that a subset of MCs develop independent of Eda signaling.”

• Throughout the manuscript, the authors use subjective language (e.g. line 106). While this reviewer does not wish to suppress or alter the authors' voices, careful consideration should be used when employing these types of descriptors. Furthermore, the authors use suggestively quantitative language inappropriately or unjustifiably. For example, in line 221, the authors use "extensive" when describing the co-labeling between atoh1a+ MCs and lineage-traced basal keratinocytes; the percentage of co-labeled cells ranged from 29-32%. Other quantitative descriptors such as "frequently" (line 171) or "uniform" (line 249) describe various features or phenomena without quantification in figures or supplements.

Thank you for this comment. We have paid careful attention to our subjective/statistical language in the revision. Regarding the usage of “uniform” - we have added the wording “relatively uniformly” to descriptions and a statement that our term “uniform” was not specifically quantified. Although the uniform appearance was not specifically quantified, we believe this provides an accurate description of the MC localization pattern in certain skin compartments.

Example word change in Results:

“For example, MCs were distributed relatively uniformly across the eye, although this spatial pattern was not specifically quantified.”

• In the lineage tracing experiments (Figure 4), the authors note that "recombination is not complete" (lines 1016-1017) to explain why not all zMCs express the basal keratinocyte lineage marker. While this idea could be supported by Figure 4-figure supplement 1, one could postulate that zMCs are derived from multiple progenitor lineages. Using the basal keratinocyte lineage tracing validation, the authors could in theory calculate a "recombination efficiency" of this transgenic line and determine approximately the percent of zMCs they 'lose' as a result. Otherwise, the authors could perform other experiments to support the claim that zMCs derive from basal keratinocytes. For example, could the authors photoconvert basal keratinocytes at 1 dpf and see how many derived MCs are still photoconverted later? Could they do this photoconversion experiment with neural crest cells? Could they ablate neural crest cells and determine if MC number is affected? These additional experiments are not necessarily required for publication, but some explanation of the unexpectedly low percentage of basal keratinocyte lineage marker-labeled MCs would suffice.

We thank the reviewer for raising this important point and the suggestion of calculating a “recombination efficiency”. We note that Cre responsive transgenes are far from a perfect technology in zebrafish as recently characterized by Lalonde et al. (2022; PMID:35582941). In response to the reviewer’s comment, we added an estimate of the recombination efficiency to Figure 4 (panels E, G, H). Importantly, a comparison between the recombination efficiency and percentage of MCs labeled by the basal keratinocyte Cre tracing was not significantly different. Our revised results section reads as follows: “After raising 4-OHT-treated animals to adulthood, we observed variable (2-81%) co-labeling between the basal keratinocyte lineage trace and MC reporters (Figure 4D’,F). We note that our lineage tracing strategy did not label all basal keratinocytes (Figure 4D; Figure 4—figure supplement 1), suggestive of incomplete Cre-ERT2 induction and/or transgene recombination. Consistent with the latter possibility, a recent analysis demonstrated Tg(actb2:LOXP-BFP-LOXP-DsRed) has a low recombination efficiency compared to other Cre reporter transgenes (Lalonde et al., 2022). To estimate the local recombination efficiency in imaged regions, we thresholded the DsRed channel and calculated the fraction of skin cells labeled (Figure 4E). Importantly, the proportion of MCs labeled by the basal keratinocyte lineage trace was not significantly different from the local recombination efficiency (Figure 4G-H). These observations support a basal keratinocyte origin of most or all zebrafish MCs.”

• The authors use appropriate statistics and have sufficient replicates when this information is presented. Yet, the presence or absence of these data is not consistent within figure captions. The authors must ensure that they provide the N of adults and scales (when appropriate), the SL range of adults, and transgenic lines used. Statistics are missing in some figures (for example: Figures 4E, 5D, 5E, 6F, 8S-1E, 9E-H) where it would be appropriate to include them. In some figures, the N changes over time (example: 5D, 5E); an explanation in the 'Methods' section would suffice.

Thank you for noting the need for additional statistics. We have added statistics to the above figures. For Figure 9E-H, we have not added additional statistics. Figure 9E-H serve to graphically visualize differences. We show statistical differences in Figure 9—figure supplement 1 for scale area, aspect ratio, and Feret’s diameter. We have added an explanation related to Figure 5D,E in the methods section: “Animals that died over the course of the experiment were excluded from further analysis.”

Minor comments

While the authors present an extensive argument for their claims, addressing these additional comments would further strengthen their story.

• Are zMC nuclei lobulated? This ultrastructure characteristic seems to be common in MCs (Chew & Leung, 1994; Tachibana & Nawa, 2002; Moll, 2005; Boulais, 2009).

We have not observed any lobulation of the MC nuclei by TEM, nor was this commented on in the TEM studies of Whitear and colleagues in other teleosts (Lane and Whitear, 1977; PMID: 198137; Whitear, 1989; PMID: 2510796). Nevertheless, we cannot rule out the possibility that serial sectioning or other high resolution analysis of the nuclear shape may reveal such features. In response to the reviewer’s comment, we have added the following paragraph to the discussion: “While our characterization revealed substantial similarities between mammalian and zebrafish MCs, we did observe anatomical differences in line with previous ultrastructural characterizations of teleost MCs (Lane and Whitear, 1977; Whitear, 1989). For example, the nuclei of mammalian MC are commonly lobulated (Boulais et al., 2009; Cheng Chew and Leung, 1994; Moll et al., 2005; Tachibana and Nawa, 2002). While we did not observe lobulation of zebrafish MC nuclei by TEM, we cannot rule out that serial sectioning or high-resolution reconstruction of nuclear shape would reveal lobulation. Mammalian MCs typically localize adjacent to basal keratinocytes (Boot et al., 1992; Cheng Chew and Leung, 1994; Fradette et al., 1995; Mihara et al., 1979; Moll et al., 1996; Smith, Jr, 1977), whereas zebrafish MCs appear in upper strata, typically beneath the periderm (Figure 1D,G’’). As the majority of the analyses completed here focus on MCs found in the trunk epidermis, it will be intriguing to determine whether all MCs in different skin compartments in the juvenile and adult zebrafish share similar properties.”

• In Figure 3C and 3", the authors show that AM1-43 labels zMCs. Yet, this technique should also stain sensory axons that associate with MCs (Meyers, 2003). Are axons also stained? Other positive controls for the stains could be useful as a supplement.

The reviewer is correct that Meyers et al., (2003; PMID: 12764092) report AM1-43 staining of neurites that innervate MCs in the whisker follicle. However, they did not report similar staining of neurites innervating touch dome MCs. In murine hairy skin, the related styryl dye FM1-43 appears to most prominently stain MCs and hair follicle-associated lanceolate endings (Banks et al., 2013 PMID: 23440964; Villarino et al., 2022 preprint DOI: 10.1101/2022.05.26.493600). Our revised legend for Figure 3 now includes the following: “AM1-43 has been reported to stain neurites innervating MCs in murine whisker vibrissae (Meyers et al., 2003). However, our AM1-43 staining regiment did not strongly label cutaneous axons, although we cannot exclude low levels of staining.”

All of the stains used in our original Figure 3 have been previously validated in zebrafish, which we have more clearly stated and cited in the corresponding results section of the revision. Because these reagents have all been previously validated and our staining patterns are consistent with the literature, we respectfully suggest that positive controls would add little value to the current manuscript. Nevertheless, in response to the reviewer’s comment, we confirmed our piezo2 FISH staining using an independent method (a piezo2 HCR probe). We have included these HCR results as the updated Figure 3D and moved the original Figure 3D to Figure 3—figure supplement 1.

• In Figure 7, the authors argue that as scales develop, MC density increases with scale area. Did the authors compare MC densities of differently-sized scales at the same age? Is fish SL/age a potential confound in the interpretation of these data?

Thank you for the suggestion. In response to the reviewer’s comment, we have replotted the data in Figure 7G,H for animals in the range 8-10 mm SL in Figure 7—figure supplement 1. We have revised the corresponding results section as follows: “The density and number of MCs positively correlated with scale area (Figure 7G,H), although this trend was less pronounced at stages less than 10 mm (Figure 7—figure supplement 1)”. As discussed above in response to Reviewer #1’s suggestion, we also now report F-statistics and P-values for the linear regressions in the figure legends.

• The authors claim that squamation begins at ~9 mm SL (line 268), prior to which MCs were "rare" in the epidermis (supported by data in Figure 7F). However, Figures 8A and 8G suggest that MCs are not rare prior to squamation/9 mm SL. Are these data in conflict?

Thank you for raising this observation. We do not believe these data are in conflict. Figure 8A and B show images of fish 8.8-8.9 mm SL, immediately prior to squamation. MCs appear about the same time as scales develop but the exact timing varies between animals. To further strengthen this section of the manuscript, we now include quantification of the density of trunk MCs at various stages prior to 9 mm SL (new data added to the developmental timeline in Figure 7F). These data are consistent with our initial interpretation. In the revised results section we clarify this as follows: “Using reporters that label MCs and scale-forming osteoblasts, we rarely observed MCs in the epidermis prior to 8 mm SL (Figure 7B, F). Between 8-10 mm SL, MCs appeared at a low density along the trunk (Figure 7F). MC density rapidly increased from 10-15 mm SL, a period of active scale growth (Figure 7C-F).”

• In Figure 6B-E, the panels are incorrectly labeled as "atoh1a:nls-Eos" (figure caption and fluorescence localization show they are atoh1a:Lifeact-EGFP).

The low magnification panels were correctly labeled as atoh1a:nls-Eos. The insets showed atoh1a:Lifeact-EGFP as described in the figure legend. We apologize for the confusion and poor data presentation. We have revised Figure 6 to eliminate the problematic labeling/display.

• Figure 9 panels E-H are not referenced in the main body of the text.

Thank you for pointing this out. Fixed in the revised manuscript.

• In Figure 6, the authors examine MC densities in the tail, but do not quantify changes here with eda mutants as they did for other regions (eye, operculum) in Figure 8. Why was this region not examined?

We have clarified this point in the revised results section as follows: “eda mutants lack fins at the stages analyzed (Harris et al., 2008) precluding analysis of these regions in the homozygous mutants.”

• The authors do a good job in detailing the current literature regarding MCs, however, two missing areas are noticeable: (1) there is no mention of mammalian MCs that reside in the oral mucosa (Hashimoto, 1972) or whether they exist in zebrafish, and (2) no mention of Merkel-like cells (Halata, 2003) and why the cells in this paper are or are likely not Merkel-like cells.

Thank you for the suggestions. Regarding the first point, we revised the introduction to reference (Hashimoto, 1972) as follows: “...vertebrates have diverse types of skin and MCs are found in both hairy and glabrous (non-hairy) skin, as well as mucocutaneous regions such as the gingiva and palate (Hashimoto, 1972; Lacour et al., 1991; Moayedi et al., 2021).” We also imaged the mucosal tissue along the roof palate of the adult mouth and identified atoh1a+ cells (see Author response image 1). Close examination of the atoh1a:Lifeact-EGFP signal revealed these cells have a spherical morphology and extend short processes similar to the MCs described across the body regions examined in Figure 6. However, as the microvillar morphology of the palatal atoh1a+ cells is not identical to those identified in other skin regions, we hesitate to call these MCs without performing additional in-depth analyses. We feel that inclusion of these data in the manuscript could distract the reader from the main focus of our study, therefore we have included them here:

Author response image 1
Download asset Open asset
atoh1a+ cells in the adult oral epithelium.

(A,B) Low- (A) and high-magnification (B) confocal micrographs of oral roof palate epithelium in an adult expressing reporters for keratinocytes (Tg(krt4:DsRed)) and atoh1a-expressing cells (Tg(atoh1a:Lifeact-EGFP)). (B’) Reconstructed cross section along the yellow line in B showing two atoh1a+ cells in the upper strata of the oral epithelium. Scale bars: 50 µm (A) and 10 µm (B,B’).

Regarding the second point, we have added the following sentence to the first paragraph of the discussion: “Second, zebrafish MCs extend numerous short, actin-rich microvilli and complex with somatosensory axons, classic morphological hallmarks of MCs (Mihara et al., 1979; Smith, Jr, 1977; Toyoshima et al., 1998). Our morphological observations support the interpretation that these cells are MCs rather than Merkel-like cells, which lack axon association and microvillar processes (reviewed by Halata et al., 2003).”

• It may help readers understand MC morphology in context if the authors include a larger picture of the TEM data that highlights the drastic difference in ultrastructure between MCs and neighboring keratinocytes.

Thank you for the suggestion. We added a new figure (Figure 1—figure supplement 1) to the revised manuscript that contains an additional TEM image that we believe illustrates the different morphologies of keratinocytes and MCs. We hope this will help the reader contextualize the morphology and position of MCs within the zebrafish epidermis. This is now referenced in the first results section as follows: “The cells appeared relatively small and spherical with a low cytoplasmic-to-nuclear ratio compared to neighboring keratinocytes (Figure 1B,C; Figure 1—figure supplement 1) …”

Reviewer #3 (Significance (Required)):

The current manuscript provides significant advancements in various biological fields and research communities. For researchers that utilize zebrafish as a model organism, these findings present a new cell type along with novel and essential genetic tools for study. These developments open the possibilities to further understand MCs, their roles in somatosensory function, mechanisms of cell type diversification, and to engage in translational research. For those already researching MCs, this manuscript shows that fundamental questions regarding MC functioning can be rigorously addressed with a new model that can fill the methodological limitations imposed by mammalian biology. Indeed, the authors do a thorough job of introducing and contextualizing our knowledge of MCs and any outstanding gaps. The authors then sit their findings comfortably alongside previous works, largely supporting those findings, and take the extra step to address MC controversies/matters of debate. This technique of supporting the current literature and then uplifting it with new findings makes this work even more impressive. Various audiences will find value in this manuscript, including but not limited to those that study epidermal cell types, the development and influence of skin appendages, somatosensation and sensory disorders, developmental biology, and Merkel cell carcinoma.

We thank the reviewer for their positive assessment of the manuscript.

https://doi.org/10.7554/eLife.85800.sa2

Article and author information

Author details

  1. Tanya L Brown

    Department of Biology, University of Washington, Seattle, United States
    Contribution
    Conceptualization, Formal analysis, Funding acquisition, Investigation, Visualization, Methodology, Writing – original draft, Writing – review and editing
    Contributed equally with
    Emma C Horton
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0001-9554-178X

  2. Emma C Horton

    Department of Biology, University of Washington, Seattle, United States
    Contribution
    Software, Formal analysis, Investigation, Visualization, Methodology, Writing – review and editing
    Contributed equally with
    Tanya L Brown
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0001-9730-7380

  3. Evan W Craig

    Department of Biology, University of Washington, Seattle, United States
    Contribution
    Formal analysis, Funding acquisition, Investigation, Visualization, Methodology
    Competing interests
    No competing interests declared

  4. Camille EA Goo

    Department of Biology, University of Washington, Seattle, United States
    Contribution
    Formal analysis, Investigation, Visualization, Methodology, Writing – review and editing
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-9118-4006

  5. Erik C Black

    1. Department of Biology, University of Washington, Seattle, United States
    2. Molecular and Cellular Biology Program, University of Washington, Seattle, United States
    Contribution
    Formal analysis, Investigation, Visualization, Methodology, Writing – review and editing
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-2333-8923

  6. Madeleine N Hewitt

    1. Molecular and Cellular Biology Program, University of Washington, Seattle, United States
    2. Department of Biological Structure, University of Washington, Seattle, United States
    Contribution
    Resources
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-4387-327X

  7. Nathaniel G Yee

    Department of Biology, University of Washington, Seattle, United States
    Contribution
    Formal analysis, Investigation, Visualization
    Competing interests
    No competing interests declared

  8. Everett T Fan

    Department of Biology, University of Washington, Seattle, United States
    Contribution
    Formal analysis, Investigation, Visualization
    Competing interests
    No competing interests declared

  9. David W Raible

    1. Department of Biological Structure, University of Washington, Seattle, United States
    2. Department of Otolaryngology - Head and Neck Surgery, University of Washington, Seattle, United States
    3. Institute for Stem Cell and Regenerative Medicine, University of Washington, Seattle, United States
    Contribution
    Supervision
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-5342-5841

  10. Jeffrey P Rasmussen

    1. Department of Biology, University of Washington, Seattle, United States
    2. Institute for Stem Cell and Regenerative Medicine, University of Washington, Seattle, United States
    Contribution
    Conceptualization, Formal analysis, Supervision, Funding acquisition, Validation, Investigation, Visualization, Writing – original draft, Project administration, Writing – review and editing
    For correspondence
    rasmuss@uw.edu
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0001-6997-3773

Funding

National Science Foundation (2011008)

  • Tanya L Brown

National Science Foundation (DGE-2140004)

  • Evan W Craig

Eunice Kennedy Shriver National Institute of Child Health and Human Development (R01HD107108)

  • Jeffrey P Rasmussen

University of Washington Research Royalty Fund (A153025)

  • Jeffrey P Rasmussen

Cancer Consortium (NIH/NCI Cancer Center Support Grant P30 CA015704)

  • Jeffrey P Rasmussen

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Acknowledgements

We thank the LSB Aquatics staff for animal care; Wai Pang Chan and Marianne Cilluffo for TEM support; the labs of Ajay Dhaka, Jacqueline Lees, and Alvaro Sagasti for sharing zebrafish stocks; the lab of Chris Joplin for sharing the piezo2 plasmid. The authors are grateful to all members of the Rasmussen lab for discussion, technical assistance, and support. This work was funded in part by a Postdoctoral Fellowship (#2011008) from the National Science Foundation to TLB, a Graduate Research Fellowship (DGE-2140004) from the National Science Foundation to EWC, R01HD107108 from the Eunice Kennedy Shriver National Institute of Child Health and Human Development to JPR, A153025 from the University of Washington Research Royalty Fund to JPR, and a New Investigator Award from the University of Washington/Fred Hutchinson Cancer Research Center Cancer Consortium, which is supported by the NIH/NCI Cancer Center Support Grant P30 CA015704, to JPR. JPR is a Washington Research Foundation Distinguished Investigator.

Ethics

All zebrafish experiments were approved by the Institutional Animal Care and Use Committee at the University of Washington (Protocol: #4439-01).

Senior and Reviewing Editor

  1. Marianne E Bronner, California Institute of Technology, United States

Publication history

  1. Preprint posted: August 11, 2022 (view preprint)
  2. Received: December 28, 2022
  3. Accepted: January 11, 2023
  4. Accepted Manuscript published: January 17, 2023 (version 1)
  5. Version of Record published: February 6, 2023 (version 2)

Copyright

© 2023, Brown, Horton 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

  • 733
    Page views
  • 125
    Downloads
  • 0
    Citations

Article citation count generated by polling the highest count across the following sources: Crossref, PubMed Central, Scopus.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

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)

  1. Tanya L Brown
  2. Emma C Horton
  3. Evan W Craig
  4. Camille EA Goo
  5. Erik C Black
  6. Madeleine N Hewitt
  7. Nathaniel G Yee
  8. Everett T Fan
  9. David W Raible
  10. Jeffrey P Rasmussen
(2023)
Dermal appendage-dependent patterning of zebrafish atoh1a+ Merkel cells
eLife 12:e85800.
https://doi.org/10.7554/eLife.85800

Categories and tags

Research organism

Further reading

    1. Developmental Biology

    Npr3 regulates neural crest and cranial placode progenitors formation through its dual function as clearance and signaling receptor

    Arun Devotta, Hugo Juraver-Geslin ... Jean-Pierre Saint-Jeannet
    Research Article Updated

    Natriuretic peptide signaling has been implicated in a broad range of physiological processes, regulating blood volume and pressure, ventricular hypertrophy, fat metabolism, and long bone growth. Here, we describe a completely novel role for natriuretic peptide signaling in the control of neural crest (NC) and cranial placode (CP) progenitors formation. Among the components of this signaling pathway, we show that natriuretic peptide receptor 3 (Npr3) plays a pivotal role by differentially regulating two developmental programs through its dual function as clearance and signaling receptor. Using a combination of MO-based knockdowns, pharmacological inhibitors and rescue assays we demonstrate that Npr3 cooperate with guanylate cyclase natriuretic peptide receptor 1 (Npr1) and natriuretic peptides (Nppa/Nppc) to regulate NC and CP formation, pointing at a broad requirement of this signaling pathway in early embryogenesis. We propose that Npr3 acts as a clearance receptor to regulate local concentrations of natriuretic peptides for optimal cGMP production through Npr1 activation, and as a signaling receptor to control cAMP levels through inhibition of adenylyl cyclase. The intracellular modulation of these second messengers therefore participates in the segregation of NC and CP cell populations.

    1. Developmental Biology
    2. Evolutionary Biology

    Mutation of vsx genes in zebrafish highlights the robustness of the retinal specification network

    Joaquín Letelier, Lorena Buono ... Juan R Martínez-Morales
    Research Article

    Genetic studies in human and mice have established a dual role for Vsx genes in retina development: an early function in progenitors’ specification, and a later requirement for bipolar-cells fate determination. Despite their conserved expression patterns, it is currently unclear to which extent Vsx functions are also conserved across vertebrates, as mutant models are available only in mammals. To gain insight into vsx function in teleosts, we have generated vsx1 and vsx2 CRISPR/Cas9 double knockouts (vsxKO) in zebrafish. Our electrophysiological and histological analyses indicate severe visual impairment and bipolar cells depletion in vsxKO larvae, with retinal precursors being rerouted toward photoreceptor or Müller glia fates. Surprisingly, neural retina is properly specified and maintained in mutant embryos, which do not display microphthalmia. We show that although important cis-regulatory remodelling occurs in vsxKO retinas during early specification, this has little impact at a transcriptomic level. Our observations point to genetic redundancy as an important mechanism sustaining the integrity of the retinal specification network, and to Vsx genes regulatory weight varying substantially among vertebrate species.

    1. Developmental Biology

    Development: The hidden depths of zebrafish skin

    Yue Rong Tan, Megan Liaw, Chen-Hui Chen
    Insight

    Single-cell transcriptome analysis of zebrafish cells clarifies the signalling pathways controlling skin formation and reveals that some cells produce proteins required for human teeth to acquire their enamel.