- Data Note
- Open Access
- Open Peer Review
Spatially transformed fluorescence image data for ERK-MAPK and selected proteins within human epidermis
GigaScience volume 4, Article number: 63 (2015)
Phosphoprotein signalling pathways have been intensively studied in vitro, yet their role in regulating tissue homeostasis is not fully understood. In the skin, interfollicular keratinocytes differentiate over approximately 2 weeks as they traverse the epidermis. The extracellular signal-regulated kinase (ERK) branch of the mitogen-activated protein kinase (MAPK) pathway has been implicated in this process. Therefore, we examined ERK-MAPK activity within human epidermal keratinocytes in situ.
We used confocal microscopy and immunofluorescence labelling to measure the relative abundances of Raf-1, MEK1/2 and ERK1/2, and their phosphorylated (active) forms within three human skin samples. Additionally, we measured the abundance of selected proteins thought to modulate ERK-MAPK activity, including calmodulin, β1 integrin and stratifin (14-3-3σ); and of transcription factors known to act as effectors of ERK1/2, including the AP-1 components Jun-B, Fra2 and c-Fos. Imaging was performed with sufficient resolution to identify the plasma membrane, cytoplasm and nucleus as distinct domains within cells across the epidermis. The image field of view was also sufficiently large to capture the entire epidermis in cross-section, and thus the full range of keratinocyte differentiation in a single observation. Image processing methods were developed to quantify image data for mathematical and statistical analysis. Here, we provide raw image data and processed outputs.
These data indicate coordinated changes in ERK-MAPK signalling activity throughout the depth of the epidermis, with changes in relative phosphorylation-mediated signalling activity occurring along the gradient of cellular differentiation. We believe these data provide unique information about intracellular signalling as they are obtained from a homeostatic human tissue, and they might be useful for investigating intercellular heterogeneity.
Dysregulated signalling is a common oncogenic driver, and a number of newer cancer drugs target components of the intracellular signalling network . Phosphoprotein signalling has been studied extensively in vitro, and this has provided detailed knowledge of the molecular interactions that propagate signals through these networks . However, to fully understand the role that signalling protein mutations have in oncogenesis, and to design treatments with minimal side effects, we need to elucidate the role that signalling pathways have in controlling cellular behaviour within normal, homeostatic tissues in situ.
In human skin, the interfollicular epidermis is a stratified epithelial tissue where keratinocytes are arranged in a gradient of cellular differentiation across the depth of the tissue . Cellular proliferation occurs within the deepest basal layer, and keratinocytes that leave this layer undergo terminal differentiation as they traverse the epidermis towards the surface layer, a process that takes approximately 2 weeks . This process establishes a spatiotemporal differentiation gradient, such that the position of a keratinocyte within the epidermis is related to its stage of differentiation. Hence, the human epidermis is a useful model system to study intracellular signalling in situ in a homeostatic tissue.
A number of regulatory mechanisms control keratinocyte behaviour to ensure epidermal tissue function [4, 5], and ERK1/2 signalling has been implicated in controlling keratinocyte differentiation both in vivo  and in vitro . Therefore, we examined ERK-MAPK activity within human epidermis to try to elucidate the role of intracellular signalling in controlling adherent cell behaviour in situ. A selection of proteins that modulate the ERK-MAPK pathway within human epidermis were also examined, as were several components of the AP-1 transcription factor family  that are regulated, in part, by ERK1/2 activity (Fig. 1).
We performed single-target labelling against a number of proteins related to the ERK-MAPK signalling cascade (Table 1 and Fig. 1) in three samples of human epidermis. A small number of image stacks with multiple-target labelling are also provided.
Image data were collected with a pixel resolution near the diffraction limit, which enabled the cytoplasm, nucleus and plasma membrane to be distinguished, while retaining a field of view large enough to capture the entire epidermis in cross-section (Fig. 2). Within the plane of the epidermis, numerous cells at similar stages of differentiation are observable, which makes it possible to perform within-sample replicate measurements.
To perform mathematical analyses with these image data, the fluorescence signal intensity of each target was sampled across the epidermis and transformed to a one-dimensional ‘normalised distance’ coordinate (for detailed information, see Additional file 1). Briefly, the morphological features of keratinocytes that were apparent within the single-target image data were used to segment the epidermis into discrete tissue layers. These features include juxtaposition of basal keratinocytes and the underlying dermis  and flattening of keratinocytes within the transitional layer owing to caspase-14-mediated release of filaggrin . The spinous and granular layers could not be distinguished without additional histological markers and were thus considered together. After sampling, the relative position of each pixel within these layers was calculated using linear interpolation. This layer-normalised coordinate mapped fluorescence signal intensity against the gradient of keratinocyte differentiation while maintaining segregation of phenotypically distinct epidermal layers, and helped to reduce variation between and within patients associated with variations in the epidermal and tissue layer thicknesses (Figure AF1.3 in Additional file 1).
Raw image data are provided as .TIFF files, together with .TXT files containing microscope settings. Processed data are provided as .TIFF files of sample locations, and .TXT files of sampled data in an x-y coordinate format, and a normalised distance format.
Human skin samples were obtained with written informed consent under a protocol approved by the New Zealand Northern Regional X Ethics Committee (project number NTX/08/09/086) and Counties-Manukau District Health Board (project number 681).
Full experimental methods are given in Additional file 2. Briefly, human skin samples were embedded in optimal cutting temperature medium, snap frozen in liquid nitrogen and stored at -80 °C. Tissue blocks were incubated at -20 °C for at least 30 min, 20 μm sections were taken in the cryotome and adhered to SuperFrost + glass slides at approximately 20 °C (room temperature), then immediately processed or stored at -20 °C for up to 24 h.
Tissue sections were fixed in phosphate buffered saline (PBS) with 4 % (v/v) formaldehyde for 30 min at 20 °C, after which sections were washed and blocked using 5 % (v/v) fetal bovine serum in PBS for 30 min at 20 °C. Sections were probed with primary antibodies (dilutions and manufacturer details in Additional file 3) for 1 h at 20 °C before being washed again. Washed sections were then probed with a fluorophore-conjugated secondary antibody for 1 h at 20 °C, in the dark. Sections underwent final washes before being mounted in Citifluor AF-1 and incubated for at least 24 h in the dark at 20 °C before imaging or storage at 4 °C.
A Leica TCS SP2 confocal microscope was used for imaging with the Leica Confocal Software (version 2.61) and data were stored in an uncompressed TIFF format. A 63× oil-immersion objective lens (numerical aperture 1.32) was used, and 1,024 × 1,024 pixel images were captured at or near the diffraction limit (approximately 160 nm/pixel). The microscope gain and offset were adjusted to utilise the full dynamic range, minimising pixel underflow and saturation at the tissue depth with the strongest apparent signal. Image stacks were collected over a z-interval of approximately 10 μm (this value varies across image data) using 4× frame-averaging.
Commercially available antibodies were used (Additional file 3), most of which had been used in reports across the literature (please refer to vendor data sheets); however, some of these studies used alternative methods such as western blotting or an enzyme-linked immunosorbent assay (ELISA) and thus, it is difficult to state with certainty that antibodies are binding the correct target in situ. Furthermore, labelling quality varied between targets and the efficacy of different primary antibodies varies. We have summarised the image data quality with a subjective metric and some notes (Table 1), and compared our results with those from previous studies of epidermal biology (Additional file 4).
Autofluorescence and non-specific signals within the stratum corneum were fairly common across our data; however, these sources of noise were excluded from the quantitative analysis by tissue segmentation applied during image sampling.
We have used some of the LOESS-smoothed  data to examine bulk changes in signalling pathway components throughout the depth of the epidermis. Dynamic Bayesian network inference methods found evidence for a number of canonical signalling interactions  within these data. Further, quantitative protein abundance data were used to examine statistical associations between active (phosphorylated) ERK-MAPK components, and as inputs to a normalised-Hill differential equation model  of ERK-MAPK signalling. Our results suggest that spatial gradients of Ca2+ and plasma membrane calmodulin across the epidermis are sufficient to drive changes in the steady-state activity of ERK-MAPK signalling, in a manner that was consistent with the gradient of ERK-MAPK activation that we observed and previously published .
These data will be useful for researchers to investigate intercellular heterogeneity of signalling pathway components in the context of a homeostatic tissue. To aid these applications, we have included some preliminary ‘whole-cell’ segmentation masks developed for various analyses, and data that can be extracted through this approach are shown in Additional file 5.
Finally, we believe that the layer-normalisation process introduced during data processing (which is detailed in Additional file 1) should make it feasible for researchers of epidermal biology to integrate new data with our dataset and extend these analyses.
Availability and requirements
Project name: Epidermal MAPK Data Scripts
Project home page: http://www.github.com/uomsystemsbiology/epidermal_data
Operating systems: platform-independent
Programming language: MATLAB
Other requirements: MATLAB (2012b or later), including the Image Processing, Curve Fitting and Neural Network Toolboxes
License: MIT license
Restrictions to use by non-academics: none
Project name: Virtual Reference Environment for Epidermal MAPK Data Scripts (further details on virtual reference environments  can be found at http://uomsystemsbiology.github.io/research/reference-environments)
Project home page: http://www.github.com/uomsystemsbiology/epidermal_data_reference_environment
Operating systems: platform-independent
Programming languages: MATLAB
Other requirements: Vagrant (version 1.7.2 or higher), VirtualBox (version 4.3.x or higher)
License: MIT license
Restrictions to use by non-academics: to replicate these results using the reference environment, users must accept the MATLAB® Compiler Runtime Libraries License which is displayed during installation.
Availability of supporting data
The dataset supporting the results of this article is available in the GigaScience Database repository .
Extracellular signal-regulated kinase
Mitogen-activated protein kinase
Phosphate buffered saline
Samatar AA, Poulikakos PI. Targeting RAS-ERK signalling in cancer: promises and challenges. Nat Rev Drug Discov. 2014;13(12):928–42. doi:10.1038/nrd4281.
Kolch W, Pitt A. Functional proteomics to dissect tyrosine kinase signalling pathways in cancer. Nat Rev Cancer. 2010;10(9):618–29. doi:10.1038/nrc2900.
Eckert RL, Crish JF, Robinson NA. The epidermal keratinocyte as a model for the study of gene regulation and cell differentiation. Physiol Rev. 1997;77:397–424.
Blanpain C, Fuchs E. Epidermal homeostasis: a balancing act of stem cells in the skin. Nat Rev Mol Cell Biol. 2009;10:207–17. doi:10.1038/nrm2636.
Candi E, Schmidt R, Melino G. The cornified envelope: a model of cell death in the skin. Nat Rev Mol Cell Biol. 2005;6:328–40. doi:10.1038/nrm1619.
Eckert RL, Efimova T, Dashti SR, Balasubramanian S, Deucher A, Crish JF, et al. Keratinocyte survival, differentiation, and death: many roads lead to mitogen-activated protein kinase. J Investig Dermatol Symp Proc. 2002;7:36–40. doi:10.1046/j.1523-1747.2002.19634.x.
Bromberg JF, Fan Z, Brown C, Mendelsohn J, Darnell JE. Epidermal growth factor-induced growth inhibition requires Stat1 activation. Cell Growth Differ. 1998;9:505–12.
Angel P, Szabowski A, Schorpp-Kistner M. Function and regulation of AP-1 subunits in skin physiology and pathology. Oncogene. 2001;20:2413–23. doi:10.1038/sj.onc.1204380.
Cleveland WS, Devlin SJ. Locally weighted regression: an approach to regression analysis by local fitting. J Am Stat Assoc. 1988;83:596–610. doi:10.2307/2289282.
Cursons J, Hurley D, Angel CE, Dunbar R, Crampin EJ, Jacobs MD. Inference of an in situ epidermal intracellular signaling cascade. Conf Proc IEEE Eng Med Biol Soc. 2010:799–802. doi:10.1109/IEMBS.2010.5626520.
Kraeutler MJ, Soltis AR, Saucerman JJ. Modeling cardiac beta-adrenergic signaling with normalized-Hill differential equations: comparison with a biochemical model. BMC Syst Biol. 2010;4:157. doi:10.1186/1752-0509-4-157.
Cursons J, Gao J, Hurley DG, Print CG, Dunbar PR, Jacobs MD, et al. Regulation of ERK-MAPK signaling in human epidermis. BMC Syst Biol. 2015;9:41. doi:10.1186/s12918-015-0187-6.
Hurley DG, Budden DM, Crampin EJ. Virtual Reference Environments: a simple way to make research reproducible. Brief Bioinform. 2014. doi:10.1093/bib/bbu043.
Cursons, J; Angel, C, E; Hurley, D, G; Print, C, G; Dunbar, P; Jacobs, M, D; Crampin, E, J (2015): Supporting data for “Spatially-transformed fluorescence image data for ERK-MAPK and selected proteins within human epidermis”. GigaScience Database. http://dx.doi.org/10.5524/100168.
Le Novere N, Hucka M, Mi H, Moodie S, Schreiber F, Sorokin A, et al. The systems biology graphical notation. Nat Biotechnol. 2009;27(8):735–41. doi:10.1038/nbt.1558.
We gratefully acknowledge the contributions of Jenni Chen (University of Auckland, New Zealand) during development of the immunofluorescence labelling methods, and the Biomedical Imaging Research Unit at the University of Auckland (New Zealand) for their assistance with confocal microscopy. JC received funding from the Tertiary Education Commission (New Zealand) Top Achiever Doctoral Scholarship. This research was funded through a New Economy Research Fund (New Zealand) grant for Advanced Skin Imaging (contract C08X0801), and funded in part by the Australian Research Council Centre of Excellence in Convergent Bio-Nano Science and Technology (project number CE140100036).
The authors declare that they have no competing interests.
JC, CGP, PRD, MDJ and EJC planned the study. JC, CEA, CGP, PRD, MDJ and EJC designed the experimental protocol. CEA collected the skin samples. JC performed the immunofluorescence labelling and collected the confocal microscopy data. JC, DGH and MDJ developed the image processing scripts. JC performed the image sampling. DGH produced the virtual reference environment. JC and EJC wrote the manuscript. All authors have read and approved the final manuscript.
Normalised distance information. (PDF 2308 kb)
Detailed experimental methods. (PDF 165 kb)
Detailed target information (including antibodies) and putative interactions around ERK-MAPK. (PDF 243 kb)
Comparison of image data to those from previous studies of epidermal biology. (PDF 12030 kb)
Intercellular heterogeneity of phospho-MEK1/2 signalling. (PDF 567 kb)
About this article
- Interfollicular keratinocytes
- Confocal microscopy
- Homeostatic tissue
- Cellular heterogeneity