Single-cell proteomics reveals changes in expression during hair-cell development

  1. Ying Zhu
  2. Mirko Scheibinger
  3. Daniel Christian Ellwanger
  4. Jocelyn F Krey
  5. Dongseok Choi
  6. Ryan T Kelly
  7. Stefan Heller
  8. Peter G Barr-Gillespie  Is a corresponding author
  1. Pacific Northwest National Laboratory, United States
  2. Stanford University, United States
  3. Amgen Inc, United States
  4. Oregon Health & Science University, United States
  5. Kyung Hee University, Republic of Korea
  6. Brigham Young University, United States
5 figures, 4 tables and 2 additional files

Figures

Figure 1 with 2 supplements
Mass spectrometry of single cells and small cell pools from E15 chick utricle.

(A) Experimental design. The E15 chick utricle’s sensory epithelium consists of sensory hair cells and supporting cells, which are also progenitor cells. FM1-43 labels hair cells more strongly than …

https://doi.org/10.7554/eLife.50777.002
Figure 1—source data 1

MaxQuant analysis of single cell proteomics data.

Excel file containing data from Experiments 1 and 2 analyzed together by MaxQuant and Andromeda.

https://doi.org/10.7554/eLife.50777.005
Figure 1—figure supplement 1
Characterization of isolated utricle cells.

(A-C) FACS gating protocol. (A) Debris was removed based on side-scatter area (SSC-A) and forward-scatter area (FSC-A). Blue box shows cells collected. (B) Doublets were discarded based on forward …

https://doi.org/10.7554/eLife.50777.003
Figure 1—figure supplement 2
Peptide coverage of TMSB4X in mass spectrometry experiments.

Characterization of isolated utricle cells. (A–C) FACS gating protocol. (A) Debris was removed based on side-scatter area (SSC-A) and forward-scatter area (FSC-A). Blue box shows cells collected. (B)…

https://doi.org/10.7554/eLife.50777.004
Figure 2 with 3 supplements
Abundant proteins in small pools of isolated E15 chick utricle cells.

(A) Heat map showing top 60 proteins or protein groups in samples of 20 cells, sorted by the average of the 20 cell FM1-43high samples. FM1-43low and FM1-43high samples from Experiment 1 (Exp1) and …

https://doi.org/10.7554/eLife.50777.006
Figure 2—figure supplement 1
Expression levels for all proteins in samples from Experiments 1 and 2 that contain pools of three cells, five cells, or 20 cells.

