Single-cell transcriptomic profiling of the zebrafish inner ear reveals molecularly distinct hair cell and supporting cell subtypes

  1. Tuo Shi
  2. Marielle O Beaulieu
  3. Lauren M Saunders
  4. Peter Fabian
  5. Cole Trapnell
  6. Neil Segil
  7. J Gage Crump  Is a corresponding author
  8. David W Raible  Is a corresponding author
  1. Department of Stem Cell Biology and Regenerative Medicine, Keck School of Medicine, University of Southern California, United States
  2. Caruso Department of Otolaryngology-Head and Neck Surgery, Keck School of Medicine, University of Southern California, United States
  3. Department of Otolaryngology-Head and Neck Surgery, University of Washington, United States
  4. Department of Genome Sciences, University of Washington, United States
  5. Department of Biological Structure, University of Washington, United States
9 figures and 13 additional files

Figures

Anatomy of zebrafish and mouse inner ears.

(A) Illustration of the lateral line system of a 5 dpf zebrafish. Blue circles represent individual neuromasts located on the body of the fish. Boxed region indicates location of the ear. (B) Enlarged diagram of the 5 dpf zebrafish ear showing cristae (red) and macular (blue) sensory organs. (C,D) Illustrations of adult zebrafish and mouse inner ears showing homologous end organs in the semicircular canal crista ampullaris (red) and macula otolith organs (blue). Light green and dark green represent unique end organs of the lagena in zebrafish and cochlea in mice. (E) Illustration of the mouse utricle showing striolar and extrastriolar regions of the sensory organ. Arrows represent hair cell planar polarity within the sensory organ and red dashed line represents the line of polarity reversal within the striola. ac: anterior crista, c: cochlea, l: lagena, lc: lateral crista, o: otolith, pc: posterior crista, s: saccule, u: utricle.

Figure 2 with 3 supplements
Molecularly distinct cell types between the zebrafish inner ear and lateral line.

Ear and lateral line cells were selected from a whole-embryo single-nucleus RNA-seq dataset from animals between 18 and 96 hpf using known marker genes for hair cells and supporting cells. (A–B) UMAP projection of inner ear and lateral line cells grouped by (A) developmental timepoint and (B) broad cell type: ear nonsensory SC (red), lateral line nonsensory SC (green), ear HC (blue), and lateral line HC (yellow). Clusters in (B) correspond to columns of following gene expression plots. Widely accepted marker genes for (C) inner ear nonsensory cells, (D) lateral line nonsensory cells, and (E) hair cells show enriched expression in the corresponding clusters from B, confirming their identity. (F) Expression of previously identified marker genes for inner ear or lateral line hair cells was used to identify hair cell origin.

Figure 2—figure supplement 1
Gene modules for embryonic to larval inner ear and lateral line dataset.

Gene modules calculated in Monocle 3 for the embryonic to larval inner ear and lateral line dataset displayed as (A) a heatmap of module gene enrichment by cluster where red indicates higher enrichment and blue indicates de-enrichment and (B) module expression across the UMAP for the dataset. Module genes with statistical values are listed in Supplementary file 2.

Figure 2—figure supplement 2
Selection of otic sensory cells from snRNA-seq dataset.

(A) Clustering of 18 hpf to 96 hpf dataset to illustrate cell subtypes. PAGA analysis of this dataset shows strong connectivity among ear nonsensory cells and among lateral line nonsensory cells, but weak interconnectivity between these two groups. (B) Feature plots show expression of the supporting cell marker lfng, and markers of structural otic cells matn4 and col2a1a. (C) UMAP of sensory patch cells from 36 to 96hpf are clustered without structural and early otic vesicle cells. PAGA analysis again shows strong connectivity within hair cells and supporting cell groups and weak connectivity between lateral line and inner ear supporting cells. PAGA connectivity scores are listed in Supplementary file 1.

