Basal keratinocytes contribute to all strata of the adult zebrafish epidermis

Phased ERK responsiveness and developmental robustness regulate teleost skin morphogenesis

Elongation of the vertebrate embryonic axis necessitates rapid expansion of the epidermis to accommodate the growth of underlying tissues. Here, we generated a toolkit to visualize and quantify signaling in entire cell populations of the periderm, the outermost layer of the epidermis, in live developing zebrafish. We find that oriented cell divisions facilitate growth of the early periderm during axial elongation rather than cell addition from the basal layer. Activity levels of Extracellular signal-regulated kinase (ERK), a downstream effector of the MAPK pathway, gauged by a live biosensor, predict cell cycle entry, and optogenetic ERK activation regulates cell cycling dynamics. As development proceeds, rates of peridermal cell proliferation decrease, and ERK activity becomes more pulsatile and functionally transitions to promote hypertrophic cell growth. Targeted genetic blockade of cell division generates animals with oversized periderm cells, yet, unexpectedly, development to adulthood is not impaired. Our findings reveal stage-dependent differential responsiveness to ERK signaling and marked developmental robustness in growing teleost skin.

Paired and solitary ionocytes in the zebrafish olfactory epithelium

The sense of smell is generated by electrical currents that are influenced by the concentration of ions in olfactory sensory neurons and mucus. In contrast to the extensive morphological and molecular characterization of sensory neurons, there has been little description of the cells that control ion concentrations in the zebrafish olfactory system. Here, we report the molecular and ultrastructural characterization of zebrafish olfactory ionocytes. Transcriptome analysis suggests that the zebrafish olfactory epithelium contains at least three different ionocyte types, which resemble Na + /K + -ATPase-rich (NaR), Na + /Cl - cotransporter (NCC), and H + -ATPase-rich (HR) cells, responsible for calcium, chloride, and pH regulation, respectively, in the zebrafish skin. NaR-like and HR-like ionocytes are usually adjacent to one another, whereas NCC-like cells are usually solitary. The distinct subtypes are differentially distributed: NaR-like/HR-like cell pairs are found broadly within the olfactory epithelium, whereas NCC-like cells reside within the peripheral non-sensory multiciliated cell zone. Comparison of gene expression and serial-section electron microscopy analysis indicates that the NaR-like cells wrap around the HR-like cells and are connected to them by shallow tight junctions. The development of olfactory ionocyte subtypes is also differentially regulated, as pharmacological Notch inhibition leads to a loss of NaR-like and HR-like cells, but does not affect NCC-like ionocyte number. These results provide a molecular and anatomical characterization of olfactory ionocytes in a stenohaline freshwater teleost. The paired ionocytes suggest that both transcellular and paracellular transport regulate ion concentrations in the olfactory epithelium, while the solitary ionocytes may enable independent regulation of multiciliated cells.

Environmental and molecular control of tissue-specific ionocyte differentiation in zebrafish

Organisms cope with environmental fluctuations and maintain fitness in part via reversible phenotypic changes (acclimation). Aquatic animals are subject to dramatic seasonal fluctuations in water salinity, which affect osmolarity of their cells and consequently cellular function. Mechanosensory lateral line hair cells detect water motion for swimming behavior and are especially susceptible to salinity changes due to their direct contact with the environment. To maintain hair cell function when salinity decreases, neuromast (Nm)-associated ionocytes differentiate and invade lateral line neuromasts. The signals that trigger the adaptive differentiation of Nm ionocytes are unknown. We demonstrate that new Nm ionocytes are rapidly specified and selectively triggered to proliferate by low Ca2+ and Na+/Cl- levels. We further show that Nm ionocyte recruitment and induction is affected by hair cell activity. Once specified, Nm ionocyte differentiation and survival are associated with sequential activation of different Notch pathway components, a process different from other tissue-specific ionocytes. In summary, we show how environmental changes activate a signaling cascade that leads to physiological adaptation. This may prove essential for survival not only in seasonal changing environments but also in changing climates.

Periderm fate and independence of tooth formation are conserved across osteichthyans

BACKGROUND

Previous studies have reported that periderm (the outer ectodermal layer) in zebrafish partially expands into the mouth and pharyngeal pouches, but does not reach the medial endoderm, where the pharyngeal teeth develop. Instead, periderm-like cells, arising independently from the outer periderm, cover prospective tooth-forming epithelia and are crucial for tooth germ initiation. Here we test the hypothesis that restricted expansion of periderm is a teleost-specific character possibly related to the derived way of early embryonic development. To this end, we performed lineage tracing of the periderm in a non-teleost actinopterygian species possessing pharyngeal teeth, the sterlet sturgeon (Acipenser ruthenus), and a sarcopterygian species lacking pharyngeal teeth, the axolotl (Ambystoma mexicanum).

RESULTS

In sturgeon, a stratified ectoderm is firmly established at the end of gastrulation, with minimally a basal ectodermal layer and a surface layer that can be homologized to a periderm. Periderm expands to a limited extent into the mouth and remains restricted to the distal parts of the pouches. It does not reach the medial pharyngeal endoderm, where pharyngeal teeth are located. Thus, periderm in sturgeon covers prospective odontogenic epithelium in the jaw region (oral teeth) but not in the pharyngeal region. In axolotl, like in sturgeon, periderm expansion in the oropharynx is restricted to the distal parts of the opening pouches. Oral teeth in axolotl develop long before mouth opening and possible expansion of the periderm into the mouth cavity.

CONCLUSIONS

Restricted periderm expansion into the oropharynx appears to be an ancestral feature for osteichthyans, as it is found in sturgeon, zebrafish and axolotl. Periderm behavior does not correlate with presence or absence of oral or pharyngeal teeth, whose induction may depend on 'ectodermalized' endoderm. It is proposed that periderm assists in lumenization of the pouches to create an open gill slit. Comparison of basal and advanced actinopterygians with sarcopterygians (axolotl) shows that different trajectories of embryonic development converge on similar dynamics of the periderm: a restricted expansion into the mouth and prospective gill slits.

Damage-induced basal epithelial cell migration modulates the spatial organization of redox signaling and sensory neuron regeneration

Epithelial damage leads to early reactive oxygen species (ROS) signaling, which regulates sensory neuron regeneration and tissue repair. How the initial type of tissue injury influences early damage signaling and regenerative growth of sensory axons remains unclear. Previously we reported that thermal injury triggers distinct early tissue responses in larval zebrafish. Here, we found that thermal but not mechanical injury impairs sensory axon regeneration and function. Real-time imaging revealed an immediate tissue response to thermal injury characterized by the rapid Arp2/3-dependent migration of keratinocytes, which was associated with tissue scale ROS production and sustained sensory axon damage. Isotonic treatment was sufficient to limit keratinocyte movement, spatially restrict ROS production, and rescue sensory neuron function. These results suggest that early keratinocyte dynamics regulate the spatial and temporal pattern of long-term signaling in the wound microenvironment during tissue repair.

Keratinocytes drive the epithelial hyperplasia key to sea lice resistance in coho salmon

BACKGROUND

Salmonid species have followed markedly divergent evolutionary trajectories in their interactions with sea lice. While sea lice parasitism poses significant economic, environmental, and animal welfare challenges for Atlantic salmon (Salmo salar) aquaculture, coho salmon (Oncorhynchus kisutch) exhibit near-complete resistance to sea lice, achieved through a potent epithelial hyperplasia response leading to rapid louse detachment. The molecular mechanisms underlying these divergent responses to sea lice are unknown.

RESULTS