Each sample (indicated by letters) under the callout for cell numbers contains that number of cells (i.e., on the left, samples a-c under labels indicating ‘Experiment 1’ and ‘FM1-43 low’ each …

https://doi.org/10.7554/eLife.50777.007
Figure 2—figure supplement 2
Expression levels for all proteins in samples from Experiments 1 and 2 from single-cell samples.

Each sample is indicated by letters. ‘Avg’ indicates the average of the 20 cell samples for that condition and experiment. Log10 riBAQ color scale is indicated at bottom. Heat map is sorted by the …

https://doi.org/10.7554/eLife.50777.008
Figure 2—figure supplement 3
Expression levels for all proteins in single-cell samples in Experiments 1 and 2.

Each sample is indicated by letters. ‘Avg’ indicates the average of the 20 cell samples for that condition and experiment. Heat map, which shows all identified proteins or protein groups in single …

https://doi.org/10.7554/eLife.50777.009
Figure 3 with 6 supplements
Immunolocalization of proteins enriched in hair cells or supporting cells of E15 chick utricle.

Confocal z-stacks of vibratome cross-sections of the whole utricle were imaged with the tiling and stitching function in Zeiss ZEN. Confocal z-stacks for magnified extrastriolar and striolar regions …

https://doi.org/10.7554/eLife.50777.010
Figure 3—figure supplement 1
Immunolocalization of AGR3 in E15 chick utricle.

The sensory epithelium and the stroma are indicated in the transmitted-light section. Nuclei were stained with DAPI and F-actin with phalloidin; OTOF and AGR3 were stained with specific antibodies. …

https://doi.org/10.7554/eLife.50777.011
Figure 3—figure supplement 2
Immunolocalization of CRABP1 in E15 chick utricle.

The sensory epithelium and the stroma are indicated in the transmitted-light section. Nuclei were stained with DAPI and F-actin with phalloidin; MYO7A and CRABP1 were stained with specific …

https://doi.org/10.7554/eLife.50777.012
Figure 3—figure supplement 3
TMSB4X localization.

The sensory epithelium and the stroma are indicated in the transmitted-light section. Nuclei were stained with DAPI and F-actin with phalloidin; OTOF and TMSB4X were stained with specific …

https://doi.org/10.7554/eLife.50777.013
Figure 3—figure supplement 4
G-actin localization with JLA20 antibody.

The sensory epithelium and the stroma are indicated in the transmitted-light section. Nuclei were stained with DAPI and F-actin with phalloidin; MYO7A and G-actin were stained with specific …

https://doi.org/10.7554/eLife.50777.014
Figure 3—figure supplement 5
OCM localization with anti-PV3 antibody.

The sensory epithelium and the stroma are indicated in the transmitted-light section. Nuclei were stained with DAPI and F-actin with phalloidin; OTOF and OCM were stained with specific antibodies. …

https://doi.org/10.7554/eLife.50777.015
Figure 3—figure supplement 6
Identification of striola and extrastriola regions.

The sensory epithelium and the stroma are indicated in the transmitted-light section. Nuclei were stained with DAPI and F-actin with phalloidin; SOX2, MYO7A and tubulin beta-3 (TUJ1) were stained …

https://doi.org/10.7554/eLife.50777.016
Pseudotemporal ordering of single utricle cells based on proteomics measurements.

(A) Relationship between variance and mean expression, distinguishes proteins or protein groups with low variance (blue) or high variance (salmon). (B) First three components of CellTrails' spectral …

https://doi.org/10.7554/eLife.50777.018
Developmental trajectory identified from single utricle cell RNA-seq measurements.

(A) CellTrails identifies nine distinct states (cellular subgroups). Shown are 328 points representing single utricle cells projected into two-dimensional space using CellTrails t-distributed …

https://doi.org/10.7554/eLife.50777.019
Figure 5—source data 1

CellTrails analysis of single-cell RNA-seq data.

Excel file including: Tab S1A: Single-cell RNA-seq measurements. Count matrix after QC and Normalization (ScNorm) of 186 genes and 328 chicken utricle single cells. Rows display single cells and columns display genes. Non-detected signals were set to 0. Tab S1B: tSNE plots and CellTrails maps. Listed are tSNE plots and CellTrails maps for each gene. The first column lists the gene symbol names, the second column shows the tSNE plots, the third column shows the CellTrails maps with the raw expression values, the fourth column shows the CellTrails maps with smoothed expression values, and the fifth column shows the CellTrails maps with the topographical expression surface. Tab S1C: Inferred expression dynamics. Shown is individual pseudotime gene expression along the extrastriolar (TrES, TrES*) and striolar (TrS) trajectory.

https://doi.org/10.7554/eLife.50777.020

Tables

Table 1
Summary statistics with confidence intervals for Figure 2C.
https://doi.org/10.7554/eLife.50777.017
95% Confidence interval
ComparisonEstimateStd. errordft-valueLower boundUpper boundp-value
ACTG1 FM-high – ACTG1 FM-low−0.0170.00922.3−1.830−0.0360.0020.081
ACTG1 FM-high – TMSB4X FM-high0.0360.00812.04.6070.0190.0540.001
ACTG1 FM-low - TMSB4X FM-low0.0040.00812.00.451−0.0140.0210.66
TMSB4X FM-low - TMSB4X FM-low−0.0500.00922.3−5.356−0.069−0.031<0.001
Key resources table
Reagent type
(species) or
resource
DesignationSource or referenceIdentifiersAdditional
information
Biological sample (Gallus gallus)Embryonic day 15 utricleeggs from Texas A and M University Poultry Science Departmentn/aFreshly isolated from Gallus gallus
Chemical compound, drugFM1-43FXThermo FisherCat# F3535510 µM
Chemical compound, drugSYTOX Red Dead Cell StainThermo FisherCat# S348591:1000 (final 5 nM)
Chemical compound, drugthermolysin from geobacillus stearothermophilusSigma-AldrichCat# T79020.5 mg/ml
Chemical compound, drugAccutaseInnovative Cell TechnologiesCat# AT104full strength
Chemical compound, drugn-dodecyl β-D-maltosideSigma-AldrichCat# D46410.1% (w/v)
Chemical compound, drugtrypsinPromegaCat# V528010 ng/µl
Chemical compound, drugLys-CPromegaCat# V167110 ng/µl
Chemical compound, drugDAPIThermo FisherCat# D13061 µg/ml
Chemical compound, drugAlexa Fluor 488-conjugated phalloidinThermo FisherCat# A123791:1000
Antibodyrabbit polyclonal anti-TMSB4XProteintechCat# 19850–1-AP, RRID: AB_106424371:250
Antibodyrabbit polyclonal anti-AGR3ProteintechCat# 11967–1-AP, RRID: n/a1:250
Antibodyrabbit polyclonal anti-CRABP1ProteintechCat# 12588–1-AP, RRID: AB_22922711:250
Antibodymouse monoclonal anti-G-actinDevelopmental Studies Hybridoma BankCat# JLA20, RRID: AB_5280681:250
Antibodymouse monoclonal anti-otoferlinDevelopmental Studies Hybridoma BankCat# HCS-1, RRID: AB_108042961:250
Antibodygoat polyclonal anti-SOX2Santa Cruz BiotechnologyCat# sc-17320, RRID: AB_22866841:100
Antibodymouse monoclonal anti-tubulin beta-3BioLegendCat# TUJ1, RRID: AB_23137731:250
Antibodyrabbit polyclonal anti-parvalbumin-3/oncomodulinHeller laboratoryn/a1:1000
Antibodyrabbit polyclonal anti-MYO7AProteus BiosciencesCat# 25–6790, RRID: AB_23148381:1000
AntibodyAlexa Fluor 546 donkey anti-rabbit polyclonalThermo FisherCat# A10040, RRID: AB_25340161:250
AntibodyAlexa Fluor 647 donkey anti-mouse polyclonalThermo FisherCat# A31571, RRID: AB_1625421:100
AntibodyAlexa Fluor 488
donkey anti-goat
polyclonal
Thermo FisherCat# A11055, RRID: AB_25341021:250
Commercial assay or kitSMARTscribeClontechCat# 639538
Commercial
assay or kit
Hifi HotStart ReadyMix (2X)Kapa BiosystemsCat# KK2602
Software, algorithmMaxQuantCox lab, Max Planck Institute of Biochemistryhttps://www.maxquant.org
Software, algorithmFIJI (ImageJ)n/ahttp://fiji.sc/
Software, algorithmRThe R Project for Statistical Computinghttps://www.r-project.org/
Software, algorithmlimma R packageBioconductorDOI: 10.18129/B9.bioc.limma
Software, algorithmsva R packageBioconductorDOI: 10.18129/B9.bioc.sva
Software, algorithmCellTrails R packageBioconductorDOI: 10.18129/B9.bioc.CellTrails
OtherMedium 199Thermo FisherCat# 12350039
OtherSPHERO Drop Delay Calibration ParticlesSpherotechCat# DDCP-70–2
OthernanoPOTS chipsCustom builtn/a
Table 2
Antibodies used.
https://doi.org/10.7554/eLife.50777.021
AntibodyDilutionSupplierCatalog #Validation
rabbit anti-TMSB4X1:250Proteintech19850–1-AP(Zhou et al., 2013; Li et al., 2018)
rabbit anti-AGR3 1:2501:250Proteintech11967–1-APFrom manufacturer’s website: correct-sized band by protein immunoblot in SKOV-3 cells, MCF-7 cells (chicken protein 88% identical to mouse)
rabbit anti-CRABP1 1:2501:250Proteintech12588–1-APFrom manufacturer’s website: correct-sized band by protein immunoblot in human spleen tissue, transfected HEK-293 cells (chicken protein 95% identical to human)
mouse anti-G-actin 1:2501:250Developmental Studies Hybridoma BankJLA20(Lin, 1981)
mouse anti-otoferlin 1:2501:250Developmental Studies Hybridoma BankHCS-1(Goodyear et al., 2010) (chicken protein 73% identical to mouse)
goat anti-SOX2 1:1001:100Santa Cruz Biotechnologysc-17320From manufacturer’s website: ‘recommended for detection of Sox-2 of mouse, rat, human and avian origin by WB, IP, IF, IHC(P) and ELISA’
mouse anti-tubulin beta-3 TUJ1. 1:2501:250BioLegend801202(Lee et al., 1990)
rabbit anti-parvalbumin3/oncomodulin1:1000Heller lab16910(Heller et al., 2002)
rabbit anti-MYO7A1:1000Proteus Biosciences25–6790(Morgan et al., 2016)
Table 3
Genes included in scRNA-seq CellTrails analysis.
https://doi.org/10.7554/eLife.50777.022
ABCA5(WHRN)(NSG2) PLSCR5TMSB4X
ACO1 DFNB59 LOC423919 PNPT1 TNNC2
ACTB(PJVK)(SHTN1) PODXL2 TOLLIP
ACTG1 DIAPH1 LOC772075 POU4F3 TPM1
ACTN1 DNM1(XIRP2) PPP1R14D TPM3
ADGRV1 DPF3 LOXHD1 PRPS1 TPRN
AGR3 DRGX MAP1A PTPRQ TRIOBP
AK1 EFCAB6 MAPK10 PTPRT TTLL12
AKAP5 EFR3A MCOLN3 PTPRZ1 TUBA3E
ANKRD24 ELMOD1 MPRIP RAB26 TUBAL3
APPL2 EML1 MSN RDX TUBB2B
ARF1 EPS8L2 MSRB3 RFX8 TUBB6
ARF4 ESPN MYH9 RPS6KA2 TWF2
ARHGAP17 ESPNL MYO15A RPS6KA5 USH1C
ARMC4 EZR MYO1C RSPH1 USH1G
ATOH1 FOXJ1 MYO1H RSPH9 USH2A
ATP2B1 FSCN1 MYO3A SCG3
ATP2B2 FSCN2 MYO3B SERPINB6
ATP2B4 GALNT9 MYO6 SGCB
ATP6V1B2 GAPDH MYO7A SGCG
ATP6V1E1 GDI2 NFATC1 SGIP1
ATP8B1 GNAI1 NMNAT2 SH3GLB2
B3GNTL1 GNAI2 NPEPPS SKOR2
BAIAP2L2 GNAI3 OCM SLC17A8
BRSK2 GNAL OSBP2 SLC8A1
CAB39L GNAS OSBPL11 SLC9A3R2
CACNA2D2 GNG4 OSBPL1A SMPX
CALB2 GPSM2 OTOA SPAG1
CAPZA1 GPX2 OTOF SPTAN1
CAPZA2 GRXCR1 PAICS STARD10
CAPZB GRXCR2 PAK1 STXBP1
CCDC50 GSTO1 PAK2 SYN3
CDH23 HSF5 PAK3 TECTA
CHRNA10 HYDIN PCDH15 TECTB
CHRNA9 IRX2 PDCD6IP TMC1
CIB2 KIAA1211L PDK4 TMC2
CKB(KIAA1211) PDZD7 TMC5
CLIC5 KIAA1549 PGM2L1 TMCC2
CORO2B KIF1A PHF21B TMEM117
CRABP1 KLHDC7A PI4KA TMEM255B
CSNK2A1 LCP1 PITPNA TMEM30A
CTH LHFPL5 PITPNB TMIE
CUL1 LHX3 PLS1 TMPRSS3
DFNB31 LOC416212 PLS3 TMPRSS7
ABCA5(WHRN)(NSG2) PLSCR5TMSB4X
ACO1 DFNB59 LOC423919 PNPT1 TNNC2
ACTB(PJVK)(SHTN1) PODXL2 TOLLIP
ACTG1 DIAPH1 LOC772075 POU4F3 TPM1
ACTN1 DNM1(XIRP2) PPP1R14D TPM3
ADGRV1 DPF3 LOXHD1 PRPS1 TPRN
AGR3 DRGX MAP1A PTPRQ TRIOBP
AK1 EFCAB6 MAPK10 PTPRT TTLL12
AKAP5 EFR3A MCOLN3 PTPRZ1 TUBA3E
ANKRD24 ELMOD1 MPRIP RAB26 TUBAL3
APPL2 EML1 MSN RDX TUBB2B
ARF1 EPS8L2 MSRB3 RFX8 TUBB6
ARF4 ESPN MYH9 RPS6KA2 TWF2
ARHGAP17 ESPNL MYO15A RPS6KA5 USH1C
ARMC4 EZR MYO1C RSPH1 USH1G
ATOH1 FOXJ1 MYO1H RSPH9 USH2A
ATP2B1 FSCN1 MYO3A SCG3
ATP2B2 FSCN2 MYO3B SERPINB6
ATP2B4 GALNT9 MYO6 SGCB
ATP6V1B2 GAPDH MYO7A SGCG
ATP6V1E1 GDI2 NFATC1 SGIP1
ATP8B1 GNAI1 NMNAT2 SH3GLB2
B3GNTL1 GNAI2 NPEPPS SKOR2
BAIAP2L2 GNAI3 OCM SLC17A8
BRSK2 GNAL OSBP2 SLC8A1
CAB39L GNAS OSBPL11 SLC9A3R2
CACNA2D2 GNG4 OSBPL1A SMPX
CALB2 GPSM2 OTOA SPAG1
CAPZA1 GPX2 OTOF SPTAN1
CAPZA2 GRXCR1 PAICS STARD10
CAPZB GRXCR2 PAK1 STXBP1
CCDC50 GSTO1 PAK2 SYN3
CDH23 HSF5 PAK3 TECTA
CHRNA10 HYDIN PCDH15 TECTB
CHRNA9 IRX2 PDCD6IP TMC1
CIB2 KIAA1211L PDK4 TMC2
CKB(KIAA1211) PDZD7 TMC5
CLIC5 KIAA1549 PGM2L1 TMCC2
CORO2B KIF1A PHF21B TMEM117
CRABP1 KLHDC7A PI4KA TMEM255B
CSNK2A1 LCP1 PITPNA TMEM30A
CTH LHFPL5 PITPNB TMIE
CUL1 LHX3 PLS1 TMPRSS3
DFNB31 LOC416212 PLS3 TMPRSS7

Additional files

Download links