Figure 2—figure supplement 3
Gene expression differences between lateral line and inner ear hair cells.

(A) UMAP of our 12 mpf hair cell dataset integrated by Scanorama with published lateral line hair cell datasets. Lateral line hair cells cluster separately from inner ear hair cells. (B) Differential gene expression analysis identifies novel marker genes specific to either lateral line or inner ear hair cells.

Figure 3 with 4 supplements
Cell subtypes in the zebrafish inner ear end organs.

(A–D) Integration and analysis of single-cell RNAseq data generated by sci-Seq (sci) or 10x Chromium sequencing (10x) for inner ear hair cells and supporting cells from embryonic (sci), larval (sci,10x), and adult (10x) stages. UMAP projection of cells are grouped by (A) dataset of origin and (B) timepoint. (C) Unsupervised clustering divides cells into 10 clusters that were grouped into 9 cell subtypes. (D) Feature plots showing hair cell marker myo6b, nascent hair cell marker dla, supporting cell marker lfng, and putative progenitor marker fgfr2 expression in the integrated dataset. (E) Differentially expressed genes across the 10 cell clusters.

Figure 3—figure supplement 1
scRNA-seq of 12 mpf zebrafish inner ear captures sensory hair cells and supporting cells as well as non-sensory supporting cells.

(A) Clustering of 12 mpf dataset to illustrate cell types in the adult zebrafish inner ear. (B–I) Feature plots of 12 mpf zebrafish scRNA-seq dataset alone showing expression of hair cell markers (B) myo6b and (C) strc, pan-supporting cell marker (D) stm, sensory supporting cell markers (E) lfng and (F) hey1, and pan-otic marker (G) otomp, and non-sensory supporting cell markers (H) matn4 and (I) col2a1a.

Figure 3—figure supplement 2
Hair cell and supporting cell marker expression in the integrated scRNA-seq dataset.

Feature plots of integrated zebrafish scRNA-seq datasets showing expression of nascent hair cell marker (A) atoh1a, inner ear hair cell markers (B) strc, (C) gfi1aa, and (D) lhfpl5a, and pan-supporting cell marker (E) stm.

Figure 3—figure supplement 3
Putative progenitor marker expression in individual progenitor and supporting cell clusters.

(A) Combined and individual UMAP projections of putative bipotent progenitor cluster (cluster 0), macular supporting cell cluster (cluster 6), and cristae supporting cell cluster (cluster 7) from the integrated zebrafish inner ear dataset. (B–E) Feature plots show expression of putative progenitor genes in the integrated dataset, as well as in individual clusters of 0, 6, and 7. (F) Violin plots showing differential gene expression of fgfr2, igsf3, fat1a, and pard3bb among clusters 0, 6, and 7. Wilcoxon rank sum test, *: p ≤ 0.05, ***: p ≤ 1e-3, ****: p ≤ 1e-4.

Figure 3—figure supplement 4
Gene modules for integrated inner ear sensory patch dataset.

Gene modules calculated in Monocle 3 for the integrated inner ear sensory patch dataset displayed as (A) a heatmap of module gene enrichment by cluster where red indicates higher enrichment and blue indicates de-enrichment and (B) module expression across the UMAP for the dataset. Module genes with statistical values are listed in Supplementary file 5.

Figure 4 with 2 supplements
Pseudotime analysis reveals developmental trajectories in the zebrafish inner ear.

(A,B) Pseudotime analysis of macular cells showing simulated developmental trajectories of a putative bipotent progenitor population into hair cell and supporting cell clusters. (C,D) Changes in putative progenitor markers along (C) hair cell and (D) supporting cell trajectories. fat1a and pard3bb only decrease along the hair cell trajectory, while fgfr2 and igsf3 decrease along both hair cell and supporting cell trajectories. (E) Transient expression of early hair cell genes dla and atoh1a along hair cell trajectories. (F) Increases in gene expression levels of gfi1aa and myo6b along hair cell trajectories. (G) Increases in stm and lfng along supporting cell trajectories.