We characterized the cellular and molecular responses of Atlantic salmon and coho salmon to sea lice using single-nuclei RNA sequencing. Juvenile fish were exposed to copepodid sea lice (Lepeophtheirus salmonis), and lice-attached pelvic fin and skin samples were collected 12 h, 24 h, 36 h, 48 h, and 60 h after exposure, along with control samples. Comparative analysis of control and treatment samples revealed an immune and wound-healing response that was common to both species, but attenuated in Atlantic salmon, potentially reflecting greater sea louse immunomodulation. Our results revealed unique but complementary roles of three layers of keratinocytes in the epithelial hyperplasia response leading to rapid sea lice rejection in coho salmon. Our results suggest that basal keratinocytes direct the expansion and mobility of intermediate and, especially, superficial keratinocytes, which eventually encapsulate the parasite.

CONCLUSIONS

Our results highlight the key role of keratinocytes in coho salmon's sea lice resistance and the diverged biological response of the two salmonid host species when interacting with this parasite. This study has identified key pathways and candidate genes that could be manipulated using various biotechnological solutions to improve Atlantic salmon sea lice resistance.

Phased ERK-responsiveness and developmental robustness regulate teleost skin morphogenesis

Elongation of the vertebrate embryonic axis necessitates rapid expansion of the epidermis to accommodate the growth of underlying tissues. Here, we generated a toolkit to visualize and quantify signaling in entire cell populations of periderm, the outermost layer of the epidermis, in live developing zebrafish. We find that oriented cell divisions facilitate growth of the early periderm during axial elongation rather than cell addition from the basal layer. Activity levels of ERK, a downstream effector of MAPK pathway, gauged by a live biosensor, predicts cell cycle entry, and optogenetic ERK activation controls proliferation dynamics. As development proceeds, rates of peridermal cell proliferation decrease, ERK activity becomes more pulsatile and functionally transitions to promote hypertrophic cell growth. Targeted genetic blockade of cell division generates animals with oversized periderm cells, yet, unexpectedly, development to adulthood is not impaired. Our findings reveal stage-dependent differential responsiveness to ERK signaling and marked developmental robustness in growing teleost skin.

Damage-induced basal epithelial cell migration modulates the spatial organization of redox signaling and sensory neuron regeneration

Epithelial damage leads to early reactive oxygen species (ROS) signaling, which regulates sensory neuron regeneration and tissue repair. How the initial type of tissue injury influences early damage signaling and regenerative growth of sensory axons remains unclear. Previously we reported that thermal injury triggers distinct early tissue responses in larval zebrafish. Here, we found that thermal but not mechanical injury impairs sensory axon regeneration and function. Real-time imaging revealed an immediate tissue response to thermal injury characterized by the rapid Arp2/3-dependent migration of keratinocytes, which was associated with tissue-scale ROS production and sustained sensory axon damage. Isotonic treatment was sufficient to limit keratinocyte movement, spatially restrict ROS production and rescue sensory neuron function. These results suggest that early keratinocyte dynamics regulate the spatial and temporal pattern of long-term signaling in the wound microenvironment during tissue repair.

Cytoneme-mediated intercellular signaling in keratinocytes essential for epidermal remodeling

The skin, the largest organ, functions as a primary defense mechanism. Epidermal stem cells supply undifferentiated keratinocytes that differentiate as they migrate toward the outermost skin layer. Although such a replenishment process is disrupted in various human skin diseases, its underlying mechanisms remain elusive. With high-resolution live imaging and in vivo manipulations, we revealed that Notch signaling between keratinocytes is mediated by signaling filopodia called cytonemes and is essential for proper keratinocyte differentiation and proliferation. Inhibiting keratinocyte cytonemes reduced Notch expression within undifferentiated keratinocytes, leading to abnormal differentiation and hyperproliferation, resembling human skin disease phenotypes. Overproduction of Interleukin (IL)-17 signal, associated with skin diseases like psoriasis, induces psoriatic phenotypes via cytonemes in zebrafish. Our study suggests that intercellular signaling between keratinocytes through cytonemes is critical for epidermal maintenance, and its misregulation could be an origin of human skin diseases.

Live tracking of basal stem cells of the epidermis during growth, homeostasis and injury response in zebrafish

Basal stem cells of the epidermis continuously differentiate into keratinocytes and replenish themselves via self-renewal to maintain skin homeostasis. Numerous studies have attempted to reveal how basal cells undergo differentiation or self-renewal; however, this has been hampered by a lack of robust basal cell markers and analytical platforms that allow single-cell tracking. Here, we report that zebrafish integrin beta 4 is a useful marker for basal cell labelling, irrespective of the body region, stage and regenerative status. We employed Cre-loxP recombination in combination with live cell tracking of single basal clones in the caudal fin and investigated the embryonic origin and behaviour of basal cells during fish growth and homeostasis. Although most basal cells, including those in fins, became quiescent in the adult stage, genetic cell ablation showed that basal cells were reactivated to either self-renew or differentiate, depending on the injured cell type. Our study provides a simple and easy-to-use platform for quantitative in vivo imaging of basal stem cells at wider stages and under various conditions.

Environmental and molecular control of tissue-specific ionocyte differentiation in zebrafish

Organisms adjust their physiology to cope with environmental fluctuations and maintain fitness. These adaptations occur via genetic changes over multiple generations or through acclimation, a set of reversible phenotypic changes that confer resilience to the individual. Aquatic organisms are subject to dramatic seasonal fluctuations in water salinity, which can affect the function of lateral line mechanosensory hair cells. To maintain hair cell function when salinity decreases, ion-regulating cells, Neuromast-associated ionocytes (Nm ionocytes), increase in number and invade lateral line neuromasts. How environmental changes trigger this adaptive differentiation of Nm ionocytes and how these cells are specified is still unknown. Here, we identify Nm ionocyte progenitors as foxi3a/foxi3b-expressing skin cells and show that their differentiation is associated with sequential activation of different Notch pathway components, which control ionocyte survival. We demonstrate that new Nm ionocytes are rapidly specified by absolute salinity levels, independently of stress response pathways. We further show that Nm ionocyte differentiation is selectively triggered by depletion of specific ions, such as Ca2+ and Na+/Cl-, but not by low K+ levels, and is independent of media osmolarity. Finally, we demonstrate that hair cell activity plays a role in Nm ionocyte recruitment and that systemic factors are not necessary for Nm ionocyte induction. In summary, we have identified how environmental changes activate a signaling cascade that triggers basal skin cell progenitors to differentiate into Nm ionocytes and invade lateral line organs. This adaptive behavior is an example of physiological plasticity that may prove essential for survival in changing climates.

Transcriptomic landscape of Atlantic salmon (Salmo salar L.) skin

In this study, we present the first spatial transcriptomic atlas of Atlantic salmon skin using the Visium Spatial Gene Expression protocol. We utilized frozen skin tissue from 4 distinct sites, namely the operculum, pectoral and caudal fins, and scaly skin at the flank of the fish close to the lateral line, obtained from 2 Atlantic salmon (150 g). High-quality frozen tissue sections were obtained by embedding tissue in optimal cutting temperature media prior to freezing and sectioning. Further, we generated libraries and spatial transcriptomic maps, achieving a minimum of 80 million reads per sample with mapping efficiencies ranging from 79.3 to 89.4%. Our analysis revealed the detection of over 80,000 transcripts and nearly 30,000 genes in each sample. Among the tissue types observed in the skin, the epithelial tissues exhibited the highest number of transcripts (unique molecular identifier counts), followed by muscle tissue, loose and fibrous connective tissue, and bone. Notably, the widest nodes in the transcriptome network were shared among the epithelial clusters, while dermal tissues showed less consistency, which is likely attributable to the presence of multiple cell types at different body locations. Additionally, we identified collagen type 1 as the most prominent gene family in the skin, while keratins were found to be abundant in the epithelial tissue. Furthermore, we successfully identified gene markers specific to epithelial tissue, bone, and mesenchyme. To validate their expression patterns, we conducted a meta-analysis of the microarray database, which confirmed high expression levels of these markers in mucosal organs, skin, gills, and the olfactory rosette.

