Single-cell transcriptomic profiling of the zebrafish inner ear reveals molecularly distinct hair cell and supporting cell subtypes
Abstract
A major cause of human deafness and vestibular dysfunction is permanent loss of the mechanosensory hair cells of the inner ear. In non-mammalian vertebrates such as zebrafish, regeneration of missing hair cells can occur throughout life. While a comparative approach has the potential to reveal the basis of such differential regenerative ability, the degree to which the inner ears of fish and mammals share common hair cells and supporting cell types remains unresolved. Here, we perform single-cell RNA sequencing of the zebrafish inner ear at embryonic through adult stages to catalog the diversity of hair cells and non-sensory supporting cells. We identify a putative progenitor population for hair cells and supporting cells, as well as distinct hair and supporting cell types in the maculae versus cristae. The hair cell and supporting cell types differ from those described for the lateral line system, a distributed mechanosensory organ in zebrafish in which most studies of hair cell regeneration have been conducted. In the maculae, we identify two subtypes of hair cells that share gene expression with mammalian striolar or extrastriolar hair cells. In situ hybridization reveals that these hair cell subtypes occupy distinct spatial domains within the three macular organs, the utricle, saccule, and lagena, consistent with the reported distinct electrophysiological properties of hair cells within these domains. These findings suggest that primitive specialization of spatially distinct striolar and extrastriolar hair cells likely arose in the last common ancestor of fish and mammals. The similarities of inner ear cell type composition between fish and mammals validate zebrafish as a relevant model for understanding inner ear-specific hair cell function and regeneration.
Editor's evaluation
This important study describes transcriptomic profiles of sensory and non-sensory cells of the zebrafish inner ear at single-cell resolution in embryonic through adult stages. These solid results catalogue transcriptomic data and show evidence that distinct cell subtypes exist between cells of the ear and the lateral line as well as within subcellular compartments in the inner ear. These findings provide information towards comparative studies of inner ear hair cell function and regeneration.
https://doi.org/10.7554/eLife.82978.sa0Introduction
Mechanosensory hair cells of the inner ear are responsible for sensing sound and head position in vertebrates. Hair cells are notoriously susceptible to damage from multiple types of insults, including noise and ototoxic drug exposure. Studies of hair cell physiology in mammals are limited by the location of the inner ear within the temporal bone, which precludes many targeted manipulations and in vivo imaging beyond the neonatal stage. As a result, non-mammalian vertebrates with analogous, more easily accessible hair cells have become useful models for studying hair cell development, death, and regeneration. Non-mammalian vertebrates such as birds and fish can regenerate hair cells of the auditory and vestibular systems that are lost due to injury (Stone and Cotanche, 2007; Monroe et al., 2015). This differs from mammals, where cochlear hair cell death leads to permanent hearing loss (Corwin and Cotanche, 1988; Yamasoba and Kondo, 2006), and limited regeneration of vestibular hair cells results in minimal recovery of function (Golub et al., 2012). Non-mammalian model systems of hair cell regeneration have the potential to reveal conserved pathways that can be targeted to promote hair cell survival and regeneration in humans. However, the extent of hair cell molecular homology across vertebrates remains unclear.
Due to its accessibility for manipulation and imaging, the zebrafish lateral line system has been widely used to study mechanisms of hair cell physiology (Pickett and Raible, 2019; Sheets et al., 2021). The lateral line is an external sensory system that allows aquatic vertebrates to detect local movement of water. Sensory organs of the lateral line, called neuromasts, contain hair cells and supporting cells that share properties with those of the inner ear. However, relative to the lateral line, cells in the zebrafish inner ear are likely more similar to their mammalian counterparts, raising the potential for it to be a more comparable system in which to study hair cell function.
Zebrafish and mammals share several inner ear sensory organs. Three semicircular canals with sensory end organs called cristae sense angular rotation of the head. Two additional sensory end organs detect linear acceleration and gravity: the utricular and saccular macula each with an associated otolith crystal (Figure 1). Fish lack a specific auditory structure such as the mammalian cochlea and instead sense sound through the saccule, utricle, and a third otolith organ, the lagena. Although historically the utricle was thought to be for vestibular function and the saccule and lagena analogous to the cochlea for sound detection, there is now substantial evidence for all three otolith end organs being used for sound detection with diverse specializations across fishes (Popper and Fay, 1993). Zebrafish exhibit behavioral responses to sound frequencies between 100 and 1200 Hz (Zeddies and Fay, 2005; Bhandiwad et al., 2013), and neural responses up to 4000 Hz (Poulsen et al., 2021). In larval zebrafish, both saccule and utricle hair cells respond to vibration stimuli, with the utricle responding to relatively lower frequencies than the saccule, as well as additive effects when both are stimulated (Yao et al., 2016; Favre-Bulle et al., 2020).
Within the mammalian utricle and saccule, there are both morphological and spatial differences between hair cells (Lysakowski and Goldberg, 2004; Eatock and Songer, 2011). Hair cells are broadly classified by their morphology and innervation, with Type I hair cells having calyx synapses surrounding the hair cell body and Type II hair cells having bouton synapses. Both Type I and Type II cells can be found within the central region of the macular organs known as the striola and in the surrounding extrastriolar zones. Although the role of spatial segregation into striolar versus extrastriolar zones has not been fully elucidated, hair cells across these regions vary in morphology, electrophysiology, and synaptic structure (Desai et al., 2005; Li et al., 2008). The striola is characterized by hair cells with taller ciliary bundles and encompasses a line of polarity reversal where hair cells change their stereocilia orientation (Figure 1E). Whereas distinct Type I and Type II hair cells, and in particular the calyx synapses typical of Type I cells, have not been identified in the maculae of fishes, afferent innervation with some calyx-like properties has been reported in goldfish cristae (Lanford and Popper, 1996). Spatial heterogeneity in the maculae, including those of zebrafish, has also been previously noted (Chang et al., 1992; Platt, 1993; Collin et al., 2000; Liu et al., 2022). However, the homologies of cells at the cellular and molecular levels have remained unknown.
Recent single-cell and single-nucleus RNA-sequencing efforts have generated a wealth of transcriptomic data from hair cells in several model systems, facilitating more direct comparison of cell types and gene regulatory networks between species. Although single-cell transcriptomic data have recently been published for the zebrafish inner ear (Jimenez et al., 2022; Qian et al., 2022), the diversity of hair cell and supporting cell subtypes has not been thoroughly analyzed. In order to better understand the diversification of cell types in the zebrafish inner ear, and their relationships to those in mammals, here we perform single-cell and single-nucleus RNA sequencing of the zebrafish inner ear from embryonic through adult stages. We find that hair and supporting cells from the zebrafish inner ear and lateral line are transcriptionally distinct, and that hair and supporting cells differ between the cristae and maculae. All of these distinct cell types are present during larval development and are maintained into adulthood. In situ hybridization reveals that these hair cell subtypes occupy distinct spatial domains within the utricle, saccule, and lagena, and computational comparison of hair cell types reveals homology with striolar and extrastriolar hair cell types in mammals. These findings point to an origin of striolar and extrastriolar hair cell types in at least the last common ancestor of fish and mammals.
Results
Inner ear hair cells and supporting cells are distinct from those of the lateral line
To assess differences between inner ear and lateral line cells, we analyzed a subset of cells from a large single-nucleus RNA-seq dataset of whole zebrafish at embryonic and larval stages (36–96 hours post-fertilization (hpf)), which was prepared by single-nucleus combinatorial indexing and sequencing (‘sci-Seq’; Saunders et al., 2022). Within an initial dataset of 1.25 million cells from 1233 embryos spanning 18 timepoints between 18 and 96 hr (see Saunders et al., 2022 for more detail), a total of 16,517 inner ear and lateral line cells were isolated, combined, and re-processed using Monocle 3 (Figure 2A–B). Initially, otic vesicle and lateral line cell clusters were identified by eya1 expression (Sahly et al., 1999) in combination with the following known marker genes. Inner ear nonsensory cells were identified by expression of the transcription factor gene sox10 (Dutton et al., 2009) in combination with inner ear supporting cell genes (stm, otog, otogl, otomp, tecta, and oc90; Figure 2C; Söllner et al., 2003; Kalka et al., 2019; Petko et al., 2008; Stooke-Vaughan et al., 2015). Lateral line nonsensory cells were identified by expression of known markers fat1b, tfap2a, tnfsf10l3, lef1, cxcr4b, fgfr1a, and hmx3a (Figure 2D; Steiner et al., 2014; Thomas and Raible, 2019; McGraw et al., 2011; Haas and Gilmour, 2006; Lee et al., 2016; Feng and Xu, 2010). We identified hair cells by expression of the pan-hair cell genes otofb, cdh23, pcdh15a, ush1c, myo7aa, slc17a8, and cacna1da (Figure 2E; Chatterjee et al., 2015; Söllner et al., 2004; Seiler et al., 2005; Phillips et al., 2011; Ernest et al., 2000; Obholzer et al., 2008; Sheets et al., 2012). To distinguish between inner ear and lateral line hair cells, we queried expression of previously described markers for inner ear (gpx2, kifl, strc, and lhfpl5a) and lateral line (strc1, lhfpl5b, and s100t) (Erickson et al., 2019; Erickson and Nicolson, 2015). Although many of these markers are at low abundance, these populations are marked distinctly by strc and s100t (Figure 2F). We used Monocle3 to identify differentially expressed genes (Supplementary file 1) and to generate modules of co-expressed genes (Figure 2—figure supplement 1, Supplementary file 2).
Both hair cells and nonsensory supporting cells from the inner ear and lateral line formed distinct clusters, with nonsensory cells from the two mechanosensory organs showing greater distinction than hair cells (Figure 2B, Figure 2—figure supplement 2A). To confirm the relative differences between inner ear and lateral line hair cells and nonsensory cells, Partition-based Graph Abstraction (PAGA) analysis was used to measure the connectivity of clusters (Wolf et al., 2019). PAGA analysis revealed strong connectivity within inner ear supporting cell clusters and within lateral line supporting cell clusters but little connectivity between them (Figure 2—figure supplement 2A, Supplementary file 3).
The inner ear nonsensory cluster includes structural cells forming the otic capsule, identified by expression of the extracellular matrix protein-encoding genes collagen type 2 a1a (col2a1a) and matrilin 4 (matn4) (Xu et al., 2018), as well as sensory supporting cells expressing lfng (Figure 3D; Figure 2—figure supplement 2B). Inner ear and lateral line supporting cells remain as distinct clusters even when structural matn4+ cells are excluded from analysis (Figure 2—figure supplement 2C). Thus, both hair cells and supporting cells have distinct gene expression profiles between the inner ear and lateral line at embryonic and larval stages.
Single-cell RNA-seq reveals distinct hair cell and supporting cell populations in the juvenile and adult inner ear of zebrafish
To identify distinct subtypes of inner ear hair cells and supporting cells from larval through adult stages, we first re-analyzed single-cell RNA sequencing (scRNA-seq) datasets from larval stages (72 and 120 hpf) (Fabian et al., 2022), in which otic placode cells and their descendants were labeled with Sox10:Cre to induce recombination of an ubiquitous ubb:LOXP-EGFP-STOP-LOXP-mCherry transgene (Kague et al., 2012). We also performed additional scRNA-seq using these transgenic lines by dissecting ears from juvenile (14 days post-fertilization (dpf)), and adult (12 months post-fertilization (mpf)) animals. Following cell dissociation and fluorescence-activated cell sorting (FACS) to purify mCherry + cells, we constructed scRNA-seq libraries using 10x Chromium technology. For all datasets, hair cells and supporting cells were identified for further analysis based on the expression of hair cell markers myo6b and strc and supporting cell markers stm and lfng; structural cells were removed from further analysis based on expression of matn4 and col2a1a (Figure 3—figure supplement 1). Using Seurat, we integrated this dataset with the sci-Seq embryonic and larval dataset (36–96 hpf) (Figure 3A and B). The combined dataset comprises 3246 inner ear cells separated into 10 groups based on unsupervised clustering, with differentially expressed genes for each cluster shown in Figure 3E and Supplementary file 4. We identified six clusters of hair cells based on shared expression of myo6b, strc, lhfpl5a, and gfi1aa (Yu et al., 2020), a nascent hair cell cluster based on expression of atoh1a (Millimaki et al., 2007) and the Notch ligand dla (Riley et al., 1999), and two clusters of supporting cells based on expression of lfng and stm (Figure 3C and D, Figure 3—figure supplement 2). An additional putative progenitor cluster (cluster 0), enriched for cells from embryonic stages, is characterized by expression of genes such as fgfr2 (Rohs et al., 2013), fat1a (Down et al., 2005), igsf3, and pard3bb (Figure 3—figure supplement 3). Although these marker genes are differentially expressed in the putative progenitor cluster, some of them (e.g. fat1a and pard3bb) retain a lower expression level in supporting cell populations (Figure 3—figure supplement 3F). This is further demonstrated by gene modules of these clusters (Figure 3—figure supplement 4, Supplementary file 5), where the progenitor signature module genes (Module 1) are expressed in lower levels in the supporting cell clusters. This transcriptional relatedness between progenitors and supporting cells may underlie the role of supporting cells as a resident stem cell population during zebrafish hair cell regeneration.
Developmental trajectories in the inner ear
To understand potential lineage relationships between clusters, we performed pseudotime trajectory analysis using Monocle3. We anchored the pseudotime projection at the putative progenitor cell cluster. Analysis revealed two major trajectories toward hair cell and supporting cell clusters for both maculae and cristae (Figure 4A and B, Figure 4—figure supplement 1), with distinct patterns of gene expression along each trajectory (Supplementary file 6). We find that average gene expression of the putative progenitor (Cluster 0) markers follow two patterns: decreasing along both hair cell and supporting cell trajectories (fgfr2 and igsf3) and decreasing only along the hair cell trajectory (fat1a and pard3bb) (Figure 4C and D, Figure 4—figure supplement 1B and C). The hair cell trajectory progresses first through a stage marked by expression of dla and then atoh1a (Cluster 2, Figure 4E, Figure 4—figure supplement 1D). Concurrent with decreasing expression of nascent hair cell genes, we observe increasing expression of mature hair cell genes gfi1aa and myo6b (Figure 4F, Figure 4—figure supplement 1E). Along the supporting cell trajectory we observed upregulation of supporting cell-specific markers, including stm and lfng (Figure 4G, Figure 4—figure supplement 1F). These bifurcating lineage trajectories from Cluster 0 (Figure 4A) to hair and supporting cell clusters are consistent with the identification of Cluster 0 as a population of bipotent progenitors regulated by Notch signaling during early development (Haddon et al., 1998; Riley et al., 1999). To localize these developmental stages in vivo, we examined dla expression by in situ hybridization (Figure 4—figure supplement 2). We find that dla is expressed in supporting cells adjacent to myo6:GFP hair cells in both cristae and maculae, consistent with peripheral addition of new cells at the margins of the sensory patches.
Distinct supporting cell types in the cristae versus maculae
Supporting cells comprise two major clusters that can be distinguished by expression of tectb and zpld1a among other genes (Figure 3C, see Supplementary file 7 for differentially expressed genes). The tectb gene encodes Tectorin beta, a component of the tectorial membrane associated with cochlear hair cells in mammals (Goodyear et al., 2017), and a component of otoliths in zebrafish (Kalka et al., 2019). The zpld1a gene, encoding Zona-pellucida-like domain containing protein 1 a, is expressed in the cristae in fish (Dernedde et al., 2014; Yang et al., 2011) and mouse (Vijayakumar et al., 2019). Using fluorescent in situ hybridization, we find that tectb is expressed in the macular organs but not cristae, and zpld1a is expressed in cristae but not maculae (Figure 5C and D). Neither were detected in lateral line neuromasts (Figure 5C and D), showing they are inner ear-specific genes. Both tectb and zpld1a are expressed primarily in supporting cells, as they show little overlap in expression with the hair cell marker myo6b:GFP, similar to expression of the supporting cell marker lfng (Figure 5B–D, Figure 5—figure supplement 1). These results demonstrate the presence of distinct supporting cell subtypes for the maculae and cristae.
Distinct types of hair cells in the zebrafish inner ear
While inner ear and lateral line hair cells share many structural and functional features, we sought to determine if these cells also have distinct molecular signatures. We compared published datasets of lateral line hair cells (Baek et al., 2022; Kozak et al., 2020; Ohta et al., 2020) to our data, restricting analysis to datasets generated by 10x Chromium preparation to avoid technical batch effects across studies. Using Scanorama for alignments (Hie et al., 2019), hair cells from the inner ear and lateral line form distinct clusters, with a number of differentially expressed genes (Figure 2—figure supplement 3), including the known markers for lateral line (s100t) and inner ear (strc) (Figure 2). This analysis suggests that inner ear hair cells of the maculae and cristae are more similar to each other than to lateral line hair cells.
Within the maculae and cristae, we find that hair cells can be subdivided into two major groups (clusters 1 and 3 versus cluster 4). These clusters are distinguished by differential expression of a number of genes including two calcium binding protein genes, cabp1b and cabp2b (Di Donato et al., 2013; Figure 3E). Hair cell cluster 5 has a mixed identity with co-expression of a number of genes shared between these two groups, including cabp1b and cabp2b.
We next tested the in vivo expression of genes in each cluster using in situ hybridization, choosing cabp1b and cabp2b as representative markers for each cluster (Figure 6A). In the larval cristae, utricle, and saccule, cabp1b and cabp2b mark myo6b+hair cells in largely non-overlapping zones (Figure 6B–D). By adult stages, complementary domains of cabp1b+and cabp2b+hair cells become clearly apparent (Figure 6E–K). In the adult utricle, a central crescent of cabp2b+; myo6b+hair cells is surrounded by a broad domain of cabp1b+; myo6b+hair cells. In the saccule and lagena, a late developing sensory organ, central cabp2b+; myo6b+hair cells are surrounded by peripheral cabp1b+; myo6b+hair cells. We also find several genes that are specific for hair cells in the cristae, utricle, or saccule (Figure 7A). These include the calcium binding protein gene cabp5b in the cristae, the transcription factor skor2 in the utricle, and the deafness gene loxhd1b in the saccule (Figure 7B–D, Figure 7—figure supplement 1).
The domain organization of hair cells in the adult macular organs resembles that of striolar and extrastriolar hair cells in the mammalian utricle. We therefore examined expression of pvalb9, the zebrafish ortholog of the mouse striolar hair cell marker Ocm (Hoffman et al., 2018; Jiang et al., 2017; Figure 8, Figure 8—figure supplement 1). In the larval utricle, we observe near complete overlap of pvalb9 with cabp2b (Figure 8B–D). In the adult utricle, there is substantial overlap of pvalb9 with cabp2b expression (except for a thin strip of pvalb9+; cabp2b- cells), and little overlap with cabp1b expression (Figure 8F and G). In addition, anti-Spectrin staining of hair bundles reveals a line of polarity reversal within the cabp2b+domain of the utricle (Figure 8H, I), consistent with polarity reversal occurring within the striolar domains of mammalian macular organs (Li et al., 2008). Cluster 1/3 (cabp1b+) and Cluster 4 (cabp2b+) populations also differentially express genes related to stereocilia tip link and mechanotransduction channel components (Figure 8—figure supplement 2, Supplementary file 8) and various calcium and potassium channels (Figure 8—figure supplement 3, Supplementary file 8). We also note that the utricle marker skor2 labels primarily extrastriolar hair cells within this end organ, with loxhd1b labeling striolar hair cells within the saccule. These findings suggest that zebrafish Cluster 4 (cabp2b+) and Cluster 1/3 (cabp1b+) hair cells largely correspond to striolar and extrastriolar hair cells, respectively, with distinct mechanotransduction and synaptic properties.
Global homology of striolar and extrastriolar hair cells between fish and mice
To further probe similarities between zebrafish Cluster 4 (cabp2b+) and Cluster 1/3 (cabp1b+) hair cells versus striolar and extrastriolar hair cells in mammals, we utilized the Self-Assembling Manifold mapping (SAMap) algorithm (Tarashansky et al., 2021; Musser et al., 2021) to compare cell types across distant species. A strength of this algorithm is that it compares not only homologous gene pairs but also close paralogs, which is especially useful considering the extensive paralog switching observed between vertebrate clades (Postlethwait, 2007), as well as the extra round of genome duplication in the teleost lineage leading to zebrafish. When comparing adult zebrafish maculae with the postnatal mouse utricle (Jan et al., 2021), we find the highest alignment score between supporting cells (Figure 9A). Consistent with the spatial domains revealed by our in situ gene expression analysis, we find that mouse striolar Type I hair cells exclusively map to zebrafish Cluster 4 (cabp2b+) hair cells, and mouse extrastriolar Type I and Type II hair cells predominantly to zebrafish Cluster 1/3 (cabp1b+) hair cells. In contrast, zebrafish lateral line hair cells (Lush et al., 2019) align exclusively to mouse extrastriolar and not striolar hair cells (Figure 9—figure supplement 1). The small degree of mapping of mouse extrastriolar Type I hair cells to zebrafish Cluster 4 (cabp2b+) hair cells suggests that zebrafish Cluster 4 (cabp2b+) hair cells may have more of a Type I identity than Cluster 1/3 (cabp1b+) cells in general. Gene pairs driving the homology alignment include striolar markers Ocm, Loxhd1, and Atp2b2 for zebrafish Cluster 4 (cabp2b+) hair cells, and mouse extrastriolar markers Tmc1, Atoh1, and Jag2 for zebrafish Cluster 1/3 (cabp1b+) hair cells (Supplementary file 9). Thus, zebrafish Cluster 4 (cabp2b+) macular hair cells are closely related to striolar cells of the mouse utricle, with zebrafish lateral line and Cluster 1/3 (cabp1b+) macular hair cells more closely related to mouse extrastriolar hair cells.
A recent single-cell study revealed distinct central versus peripheral hair cell subpopulations in postnatal mouse cristae, reminiscent of the striolar and extrastriolar populations in the maculae (Wilkerson et al., 2021). As our zebrafish cristae hair cells also separate into distinct clusters, Cluster 9 (cabp1b+) and Cluster 8 (cabp2b+) (Figure 6A and B), we performed SAMap analysis between the crista cell populations of the two species to investigate cell type homology. Similar to what we observed for the utricle, zebrafish centrally located Cluster 8 crista hair cells predominantly map to mouse central crista hair cells, and zebrafish peripherally located Cluster 9 crista hair cells exclusively map to mouse peripheral crista hair cells (Figure 9B, see Supplementary file 10 for differentially expressed genes in Cluster 8 and Cluster 9 hair cells and Supplementary file 11 for gene pairs driving homology). Conserved types of spatially segregated HCs therefore exist in both the maculae and cristae of zebrafish and mouse.
Discussion
Our single-cell transcriptomic profiling of the embryonic to adult zebrafish inner ear reveals a diversity of hair cell and supporting cell subtypes that differ from those of the lateral line. As much of our knowledge about zebrafish hair cell regeneration comes from studies of the lateral line, understanding similarities and differences between the lateral line and inner ear has the potential to uncover mechanisms underlying the distinct regenerative capacity of inner ear hair cell subtypes. Recent tools to systematically damage inner ear hair cells in zebrafish (Jimenez et al., 2021) should enable such types of comparative studies.
We identify hair cells and supporting cells specific for maculae versus cristae, as well as two spatially segregated types of zebrafish inner ear hair cells with similarities to mammalian striolar and extrastriolar hair cells. These molecular signatures are conserved across larval and adult stages. However, consistent with other recent work (Jimenez et al., 2022; Qian et al., 2022), we were not able to resolve distinct clusters of hair cells or supporting cells corresponding to the distinct types of maculae: i.e. utricle, saccule, and lagena.
The division of auditory and vestibular function across the otolith organs in zebrafish remains somewhat unclear. The saccule is thought to act as the primary auditory organ of larval zebrafish, as the utricle is not necessary for sound detection above low frequencies (Yao et al., 2016). In the zebrafish adult, excess sound exposure can damage the saccule, while damage to the utricle is unknown (Schuck and Smith, 2009). Conversely, the utricle is critical for larval vestibular function, while input from the saccule is unnecessary (Riley and Moorman, 2000). However, there is contrasting evidence for overlap in function of both saccule and utricle for sound detection in larvae (Favre-Bulle et al., 2020; Poulsen et al., 2021). Currently we are not able to identify clearly distinct hair cell types in the utricle compared to the saccule that might reflect functional differences; whether such genetic signatures exist remains an important question that will require further in-depth analysis. It is interesting to note that mammalian vestibular end organs are also capable of responding to high-frequency sound stimuli (reviewed in Curthoys, 2017), suggesting that sound detection by hair cells may not be linked to a distinct end organ-specific molecular signature.
Our study supports zebrafish possessing distinct types of striolar and extrastriolar hair cells in the maculae and cristae, with molecular differences between these subtypes implying different physiological properties. In the zebrafish utricle, vibration is preferentially transduced by striolar cells while static tilt is received by extrastriolar cells (Tanimoto et al., 2022). Consistent with use of a s100s-hs:tdTomato transgene to mark striolar cells in this previous study, s100s is a highly specific marker for our striolar hair cell cluster (Figure 3E). We also find zebrafish striolar and extrastriolar hair cell subtypes express distinct combinations of ion channel genes and mechanotransduction components, consistent with previous reports of distinct current profiles in central versus peripheral hair cells in the zebrafish utricle, saccule, and lagena (Haden et al., 2013; Olt et al., 2014), as well as spatial differences in ciliary bundle morphology and synaptic innervation in the larval zebrafish utricle (Liu et al., 2022). The distinct spatial distribution, channel expression, and hair bundle morphologies in these hair cells resembles the known spatial, electrophysiological, and hair bundle compositional differences seen in the striolar versus extrastriolar hair cells in the amniote vestibular end organs (Holt et al., 2007; Kharkovets et al., 2000; Lapeyre et al., 1992; Meredith and Rennie, 2016; Moravec and Peterson, 2004; Rüsch et al., 1998; Xue and Peterson, 2006).
In each of the zebrafish end organs, striolar and extrastriolar hair cells can be defined by differential expression of calcium binding proteins, in particular cabp1b versus cabp2b. As these calcium binding proteins closely interact with synaptic calcium channels (Cui et al., 2007; Picher et al., 2017) with potential functionally different consequences (Yang et al., 2018), their differential expression may confer unique electrophysiological properties to each cell type. Mutations in human CABP2 associated with the autosomal recessive locus DFNB93 result in hearing loss (Schrauwen et al., 2012; Picher et al., 2017), underlining its functional importance. Even though we chose cabp1b and cabp2b as characteristic markers for zebrafish extrastriolar and striolar regions, it is worth noting that Cabp2, but not Cabp1, is expressed in all mouse postnatal utricular hair cells with differentially higher expression in the striola (Jan et al., 2021). Of note, lateral line hair cells express higher levels of cabp2b than cabp1b (Lush et al., 2019), despite our analysis suggesting that they are more closely related to extrastriolar hair cells. These observations emphasize the importance of examining global patterns of gene expression rather than individual markers when assigning homology of cell types.
By contrast, we found no clear homology of zebrafish inner ear hair cells with mammalian Type I and Type II hair cells. The lack of molecular signatures corresponding to Type I hair cells is consistent with previous reports that one of their major features, calyx synapses, are absent from macular organs in fishes (Lysakowski and Goldberg, 2004, but see Lanford and Popper, 1996 for evidence of calyx synapses in goldfish cristae). These findings suggest that the diversification of inner ear hair cells into Type I and Type II cells may have largely emerged after the evolutionary split of ray-finned fishes from the lineage leading to mammals.
We recognize that identifying cell type homology across tissues and species through molecular analysis has several potential caveats. Although we have collected transcriptomic data from the zebrafish inner ear from a wide range of developmental stages, we are limited by the fact that the publicly available datasets for zebrafish lateral line and mouse utricle and cristae are restricted to immature stages. Thus, cell maturity could be a confounder in our analyses. However, when we limited the comparison of lateral line hair cells and postnatal mouse vestibular hair cells to 3–5 dpf inner ear hair cells, we see similar alignments as when we used our 12 mpf data (Figure 9—figure supplement 1). In addition, we collected fewer supporting cells from adult zebrafish than expected, skewing cell type representation towards hair cells (Figure 3C). Thus, additional optimization may be needed to further interrogate the cell subtypes within zebrafish inner ear supporting cell populations.
Nonetheless, our integrated dataset reveals distinct molecular characteristics of hair cells and supporting cells in the zebrafish inner ear sensory organs, with conservation of these patterns from larval stages to adults. Although not discussed in detail here, our data include additional cell populations of the zebrafish inner ear that express extracellular matrix-associated genes important for otic capsule structure and ion channel-associated genes associated with fluid regulation. These data form a resource that can be further explored to inform molecular aspects of hair cell electrophysiology, mechanotransduction, sound versus motion detection, maintenance of inner ear structure and ionic balance, and inner ear-specific hair cell regeneration.
Materials and methods
Zebrafish lines
Request a detailed protocolThis study was performed in strict accordance with the recommendations in the Guide for the Care and Use of Laboratory Animals of the National Institutes of Health. The Institutional Animal Care and Use Committees of the University of Southern California (Protocol 20771) and University of Washington (Protocol 2997–01) approved all animal experiments. Experiments were performed on zebrafish (Danio rerio) of AB or mixed AB/Tubingen background. For adult stages, mixed sexes of animals were used for constructing single-cell libraries, as well as RNAScope experiments. Published lines include Tg(Mmu.Sox10-Mmu.Fos:Cre)zf384 (Kague et al., 2012); Tg(–3.5ubb:LOXP-EGFP-STOP-LOXP-mCherry)cz1701Tg (Mosimann et al., 2011); and Tg(myosin 6b:GFP)w186 (Hailey et al., 2017).
In situ hybridization and RNAScope
Request a detailed protocolHybridization chain reaction in situ hybridizations (Molecular Instruments, HCR v3.0) were performed on 5 dpf myo6b:GFP larvae as directed for whole-mount zebrafish embryos and larvae (Choi et al., 2016; Choi et al., 2018). Briefly, embryos were treated with 1-phenyl 2-thiourea (PTU) beginning at 24 hpf. At 5 dpf, larvae were fixed in 4% PFA overnight at 4 °C. Larvae were washed with PBS and then stored in MeOH at –20 °C until use. Larvae were rehydrated using a gradation of MeOH and PBST washes, treated with proteinase K for 25 min and post-fixed with 4% PFA for 20 min at room temperature. For the detection phase, larvae were pre-hybridized with a probe hybridization buffer for 30 min at 37 °C, then incubated with probes overnight at 37 °C. Larvae were washed with 5 X SSCT to remove excess probes. For the amplification stage, larvae were pre-incubated with an amplification buffer for 30 min at room temperature and incubated with hairpins overnight in the dark at room temperature. Excess hair pins were removed by washing with 5 X SSCT. Larvae were treated with DAPI and stored at 4 °C until imaging. All HCR in situ patterns were confirmed in at least three independent animals. Transcript sequences submitted to Molecular Instruments for probe generation are listed in Supplementary file 12. The cabp1b probes were tested on 3 separate occasions and imaged in at least 6 animals; cabp2b probes were tested on 5 separate occasions and imaged in at least 20 different animals; cabp5b probes were tested on 3 separate occasions and imaged in at least 9 different animals; lfng probes were tested on 2 separate occasions and imaged in at least 5 different animals; loxhd1b probes were tested on 2 separate occasions and imaged in at least 7 animals; pvalb9 probes were tested on 2 separate occasions and imaged in at least 6 different animals; skor2 probes were tested on 2 separate occasions and imaged in at least 6 different animals; tectb probes were tested on 4 separate occasions and imaged in at least 10 different animals; zpld1a probes were tested on 3 separate occasions and imaged in at least 9 different animals.
RNAScope samples were prepared by fixation in 4% paraformaldehyde either at room temperature for 2 hr or at 4 °C overnight. Adult (28–33 mm) inner ears were dissected and dehydrated in methanol for storage. RNAScope probes were synthesized by Advanced Cell Diagnostics (ACD): Channel 1 probe myo6b (1045111-C1), Channel 2 probe pvalb9 (1174621-C2), and Channel 3 probes cabp1b (1137731-C3) and cabp2b (1137741-C3). Whole inner ear tissues were processed through the RNAScope Fluorescent Multiplex V2 Assay (ACD Cat. No. 323100) according to manufacturer’s protocols with the ACD HybEZ Hybridization oven. cabp1b probe was tested on 4 separate occasions with 6 animals or 12 ears total; cabp2b probe was tested on 4 separate occasions with 7 animals or 14 ears total; pvalb9 probe was tested on 2 separate occasions with 6 animals or 12 ears total. myo6b probe was used with each of the above probes.
Immunofluorescence staining
Request a detailed protocolImmediately following the RNAScope protocol, samples were prepared for immunofluorescence staining using mouse anti-β-Spectrin II antibody (BD Bioscience Cat. No. 612562, RRID: AB_399853). Briefly, RNAScope probed zebrafish ears were rehydrated in PBS for 5 min and rinsed in PBDTx (0.5 g bovine serum albumin, 500 μL DMSO, 250 μL 20% Triton-X in 50 mL PBS, pH = 7.4) for 15 min at room temperature. They were then blocked in 2% normal goat serum (NGS) in PBDTx for 3 hr at room temperature, and incubated with 1:500 dilution of mouse anti-β-Spectrin II antibody in PBDTx containing 2% NGS overnight at 4 °C. After three washes in PBDTx for 20 min each at room temperature, samples were incubated with 1:1000 dilution of Alexa 647 goat-anti-mouse IgG1 secondary antibody (Invitrogen Cat. No. A-21240, RRID: AB_2535809) for 5 hr at room temperature. They were then washed 2 times in PBSTx (250 μL 20% Triton-X in 50 mL PBS) for 5 min each before imaging. Three animals or 6 ears total were subjected to Spectrin detection on 2 separate occasions.
Imaging
Request a detailed protocolConfocal images of whole-mount RNAScope samples were captured on a Zeiss LSM800 microscope (Zeiss, Oberkochen, Germany) using ZEN software. HCR-FISH imaging was performed on a Zeiss LSM880 microscope (Zeiss, Oberkochen, Germany) with Airyscan capability. Whole larvae were mounted between coverslips sealed with high vacuum silicone grease (Dow Corning) to prevent evaporation. Z-stacks were taken through the ear at intervals of 1.23 μm using a 10 X objective or through individual inner ear organs at an interval of 0.32 μm using a 20 X objective. 3D Airyscan processing was performed at standard strength settings using Zen Blue software.
Single-cell preparation and analysis
scRNA-seq library preparation and alignment
Request a detailed protocolFor 14 dpf animals (n=35), heads from converted Sox10:Cre; ubb:LOXP-EGFP-STOP-LOXP-mCherry fish were decapitated at the level of the pectoral fin with eyes and brains removed. For 12 mpf animals (n=6, 27–31 mm), utricle, saccule, and lagena were extracted from converted Sox10:Cre; ubb:LOXP-EGFP-STOP-LOXP-mCherry fish after brains and otolith crystals were removed. Dissected heads and otic sensory patches were then incubated in fresh Ringer’s solution for 5–10 min, followed by mechanical and enzymatic dissociation by pipetting every 5 min in protease solution (0.25% trypsin (Life Technologies, 15090–046), 1 mM EDTA, and 400 mg/mL Collagenase D (Sigma, 11088882001) in PBS) and incubated at 28.5 °C for 20–30 min or until full dissociation. Reaction was stopped by adding 6×stop solution (6 mM CaCl2 and 30% fetal bovine serum (FBS) in PBS). Cells were pelleted (376 × g, 5 min, 4 °C) and resuspended in suspension media (1% FBS, 0.8 mM CaCl2, 50 U/mL penicillin, and 0.05 mg/mL streptomycin (Sigma-Aldrich, St. Louis, MO) in phenol red-free Leibovitz’s L15 medium (Life Technologies)) twice. Final volumes of 500 μL resuspended cells were placed on ice and fluorescence-activated cell sorted (FACS) to isolate live cells that excluded the nuclear stain DAPI. For scRNAseq library construction, barcoded single-cell cDNA libraries were synthesized using 10 X Genomics Chromium Single Cell 3′ Library and Gel Bead Kit v.3.1 (14 dpf) or Single Cell Multiome ATAC +Gene Expression kit (12 mpf, single library built with all three sensory patches combined prior to library preparation, ATAC data not shown) per the manufacturer’s instructions. Libraries were sequenced on Illumina NextSeq or HiSeq machines at a depth of at least 1,000,000 reads per cell for each library. Read2 was extended from 98 cycles, per the manufacturer’s instructions, to 126 cycles for higher coverage. Cellranger v6.0.0 (10X Genomics) was used for alignment against GRCz11 (built with GRCz11.fa and GRCz11.104.gtf) and gene-by-cell count matrices were generated with default parameters.
Data processing of scRNA-seq
Request a detailed protocolCount matrices of inner ear and lateral line cells from embryonic and larval timepoints (18–96 hpf) were analyzed using the R package Monocle3 (v1.0.0) (Cao et al., 2019). Matrices were processed using the standard Monocle3 workflow (preprocess_cds, detect_genes, estimate_size_factors, reduce_dimension(umap.min_dist = 0.2, umap.n_neighbors = 25 L)). This cell data set was converted to a Seurat object for integration with 10 X Chromium sequencing data using SeuratWrappers. The count matrices of scRNA-seq data (14 dpf and 12 mpf) were analyzed by R package Seurat (v4.1.0) (Hao et al., 2021). Cells of neural crest origins were removed bioinformatically based on our previous study (Fabian et al., 2022). The matrices were normalized (NormalizeData) and integrated with normalized scRNA-seq data from the embryonic and larval time points according to package instruction (FindVariableFeatures, SelectIntegrationFeatures, FindIntegrationAnchors, IntegrateData; features = 3000). The integrated matrices were then scaled (ScaleData) and dimensionally reduced to 30 principal components. The data were then subjected to neighbor finding (FindNeighbors, k = 20) and clustering (FindClusters, resolution = 0.5), and then visualized through UMAP with 30 principal components as input. After data integration and processing, RNA raw counts from all matrices were normalized and scaled according to package instructions to determine gene expression for all sequenced genes, as the integrated dataset only contained selected features for data integration.
Mouse utricle scRNA-seq data (Jan et al., 2021) was downloaded from NCBI Gene Expression Omnibus (GSE155966). The count matrix was analyzed by R package Seurat (v4.1.0). Matrices were normalized (NormalizeData) and scaled for the top 2000 variable genes (FindVariableFeatures and ScaleData). The scaled matrices were dimensionally reduced to 15 principal components. The data were then subjected to neighbor finding (FindNeighbors, k = 20) and clustering (FindClusters, resolution = 1) and visualized through UMAP with 15 principal components as input. Hair cells and supporting cells were bioinformatically selected based on expression of hair cells and supporting cell markers Myo6 and Lfng, respectively. Hair cells were further subcategorized into striola type I hair cells by co-expression of striola marker Ocm and type I marker Spp, extrastriola type I hair cells by expression of Spp without Ocm, and extrastriola type II hair cells by expression of Anxa4 without Ocm.
Mouse crista scRNA-seq data (Wilkerson et al., 2021) was downloaded from NCBI Gene Expression Omnibus (GSE168901). The count matrix was analyzed by R package Seurat (v4.1.0). Matrices were normalized (NormalizeData) and scaled for the top 2000 variable genes (FindVariableFeatures and ScaleData). The scaled matrices were dimensionally reduced to 15 principal components. The data were then subjected to neighbor finding (FindNeighbors, k = 20) and clustering (FindClusters, resolution = 1) and visualized through UMAP with 15 principal components as input. Hair cells and supporting cells were bioinformatically selected based on expression of hair cell and supporting cell markers Pou4f3 and Sparcl1, respectively. Hair cells were further subcategorized into central hair cells by expression of Ocm and peripheral hair cells by expression of Anxa4.
Pseudotime analysis
Request a detailed protocolWe used the R package Monocle3 (v1.0.1) to predict the pseudo temporal relationships within the integrated scRNA-seq dataset of sensory patches from 36 hpf to 12 mpf. Cell paths were predicted by the learn_graph function of Monocle3. We set the origin of the cell paths based on the enriched distribution of 36–48 hpf cells. Hair (all macular hair cells, clusters 0–5) and supporting (macular supporting cells clusters 0 and 6) cell paths were selected separately (choose_cells) to plot hair cells and supporting cell marker expression along pseudotime (plot_genes_in_pseudotime).
Differential gene expression
Request a detailed protocolWe utilized presto package’s differential gene expression function to identify differentially expressed genes among the different cell types. Wilcox rank sum test was performed by the function wilcox usc. We then filtered for genes with log2 fold change greater than 0.5 and adjusted p-value less than 0.01. To compare inner ear hair cells to lateral line hair cells, we used the following datasets from GEO: 6–7 dpf lateral line hair cells (GSE144827, Kozak et al., 2020), 4 dpf lateral line hair cells (GSE152859, Ohta et al., 2020), and 5 dpf lateral line hair cells and supporting cells (GSE196211, Baek et al., 2022). Hair cells were selected from datasets by expression of otofb and integrated along with our 10 x Chromium dataset with Scanorama (Hie et al., 2019). Gene modules were computed in Monocle3 (v1.0.1) with a q-value cutoff of 1 x e-50.
SAMap analysis for cell type homology
Request a detailed protocolWe used the python package SAMap (v1.0.2) (Tarashansky et al., 2021) to correlate gene expression patterns and determine cell type homology between mouse utricle (GSE155966) (Jan et al., 2021) or crista (GSE168901) (Wilkerson et al., 2021) hair cells and supporting cells and our 12 mpf zebrafish inner ear scRNA-seq data. Zebrafish lateral line hair cell sc-RNA data (GSE123241) (Lush et al., 2019) was integrated with our 12 mpf inner ear data using Seurat in order to compare to mice. First, a reciprocal BLAST result of the mouse and zebrafish proteomes was obtained by performing blastp (protein-protein BLAST, NCBI) in both directions using in-frame translated peptide sequences of zebrafish and mouse transcriptome, available from Ensembl (Danio_rerio.GRCz11.pep.all.fa and Mus_musculus.GRCm38.pep.all.fa). The generated maps were then used for the SAMap algorithm. Raw count matrices of zebrafish and mouse scRNA-seq Seurat objects with annotated cell types were converted to h5ad format using SeuratDisk package (v0.0.0.9020) and loaded into Python 3.8.3. Raw data were then processed and integrated by SAMap. Mapping scores between cell types of different species were then calculated by get_mapping_scores and visualized by sankey_plot. Gene pairs driving cell type homology were identified by GenePairFinder.
Single-cell RNA seq datasets are available from the NCBI Gene Expression Omnibus with Gene Set Accession number GSE211728.
Data availability
Sequencing data have been deposited in GEO under accession code GSE211728.
-
NCBI Gene Expression OmnibusID GSE211728. Single-Cell Transcriptomic Profiling of the Zebrafish Inner Ear Reveals Molecularly Distinct Hair and Supporting Cell Subtypes.
-
NCBI Gene Expression OmnibusID GSE144827. Epithelial planar bipolarity emerges from Notch-mediated asymmetric inhibition of Emx2.
-
NCBI Gene Expression OmnibusID GSE152859. Emx2 defines bidirectional polarity of neuromasts by changing hair-bundle orientation and not hair-cell positions.
-
NCBI Gene Expression OmnibusID GSE196211. High-resolution single cell transcriptome analysis of zebrafish sensory hair cell regeneration.
-
NCBI Gene Expression OmnibusID GSE155966. Single cell RNA seq analysis of postnatal mouse utricle.
-
NCBI Gene Expression OmnibusID GSE168901. Single-cell Transcriptomic Analysis of the Mouse Crista Ampullaris.
-
NCBI Gene Expression OmnibusID GSE123241. Single cell RNA-Seq reveals Fgf signaling dynamics during sensory hair cell regeneration.
-
NCBI Gene Expression OmnibusID GSE178969. Single-cell profiling of cranial neural crest diversification across a vertebrate lifetime.
References
-
Auditory sensitivity of larval zebrafish (Danio rerio) measured using a behavioral prepulse inhibition assayThe Journal of Experimental Biology 216:3504–3513.https://doi.org/10.1242/jeb.087635
-
Heterogeneity of sensory hair cells in a fish earThe Journal of Comparative Neurology 324:621–640.https://doi.org/10.1002/cne.903240413
-
Mapping a multiplexed zoo of mrna expressionDevelopment 143:3632–3637.https://doi.org/10.1242/dev.140137
-
Hair cell heterogeneity and ultrasonic hearing: recent advances in understanding fish hearingPhilosophical Transactions of the Royal Society of London. Series B 355:1277–1280.https://doi.org/10.1098/rstb.2000.0683
-
Ca2+-Binding proteins tune ca2+-feedback to Cav1.3 channels in mouse auditory hair cellsThe Journal of Physiology 585:791–803.https://doi.org/10.1113/jphysiol.2007.142307
-
The new vestibular stimuli: sound and vibration-anatomical, physiological and clinical evidenceExperimental Brain Research 235:957–972.https://doi.org/10.1007/s00221-017-4874-y
-
Comparative morphology of rodent vestibular peripheryI. Saccular and Utricular Maculae. J. Neurophysiol 93:251–266.https://doi.org/10.1152/jn.00746.2003
-
Cloning and expression of the large zebrafish protocadherin gene, fatGene Expression Patterns 5:483–490.https://doi.org/10.1016/j.modgep.2004.12.005
-
A zebrafish model for waardenburg syndrome type IV reveals diverse roles for SOX10 in the otic vesicleDisease Models & Mechanisms 2:68–83.https://doi.org/10.1242/dmm.001164
-
Vestibular hair cells and afferents: two channels for head motion signalsAnnual Review of Neuroscience 34:501–534.https://doi.org/10.1146/annurev-neuro-061010-113710
-
Mariner is defective in myosin VIIA: a zebrafish model for human hereditary deafnessHuman Molecular Genetics 9:2189–2196.https://doi.org/10.1093/hmg/9.14.2189
-
Sound generation in zebrafish with bio-opto-acousticsNature Communications 11:6120.https://doi.org/10.1038/s41467-020-19982-5
-
Pivotal role of hmx2 and hmx3 in zebrafish inner ear and lateral line developmentDevelopmental Biology 339:507–518.https://doi.org/10.1016/j.ydbio.2009.12.028
-
Hair cell replacement in adult mouse utricles after targeted ablation of hair cells with diphtheria toxinThe Journal of Neuroscience 32:15093–15105.https://doi.org/10.1523/JNEUROSCI.1709-12.2012
-
Fluorescent aminoglycosides reveal intracellular trafficking routes in mechanosensory hair cellsThe Journal of Clinical Investigation 127:472–486.https://doi.org/10.1172/JCI85052
-
Efficient integration of heterogeneous single-cell transcriptomes using scanoramaNature Biotechnology 37:685–691.https://doi.org/10.1038/s41587-019-0113-3
-
Oncomodulin expression reveals new insights into the cellular organization of the murine utricle striolaJournal of the Association for Research in Otolaryngology 19:33–51.https://doi.org/10.1007/s10162-017-0652-6
-
Dominant-negative inhibition of M-like potassium conductances in hair cells of the mouse inner earThe Journal of Neuroscience 27:8940–8951.https://doi.org/10.1523/JNEUROSCI.2085-07.2007
-
Aminoglycoside-induced cochleotoxicity: a reviewFrontiers in Cellular Neuroscience 11:308.https://doi.org/10.3389/fncel.2017.00308
-
Vestibular and auditory hair cell regeneration following targeted ablation of hair cells with diphtheria toxin in zebrafishFrontiers in Cellular Neuroscience 15:721950.https://doi.org/10.3389/fncel.2021.721950
-
In vivo and in vitro analysis of starmaker activity in zebrafish otolith biomineralizationFASEB J. Off. Publ. Fed. Am. Soc. Exp. Biol 33:6877–6886.https://doi.org/10.1096/fj.201802268R
-
Novel afferent terminal structure in the crista ampullaris of the goldfish, carassius auratusThe Journal of Comparative Neurology 366:572–579.https://doi.org/10.1002/(SICI)1096-9861(19960318)366:4<572::AID-CNE2>3.0.CO;2-1
-
Architecture of the mouse utricle: macular organization and hair bundle heightsJournal of Neurophysiology 99:718–733.https://doi.org/10.1152/jn.00831.2007
-
Organization of the gravity-sensing system in zebrafishNature Communications 13:5060.https://doi.org/10.1038/s41467-022-32824-w
-
BookMorphophysiology of the vestibular peripheryIn: Lysakowski A, editors. The Vestibular System. Springer. pp. 57–152.https://doi.org/10.1007/0-387-21567-0_3
-
Sensory hair cell death and regeneration in fishesFrontiers in Cellular Neuroscience 9:131.https://doi.org/10.3389/fncel.2015.00131
-
Differences between stereocilia numbers on type I and type II vestibular hair cellsJournal of Neurophysiology 92:3153–3160.https://doi.org/10.1152/jn.00428.2004
-
Vesicular glutamate transporter 3 is required for synaptic transmission in zebrafish hair cellsThe Journal of Neuroscience 28:2110–2118.https://doi.org/10.1523/JNEUROSCI.5230-07.2008
-
Otoc1: a novel otoconin-90 ortholog required for otolith mineralization in zebrafishDevelopmental Neurobiology 68:209–222.https://doi.org/10.1002/dneu.20587
-
Harmonin (USH1C) is required in zebrafish Müller glial cells for photoreceptor synaptic development and functionDisease Models & Mechanisms 4:786–800.https://doi.org/10.1242/dmm.006429
-
Water waves to sound waves: using zebrafish to explore hair cell biologyJournal of the Association for Research in Otolaryngology 20:1–19.https://doi.org/10.1007/s10162-018-00711-1
-
Sound detection and processing by fish: critical review and major research questionsBrain, Behavior and Evolution 41:14–38.https://doi.org/10.1159/000113821
-
The zebrafish genome in context: ohnologs gone missingJournal of Experimental Zoology. Part B, Molecular and Developmental Evolution 308:563–577.https://doi.org/10.1002/jez.b.21137
-
Single-cell RNA-sequencing of zebrafish hair cells reveals novel genes potentially involved in hearing lossCellular and Molecular Life Sciences 79:385.https://doi.org/10.1007/s00018-022-04410-2
-
Neuronal expression of fibroblast growth factor receptors in zebrafishGene Expression Patterns 13:354–361.https://doi.org/10.1016/j.gep.2013.06.006
-
The zebrafish Eya1 gene and its expression pattern during embryogenesisDevelopment Genes and Evolution 209:399–410.https://doi.org/10.1007/s004270050270
-
A mutation in CaBP2, expressed in cochlear hair cells, causes autosomal-recessive hearing impairmentAmerican Journal of Human Genetics 91:636–645.https://doi.org/10.1016/j.ajhg.2012.08.018
-
Presynaptic cav1.3 channels regulate synaptic ribbon size and are required for synaptic maintenance in sensory hair cellsThe Journal of Neuroscience 32:17273–17286.https://doi.org/10.1523/JNEUROSCI.3005-12.2012
-
How zebrafish can drive the future of genetic-based hearing and balance researchJournal of the Association for Research in Otolaryngology 22:215–235.https://doi.org/10.1007/s10162-021-00798-z
-
Hair cell regeneration in the avian auditory epitheliumThe International Journal of Developmental Biology 51:633–647.https://doi.org/10.1387/ijdb.072408js
-
Hair bundle heights in the utricle: differences between macular locations and hair cell typesJournal of Neurophysiology 95:171–186.https://doi.org/10.1152/jn.00800.2005
-
Supporting cell proliferation after hair cell injury in mature guinea pig cochlea in vivoCell and Tissue Research 325:23–31.https://doi.org/10.1007/s00441-006-0157-9
-
The role of gfi1.2 in the development of zebrafish inner earHearing Research 396:108055.https://doi.org/10.1016/j.heares.2020.108055
-
Development of the acoustically evoked behavioral response in zebrafish to pure tonesThe Journal of Experimental Biology 208:1363–1372.https://doi.org/10.1242/jeb.01534
Article and author information
Author details
Funding
National Institute on Deafness and Other Communication Disorders (R21DC019948)
- David W Raible
National Institute on Deafness and Other Communication Disorders (F31DC020898)
- Marielle O Beaulieu
Hamilton and Mildred Kellogg Trust
- David W Raible
The Whitcraft Family Gift
- David W Raible
Hearing Health Foundation
- David W Raible
Paul G. Allen Frontiers Group (Allen Discovery Center for Cell Lineage Tracing)
- Cole Trapnell
National Human Genome Research Institute (UM1HG011586)
- Cole Trapnell
National Human Genome Research Institute (1R01HG010632)
- Cole Trapnell
National Institute on Deafness and Other Communication Disorders (F31DC020633)
- Tuo Shi
National Institute of Dental and Craniofacial Research (R35DE027550)
- J Gage Crump
National Institute on Deafness and Other Communication Disorders (R01DC015829)
- Neil Segil
National Institute on Deafness and Other Communication Disorders (T32DC009975)
- Tuo Shi
- Neil Segil
National Institute on Deafness and Other Communication Disorders (T32DC005361)
- Marielle O Beaulieu
- David W Raible
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
This manuscript is dedicated to Neil Segil, who was a wonderful colleague, friend, and mentor to many. We thank Megan Matsutani for fish care, the USC Stem Cell Flow Cytometry Core, and the CHLA Next-Generation Sequencing Core. We also thank David White and the UW Zebrafish Facility staff for fish care.
Ethics
This study was performed in strict accordance with the recommendations in the Guide for the Care and Use of Laboratory Animals of the National Institutes of Health. The Institutional Animal Care and Use Committees of the University of Southern California (Protocol 20771) and University of Washington (Protocol 2997-01) approved all animal experiments.
Copyright
© 2023, Shi, Beaulieu 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
-
- 3,920
- views
-
- 528
- downloads
-
- 80
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Developmental Biology
Premature infants with bronchopulmonary dysplasia (BPD) have impaired alveolar gas exchange due to alveolar simplification and dysmorphic pulmonary vasculature. Advances in clinical care have improved survival for infants with BPD, but the overall incidence of BPD remains unchanged because we lack specific therapies to prevent this disease. Recent work has suggested a role for increased transforming growth factor-beta (TGFβ) signaling and myofibroblast populations in BPD pathogenesis, but the functional significance of each remains unclear. Here, we utilize multiple murine models of alveolar simplification and comparative single-cell RNA sequencing to identify shared mechanisms that could contribute to BPD pathogenesis. Single-cell RNA sequencing reveals a profound loss of myofibroblasts in two models of BPD and identifies gene expression signatures of increased TGFβ signaling, cell cycle arrest, and impaired proliferation in myofibroblasts. Using pharmacologic and genetic approaches, we find no evidence that increased TGFβ signaling in the lung mesenchyme contributes to alveolar simplification. In contrast, this is likely a failed compensatory response, since none of our approaches to inhibit TGFβ signaling protect mice from alveolar simplification due to hyperoxia while several make simplification worse. In contrast, we find that impaired myofibroblast proliferation is a central feature in several murine models of BPD, and we show that inhibiting myofibroblast proliferation is sufficient to cause pathologic alveolar simplification. Our results underscore the importance of impaired myofibroblast proliferation as a central feature of alveolar simplification and suggest that efforts to reverse this process could have therapeutic value in BPD.
-
- Developmental Biology
Cells called alveolar myofibroblasts, which have a central role in the development of the lung after birth, receive an orchestrated input from a range of different signaling pathways.