Figure 4—figure supplement 1
Pseudotime analysis of cristae hair and supporting cells in the zebrafish inner ear.

(A) Pseudotime analysis showing simulated developmental trajectories of a putative bipotent progenitor population into both cristae hair and supporting cell clusters. (B,C) Changes in putative progenitor markers along (B) hair cell and (C) supporting cell trajectories. fat1a and pard3bb only decrease along the hair cell trajectory, while fgfr2 and igsf3 decrease along both hair cell and supporting cell trajectories. (D) Transient expression of early hair cell genes dla and atoh1a along hair cell trajectories. (E) Increases in gene expression levels of gfi1aa and myo6b along hair cell trajectories. (F) Increases in stm and lfng along supporting cell trajectories.

Figure 4—figure supplement 2
dla labels putative hair cell progenitors in the cristae and maculae.

HCR in situ hybridization of 5 dpf zebrafish. Maximum intensity projections of (A) posterior crista (lateral view), (B) utricle (dorsal view), and (C) saccule (lateral view) showing dla expression in a subset of support cells (arrowheads) peripheral to myo6b+hair cells. Scale bars = 10 μm.

Figure 5 with 1 supplement
Distinct markers separate macula and crista supporting cells.

(A) Feature plots showing expression of macula supporting cell marker tectb and crista supporting cell marker zpld1a. (B–D) HCR in situ hybridization in myo6b:GFP transgenic animals. Each set of images shown represents a projection of one z-stack split into cristae (lateral) and macula (medial) slices. Lateral line neuromasts positioned over the ear are visible in lateral slices. Expression pattern for (B) the pan-supporting cell marker lfng, (C) macula-specific marker tectb, and (D) crista-specific marker zpld1a in 5 dpf myo6b:GFP fish. Each set of images shown represents a projection of one z-stack split into cristae (lateral) and macula (medial) slices. ac: anterior crista, lc: lateral crista, nm: neuromast, pc: posterior crista, u: utricle, s: saccule. Scale bars = 20 μm.

Figure 5—figure supplement 1
zpld1a and tectb are primarily expressed in supporting cells.

HCR in situ hybridization of 5 dpf myo6b:GFP zebrafish. (A–B) Confocal slices through (A) anterior crista and (B) lateral crista (lateral view) show localization of cabp5b in hair cells and zpld1a in supporting cells. (C) Slice through utricle (dorsal view) shows cabp2b expression in hair cells and tectb expression primarily in the surrounding supporting cells. (D) Slices through saccule (lateral view) at the level of hair cell bodies (top row) and supporting cell bodies (bottom row). cabp2b is primarily expressed in hair cells and tectb is primarily expressed in supporting cells. Scale bars = 10 μm.

cabp1b+and cabp2b+label hair cells in distinct regions of sensory end organs.

(A) Feature plots showing differential expression of cabp1b and cabp2b among crista and macula hair cells. (B–D) HCR in situ projections of individual sensory patches from 5 dpf myo6b:GFP fish showing differential spatial expression patterns of cabp1b and cabp2b. (B) cabp1b is expressed at the ends of the cristae, while cabp2b is expressed centrally. Anterior crista is shown. (C) In the utricle, cabp1b is expressed medially and cabp2b is expressed laterally. (D) In the saccule, cabp1b is expressed in peripheral cells at the dorsal and ventral edges of the organ. cabp2b is expressed centrally. Scale bars for HCR images = 10 μm. (E) Cartoon illustrations of the zebrafish utricle, saccule, and lagena, and the expression patterns of cabp1b (yellow) and cabp2b (magenta) within each sensory patch. (F–H) Wholemount RNAScope confocal images of adult inner ear organs showing peripheral expression pattern of cabp1b (n=3) in the adult zebrafish (F) utricle, (G) saccule, and (H) lagena. (I–K) Whole-mount RNAScope confocal images showing central expression pattern of cabp2b (n=4) in the adult zebrafish (I) utricle, (J) saccule, and (K) lagena. Scale bars for RNAScope images = 25 μm.