Transcriptomic profiling of tissue environments critical for post-embryonic patterning and morphogenesis of zebrafish skin

Pigment patterns and skin appendages are prominent features of vertebrate skin. In zebrafish, regularly patterned pigment stripes and an array of calcified scales form simultaneously in the skin during post-embryonic development. Understanding the mechanisms that regulate stripe patterning and scale morphogenesis may lead to the discovery of fundamental mechanisms that govern the development of animal form. To learn about cell types and signaling interactions that govern skin patterning and morphogenesis, we generated and analyzed single-cell transcriptomes of skin from wild-type fish as well as fish having genetic or transgenically induced defects in squamation or pigmentation. These data reveal a previously undescribed population of epidermal cells that express transcripts encoding enamel matrix proteins, suggest hormonal control of epithelial-mesenchymal signaling, clarify the signaling network that governs scale papillae development, and identify a critical role for the hypodermis in supporting pigment cell development. Additionally, these comprehensive single-cell transcriptomic data representing skin phenotypes of biomedical relevance should provide a useful resource for accelerating the discovery of mechanisms that govern skin development and homeostasis.

A genetically encoded sensor for visualizing leukotriene B4 gradients in vivo

Leukotriene B₄ (LTB₄) is a potent lipid chemoattractant driving inflammatory responses during host defense, allergy, autoimmune and metabolic diseases. Gradients of LTB₄ orchestrate leukocyte recruitment and swarming to sites of tissue damage and infection. How LTB₄ gradients form and spread in live tissues to regulate these processes remains largely elusive due to the lack of suitable tools for monitoring LTB₄ levels in vivo. Here, we develop GEM-LTB₄, a genetically encoded green fluorescent LTB₄ biosensor based on the human G-protein-coupled receptor BLT1. GEM-LTB₄ shows high sensitivity, specificity and a robust fluorescence increase in response to LTB₄ without affecting downstream signaling pathways. We use GEM-LTB₄ to measure ex vivo LTB₄ production of murine neutrophils. Transgenic expression of GEM-LTB₄ in zebrafish allows the real-time visualization of both exogenously applied and endogenously produced LTB₄ gradients. GEM-LTB₄ thus serves as a broadly applicable tool for analyzing LTB₄ dynamics in various experimental systems and model organisms.

Macrophage Migration Inhibitory Factor on Apoptotic Extracellular Vesicles Regulates Compensatory Proliferation

Apoptotic cells can signal to neighboring cells to stimulate proliferation and compensate for cell loss to maintain tissue homeostasis. While apoptotic cell-derived extracellular vesicles (AEVs) can transmit instructional cues to mediate communication with neighboring cells, the molecular mechanisms that induce cell division are not well understood. Here we show that macrophage migration inhibitory factor (MIF)-containing AEVs regulate compensatory proliferation via ERK signaling in epithelial stem cells of larval zebrafish. Time-lapse imaging showed efferocytosis of AEVs from dying epithelial stem cells by healthy neighboring stem cells. Proteomic and ultrastructure analysis of purified AEVs identified MIF localization on the AEV surface. Pharmacological inhibition or genetic mutation of MIF, or its cognate receptor CD74, decreased levels of phosphorylated ERK and compensatory proliferation in the neighboring epithelial stem cells. Disruption of MIF activity also caused decreased numbers of macrophages patrolling near AEVs, while depletion of the macrophage lineage resulted in a reduced proliferative response by the epithelial stem cells. We propose that AEVs carrying MIF directly stimulate epithelial stem cell repopulation and guide macrophages to cell non-autonomously induce localized proliferation to sustain overall cell numbers during tissue maintenance.

Cytoskeletal Keratins Are Overexpressed in a Zebrafish Model of Idiopathic Scoliosis

Idiopathic scoliosis (IS) is a three-dimensional rotation of the spine >10 degrees with an unknown etiology. Our laboratory established a late-onset IS model in zebrafish (Danio rerio) containing a deletion in kif7. A total of 25% of kif7^co63/co63 zebrafish develop spinal curvatures and are otherwise developmentally normal, although the molecular mechanisms underlying the scoliosis are unknown. To define transcripts associated with scoliosis in this model, we performed bulk mRNA sequencing on 6 weeks past fertilization (wpf) kif7^co63/co63 zebrafish with and without scoliosis. Additionally, we sequenced kif7^co63/co63, kif7^co63/+, and AB zebrafish (n = 3 per genotype). Sequencing reads were aligned to the GRCz11 genome and FPKM values were calculated. Differences between groups were calculated for each transcript by the t-test. Principal component analysis showed that transcriptomes clustered by sample age and genotype. kif7 mRNA was mildly reduced in both homozygous and heterozygous zebrafish compared to AB. Sonic hedgehog target genes were upregulated in kif7^co63/co63 zebrafish over AB, but no difference was detected between scoliotic and non-scoliotic mutants. The top upregulated genes in scoliotic zebrafish were cytoskeletal keratins. Pankeratin staining of 6 wpf scoliotic and non-scoliotic kif7^co63/co63 zebrafish showed increased keratin levels within the zebrafish musculature and intervertebral disc (IVD). Keratins are major components of the embryonic notochord, and aberrant keratin expression has been associated with intervertebral disc degeneration (IVDD) in both zebrafish and humans. The role of increased keratin accumulation as a molecular mechanism associated with the onset of scoliosis warrants further study.

A single-shot hyperspectral phasor camera for fast, multi-color fluorescence microscopy

Hyperspectral fluorescence imaging improves multiplexed observations of biological samples by utilizing multiple color channels across the spectral range to compensate for spectral overlap between labels. Typically, spectral resolution comes at a cost of decreased detection efficiency, which both hampers imaging speed and increases photo-toxicity to the samples. Here, we present a high-speed, high-efficiency snapshot spectral acquisition method, based on optical compression of the fluorescence spectra via Fourier transform, that overcomes the challenges of discrete spectral sampling: single-shot hyperspectral phasor camera (SHy-Cam). SHy-Cam captures fluorescence spatial and spectral information in a single exposure with a standard scientific CMOS camera, with photon efficiency of over 80%, easily and with acquisition rates exceeding 30 datasets per second, making it a powerful tool for multi-color in vivo imaging. Its simple design, using readily available optical components, and its easy integration provide a low-cost solution for multi-color fluorescence imaging with increased efficiency and speed.

Actomyosin contractility in olfactory placode neurons opens the skin epithelium to form the zebrafish nostril

Despite their barrier function, epithelia can locally lose their integrity to create physiological openings during morphogenesis. The mechanisms driving the formation of these epithelial breaks are only starting to be investigated. Here, we study the formation of the zebrafish nostril (the olfactory orifice), which opens in the skin epithelium to expose the olfactory neurons to external odorant cues. Combining live imaging, drug treatments, laser ablation, and tissue-specific functional perturbations, we characterize a mechanical interplay between olfactory placode neurons and the skin, which plays a crucial role in the formation of the orifice: the neurons pull on the overlying skin cells in an actomyosin-dependent manner which, in combination with a local reorganization of the skin epithelium, triggers the opening of the orifice. This work identifies an original mechanism to break an epithelial sheet, in which an adjacent group of cells mechanically assists the epithelium to induce its local rupture.

Dermal appendage-dependent patterning of zebrafish atoh1a+ Merkel cells

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.

Confined keratocytes mimic in vivo migration and reveal volume-speed relationship

