henhaobio
Sign in

eLifeDevelopmental BiologyPeer-reviewed journal article

A protein phosphatase network controls temporal and spatial dynamics of differentiation commitment in human epidermis

Ajay MishraiD, Bénédicte OulèsiD, Angela Oliveira PiscoiD, Tony LyiD, Kifayathullah Liakath‐AliiD, Gernot WalkoiD, Priyalakshmi ViswanathaniD, Matthieu TihyiD, Jagdeesh Nijjher, Sara-Jane DunniD, Angus I. LamondiD, Fiona M. WattiD

PDFDOIPreprint

In-depth summary

Claude Opus 5.5 · xhigh effort · 8 Nov 2026 · from the full text of version 1 (8 Apr 2017) · CC BY 4.0

In suspended human keratinocytes, five induced phosphatases dampen ERK to drive differentiation commitment, opposed by DUSP10, in a time-rewired network that renders commitment an unstable switch.

Commitment, the transient step between epidermal stem cell renewal and terminal differentiation, has no markers and is missed by end-point analyses, partly because fast state changes can depend on post-translational events such as phosphorylation. Time-resolved transcript, protein and phosphosite profiling plus functional screens identify a phosphatase network that sets when commitment starts and how long it lasts.

Key results

  1. ShownThe 4 h commitment window is a distinct cell state. Protein levels track transcripts poorly and dephosphorylation dominates, which points to post-transcriptional control of commitment.

    Clonogenicity fell sharply between 4 and 8 h. t-SNE and clustering separated 4 h from 8 and 12 h. Protein-transcript Pearson r was 0.19 at 4 h vs 0.51 (8 h) and 0.68 (12 h). Over 3,500 phosphosites (Andromeda score >= 30) were quantified; about two thirds of 4 h changes were dephosphorylation not explained by protein abundance. · Results and Discussion; Figure legends

  2. ShownAn siRNA screen of phosphatases induced at commitment identified DUSP6, PPTC7, PTPN1, PTPN13 and PPP3CA as pro-commitment factors. Depleting them increases self-renewal, and depleting DUSP6, PTPN1, PPP3CA or PTPN13 blurs basal versus suprabasal identity in reconstituted epidermis.

    45 phosphatases changed from 0 to 4 h (20 up). Silencing seven of the 20 raised clonal growth, the five named most strongly (p-value <0.001; n=3 screens, 3 cultures each). Knockdowns delayed TP63 loss and TGM1 gain. In organotypic cultures (n=2 transfections) TP63+ fractions rose and TP63/Ki67+ cells appeared suprabasally. · Results and Discussion; Figure legends

  3. ShownPro-commitment phosphatases lower ERK1/2 activity and promote induction of MAF-family AP1 factors, since their knockdown blunts it. DUSP10, also induced at 4 h, opposes commitment, and its knockdown does not raise ERK phosphorylation.

    ERK2 (MAPK1) and ERK1 (MAPK3) sites ranked in the top 15 most dephosphorylated at 4 h. Knockdowns raised phospho-ERK1/2 vs scrambled control (n = 2 blots) and blunted MAF, MAFB, MAFG induction. DUSP10 knockdown cut clonogenicity (p-value<0.001) without raising phospho-ERK. Overexpressed DUSP10 raised and DUSP6 lowered colony formation; DUSP6 C293S had no effect. · Results and Discussion; Figure legends

  4. SuggestedPhosphatase cross-regulation rewires over time: mostly negative at 0, 8 and 12 h, but almost entirely positive around a self-amplifying DUSP6 at 4 h. Boolean reasoning indicates the network must switch, framing commitment as an unstable switch.

    Knockdown-by-time RT-qPCR networks (2-way ANOVA, Holm-Sidak) show all but one interaction positive at 4 h. RE:SIN needed added possible edges to satisfy constraints. RE:IN found that no single Abstract Boolean Network reproduces the dynamics. PKCi 12 h expression fit the 0 h network; neither PKCi nor TSA reached the 12 h network. · Results and Discussion; Network analysis; Figure legends

  5. SuggestedCommitment phosphatases are spatially patterned in the basal layer of human epidermis, and the pattern is reproduced on topographic substrates in vitro. This implies that dermal cells are not required for the patterning.

    Wholemount confocal imaging: DUSP6, PTPN1 and PPP3CA were highest in beta1-integrin-high stem cell clusters, PPTC7 was enriched in integrin-low regions, and PTPN13 and DUSP10 were uniform. On collagen-coated PDMS substrates mimicking the epidermal-dermal junction, DUSP6-high clusters formed on feature tops while DUSP10 stayed uniform. · Results and Discussion; Figure legends

Methods, data and reproducibility