Figure 7 with 1 supplement
Distinct markers separate macula and crista hair cells.

(A) Feature plots showing marker genes enriched in organ-specific subsets of inner ear hair cells: cabp5b, skor2, and loxhd1b. (B–D) HCR in situs in 5 dpf myo6b:GFP fish show expression of (B) cabp5b in crista but not macula hair cells, (C) skor2 in the utricle only, and (D) loxhd1b in the saccule, as well as lateral line neuromast hair cells. Each set of images represents an orthogonal projection of one z-stack split into cristae (lateral) and macular (medial) slices. ac: anterior crista, lc: lateral crista, nm: neuromast, pc: posterior crista, s: saccule, u: utricle. Scale bar = 20 μm.

Figure 7—figure supplement 1
skor2 and loxhd1b label subsets of hair cells in utricle or saccule.

HCR in situ hybridization of 5 dpf zebrafish. (A) Maximum intensity projection of utricle (dorsal view) showing skor2 expression in medially located myo6b:GFP+ hair cells. (B) Maximum intensity projection of saccule (lateral view) showing loxhd1b expression in a peripheral subset of hair cells. Scale bars = 10 μm.

Figure 8 with 3 supplements
Zebrafish cabp2b+domain shares features with the mouse striolar region.

(A) Feature plot shows enrichment for the striola marker pvalb9 in cabp2b-expressing striolar cells. (B–D) HCR in situs in 5 dpf myo6b:GFP fish shows pvalb9 and cabp2b co-expression in the utricle. Scale bar = 10 μm. (E) Cartoon illustration of overlapping expression of pvalb9 (white) and cabp2b (magenta) that coincides with the line of hair cell polarity reversal. (F, G) Whole-mount RNAScope confocal images of adult zebrafish utricles showing expression of pvalb9 relative to (F) cabp1b (n=3) and (G) cabp2b (n=4). Scale bar = 25 μm. (H,I) Whole-mount RNAScope RNA and protein co-detection assay showing co-localization of cabp2b expression (RNA) and the hair cell line of polarity reversal indicated by Spectrin (protein) staining (n=3). Scale bar = 25 μm. Arrows denote hair cell polarity and dotted line outlines line of polarity reversal.

Figure 8—figure supplement 1
Striola marker pvalb9 is expressed in all inner ear sensory end organs.

HCR in situ hybridization of 5 dpf zebrafish. (A) Maximum intensity projection of saccule (lateral view) shows pvalb9 expression in centrally located myo6b:GFP+ hair cells. (B) Slice through the anterior crista shows pvalb9 expression in a subset of crista hair cells. Scale bars = 10 μm.

Figure 8—figure supplement 2
Inner ear hair cell subtypes differentially express mechanosensory apparatus genes.

Feature plots for mechanosensory transduction genes from the integrated zebrafish scRNA-seq dataset of Figure 3.

Figure 8—figure supplement 3
Inner ear hair cell subtypes differentially express voltage-gated calcium and potassium channel genes.

Feature plots for ion channel genes from the integrated zebrafish scRNA-seq dataset of Figure 3.

Figure 9 with 1 supplement
SAMap analysis reveals conserved gene expression patterns between mouse and zebrafish hair cell types.

(A–B) Sankey plot showing the SAMap mapping scores (0–1) that indicate transcriptome relatedness between (A) mouse utricular and zebrafish macular single-cell clusters and (B) mouse and zebrafish cristae single-cell clusters. A mapping score of 0 indicates no evolutionary correlation in transcriptome while a mapping score of 1 indicates perfect correlation. Correlations below 0.15 were not plotted.

Figure 9—source data 1

Mapping Scores between mouse utricle and zebrafish maculae hair and supporting cells.