Fish basal epidermal cells, known as keratocytes, are well-suited for cell migration studies. In vitro, isolated keratocytes adopt a stereotyped shape with a large fan-shaped lamellipodium and a nearly spherical cell body. However, in their native in vivo environment, these cells adopt a significantly different shape during their rapid migration toward wounds. Within the epidermis, keratocytes experience two-dimensional (2D) confinement between the outer epidermal cell layer and the basement membrane; these two deformable surfaces constrain keratocyte cell bodies to be flatter in vivo than in isolation. In vivo keratocytes also exhibit a relative elongation of the front-to-back axis and substantially more lamellipodial ruffling, as compared to isolated cells. We have explored the effects of 2D confinement, separated from other in vivo environmental cues, by overlaying isolated cells with an agarose hydrogel with occasional spacers, or with a ceiling made of polydimethylsiloxane (PDMS) elastomer. Under these conditions, isolated keratocytes more closely resemble the in vivo migratory shape phenotype, displaying a flatter apical-basal axis and a longer front-to-back axis than unconfined keratocytes. We propose that 2D confinement contributes to multiple dimensions of in vivo keratocyte shape determination. Further analysis demonstrates that confinement causes a synchronous 20% decrease in both cell speed and volume. Interestingly, we were able to replicate the 20% decrease in speed using a sorbitol hypertonic shock to shrink the cell volume, which did not affect other aspects of cell shape. Collectively, our results suggest that environmentally imposed changes in cell volume may influence cell migration speed, potentially by perturbing physical properties of the cytoplasm.

Zebrafish cornea formation and homeostasis reveal a slow maturation process, similarly to terrestrial vertebrates' corneas

Corneal blindness is the fourth leading cause of blindness worldwide. The superficial position of cornea on the eye makes this tissue prone to environmental aggressions, which can have a strong impact on sight. While most corneal pathology studies utilize terrestrial models, the knowledge on zebrafish cornea is too scarce to comprehend its strategy for the maintenance of a clear sight in aquatic environment. In this study, we deciphered the cellular and molecular events during corneal formation and maturation in zebrafish. After describing the morphological changes taking place from 3 days post fertilization (dpf) to adulthood, we analyzed cell proliferation. We showed that label retaining cells appear around 14 to 21dpf. Our cell proliferation study, combined to the study of Pax6a and krtt1c19e expression, demonstrate a long maturation process, ending after 45dpf. This maturation ends with a solid patterning of corneal innervation. Finally, we demonstrated that corneal wounding leads to an intense dedifferentiation, leading to the recapitulation of corneal formation and maturation, via a plasticity period. Altogether, our study deciphers the maturation steps of an aquatic cornea. These findings demonstrate the conservation of corneal formation, maturation and wound healing process in aquatic and terrestrial organisms, and they will enhance the use of zebrafish as model for corneal physiology studies.

DeepProjection: specific and robust projection of curved 2D tissue sheets from 3D microscopy using deep learning

The efficient extraction of image data from curved tissue sheets embedded in volumetric imaging data remains a serious and unsolved problem in quantitative studies of embryogenesis. Here, we present DeepProjection (DP), a trainable projection algorithm based on deep learning. This algorithm is trained on user-generated training data to locally classify 3D stack content, and to rapidly and robustly predict binary masks containing the target content, e.g. tissue boundaries, while masking highly fluorescent out-of-plane artifacts. A projection of the masked 3D stack then yields background-free 2D images with undistorted fluorescence intensity values. The binary masks can further be applied to other fluorescent channels or to extract local tissue curvature. DP is designed as a first processing step than can be followed, for example, by segmentation to track cell fate. We apply DP to follow the dynamic movements of 2D-tissue sheets during dorsal closure in Drosophila embryos and of the periderm layer in the elongating Danio embryo. DeepProjection is available as a fully documented Python package.

A non-disruptive and efficient knock-in approach allows fate tracing of resident osteoblast progenitors during repair of vertebral lesions in medaka

During bone development and repair, osteoblasts are recruited to bone deposition sites. To identify the origin of recruited osteoblasts, cell lineage tracing using Cre/loxP recombination is commonly used. However, a confounding factor is the use of transgenic Cre drivers that do not accurately recapitulate endogenous gene expression or the use of knock-in Cre drivers that alter endogenous protein activity or levels. Here, we describe a CRISPR/Cas9 homology-directed repair knock-in approach that allows efficient generation of Cre drivers controlled by the endogenous gene promoter. In addition, a self-cleaving peptide preserves the reading frame of the endogenous protein. Using this approach, we generated col10a1p2a-CreERT2 knock-in medaka and show that tamoxifen-inducible CreERT2 efficiently recombined loxP sites in col10a1 cells. Similar knock-in efficiencies were obtained when two unrelated loci (osr1 and col2a1a) were targeted. Using live imaging, we traced the fate of col10a1 osteoblast progenitors during bone lesion repair in the medaka vertebral column. We show that col10a1 cells at neural arches represent a mobilizable cellular source for bone repair. Together, our study describes a previously unreported strategy for precise cell lineage tracing via efficient and non-disruptive knock-in of Cre.

Skin cells undergo asynthetic fission to expand body surfaces in zebrafish

As an animal's surface area expands during development, skin cell populations must quickly respond to maintain sufficient epithelial coverage. Despite much progress in understanding of skin cell behaviours in vivo^1,2, it remains unclear how cells collectively act to satisfy coverage demands at an organismic level. Here we created a multicolour cell membrane tagging system, palmskin, to monitor the entire population of superficial epithelial cells (SECs) in developing zebrafish larvae. Using time-lapse imaging, we found that many SECs readily divide on the animal body surface; during a specific developmental window, a single SEC can produce a maximum of four progeny cells over its lifetime on the surface of the animal. Remarkably, EdU assays, DNA staining and hydroxyurea treatment showed that these terminally differentiated skin cells continue splitting despite an absence of DNA replication, causing up to 50% of SECs to exhibit reduced genome size. On the basis of a simple mathematical model and quantitative analyses of cell volumes and apical surface areas, we propose that 'asynthetic fission' is used as an efficient mechanism for expanding epithelial coverage during rapid growth. Furthermore, global or local manipulation of body surface growth affects the extent and mode of SEC division, presumably through tension-mediated activation of stretch-activated ion channels. We speculate that this frugal yet flexible mode of cell proliferation might also occur in contexts other than zebrafish skin expansion.

Embryologic and Fetal Development of the Eyelid and the Lacrimal Drainage System

The embryological and fetal morphogenesis of the eyelids and lacrimal drainage system is multifactorial and complex. This protracted process begins in the 5th week of prenatal life and involves a successive series of subtle and tightly regulated morphogenetic events. Major milestones of the embryological and fetal development of the eyelids include the beginning of eyelid formation during the 6th week, eyelid fusion by the 8th week, and the development of eyelid structures beginning in the 9th week (immediately following eyelid fusion), with progression until eyelid separation by the 24th week. After eyelid separation, the eyelids begin to assume their newly developed shape. Around the 32nd week, eyelids are almost fully developed and fully separated but still visibly closed. Key development steps of the lacrimal drainage system include formation of the lacrimal lamina (the primordium of the future lacrimal system) and the lacrimal cord (the primordium of the lacrimal canaliculi) in the 7th week, with canalization starting from the 10th week. During the 10th week, the excretory lacrimal system displays a lumen with a true lacrimal duct that can be distinguished. The epithelium of the lacrimal canaliculi is in contact with the palpebral conjunctival epithelium. The two epithelia form a continuous epithelial lamina. The caudal extreme of the lacrimal duct and the inferior meatal lamina join and the latter begins to cavitate. Understanding this multidimensional process of development in prenatal life, as well as identifying and linking signaling cascades and regulatory genes to existing diseases, may pave the way for developing minimally invasive interventions and scar reducing surgical methods, controlling the spread of malignancies, and the use of progenitor/stem cell and even regenerative therapy.

