Single-cell transcriptional atlas of cochlear and vestibular cells.

A: Schematic drawing of the organ of Corti (top panel) and representative images of IHCs and OHCs from adult mouse cochleae. B: Schematic drawing of the utricle (top panel) and confocal images of the utricle prepared from an adult mouse. HCs are stained with MYO6, SOX2 and DAPI. SOX2-positive cells include type II HCs and SCs underneath HCs. C: Representative images of type I and type II HCs from utricular and saccular maculae as well as crista ampullaris from adult mice. D, E: tSNE plots of distinct cell types detected in the adult CBA mouse cochlea (D) and utricle, saccule and crista. Different cell types are color-coded. F: Feature plots of the expression 6 marker genes in different HC populations. G: Dot plot heatmap of average expression and cellular detection rate of 26 representative marker genes in different HC types in cochleae and vestibular end organs. Abbreviations: IHC (inner HC); OHC (outer HC); SGN (spiral ganglion neuron); SC (Schwann cell)/SGC (satellite glial cell); DC (Deiters’ cell)/PC (pillar cell); IPhC (inner phalangeal cell)/IB (inner border cell); ISC (inner sulcus cell); HeC (Hensen’s cell); SpC (spindle cell)/RC (root cell); MC (marginal cell)/IC (intermediate cell); BaC (basal cell); MP (macrophage); RM (Reissner’s membrane); BC (B cell); TC (T cell); Gran (granulocyte); Mono (monocyte); RBC (red blood cell); OSC (outer sulcus cell); TBC (tympanic border cell). Type I HC (type I HC); Type II HC (type II HC); SESC (sensory epithelial SC); TEC (transitional epithelial cell); PTEC (peripheral transitional epithelial cell); ERC (epithelial roof cell). Neurons (vestibular neurons). For those cells whose definite identities cannot be annotated, the top expressed genes were used for identification and annotation. These cells include Tectb+, Mgp+, Coch+, Prl+, Otos+, Emilin2+/Cavin2+, Fxyd2+/Kcnk2+, Prtn3+/Ccl5+, and Ptgds+/Coch+.

Similarity and difference among different HC types and biological processes enriched in cochlear and vestibular HCs.

A: PCA plot showing similarity based on pseudo bulk RNA-seq data from four HC types. B: PCA plot showing similarity based on individual HC gene expression among the four HC types. C: Venn diagram depicting the number of expressed genes (RPKM > 0) in four HC types. D: The number and percentage of the total genes shared among two or more cell types or those uniquely expressed by a single cell type. E: Volcano plot showing differentially expressed genes between different HC types. Red dots indicate differentially expressed genes with P-value < 10e-5 and log2 fold change > 1. Only the top 20 differentially expressed genes are labeled. F: Biological processes enriched in vestibular HCs compared to cochlear HCs. Biological processes related to motile cilia are highlighted by red asterisks. G: Biological processes enriched in type I HCs compared to type II HCs. H: Biological processes enriched in type II HCs compared to type I HCs.

Shared and unique genes expressed in cochlear and vestibular HCs.

A: Violin plots of the expression 72 genes in four different HC types. B: Validation of differential expression of 9 genes (with underline in A) in cochlear and vestibular HCs in thin section. Bar: 5 μm for all images in B. C: Confocal images of expression of DNM1, SLC7A14, and TJAP1 in cochlear and vestibular HCs. Bar: 10 μm for all images in C.

Expression of genes related to HC specialization.

A: Heatmap showing expression of genes related to stereocilia bundles, mechanotransduction, ion channels and synaptic structure. B: Validation of gene expression using smFISH. Bar: 10 μm.

Cilia-related genes detected in cochlear and vestibular HCs.

A: Schematic illustration of HC hair bundle (adapted from Schwander et al., 2010). B: Venn diagram illustrating the number of genes in each database and the cilia-related genes detected in HC transcriptomes. C: Schematic illustration of primary and motile cilia, highlighting 9+0 or 9+2 arrangement of microtubules for primary and motile cilia, respectively: radial spokes (RS), central pair complex (CPC), nexin–dynein regulatory complex (N-DRC), microtubule inner proteins (MIPs), inner and outer dynein arms (IDA and ODA). This panel was created using BioRender.com. D: Expression of top 50 cilia-related genes and genes related to IFT in the four types of HCs. E: Immunostaining of IFT172 and CLUAP1 expression in vestibular kinocilia. Bar: 5 µm. F: Violin plots showing aggregated expression of genes associated with 96-nm repeat. G: Heatmaps showing the comparison of gene expression related to motile-cilia machinery: radial spokes (RS), central pair complex (CPC), nexin–dynein regulatory complex (N-DRC), microtubule inner proteins (MIPs), inner and outer dynein arms (IDA and ODA), external coiled coils, and motile cilia transcriptional regulators in cochlear and vestibular HCs.

Expression of motile cilia-related genes/proteins in the vestibular HCs.

A: Expression of motile cilia-related genes in mouse, zebrafish, and human vestibular HCs. Mouse gene nomenclature is used in the heatmaps. B: Confocal images of the expression of key motile cilia-related proteins. Scale bars represent 5 µm. C: SEM micrograph of hair bundles of OHCs from P2 cochlea. Kinocilia (in magenta) are still present at this age. D: Comparison of expression of motile cilia-related genes between P2 cochlear and vestibular HCs. E: Confocal images of expression of CCDC39, CCDC40 and DNAH6 in P2 vestibular and cochlear HCs. CCDC39, CCDC40, and DNAH6 were not expressed in cochlear HCs at P2. Bar: 5 μm.

Kinocilia morphology and motility.

A: TEM images of stereocilia and kinocilium from bullfrog crista HCs. Different regions of the kinocilium in higher magnification are also shown. Long black arrows indicate where the magnified images were taken. Bar: 250 nm. Red arrow indicates two central microtubule singlets. Short black arrows mark absence of central microtubule singlets in the distal regions near the tip of kinocilium and transition zone. B: Images captured from in vitro live imaging of kinocilium and bundle motion of a bullfrog crista HC. The images were captured with speed of 15 frames per second. Black arrows indicate kinocilium. C: DIC image of cilia from mouse middle ear. Black frame marks where the motion was measured. Bar: 5 μm. D: Representative waveforms of spontaneous ciliary motion from middle ear tissue. The FFT analysis of cilia motion is also shown. E: DIC image of hair bundles from mouse crista HCs. Bar: 10 μm. F: Three representative waveforms of spontaneous motion of hair bundles. The response in blue shows measurements taken from a hair bundle with no spontaneous motility. Time scale in F also applies to D. FFT analysis of bundle motion is also shown. Response waveforms and spectra are color-coded and -matched.

Predicted models of the molecular architecture of 96 nm axonemal repeat of vestibular kinocilia.

A-B: Longitudinal and cross-sectional views of the doublet microtubule and associated structure in 96-nm repeat, derived from combining cryo-EM data and single-cell transcriptomic analysis from human respiratory (A) and bovine sperm flagella (B). Key axonemal motile-machinery components are color-coded: ODA (Indian red), IDA (cyan), N-DRC (green), MIPs (orchid), RS (purple), and external coiled-coils (blue). Radial spoke 3 (RS3) has not been resolved to atomic resolution, but its shorter form (RS3s) is depicted. Doublet microtubules (DMT) are represented in gray. Regions highlighted in gold indicate the absence of corresponding transcripts in our mouse transcriptomic data. C: Genes which are not detected in mouse and human vestibular HC transcriptomes and related to motility-relevant compartments are listed in the table. The roles of these sgenes in the 96-nm repeat module and cilia motility and ciliopathy are also included.