https://cdn.elifesciences.org/articles/82978/elife-82978-fig9-data1-v2.zip
Figure 9—source data 2

Mapping Scores between mouse cristae and zebrafish cristae hair and supporting cells.

https://cdn.elifesciences.org/articles/82978/elife-82978-fig9-data2-v2.zip
Figure 9—figure supplement 1
SAMap analysis of mouse utricle versus zebrafish macular and lateral line cells.

(A–B) Sankey plot showing the SAMap mapping scores (0–1) that indicate transcriptome relatedness between mouse utricular and integrated zebrafish macular and lateral line single-cell clusters. (A) Zebrafish 12 mpf macular HCs integrated with 5 dpf lateral line HCs. (B) Zebrafish 3–5 dpf macular HCs integrated with 5 dpf lateral line HCs. A mapping score of 0 indicates no evolutionary correlation in transcriptome, and a mapping score of 1 indicates perfect correlation. Correlations below 0.2 were not plotted.

Figure 9—figure supplement 1—source data 1

Mapping Scores between mouse utricle and zebrafish maculae and lateral line hair cells.

https://cdn.elifesciences.org/articles/82978/elife-82978-fig9-figsupp1-data1-v2.txt

Additional files

Supplementary file 1

Differentially expressed genes across inner ear and lateral line clusters.

https://cdn.elifesciences.org/articles/82978/elife-82978-supp1-v2.csv
Supplementary file 2

PAGA scores for relative connectivity between clusters (related to Figure 2—figure supplement 2).

https://cdn.elifesciences.org/articles/82978/elife-82978-supp2-v2.csv
Supplementary file 3

Gene modules for embryonic to larval inner ear and lateral line dataset.

https://cdn.elifesciences.org/articles/82978/elife-82978-supp3-v2.xlsx
Supplementary file 4

Differentially expressed genes in inner ear cell clusters.

https://cdn.elifesciences.org/articles/82978/elife-82978-supp4-v2.xlsx
Supplementary file 5

Gene modules for inner ear sensory patch dataset.

https://cdn.elifesciences.org/articles/82978/elife-82978-supp5-v2.xlsx
Supplementary file 6

Genes enriched along pseudotime trajectories.

https://cdn.elifesciences.org/articles/82978/elife-82978-supp6-v2.csv
Supplementary file 7

Genes enriched in supporting cell clusters.

https://cdn.elifesciences.org/articles/82978/elife-82978-supp7-v2.xlsx
Supplementary file 8

Genes enriched in cabp1b+and cabp2b+macular cells.

https://cdn.elifesciences.org/articles/82978/elife-82978-supp8-v2.xlsx
Supplementary file 9

Genes driving macular SAMap alignment.

https://cdn.elifesciences.org/articles/82978/elife-82978-supp9-v2.txt
Supplementary file 10

Genes enriched in cabp1b+and cabp2b+crista cells.

https://cdn.elifesciences.org/articles/82978/elife-82978-supp10-v2.xlsx
Supplementary file 11

Genes driving crista SAMap alignment.

https://cdn.elifesciences.org/articles/82978/elife-82978-supp11-v2.txt
Supplementary file 12

cDNA sequences used for HCR in situ hybridization probes.

https://cdn.elifesciences.org/articles/82978/elife-82978-supp12-v2.xlsx
MDAR checklist
https://cdn.elifesciences.org/articles/82978/elife-82978-mdarchecklist1-v2.docx

Download links

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

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Tuo Shi
  2. Marielle O Beaulieu
  3. Lauren M Saunders
  4. Peter Fabian
  5. Cole Trapnell
  6. Neil Segil
  7. J Gage Crump
  8. David W Raible
(2023)
Single-cell transcriptomic profiling of the zebrafish inner ear reveals molecularly distinct hair cell and supporting cell subtypes
eLife 12:e82978.
https://doi.org/10.7554/eLife.82978