Filtered primary human keratinocytes were held in methylcellulose suspension for 0 to 12 h and profiled by microarray, SILAC proteomics and phosphoproteomics. Upregulated phosphatases were tested by siRNA knockdown, overexpression, clonogenicity, organotypic culture and RT-qPCR, and their cross-regulation was modelled with Abstract Boolean Networks.

  • Illumina BeadArray time course of filtered neonatal foreskin keratinocytes (strain km) at 0, 4, 8 and 12 h in methylcellulose suspension · human · n = n=3 independent replicates per condition · GSE73147
  • SILAC LC-MS/MS total proteome (4, 8, 12 h vs 0 h) and HILIC/TiO2-enriched phosphoproteome (4 and 8 h vs 0 h) from the same suspension model · human · n = n=3 replicates (proteome); over 3,500 phosphosites · PXD003281
  • Published dataset of calcium-induced stratification of primary keratinocytes (Hopkin et al., 2012), reanalysed for phosphatase expression · human · GSE38628
  • Published dataset of ex vivo epidermal reconstitution, reanalysed for phosphatase expression · human · GSE52651
  • siRNA knockdown screen of 20 commitment-upregulated phosphatases plus INPP5J control, scored by clonogenicity and live-cell growth · human · n = n=3 independent screens, 3 cultures per screen
  • Organotypic epidermis on de-epidermised human dermis after knockdown of DUSP6, PTPN1, PPP3CA or PTPN13, cultured three weeks at the air-liquid interface · human · n = n=2 independent transfections
  • Epidermal wholemounts from breast or abdominal surgical-waste skin, immunostained for phosphatases and beta1 integrin · human
  • Keratinocytes cultured on collagen-coated PDMS substrates mimicking the topography of the epidermal-dermal junction · human
Methods & tools
Illumina BeadArray microarrays with BeadStudio, GeneSpring GX13.1, SILAC LC-MS/MS (Q-Exactive Orbitrap) with HILIC/TiO2 phosphopeptide enrichment, MaxQuant with Andromeda search engine, Revigo, t-SNE and hierarchical clustering, RE:IN and RE:SIN (Abstract Boolean Network reasoning), Cytoscape V3.2.1, SMART pool ON-TARGET plus siRNAs (Amaxa nucleofection), Doxycycline- and Cumate-inducible lentiviral vectors (pCW57-GFP-2A-MCS, QM812B-1), GraphPad Prism 7 and MeV_4_8, ImageJ, Incucyte live-cell imaging, NIS Elements
Reagents & environment
GeneSpring GX13.1 (v13), Cytoscape V3.2.1, GraphPad Prism 7, MeV_4_8, NIS Elements version 4.00.04; R and MaxQuant used without stated versions.
Reproducibility gaps
Data are in GEO GSE73147 and PRIDE PXD003281, but no repository is given for custom analysis or network-inference scripts. The RE:IN/RE:SIN software is accessible at research.microsoft.com/rein.
Limitations (5)
  • Phosphatase networks inferred from single knockdowns could not satisfy the Boolean experimental constraints, which implies undetected interactions. Extra possible edges, taken from overexpression and protein-level data, were added to obtain consistent models. · stated by the authors
  • No molecular markers of commitment exist, so its onset (approximately 4 h) and progression are defined indirectly, by loss of clonal growth after cells are replated from suspension. · stated by the authors
  • The proposed link to psoriasis, based on basal markers appearing suprabasally after knockdown, is framed by the authors as speculation; no diseased skin was examined. · stated by the authors
  • Counts differ between the main text and the figure legends (45 vs 47 phosphatases differentially expressed at 4 h; 3 vs 8 Boolean experimental constraints), which makes it harder to reconstruct the screen and the model exactly. · noted in this summary
  • Functional work uses one neonatal foreskin keratinocyte strain (km), and several key readouts rest on few replicates (n=2 organotypic transfections, n = 2 blots for phospho-ERK quantification). · noted in this summary

Read this if you study epidermal stem cell fate decisions or phosphatase control of ERK/AP1 signalling, or want an example that combines time-resolved multi-omics, siRNA screening and Boolean network modelling of a cell-state switch.

Citations
69
FWCI
3.46
Field percentile
top 6%

Abstract

Epidermal homeostasis depends on a balance between stem cell renewal and terminal differentiation. The transition between the two cell states, termed commitment, is poorly understood. Here we characterise commitment by integrating transcriptomic and proteomic data from disaggregated primary human keratinocytes held in suspension for up to 12h to induce differentiation. We find that cell detachment induces a network of protein phosphatases. The pro-commitment phosphatases - including DUSP6, PPTC7, PTPN1, PTPN13 and PPP3CA - promote differentiation by negatively regulating ERK MAPK and positively regulating AP1 transcription factors. Their activity is antagonised by concomitant upregulation of DUSP10. Boolean network modelling of phosphatase interactions identifies commitment as an inherently unstable biological switch between the stem and differentiated cell states. Furthermore, phosphatase expression is spatially regulated both in vivo and in vitro. We conclude that an auto-regulatory phosphatase network maintains epidermal homeostasis by controlling the onset and duration of commitment.

Reproducibility ledger

Nobody has reported trying to reproduce this paper yet. Did you replicate an experiment or run the code? Your report saves the next lab a week.

Sign in to file a reproducibility report.