Shifting the focus of zebrafish toward a model of the tumor microenvironment

Cancer cells exist in a complex ecosystem with numerous other cell types in the tumor microenvironment (TME). The composition of this tumor/TME ecosystem will vary at each anatomic site and affects phenotypes such as initiation, metastasis, and drug resistance. A mechanistic understanding of the large number of cell-cell interactions between tumor and TME requires models that allow us to both characterize as well as genetically perturb this complexity. Zebrafish are a model system optimized for this problem, because of the large number of existing cell-type-specific drivers that can label nearly any cell in the TME. These include stromal cells, immune cells, and tissue resident normal cells. These cell-type-specific promoters/enhancers can be used to drive fluorophores to facilitate imaging and also CRISPR cassettes to facilitate perturbations. A major advantage of the zebrafish is the ease by which large numbers of TME cell types can be studied at once, within the same animal. While these features make the zebrafish well suited to investigate the TME, the model has important limitations, which we also discuss. In this review, we describe the existing toolset for studying the TME using zebrafish models of cancer and highlight unique biological insights that can be gained by leveraging this powerful resource.

The conundrum of pharyngeal teeth origin: the role of germ layers, pouches, and gill slits

There are several competing hypotheses on tooth origins, with discussions eventually settling in favour of an 'outside-in' scenario, in which internal odontodes (teeth) derived from external odontodes (skin denticles) in jawless vertebrates. The evolution of oral teeth from skin denticles can be intuitively understood from their location at the mouth entrance. However, the basal condition for jawed vertebrates is arguably to possess teeth distributed throughout the oropharynx (i.e. oral and pharyngeal teeth). As skin denticle development requires the presence of ectoderm-derived epithelium and of mesenchyme, it remains to be answered how odontode-forming skin epithelium, or its competence, were 'transferred' deep into the endoderm-covered oropharynx. The 'modified outside-in' hypothesis for tooth origins proposed that this transfer was accomplished through displacement of odontogenic epithelium, that is ectoderm, not only through the mouth, but also via any opening (e.g. gill slits) that connects the ectoderm to the epithelial lining of the pharynx (endoderm). This review explores from an evolutionary and from a developmental perspective whether ectoderm plays a role in (pharyngeal) tooth and denticle formation. Historic and recent studies on tooth development show that the odontogenic epithelium (enamel organ) of oral or pharyngeal teeth can be of ectodermal, endodermal, or of mixed ecto-endodermal origin. Comprehensive data are, however, only available for a few taxa. Interestingly, in these taxa, the enamel organ always develops from the basal layer of a stratified epithelium that is at least bilayered. In zebrafish, a miniaturised teleost that only retains pharyngeal teeth, an epithelial surface layer with ectoderm-like characters is required to initiate the formation of an enamel organ from the basal, endodermal epithelium. In urodele amphibians, the bilayered epithelium is endodermal, but the surface layer acquires ectodermal characters, here termed 'epidermalised endoderm'. Furthermore, ectoderm-endoderm contacts at pouch-cleft boundaries (i.e. the prospective gill slits) are important for pharyngeal tooth initiation, even if the influx of ectoderm via these routes is limited. A balance between sonic hedgehog and retinoic acid signalling could operate to assign tooth-initiating competence to the endoderm at the level of any particular pouch. In summary, three characters are identified as being required for pharyngeal tooth formation: (i) pouch-cleft contact, (ii) a stratified epithelium, of which (iii) the apical layer adopts ectodermal features. These characters delimit the area in which teeth can form, yet cannot alone explain the distribution of teeth over the different pharyngeal arches. The review concludes with a hypothetical evolutionary scenario regarding the persisting influence of ectoderm on pharyngeal tooth formation. Studies on basal osteichthyans with less-specialised types of early embryonic development will provide a crucial test for the potential role of ectoderm in pharyngeal tooth formation and for the 'modified outside-in' hypothesis of tooth origins.

Whole-body clonal mapping identifies giant dominant clones in zebrafish skin epidermis

Skin expansion during development is predominantly driven by growth of basal epithelial cell (BEC)-derived clonal populations, which often display varied sizes and shapes. However, little is known about the causes of clonal heterogeneity and the maximum size to which a single clone can grow. Here, we created a zebrafish model, basebow, for capturing clonal growth behavior in the BEC population on a whole-body, centimeter scale. By tracking 222 BECs over the course of a 28-fold expansion of body surface area, we determined that most BECs survive and grow clonal populations with an average size of 0.013 mm2. An extensive survey of 742 sparsely labeled BECs further revealed that giant dominant clones occasionally arise on specific body regions, covering up to 0.6% of the surface area. Additionally, a growth-induced extracellular matrix component, Lamb1a, mediates clonal growth in a cell-autonomous manner. Altogether, our findings demonstrate how clonal heterogeneity and clonal dominance may emerge to enable post-embryonic growth of a vertebrate organ, highlighting key cellular mechanisms that may only become evident when visualizing single cell behavior at the whole-animal level.

Metabolic Alterations in Preneoplastic Development Revealed by Untargeted Metabolomic Analysis

Metabolic rewiring is a critical hallmark of tumorigenesis and is essential for the development of cancer. Although many key features of metabolic alteration that are crucial for tumor cell survival, proliferation and progression have been identified, these are obtained from studies with established tumors and cancer cell lines. However, information on the essential metabolic changes that occur during pre-neoplastic cell (PNC) development that enables its progression to full blown tumor is still lacking. Here, we present an untargeted metabolomics analysis of human oncogene HRASG12V induced PNC development, using a transgenic inducible zebrafish larval skin development model. By comparison with normal sibling controls, we identified six metabolic pathways that are significantly altered during PNC development in the skin. Amongst these altered pathways are pyrimidine, purine and amino acid metabolism that are common to the cancer metabolic changes that support rapid cell proliferation and growth. Our data also suggest alterations in post transcriptional modification of RNAs that might play a role in PNC development. Our study provides a proof of principle work flow for identifying metabolic alterations during PNC development driven by an oncogenic mutation. In the future, this approach could be combined with transcriptomic or proteomic approaches to establish the detailed interaction between signaling networks and cellular metabolic pathways that occur at the onset of tumor progression.

Fxr signaling and microbial metabolism of bile salts in the zebrafish intestine

Bile salt synthesis, secretion into the intestinal lumen, and resorption in the ileum occur in all vertebrate classes. In mammals, bile salt composition is determined by host and microbial enzymes, affecting signaling through the bile salt-binding transcription factor farnesoid X receptor (Fxr). However, these processes in other vertebrate classes remain poorly understood. We show that key components of hepatic bile salt synthesis and ileal transport pathways are conserved and under control of Fxr in zebrafish. Zebrafish bile salts consist primarily of a C27 bile alcohol and a C24 bile acid that undergo multiple microbial modifications including bile acid deconjugation that augments Fxr activity. Using single-cell RNA sequencing, we provide a cellular atlas of the zebrafish intestinal epithelium and uncover roles for Fxr in transcriptional and differentiation programs in ileal and other cell types. These results establish zebrafish as a nonmammalian vertebrate model for studying bile salt metabolism and Fxr signaling.

Matriptase activation of Gq drives epithelial disruption and inflammation via RSK and DUOX

