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.

Skeletal muscles support the stability and mobility of the skeleton but differ in biomechanical properties and physiological functions. The intrinsic factors that regulate muscle-specific characteristics are poorly understood. To study these, we constructed a large atlas of RNA-seq profiles from six leg muscles and two locations from one muscle, using biopsies from 20 healthy young males. We identified differential expression patterns and cellular composition across the seven tissues using three bioinformatics approaches confirmed by large-scale newly developed quantitative immune-histology procedures. With all three procedures, the muscle samples clustered into three groups congruent with their anatomical location. Concomitant with genes marking oxidative metabolism, genes marking fast- or slow-twitch myofibers differed between the three groups. The groups of muscles with higher expression of slow-twitch genes were enriched in endothelial cells and showed higher capillary content. In addition, expression profiles of Homeobox (HOX) transcription factors differed between the three groups and were confirmed by spatial RNA hybridization. We created an open-source graphical interface to explore and visualize the leg muscle atlas (https://tabbassidaloii.shinyapps.io/muscleAtlasShinyApp/). Our study reveals the molecular specialization of human leg muscles, and provides a novel resource to study muscle-specific molecular features, which could be linked with (patho)physiological processes.