Epithelial tissues are primed to respond to insults by activating epithelial cell motility and rapid inflammation. Such responses are also elicited upon overexpression of the membrane-bound protease, Matriptase, or mutation of its inhibitor, Hai1. Unrestricted Matriptase activity also predisposes to carcinoma. How Matriptase leads to these cellular outcomes is unknown. We demonstrate that zebrafish hai1a mutants show increased H₂O₂, NfκB signalling, and IP₃R -mediated calcium flashes, and that these promote inflammation, but do not generate epithelial cell motility. In contrast, inhibition of the Gq subunit in hai1a mutants rescues both the inflammation and epithelial phenotypes, with the latter recapitulated by the DAG analogue, PMA. We demonstrate that hai1a has elevated MAPK pathway activity, inhibition of which rescues the epidermal defects. Finally, we identify RSK kinases as MAPK targets disrupting adherens junctions in hai1a mutants. Our work maps novel signalling cascades mediating the potent effects of Matriptase on epithelia, with implications for tissue damage response and carcinoma progression.

Cancer occurs when normal processes in the cell become corrupted or unregulated. Many proteins can contribute, including one enzyme called Matriptase that cuts other proteins at specific sites. Matriptase activity is tightly controlled by a protein called Hai1. In mice and zebrafish, when Hai1 cannot adequately control Matriptase activity, invasive cancers with severe inflammation develop. However, it is unclear how unregulated Matriptase leads to both inflammation and cancer invasion.

One outcome of Matriptase activity is removal of proteins called Cadherins from the cell surface. These proteins have a role in cell adhesion: they act like glue to stick cells together. Without them, cells can dissociate from a tissue and move away, a critical step in cancer cells invading other organs. However, it is unknown exactly how Matriptase triggers the removal of Cadherins from the cell surface to promote invasion.

Previous work has shown that Matriptase switches on a receptor called Proteinase-activated receptor 2, or Par2 for short, which is known to activate many enzymes, including one called phospholipase C. When activated, this enzyme releases two signals into the cell: a sugar called inositol triphosphate, IP3; and a lipid or fat called diacylglycerol, DAG. It is possible that these two signals have a role to play in how Matriptase removes Cadherins from the cell surface.

To find out, Ma et al. mapped the effects of Matriptase in zebrafish lacking the Hai1 protein. This revealed that Matriptase increases IP3 and DAG levels, which initiate both inflammation and invasion. IP3 promotes inflammation by switching on pro-inflammatory signals inside the cell such as the chemical hydrogen peroxide. At the same time, DAG promotes cell invasion by activating a well-known cancer signalling pathway called MAPK. This pathway activates a protein called RSK. Ma et al. show that this protein is required to remove Cadherins from the surface of cells, thus connecting Matriptase’s activation of phospholipase C with its role in disrupting cell adhesion.

An increase in the ratio of Matriptase to HAI-1 (the human equivalent of Hai1) is present in many cancers. For this reason, the signal cascades described by Ma et al. may be of interest in developing treatments for these cancers. Understanding how these signals work together could lead to more direct targeted anti-cancer approaches in the future.

Transparent Touch: Insights From Model Systems on Epidermal Control of Somatosensory Innervation

Somatosensory neurons (SSNs) densely innervate our largest organ, the skin, and shape our experience of the world, mediating responses to sensory stimuli including touch, pressure, and temperature. Historically, epidermal contributions to somatosensation, including roles in shaping innervation patterns and responses to sensory stimuli, have been understudied. However, recent work demonstrates that epidermal signals dictate patterns of SSN skin innervation through a variety of mechanisms including targeting afferents to the epidermis, providing instructive cues for branching morphogenesis, growth control and structural stability of neurites, and facilitating neurite-neurite interactions. Here, we focus onstudies conducted in worms (Caenorhabditis elegans), fruit flies (Drosophila melanogaster), and zebrafish (Danio rerio): prominent model systems in which anatomical and genetic analyses have defined fundamental principles by which epidermal cells govern SSN development.

Genome of the butterfly hillstream loach provides insights into adaptations to torrential mountain stream life

Butterfly hillstream loach (Beaufortia kweichowensis), a benthic fish in the torrential mountain streams, possesses a totally flat ventrum, flattened craniofacial and body skeletons, and enlarged paired fins covered by substantially small keratinous structures. However, little is known about the genetic basis of these specialized morphological adaptations. Here we present a 448.52-Mb genome assembly with contig N50 length of 5.53 Mb by integrating Illumina short-read sequencing, Nanopore long-read sequencing and HiC-based chromatin map. Demographic history reconstruction of the butterfly hillstream loach reveals that the population dynamics is correlated with the different stages of uplifting of the Tibetan Plateau. Comparative genomic analysis finds evidence of six keratin genes in butterfly hillstream loach evolving under positive selection. Within these genes, two keratin genes exhibit species-specific and divergent amino acid changes, suggesting a role in the formation of the unculi. Additionally, a series of positively selected genes, rapid evolving genes, specific variant genes and expanded gene families are found, including genes related to Hedgehog, Notch and BMP pathways, which may be involved in craniofacial development. These findings may have important implications for understanding the genetic basis of phenotypic adaptation to torrential mountain stream life.

Adaptive cell invasion maintains lateral line organ homeostasis in response to environmental changes

Mammalian inner ear and fish lateral line sensory hair cells (HCs) detect fluid motion to transduce environmental signals. Actively maintained ionic homeostasis of the mammalian inner ear endolymph is essential for HC function. In contrast, fish lateral line HCs are exposed to the fluctuating ionic composition of the aqueous environment. Using lineage labeling, in vivo time-lapse imaging and scRNA-seq, we discovered highly motile skin-derived cells that invade mature mechanosensory organs of the zebrafish lateral line and differentiate into Neuromast-associated (Nm) ionocytes. This invasion is adaptive as it is triggered by environmental fluctuations. Our discovery of Nm ionocytes challenges the notion of an entirely placodally derived lateral line and identifies Nm ionocytes as likely regulators of HC function possibly by modulating the ionic microenvironment. Nm ionocytes provide an experimentally accessible in vivo system to study cell invasion and migration, as well as the physiological adaptation of vertebrate organs to changing environmental conditions.

BICRA, a SWI/SNF Complex Member, Is Associated with BAF-Disorder Related Phenotypes in Humans and Model Organisms

SWI/SNF-related intellectual disability disorders (SSRIDDs) are rare neurodevelopmental disorders characterized by developmental disability, coarse facial features, and fifth digit/nail hypoplasia that are caused by pathogenic variants in genes that encode for members of the SWI/SNF (or BAF) family of chromatin remodeling complexes. We have identified 12 individuals with rare variants (10 loss-of-function, 2 missense) in the BICRA (BRD4 interacting chromatin remodeling complex-associated protein) gene, also known as GLTSCR1, which encodes a subunit of the non-canonical BAF (ncBAF) complex. These individuals exhibited neurodevelopmental phenotypes that include developmental delay, intellectual disability, autism spectrum disorder, and behavioral abnormalities as well as dysmorphic features. Notably, the majority of individuals lack the fifth digit/nail hypoplasia phenotype, a hallmark of most SSRIDDs. To confirm the role of BICRA in the development of these phenotypes, we performed functional characterization of the zebrafish and Drosophila orthologs of BICRA. In zebrafish, a mutation of bicra that mimics one of the loss-of-function variants leads to craniofacial defects possibly akin to the dysmorphic facial features seen in individuals harboring putatively pathogenic BICRA variants. We further show that Bicra physically binds to other non-canonical ncBAF complex members, including the BRD9/7 ortholog, CG7154, and is the defining member of the ncBAF complex in flies. Like other SWI/SNF complex members, loss of Bicra function in flies acts as a dominant enhancer of position effect variegation but in a more context-specific manner. We conclude that haploinsufficiency of BICRA leads to a unique SSRIDD in humans whose phenotypes overlap with those previously reported.

Endoderm-Derived Myeloid-like Metaphocytes in Zebrafish Gill Mediate Soluble Antigen-Induced Immunity

Immune cells in the mucosal barriers of vertebrates are highly heterogeneous in their origin and function. This heterogeneity is further exemplified by the recent discovery of ectoderm-derived immune cells-metaphocytes in zebrafish epidermis. Yet, whether non-hematopoiesis-derived immune cells generally exist in barrier tissues remains obscured. Here, we report the identification and characterization of an endoderm-derived immune cell population in the gill and intestine of zebrafish. Transcriptome analysis reveals that the endoderm-derived immune cells are myeloid-like cells with high similarities to the ectoderm-derived metaphocytes in epidermis. Like metaphocytes in epidermis, the endoderm-derived immune cells are non-phagocytic but professional in external soluble antigen uptake. Depletion of the endoderm-derived immune cells in gill hinder the local immune response to external soluble stimulants. This study demonstrates a general existence of non-hematopoiesis-derived immune cells in zebrafish mucosal barriers and challenges the prevalent view that resident immune cells in mucosal barriers arise exclusively from hematopoiesis.

The dynamics of epidermal stratification during post-larval development in zebrafish

BACKGROUND

The epidermis, as a defensive barrier, is a consistent trait throughout animal evolution. During post-larval development, the zebrafish epidermis thickens by stratification or addition of new cell layers. Epidermal basal stem cells, expressing the transcription factor p63, are known to be involved in this process. Zebrafish post-larval epidermal stratification is a tractable system to study how stem cells participate in organ growth.

METHODS

We used immunohistochemistry, in combination with EdU cell proliferation detection, to study zebrafish epidermal stratification. For this procedure, we selected a window of post-larval stages (5-8 mm of standard length or SL, which normalizes age by size). Simultaneously, we used markers for asymmetric cell division and the Notch signaling pathway.

RESULTS

We found that epidermal stratification is the consequence of several events, including changes in cell shape, active cell proliferation and asymmetrical cell divisions. We identified a subset of highly proliferative epidermal cells with reduced levels of p63, which differed from the basal stem cells with high levels of p63. Additionally, we described different mechanisms that participate in the stratification process, including the phosphorylation of p63, asymmetric cell division regulated by the Par3 and LGN proteins, and expression of Notch genes.

Entosis and apical cell extrusion constitute a tumor-suppressive mechanism downstream of Matriptase

Armistead et al. show that in a bilayered epithelium in vivo, apical cell extrusion of basal cells is achieved via their engulfment by surface cells. In zebrafish hai1a mutants, this constitutes a tumor-suppressive mechanism, revealing a double face of Matriptase.

The type II transmembrane serine protease Matriptase 1 (ST14) is commonly known as an oncogene, yet it also plays an understudied role in suppressing carcinogenesis. This double face is evident in the embryonic epidermis of zebrafish loss-of-function mutants in the cognate Matriptase inhibitor Hai1a (Spint1a). Mutant embryos display epidermal hyperplasia, but also apical cell extrusions, during which extruding outer keratinocytes carry out an entosis-like engulfment and entrainment of underlying basal cells, constituting a tumor-suppressive effect. These counteracting Matriptase effects depend on EGFR and the newly identified mediator phospholipase D (PLD), which promotes both mTORC1-dependent cell proliferation and sphingosine-1-phosphate (S1P)–dependent entosis and apical cell extrusion. Accordingly, hypomorphic hai1a mutants heal spontaneously, while otherwise lethal hai1a amorphs are efficiently rescued upon cotreatment with PLD inhibitors and S1P. Together, our data elucidate the mechanisms underlying the double face of Matriptase function in vivo and reveal the potential use of combinatorial carcinoma treatments when such double-face mechanisms are involved.

Cellular diversity of the regenerating caudal fin

Zebrafish faithfully regenerate their caudal fin after amputation. During this process, both differentiated cells and resident progenitors migrate to the wound site and undergo lineage-restricted, programmed cellular state transitions to populate the new regenerate. Until now, systematic characterizations of cells comprising the new regenerate and molecular definitions of their state transitions have been lacking. We hereby characterize the dynamics of gene regulatory programs during fin regeneration by creating single-cell transcriptome maps of both preinjury and regenerating fin tissues at 1/2/4 days post-amputation. We consistently identified epithelial, mesenchymal, and hematopoietic populations across all stages. We found common and cell type-specific cell cycle programs associated with proliferation. In addition to defining the processes of epithelial replenishment and mesenchymal differentiation, we also identified molecular signatures that could better distinguish epithelial and mesenchymal subpopulations in fish. The insights for natural cell state transitions during regeneration point to new directions for studying this regeneration model.

Analysis of zebrafish periderm enhancers facilitates identification of a regulatory variant near human KRT8/18

Genome-wide association studies for non-syndromic orofacial clefting (OFC) have identified single nucleotide polymorphisms (SNPs) at loci where the presumed risk-relevant gene is expressed in oral periderm. The functional subsets of such SNPs are difficult to predict because the sequence underpinnings of periderm enhancers are unknown. We applied ATAC-seq to models of human palate periderm, including zebrafish periderm, mouse embryonic palate epithelia, and a human oral epithelium cell line, and to complementary mesenchymal cell types. We identified sets of enhancers specific to the epithelial cells and trained gapped-kmer support-vector-machine classifiers on these sets. We used the classifiers to predict the effects of 14 OFC-associated SNPs at 12q13 near KRT18. All the classifiers picked the same SNP as having the strongest effect, but the significance was highest with the classifier trained on zebrafish periderm. Reporter and deletion analyses support this SNP as lying within a periderm enhancer regulating KRT18/KRT8 expression.

Stepwise polarisation of developing bilayered epidermis is mediated by aPKC and E-cadherin in zebrafish

The epidermis, a multilayered epithelium, surrounds and protects the vertebrate body. It develops from a bilayered epithelium formed of the outer periderm and underlying basal epidermis. How apicobasal polarity is established in the developing epidermis has remained poorly understood. We show that both the periderm and the basal epidermis exhibit polarised distribution of adherens junctions in zebrafish. aPKC, an apical polarity regulator, maintains the robustness of polarisation of E-cadherin- an adherens junction component- in the periderm. E-cadherin in one layer controls the localisation of E-cadherin in the second layer in a layer non-autonomous manner. Importantly, E-cadherin controls the localisation and levels of Lgl, a basolateral polarity regulator, in a layer autonomous as well non-autonomous manner. Since periderm formation from the enveloping layer precedes the formation of the basal epidermis, our analyses suggest that peridermal polarity, initiated by aPKC, is transduced in a stepwise manner by E-cadherin to the basal layer.

Tissue-Specific Transcriptomes Reveal Gene Expression Trajectories in Two Maturing Skin Epithelial Layers in Zebrafish Embryos

Epithelial cells are the building blocks of many organs, including skin. The vertebrate skin initially consists of two epithelial layers, the outer periderm and inner basal cell layers, which have distinct properties, functions, and fates. The embryonic periderm ultimately disappears during development, whereas basal cells proliferate to form the mature, stratified epidermis. Although much is known about mechanisms of homeostasis in mature skin, relatively little is known about the two cell types in pre-stratification skin. To define the similarities and distinctions between periderm and basal skin epithelial cells, we purified them from zebrafish at early development stages and deeply profiled their gene expression. These analyses identified groups of genes whose tissue enrichment changed at each stage, defining gene flow dynamics of maturing vertebrate epithelia. At each of 52 and 72 hr post-fertilization (hpf), more than 60% of genes enriched in skin cells were similarly expressed in both layers, indicating that they were common epithelial genes, but many others were enriched in one layer or the other. Both expected and novel genes were enriched in periderm and basal cell layers. Genes encoding extracellular matrix, junctional, cytoskeletal, and signaling proteins were prominent among those distinguishing the two epithelial cell types. In situ hybridization and BAC transgenes confirmed our expression data and provided new tools to study zebrafish skin. Collectively, these data provide a resource for studying common and distinguishing features of maturing epithelia.

Bacteria evoke alarm behaviour in zebrafish

When injured, fish release an alarm substance (Schreckstoff) that elicits fear in members of their shoal. Although Schreckstoff has been proposed to be produced by club cells in the skin, several observations indicate that these giant cells function primarily in immunity. Previous data indicate that the alarm substance can be isolated from mucus. Here we show that mucus, as well as bacteria, are transported from the external surface into club cells, by cytoplasmic transfer or invasion of cells, including neutrophils. The presence of bacteria inside club cells raises the possibility that the alarm substance may contain a bacterial component. Indeed, lysate from a zebrafish Staphylococcus isolate is sufficient to elicit alarm behaviour, acting in concert with a substance from fish. These results suggest that Schreckstoff, which allows one individual to unwittingly change the emotional state of the surrounding population, derives from two kingdoms and is associated with processes that protect the host from bacteria.

Periderm invasion contributes to epithelial formation in the teleost pharynx

The gnathostome pharyngeal cavity functions in food transport and respiration. In amniotes the mouth and nares are the only channels allowing direct contact between internal and external epithelia. In teleost fish, gill slits arise through opening of endodermal pouches and connect the pharynx to the exterior. Using transgenic zebrafish lines, cell tracing, live imaging and different markers, we investigated if pharyngeal openings enable epithelial invasion and how this modifies the pharyngeal epithelium. We conclude that in zebrafish the pharyngeal endoderm becomes overlain by cells with a peridermal phenotype. In a wave starting from pouch 2, peridermal cells from the outer skin layer invade the successive pouches until halfway their depth. Here the peridermal cells connect to a population of cells inside the pharyngeal cavity that express periderm markers, yet do not invade from outside. The latter population expands along the midline from anterior to posterior until the esophagus-gut boundary. Together, our results show a novel role for the periderm as an internal epithelium becomes adapted to function as an external surface.

Proteolytic and Opportunistic Breaching of the Basement Membrane Zone by Immune Cells during Tumor Initiation

Cancer-related inflammation impacts significantly on cancer development and progression. From early stages, neutrophils and macrophages are drawn to pre-neoplastic cells in the epidermis, but before directly interacting, they must first breach the underlying extracellular matrix barrier layer that includes the basement membrane. Using several different skin cancer models and a collagen I-GFP transgenic zebrafish line, we have undertaken correlative light and electron microscopy (CLEM) to capture the moments when immune cells traverse the basement membrane. We show evidence both for active proteolytic burrowing and for the opportunistic use of pre-existing weak spots in the matrix layer. We show that these small holes, as well as much larger, cancer cell-generated or wound-triggered gaps in the matrix barrier, provide portals for immune cells to access cancer cells in the epidermis and thus are rate limiting in cancer progression.

E-cadherin expression pattern during zebrafish embryonic epidermis development

Background: E-cadherin is the major adhesion receptor in epithelial adherens junctions (AJs). On established epidermis, E-cadherin performs fine-tuned cell-cell contact remodeling to maintain tissue integrity, which is characterized by modulation of cell shape, size and packing density. In zebrafish, the organization and distribution of E-cadherin in AJs during embryonic epidermis development remain scarcely described.

Methods: Combining classical immunofluorescence, deconvolution microscopy and 3D-segmentation of AJs in epithelial cells, a quantitative approach was implemented to assess the spatial and temporal distribution of E-cadherin across zebrafish epidermis between 24 and 72 hpf.

Results: increasing levels of E-cadh protein parallel higher cell density and the appearance of hexagonal cells in the enveloping layer (EVL) as well as the establishments of new cell-cell contacts in the epidermal basal layer (EBL), being significantly between 31 and 48 hpf .

Conclusions: Increasing levels of E-cadherin in AJs correlates with extensive changes in cell morphology towards hexagonal packing during the epidermis morphogenesis.

Ectopically Expressed IL-34 Can Efficiently Induce Macrophage Migration to the Liver in Zebrafish

Interleukin-34 (IL-34) is a known cytokine that plays an important role in the survival, proliferation, and differentiation of macrophages. In previous studies, IL-34 can induce macrophage migration through syndecan-1 or focal adhesion kinase and extracellular signal-related kinase 1 and 2 pathway. These studies mainly focused on in vitro experiments, but the effect of IL-34 on macrophage migration in vivo is less understood. In our study, we artificially induced macrophage, but not neutrophil, enrichment in the skin or liver by overexpressing IL-34 in epidermal cells or hepatocytes in zebrafish. Live imaging showed that the enrichment of macrophages in the liver is due to the direct attraction of macrophages by IL-34. Our results demonstrated that ectopically expressed IL-34 can induce macrophage migration to liver in vivo.

Vitamin D Stimulates Cardiomyocyte Proliferation and Controls Organ Size and Regeneration in Zebrafish

Attaining proper organ size during development and regeneration hinges on the activity of mitogenic factors. Here, we performed a large-scale chemical screen in embryonic zebrafish to identify cardiomyocyte mitogens. Although commonly considered anti-proliferative, vitamin D analogs like alfacalcidol had rapid, potent mitogenic effects on embryonic and adult cardiomyocytes in vivo. Moreover, pharmacologic or genetic manipulation of vitamin D signaling controlled proliferation in multiple adult cell types and dictated growth rates in embryonic and juvenile zebrafish. Tissue-specific modulation of vitamin D receptor (VDR) signaling had organ-restricted effects, with cardiac VDR activation causing cardiomegaly. Alfacalcidol enhanced the regenerative response of injured zebrafish hearts, whereas VDR blockade inhibited regeneration. Alfacalcidol activated cardiac expression of genes associated with ErbB2 signaling, while ErbB2 inhibition blunted its effects on cell proliferation. Our findings identify vitamin D as mitogenic for cardiomyocytes and other cell types in zebrafish and indicate a mechanism to regulate organ size and regeneration.

Generation and characterization of a zebrafish muscle specific inducible Cre line

Zebrafish transgenic lines provide valuable insights into gene functions, cell lineages and cell behaviors during development. Spatiotemporal control over transgene expression is a critical need in many experimental approaches, with applications in loss- and gain-of-function expression, ectopic expression and lineage tracing experiments. The Cre/loxP recombination system is a powerful tool to provide this control and the demand for validated Cre and loxP zebrafish transgenics is high. One of the major challenges to widespread application of Cre/loxP technology in zebrafish is comparatively small numbers of established tissue-specific Cre or CreERT2 lines. We used Tol2-mediated transgenesis to generate Tg(CrymCherry;-1.9mylz2:CreERT2) which provides an inducible CreERT2 source driven by muscle-specific mylz2 promoter. The transgenic specifically labels the trunk and tail skeletal muscles. We assessed the temporal responsiveness of the transgenic by screening with a validated loxP reporter transgenic ubi:Switch. Further, we evaluated the recombination efficiency in the transgenic with varying concentrations of 4-OHT, for different induction time periods and at different stages of embryogenesis and observed that higher recombination efficiency is achieved when embryos are induced with 10 μM 4-OHT from 10-somites or 24 hpf till 48 or 72 hpf. The transgenic is an addition to currently available zebrafish transgenesis toolbox and a significant tool to advance muscle biology studies in zebrafish.