Addressing Differentiation in Live Human Keratinocytes by Assessment of Membrane Packing Order

Introduction

The epidermis, in the form of a stratified and desquamating tissue is a fascinating feature of many organisms; it retains very distinct characteristics which enable it to exert function, shielding us from external threats and preventing fluid evaporation. An understanding is emerging on how these unique features contribute to the important protective role of the epidermis; these lessons are learned both from a physiological setting and through the perspective of skin diseases. For example, keratinocyte differentiation and cornification are essential for prevention of allergic diathesis by forming a functional skin barrier (Kubo et al., 2012; Flohr et al., 2014; Horimukai et al., 2014; Venkataraman et al., 2014). Cornification encompasses a series of processes leading to the programmed keratinocyte death, which results in formation of a cellular remnant known as “corneocyte.” Corneocytes, embedded in a lipid-enriched matrix (Eckhart et al., 2013) provide mechanical hardness, while the lipids support elastic properties of the stratum corneum. In addition, at the cellular level, the lamellar organization of lipids and their lateral packing properties prevent evaporation (Iwai et al., 2012; Janssens et al., 2012; Narangifard et al., 2018), forming a seal over the moisture trapped in the stratum corneum. This phenomenon is achieved by high abundance of hydrophilic compounds, known as Natural Moisturizing Factor (NMF).

Multiple processes take place during the cornification advancement in the epidermis; these occur simultaneously both at the organellar and molecular levels (Figure 1). The characteristics of progression in keratinocyte differentiation in vitro are reflected by multiple morphological changes (i.e., related to the cell shape and adherence) and distinct organellar transformations). For instance, the appearance and maturation of keratohyalin granules (KHGs) (Gutowska-Owsiak et al., 2018), nuclear condensation and extrusion (Gdula et al., 2013; Rogerson et al., 2018), cytoskeleton collapse (Gutowska-Owsiak et al., 2018), and alterations in mitochondria (Ipponjima et al., 2020); in some of those events, autophagy has been shown to be involved (Akinduro et al., 2016). At the molecular level, shifts in cytokeratin expression have been documented (keratin 5, keratin 14 are progressively replaced with keratin 1 and keratin 10) while additional proteins are upregulated, e.g., involucrin and transglutaminase as well as the markers of late differentiation, defined by the expression of late envelope proteins (e.g., involucrin, filaggrin, loricrin) (Marshall et al., 2001; Candi et al., 2005, 2016; Sandilands et al., 2009; Eckhart et al., 2013). These differentiation-dependent proteins contribute to barrier strengthening by incorporation into the insoluble cornified envelope (CE; also known as cornified cell envelope, CCE) or by supporting the hydration of the stratum corneum (Kezic et al., 2009).

FIGURE 1

Figure 1. Organellar and molecular changes in keratinocytes during the differentiation process. Stratified epidermis consists of multiple layers of epidermal keratinocytes characterized by increasing advancement of cell differentiation. Cells undergo extensive organellar and molecular alterations as they gradually move to the upper layers. Where SC, Stratum Corneum; SL, Stratum Lucidum; SG, Stratum Granulosum; SS, Stratum Spinosum; SB, Stratum Basale; FLG, filaggrin; KRT5, keratin 5; KRT14, keratin 14; KRT1, keratin 1; KRT10, keratin 10; LOR, loricrin; INV, involucrin; TG1, transglutaminase 1.

Along with the expression of protein markers and organellar changes, an overall increase in cell stiffness can be observed (Plasencia et al., 2007; Janssens et al., 2012; Wennberg et al., 2018). This stiffness is a consequence of two main processes; the formation of a rigid, protein-based CE and extensive lipid remodeling leading to specialized lipid profile adaptation (Ponec et al., 1987; Schmidt et al., 1991; López et al., 2007). Specifically, the high phospholipid and cholesterol content is gradually replaced by a predominance of ceramides and saturated triglycerides (Ponec et al., 1988), especially in conditions where stratification is supported. These lipid composition variations lead to more ordered membrane packing (Janssens et al., 2012), higher stiffness and consequently, changes in the relative lateral heterogeneity properties of the plasma membrane (Bernardino de la Serna et al., 2016). This selection has specific functional properties; the ceramides and free fatty acids provide attachment points for the CE proteins (Behne et al., 2000; López et al., 2007; Akiyama, 2017) and a scaffold to replace the plasma membrane by a corneocyte lipid envelope (CLE), which links CE with the lipids of the intercellular matrix (Elias et al., 2014). It is believed that ceramide enrichment displaces the balance between ordered and disordered regions to support membrane heterogeneity, where higher order regions work as docking molecular anchors supporting the CLE function and ultimately, stabilizing the CE structure. Regulation of ceramide synthesis is mediated via peroxisome proliferator-activated receptor (PPAR) pathways (Chon et al., 2015) associated with keratinocyte differentiation (Mao-Qiang et al., 2004); disturbances in these pathways result in clinical manifestations (Akiyama, 2017).

To date, no robust experimental readout for the assessment of differentiation advancement in live primary keratinocytes in real-time is available. Protein-targeted detection methods, such as the expression quantification by proteomics-based approaches, as well as other commonly used, e.g., immunostaining or western blot, require cell permeabilization or lysis and are only suitable for end-point measurements.

Previous studies have revealed the gel nature of the membranes in the stratum corneum (Plasencia et al., 2007). Here, we built upon the known lipid profile variation and ceramide enrichment during keratinocyte differentiation; we hypothesized that lipid lateral packing and order increase can be used as a hallmark for quantitative assessment of live-cell cornification advancement. For this purpose, we took advantage of the lipophilic polarity-sensitive membrane dye Laurdan, purposely designed by Weber and Farris (1979) to encompass an electron-donor and electron-acceptor, which displays large solvent-dependent fluorescence shifts. Laurdan ubiquitously distributes at the plane of the plasma membrane, regardless of its lipid lateral packing properties, and is oriented parallel with the phospholipid's hydrophobic tails in membranes. Moreover, its emission spectra and location are independent of phospholipid head groups and notably, indirectly report about the relative packing order by sensing the degree of water penetration in its vicinity (Parasassi et al., 1990, 1994, 1997; Bagatolli et al., 1998, 1999); the studies showed that Laurdan could sense the well-known phase transition from fluid to gel at 41C in in 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) bilayers. Most of these environmentally-sensitive probes display an increase in charge separation when excited in polar solvents resulting in a larger dipolar moment (Lakowicz, 1999); the penetration of water molecules into a bilayer formed of loosely packed lipids allow more rotational and translational freedom to the probe, yielding an emission shift toward longer wavelengths. On the contrary, the probe emission shifts toward the blue when its motion is more restricted due to sensing a milieu with lesser water content. Ratiometrically measuring the emission intensities at 440 and 490 nm renders the Laurdan “general polarization index” (GP), which allows relative quantification of the membrane order or degree of lipid lateral packing (Parasassi et al., 1990, 1997; Sezgin et al., 2015). This method has been used in visualizing native membrane microdomains in planar supported bilayers, giant unilamellar vesicles (de la Serna et al., 2009; Bernardino de la Serna et al., 2013) and in living cells (Gaus et al., 2003; Owen et al., 2011; Carugo et al., 2017; Gutowska-Owsiak et al., 2018; Bernabé-Rubio et al., 2019).

Even though the membrane order can only be measured indirectly, and being aware that the stiffness of a membrane influences can be in principle characterized by fluorescence polarization or fluorescence anisotropy, in this paper we introduced membrane lateral packing order as a readout for keratinocyte differentiation advancement. Using single-cell based confocal spectral imaging or a fluorescence plate reader for high-throughput screening of cellular ensembles we utilize the Laurdan general polarization index assessment of membrane order as an efficient readout for determination of the progression of keratinocyte differentiation in living cells.

Materials and Methods

Human Keratinocyte Culture and Calcium Switch

Normal human epidermal keratinocytes (NHEKs) were purchased from Lonza (Lonza, Basel, Switzerland, neonatal, pooled) and cultured in keratinocyte KBM-2 media (Lonza, Basel, Switzerland) at conditions supporting proliferation ([Ca2+] = 0.06 mM) and subcultured by accutase treatment (Sigma Aldrich, Dorset, United Kingdom). To promote terminal differentiation, calcium switch was performed by adjusting calcium level do the final concentration of [Ca2+] = 1.5 mM for a period of 24 h before immunofluorescent labeling. To adjust to the desired calcium concentration a calcium switch was conducted over a period of 24 h by replacing the culture media with media adjusted to the desired calcium concentration by adding CaCl2 (Sigma Aldrich, Dorset, United Kingdom sigma).

Recombinant Filaggrin

The expression plasmid was constructed using a SLIC method. Briefly, the nucleotide sequence encoding 7th FLG domain with N-terminal 6-His-tagged SUMO protein sequence was cloned into pET28 vector using BamHI i XhoI restrictions sites. The construct was transformed into E. coli BL21-CodonPlus-RIL and propagated in LB liquid media with antibiotics, followed by autoinduction media (Formedium AIM- Super Broth) for 48 h at 18°C. The cells were sonicated in lysis buffer (10 mM Tris pH 8, 150 mM NaCl, 10 mM imidazole) with protease inhibitors cocktail. The purification procedure was performed on an ÄKTA™Xpress chromatography system. Lysate was loaded on Ni-NTA Agarose column (Qiagen). On-column cleavage by SUMO protease (4 μg/1 ml) was carried out in elution buffer (10 mM Tris pH 8, 150 mM NaCl, 300 mM imidazole) for 8 h at 10°C. The protein was purified by desalting and on Ni-NTA Agarose columns. Flow though fractions were loaded into a Superdex 200 column (GE Healthcare) and analyzed by SDS-PAGE; protein identification was carried out by MALDI-TOF MS.

Fluorescent Antibody Staining and Confocal Microscopy Imaging

NHEKs were grown and calcium-switched in eight-well cell culture slides (Beckton Dickinson), fixed and permeabilized by neat acetone. Blocking was carried out in freshly made blocking buffer (5% FCS, Sigma-Aldrich, Gillingham, Dorset, United Kingdom; 2% BSA, Sigma-Aldrich, Gillingham, Dorset, United Kingdom in PBS). For the primary antibody anti-filaggrin antibody (mouse monoclonal 15C10 from Leica Biosystems, Milton Keynes, United Kingdom) and anti-keratin 14 (clone EPR1612 from GeneTex, Irvine, CA, United States) were used; specificity of the antibody was verified prior with both recombinant human FLG and on epidermal lysates (Supplementary Figure 1). The staining was followed by secondary antibody labeling (anti-mouse Alexa 488 and anti-rabitt Alexa 568 Life Technologies/ThermoFisher Scientific, Waltham, MA, United States), all carried out in PBS. To visualize nuclei NucBlue reagent (Hoechst, Life Technologies/ThermoFisher Scientific, Waltham, MA, United States) was used. The cover-slides were mounted with Mowiol-488 (Sigma-Aldrich, Gillingham, Dorset, United Kingdom).

Imaging acquisition was carried out on a Zeiss 780 inverted confocal microscope (Zeiss, Jena, Germany), by recording 2D images or 3D z-stacks. We used 488 nm excitation line and detected from 500–550nm. Images were postprocessed using Zen Software (Zeiss, Jena, Germany) and ImageJ National Institutes of Health, Bethesda, MD, United States).

Immunocytochemistry

Cells were grown in culture slides (BD Biosciences, San Jose, CA, United States) until 80% confluent, and were subjected to calcium switch for 24 h. The staining was conducted with anti-E-cadherin antibody (Biolegend, San Diego, CA, United States). An EnVision+ DAB polymer immunohistochemistry system (Dako) was used for both for visualization.

Intracellular Western Blot

Intracellular filaggrin and keratin-14 detection was carried out on cells cultured on 96 well plates that were used for the high throughput Laurdan screen, following data acquisition on the plate reader. Briefly, the cells were washed in PBS and fixed in 3.7% formaldehyde; 0.1% Triton X-100 was used for permeabilization and 2% BSA was used for blocking. Anti-filaggrin clone FLG01 (Abcam, Cambridge, United Kingdom) was verified with rhFLG and on epidermal lysate (by WB) as above. Anti-filaggrin FLG01 and anti-keratin 14 (clone EPR1612 from GeneTex, Irvine, California, United States) were diluted in Odyssey blocking buffer (Lincoln, NE, United States) and incubated for 2 h. After extensive wash (0.1% Tween 20) secondary antibodies (IRDye goat anti-mouse 800 and goat anti-rabbit 680; Li-Cor, Lincoln, NE, United States) were used for 1 h incubation. The data were acquired on Li-Cor Odyssey Scanner and analyzed in Image Studio software and statistical significance was determined by one-way ANOVA test with Tukey's multiple comparison.

Quantification of Membrane Order by Spectral Imaging

Spectral imaging of the different membrane samples was performed on a Zeiss LSM 780 confocal microscope equipped with a 32-channel GaAsP detector array. Excitation laser at 405 nm was used and the lambda detection range was set between 415 and 691 nm, and intervals set at 8.9 nm for the individual detection channels. This permitted the coverage of the whole spectrum with the 32 detection channels. The images were saved in .lsm file format and analyzed with a custom plug-in compatible with Fiji/ImageJ, as described later. Selection of regions of interest was done using ImageJ and the quantification of the GP index from these regions was done as explained below.

Quantification of Plasma Membrane Order in Plate Reader

Measurements of the plasma membrane stiffness were carried out on live NHEKs monolayers, employing either a fluorometer (CLARIOstar, BMG LABTECH, Ortenberg, Germany) or optical microscopy using spectral imaging on a Zeiss 780 inverted microscope. We used an environmental polarity sensitive probe, 6-dodecanoyl-2-dimethylamino naphthalene (Laurdan; Sigma Aldrich). Briefly, for the fluorimetric assessment by the plate reader, cells were seeded out on a 96-well glass-bottom plate (Greiner, Stonehouse, UK). Cell membrane labeling with Laurdan was obtained by 5 to 10 min incubation with 5 μL of a solution in dimethyl sulfoxide (DMSO, Sigma Aldrich) and Phosphate Buffer Saline pH 7.4 (PBS, Sigma Aldrich), DMSO/PBS 1:3 v/v, at a final Laurdan concentration of 0.5 μM at room temperature; afterwards, cells were washed with PBS.

Fluorescence emission of Laurdan was exciting at 374 nm and recorded over its whole spectrum from 405 to 600 nm. The intensity of emission wavelengths at 440 ± 10 nm and 490 ± 10 nm was used to obtain the GP values. The experiments were done in triplicates employing 3 different cell batches; the fluorescence values observed were an average out of minimum 25 flashes. Additionally, for imaging purposes we added phenol-free L-15 cell media (Leibovitz, Lonza). Samples were prepared and imaged on an 8-well glass-bottom chamber #1.5 (Ibidi, Planegg/Martinsried, Germany).

Generalized Polarization Index (GP)

Calculation of GP value was carried out as:

$\begin{array}{l}GP=\frac{I_{440}-I_{490}}{I_{440}+I_{490}}\end{array}$

where I₄₄₀ and I₄₉₀ correspond to the emission intensities of Laurdan at 440 and 490 nm, respectively, using 380 nm excitation wavelength. Values of GP vary from 1 to −1, where higher numbers reflect lower fluidity or higher lateral lipid order, whereas lower numbers indicate increasing fluidity. The GP measurements on live keratinocyte monolayers at different calcium concentrations (from [Ca2+] = 0.06 to 5 mM) were on one hand performed on a micro plate fluorescence reader (CLARIOstar, BMG LABTECH), with Laurdan fluorescence emission excited at 374 nm and recorded from 405 to 600 nm. The emission intensities at 440 and 490 ± 10 nm were used to obtain the GP values according to the above equation. On the other hand, confocal spectral imaging on live NHEKs at calcium concentrations [Ca2+] = 0.06 and 1.5 mM was performed on a Zeiss LSM 780 confocal microscope equipped with a 32-channel GaAsP detector array. Fluorescence of Laurdan was excited at 405 nm and detected between 415 and 691 nm. The images were then analyzed using a custom plug-in compatible with Fiji/ImageJ, as previously described 2 using a gamma variate fit of the spectra. A frequency histogram of the GP values was generated in Origin Pro (Oregon, United States), which disclosed two populations, one with high GP values representing the plasma membrane fraction and one with low GP values revealing the cytosolic fraction (e.g., from organelles). A fit of a double Gaussian distribution to the distribution allowed determining average GP values as well as standard deviations of both populations. In contrast, the plate reader measurements gave an average value over both populations. For both plate reader and spectral imaging n = 3 biological replicate were performed, while from each keratinocyte monolayer a minimum of 25 plate reader recordings or 15 confocal images (spectral images) were acquired and analyzed.

For the calculation of the GP index employing spectral imaging and rendering a GP map, we used a purposed-coded plug-in published elsewhere (Sezgin et al., 2015) and available at https://github.com/dwaithe/GP-plugin). Briefly, the plugin has in-built the GP general polarization formula and the discrete intensity values around the peak maxima (i.e., 440, and 490 nm) are obtained from a fitting to the spectrum obtained per pixel and extracting the values around these peak maxima from the nearest wavelength intervals. This methodology produces a spatial GP map representing the GP value for each pixel of the image. The outputs from these pixels are then exported into OriginPro (OriginLab, Oregon, United States) to produce the normalized frequency counts histograms per GP value. Every image pseudo-color representation has a corresponding custom look-up-table (LUT) matching the highest and lowest GP value obtained from the data.

Results

Assessment of Membrane Order in Live Human Epidermal Keratinocytes at a Single-Cell Level

Previously, we extensively characterized the coordinated role of actin scaffolds and filaggrin granule formation during differentiation of normal human epidermal keratinocytes (NHEKs) (Gutowska-Owsiak et al., 2018). Specifically, we monitored filaggrin and actin at different calcium concentrations and described and quantified the distribution of these supramolecular assemblies in cornification. We noted that membrane molecular packing was following an interesting condensation pattern. Therefore, we decided to further investigate this and develop a fast, robust, and easy method to assess cornification development by quantifying lipid order in NHEKs. For this purpose, we used a “calcium switch” model of keratinocyte differentiation, whereby differentiation is induced by addition of calcium (Bikle et al., 2012), reflecting physiological changes observed during epidermal stratification (Menon et al., 1985). In this model a dramatic increase of filaggrin-positive granule accumulation in the cytosol and around the nucleus with higher calcium concentrations can be observed (Figure 2A). Furthermore, we observe expected changes in additional markers of keratinocyte differentiation (keratin 14 expression and formation of cellular junction), assessed by E-cadherin staining localization (Owens et al., 2000; Charest et al., 2009) as shown in Figures 2B,C, respectively. Combined, all these hallmark features span through the entire thickness of a live epidermis, i.e., represent early, intermediate and late differentiation stages of keratinocyte differentiation (keratin-14, adherens junction formation, and filaggrin, respectively) validate the assay as a suitable model to investigate keratinocyte differentiation. In order to characterize the lipid lateral packing properties in NHEK membranes, we directly labeled undifferentiated cells grown at low calcium level ([Ca2+] = 0.06 mM) with Laurdan. Next, we carried out spectral imaging by confocal microscopy to address changes in membrane stiffness reported by the quantification of the Laurdan general polarization (GP) index. This solvatochromic probe reports on the extent of water relaxation processes in its direct proximity in local molecular environment (Figure 2D). The additional advantage of Laurdan is that this dye only fluoresces when incorporated into a membrane environment, hence the GP index is unaffected by the minimal amount of Laurdan which is not associated to membranes. In cells, it is known that the overall internal organelle membrane composition differs from the plasma membrane (Bernardino de la Serna et al., 2016). Plasma membrane lipids have been reported to be more tightly packed and this is clearly demonstrated in Figure 2E GP pseudo-color map, where the plasma membrane display a redder color, indicative of a higher degree of condensation or more densely packed membranes in comparison to the intracellular membranes. Since our aim for this study was to develop a high throughput method to report on keratinocyte differentiation, we wanted to ensure that the overall quantification of the GP index would unambiguously report on the degree of cellular lipid lateral condensation, in a similar fashion to the quantification performed at the plasma membrane only. For this purpose, we compared the GP index population distribution and gaussian fit of a whole cell (Figure 2E-1) with the GP of the internal membranes (Figure 2E-2) and with the GP of the plasma membrane (Figure 2E-3). The distribution of GP values of the whole cell yielded two distinct gaussians: one with lower GP values, mostly distributed in the internal organelles and reflecting more fluid membranes (see Figure 2E in blue-green color; GP = −0.09); and the other with higher values, mainly localized at the plasma membrane, indicating higher degree of lateral packing (see Figure 2E in yellow-orange color; GP = 0.26). The values for the lower and higher GPs obtained from the whole cell, did not substantially differ from the lower GP values of the cytosolic membrane component and the higher of the plasma membrane component. This indicates that resolving the GP index distribution in cell membranes is sufficient to successfully characterize the lipid order heterogeneity.

FIGURE 2

Figure 2. Assessment of the membrane order and stiffness in live normal human epithelial keratinocytes (NHEKs) by LAURDAN polarization sensing dye. (A) Detection of filaggrin-containing keratohyalin granules (KHGs) induced in primary normal human epidermal keratinocytes upon calcium switch by immunofluorescence ([Ca2+] = 5 mM to [Ca2+] = 2.5 mM); where Chromatin stained with DAPI is shown in Blue, and Filaggrin in green; scale bar 10 μm; (B) expression of keratin-14 and (C) E-cadherin staining to assess the formation of cellular junctions during calcium-induced differentiation; (D) A cartoon representing Laurdan chemical structure, its positioning within membranes and a description of its biophysical and fluorescence spectral properties. Loosely packed lipids and high water content promote unrestricted Laurdan mobility, resulting in dye emission shift representing fluid membrane state while dense lipid packing and lower water penetration restrict free movement of the dye, resulting in reflecting of shorter wavelengths; differences in LAURDAN emission reflects increased order of the lipids and is detected by a shift in emission wavelength (right panel). This results in an increase in the value of General Polarization (GP). (E) Pseudo-color scale represents the increase of a gradual shift from fluid (blue) to stiff (red) membrane. Representative GP images as determined from 2D scanning confocal spectral imaging on live undifferentiated NHEKS [Ca2+] = 0.06 mM; overview (upper panel, scale bar 5 μm) and magnified (lower panel, gray-shaded area: selection for plasma membrane part) of area marked by a white box. Pseudo-color scale: blue to red shift represent increase from low to high GP values. Frequency histograms of GP pixel values from areas of images marked by the respective number 1–3: (1) averaging over a whole cell; (2) cytosolic membranes fraction; and (3) plasma membrane fraction. The histograms were fit by a double-Gaussian (black line, with peak values, and their accuracy from the fit stated) to outline the plasma membrane (right peaks) and cytosolic (left peaks) fraction. Analysis from 3 different cell batches, and a minimum of 5 cells per batch; Statistical GP values obtained were −0.15+/−0.03 and 0.29+/−0.03 at cytosolic membranes and plasma membrane, respectively; the analysis of the whole cell yielded a double peak pit centered at a low GP −0.1+/−0.04 of and at a high GP of 0.21+/−0.05.

Laurdan General Polarization Index Increases During in vitro Keratinocyte Differentiation

Next, we tested whether confocal spectral imaging would show membrane order changes in live NHEKs during keratinocyte differentiation at a single cell level in vitro. Since a very steep calcium gradient in the epidermis (Menon et al., 1985, 1992; Elias et al., 1998, 2002; Ahn et al., 1999) is thought to be the main factor promoting differentiation and stratification in vivo (Hennings et al., 1980; Dale et al., 1983; Pillai et al., 1990), this model directly relates to skin physiology. Initially, we validated our assay by simultaneous determination of the differentiation progression in the cells using proteomics (end point assessment). Here, we used the ratio between the expression of a late differentiation marker (filaggrin) and keratin 14, which is predominantly expressed by undifferentiated keratinocytes of lower epidermal layers, including in the basal layer characterized by high proliferation rates (Figure 3A). The increase of the filaggrin/keratin 14 ratio indicates keratinocytes advancing in their differentiation process.

FIGURE 3

Figure 3. Changes in the membrane stiffness during in vitro keratinocyte differentiation assessed by LAURDAN live imaging at the level of a single cell. (A) Representation of the keratinocyte differentiation advancement by assessment of filaggrin/keratin 14 ratio in intracellular western blot assay in the cells used for LAURDAN assay; representative of 2 separate NHEK batches, assay carried out in duplicates; (B) GP images (upper panels) and histograms of GP pixel values (lower panels, from whole image) for increasing calcium levels ([Ca²⁺] = 0.06 mM (left), 1.5 mM (right)). The histograms were fit by a double-Gaussian (black line, with peak values and their accuracy from the fit stated) to outline the increase in GP and thus stiffness for the plasma membrane fraction (right peaks) and an indifferent change for the cytosolic fraction (left peaks). Representative example from 3 different cell batches, and a minimum of 5 cells per batch. Statistical membrane GP values obtained were 0.21+/−0.03 and 0.34+/−0.05 for [Ca²⁺] = 0.06 and 1.5 mM, respectively; scale bar 5 μm.

Next, using spectral imaging we investigated whether quantifying the GP index in live NHEKs in our validated differentiation model would yield differences that could be confidently used to assess cornification by the means of membrane order quantification. Specifically, we exposed NHEKs to increasing calcium level for 24 h, labeled them with Laurdan and determined GP imaging maps and GP index population distributions. The GP maps in Figure 3B confirmed the overall cellular membrane state, showing a decrease in fluid membranes overall and an increase in more ordered membranes. Moreover, the histograms showed a clear and significant response of the membrane packing readout in both the low and high GP value with increased calcium concentration. For instance, from 0.06 mM to 1.5 mM Ca²⁺ the whole population distribution is shifted toward higher GP values: from −0.06 to 0.05 and from 0.26 to 0.45 for the low and high GP, respectively. Interestingly, the highest calcium concentration (5 mM) shows complete disruption of the cells and their membranes. Plasma membrane stiffness and thus GP values increased with calcium concentration throughout the tested range (0.06–5 mM). Given this, we concluded that GP values hold as a robust readout for determination of differentiation stages of keratinocytes at the single live-cell level.

High Throughput GP Index Quantification Assay to Determine Progression of Differentiation in Live Keratinocytes

Once we demonstrated that the GP index is a suitable reporter of cell membrane order without the need to obtain values exclusively from the plasma membrane, we progressed further in our aim of developing an inexpensive and robust methodology to asses keratinocyte differentiation with a high throughput readout and analysis. For this, we tested whether calculating the GP index at the cell-ensemble level, rather than at the single-cell, would yield corresponding results. Hence, we reproduced the previous experimental set up, but employing a plate reader instead of a confocal microscope. This assay provided average GP values from whole cell populations within sub-second time frames. In these experiments we determined GP values representative of all cell membranes, i.e., calculating the average over both the cytoplasmic and cytosolic membrane fractions. This new approach allowed us to test multiple calcium concentrations faster and in a larger field of view; ultimately allowing access to a more robust statistical analysis. A caveat to this approach is the fact that the readout is an average of the fluorescence spectra detected from an ensemble of keratinocytes within a cell monolayer. To further test whether this could be a set-back to our method, we ran the samples in parallel in a plate under a confocal microscope and in a plate reader (Figure 4A). Notably, using the same cell populations as imaged on the confocal microscope (the plate format fitted both the microscopy setup and the fluorescent reader) we could still observe a gradual increase of the Laurdan GP values along with the progression of keratinocyte differentiation in our calcium switch model (Figure 4B). The low error bars along with the short measurement times highlight the plate reader assay as a robust high-throughput readout to assess keratinocyte differentiation.

FIGURE 4

Figure 4. High throughput assessment of Laurdan General Polarization (GP) in living keratinocytes at the population level by a fluorescent plate readout. (A) Workflow diagram: primary keratinocytes were seeded and calcium-switched in a range of calcium concentrations ([Ca2+] = 0.06 mM to [Ca2+] = 2.5 mM) at confluence. Twenty four hour after cells were LAURDAN-stained live and GP values were obtained by confocal immunofluorescence readout and fluorescence spectroscopy on a plate reader. (B) Whiskers-box plot of the distribution of GP values, with average (middle vertical black lines), the standard deviations (error bars), and the outlier ranges (gray crosses).

Discussion

Epidermis consists of multiple layers of keratinocytes which represent different stages of progression of terminal differentiation. Successful and complete differentiation is critical for the functional characteristics of the skin, and affects their permeability, antimicrobial, and mechanical properties. While the cells gradually progress, their characteristics change dramatically; this results in differential cell reactivity and outcomes in functionally distinct keratinocyte populations. To date there are no specific tools available to trace the progression of keratinocyte differentiation in vitro in real-time or to precisely determine differentiation advancement stage of these cells while preserving cell viability; at neither single-cell nor population level. Such an assessment modality would be beneficial, however, as it would enable the determination of the dependence between the cell reactivity and functional outcomes of the advancement in their differentiation program. Furthermore, a plate reader format could provide an efficient tool for scaling up into a high-throughput assay. Similarly, given that several skin diseases have been found to show a degree of delayed or abnormal epidermal differentiation and resulting insufficiency in the skin barrier function (Kubo et al., 2012; Eyerich et al., 2018), there is an unmet need for robust and validated techniques which could be used to propose new diagnostic and therapeutic options. Our study provides an inexpensive and effective methodology to assess keratinocyte differentiation and has implications to both basic and applied research, e.g., in drug testing. However, it is important to note that while our assay may be used to determine the GP value and make assumptions with regards to the membrane stiffness across a spectrum of calcium concentrations, the cells grown at calcium levels of 5 mM or above may no longer be viable, likely due to differentiation-dependent cell death and not suitable for downstream analysis. Another important aspect to consider is the fact the Laurdan GP is not an absolute, but relative value; the values are indicative of a process of differentiation, where clearly the order in the cellular membranes increases upon higher calcium levels. As shown, absolute value differences in the analyses can be attributed to (i) the lower sensitivity of the plate reader compared to the microscope detectors; (ii) the excitation power and wavelength (in a LED vs. laser-based for the microscope); and (iii) the difference in the cellular organization (ensemble of cells, where cell-cell attachments influence the membrane fluidity, vs. single-cell visualization under the microscope).

Lipid membrane remodeling and resulting changes in the polarity can be employed to measure the degree of membrane condensation by utilizing fluorescent lipophilic probes such as Laurdan, where emission spectra and solvent-polarity dependence report on lipid lateral packing. The readout of the Laurdan fluorescence and GP index can be used to provide information regarding membrane heterogeneity, which allows for the estimation of the membrane stiffness and fluidity. Furthermore, lack of the apparent Laurdan toxicity and the sub-second readout time of the plate reader allowed us to study changes in keratinocyte stiffness and progression of differentiation over time. It is important to note that the GP index is not an absolute measure of fluidity, but an indirect and relative estimate. Therefore, GP values may differ from study to study. For instance, very precise measurements at room temperature of GP in gradually laterally compressed monolayers of the di-saturated lipid DPPC (1,2-dipalmitoyl-sn-glycero-3-phosphocholine), a lipid with a transition temperature of 41°C, lead to GP values between 0.25 and 0.4 (Brewer et al., 2010) at lateral pressures, similar to those experienced in a lipid bilayer in the liquid order phase (de la Serna et al., 2009; Bernardino de la Serna et al., 2013). However, when DPPC is compressed further to reach a gel state (e.g., a phase that has only been reported in skin membranes or model membranes) the GP value still reaches values below 0.5. Therefore, the values reported here for highly densely packed membranes of are in good agreement with the increase of the membrane order (de la Serna et al., 2009; Gutowska-Owsiak et al., 2018).

In this study we demonstrated that Laurdan general polarization (GP) index can be used as an effective tool to determine keratinocyte differentiation advancement and progression in live keratinocytes at both single-cell and population level, within a cell monolayer and with an option for high-throughput assessment. Laurdan has been previously used in experiments with porcine skin to determine penetration of liposomes (Carrer et al., 2008). Assessment of Laurdan emission in these studies showed a decrease of the GP value with the increased skin depth, from stratum corneum through to the dermis, including a progressive decrease at the level of the epidermis; this is in line with data presented here. While we utilized Laurdan to assess differentiation process in epidermal keratinocytes specifically, we anticipate that this methodology may be also potentially successfully used with other cell types, which undergo lipid remodeling during differentiation or maturation (Reynier et al., 1991; Wong et al., 2018).

In summary, we designed and validated a Laurdan-based assay which enables for a precise determination of differentiation advancement stage in living keratinocytes, by the assessment of changes in membrane stiffness. Due to the lack of apparent cytotoxicity, this in vitro assay may also be used to trace differentiation over time, either at the level of a single cell or scaled up for high throughput measurements.

Data Availability Statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Author Contributions

DG-O and JBS designed and performed the experiments, carried out the data processing and analysis, discussed, prepared, and wrote the manuscript. EIP carried out experiments. CE and GO supervised the work, helped with analysis, discussed, prepared, and wrote the manuscript. All authors contributed to the article and approved the submitted version.

Funding

DG-O and GO are grateful for the support from the Medical Research Council (UK) and the NIHR Oxford Biomedical Research Centre. DG-O also acknowledged funding from the European Union's Horizon 2020 research and innovation programme under the Marie Skłodowska-Curie grant agreement no. 665778 coordinated by the National Science Centre, under the Polonez project number 2016/23/P/NZ6/04056 and Foundation for Polish Science (FNP) under the First TEAM project number POIR.04.04.00-00-21FA/16-00, and the British Skin Foundation. EIP was supported by the Preludium (2013/09/N/NZ6/02405) and the Etiuda Programs funded by the National Science Centre of Poland. JBS acknowledges support from a Marie Curie Career Integration Grant (PCIG13-GA-2013-618914). CE acknowledges support from the Wolfson Imaging Centre Oxford (Christoffer Lagerholm), the Wolfson Foundation (18272), the Medical Research Council (MC_UU_12010/unit programmes G0902418 and MC_UU_12025), MRC/BBSRC/ESPRC (Grant No. MR/K01577X/1), the Wellcome Trust [104924/14/Z/14, 100262Z/12/Z, 098274/Z/12/Z, and Strategic Award 091911 (107457/Z/15/Z, Advanced Micron Bioimaging Unit)], the Deutsche Forschungsgemeinschaft (Collaborative Research Center 1278, Research unit 1905, Jena Excellence Cluster Balance of the Microverse), and internal University of Oxford funding (EPA Cephalosporin Fund and John Fell Fund).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

We are grateful to Prof. Wojciech Bal for his support to EIP and the study.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcell.2020.573230/full#supplementary-material

Supplementary Figure 1. Validation of anti-filaggrin antibodies used in the study. Western blot assessment of the antibody specificity for FLG01 and 15C10 antibodies used in the study. (A) Blotting against recombinant human FLG protein; (B) Blotting against epidermal lysates.

References

Ahn, S. K., Hwang, S. M., Jiang, S. J., Choi, E. H., and Lee, S. H. (1999). The changes of epidermal calcium gradient and transitional cells after prolonged occlusion following tape stripping in the murine epidermis. J. Invest. Dermatol. 113, 189–195. doi: 10.1046/j.1523-1747.1999.00650.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Akinduro, O., Sully, K., Patel, A., Robinson, D. J., Chikh, A., McPhail, G., et al. (2016). Constitutive autophagy and nucleophagy during epidermal differentiation. J. Invest. Dermatol. 136, 1460–1470. doi: 10.1016/j.jid.2016.03.016

PubMed Abstract | CrossRef Full Text | Google Scholar

Akiyama, M. (2017). Corneocyte lipid envelope (CLE), the key structure for skin barrier function and ichthyosis pathogenesis. J. Dermatol. Sci. 88, 3–9. doi: 10.1016/j.jdermsci.2017.06.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Bagatolli, L. A., Gratton, E., and Fidelio, G. D. (1998). Water dynamics in glycosphingolipid aggregates studied by LAURDAN fluorescence. Biophys. J. 75, 331–341. doi: 10.1016/S0006-3495(98)77517-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Bagatolli, L. A., Parasassi, T., Fidelio, G. D., and Gratton, E. (1999). A model for the interaction of 6-lauroyl-2-(N,N-dimethylamino)naphthalene with lipid environments: implications for spectral properties. Photochem. Photobiol. 70, 557–564. doi: 10.1562/0031-8655(1999)070<0557:AMFTIO>2.3.CO;2

PubMed Abstract | CrossRef Full Text | Google Scholar

Behne, M., Uchida, Y., Seki, T., de Montellano, P. O., Elias, P. M., and Holleran, W. M. (2000). Omega-hydroxyceramides are required for corneocyte lipid envelope (CLE) formation and normal epidermal permeability barrier function. J. Invest. Dermatol. 114, 185–192. doi: 10.1046/j.1523-1747.2000.00846.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Bernabé-Rubio, M., Bosch-Fortea, M., García, E., de la Serna, J. B., and Alonso, M. A. (2019). The ciliary membrane of polarized epithelial cells stems from a midbody remnant-associated membrane patch with condensed nanodomains. bioRxiv 667642. doi: 10.1101/667642

CrossRef Full Text | Google Scholar

Bernardino de la Serna, J., Hansen, S., Berzina, Z., Simonsen, A. C., Hannibal-Bach, H. K., Knudsen, J., et al. (2013). Compositional and structural characterization of monolayers and bilayers composed of native pulmonary surfactant from wild type mice. Biochim. Biophys. Acta 1828, 2450–2459. doi: 10.1016/j.bbamem.2013.07.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Bernardino de la Serna, J., Schütz, G. J., Eggeling, C., and Cebecauer, M. (2016). There is no simple model of the plasma membrane organization. Front Cell Dev. Biol. 4:106. doi: 10.3389/fcell.2016.00106

PubMed Abstract | CrossRef Full Text

Bikle, D. D., Xie, Z., and Tu, C.-L. (2012). Calcium regulation of keratinocyte differentiation. Expert Rev. Endocrinol. Metab. 7, 461–472. doi: 10.1586/eem.12.34

PubMed Abstract | CrossRef Full Text | Google Scholar

Brewer, J., de la Serna, J. B., Wagner, K., and Bagatolli, L. A. (2010). Multiphoton excitation fluorescence microscopy in planar membrane systems. Biochim. Biophys. Acta 1798, 1301–1308. doi: 10.1016/j.bbamem.2010.02.024

PubMed Abstract | CrossRef Full Text | Google Scholar

Candi, E., Knight, R. A., Panatta, E., Smirnov, A., and Melino, G. (2016). Cornification of the skin: a non-apoptotic cell death mechanism. eLS 1–10. doi: 10.1002/9780470015902.a0021583.pub2

CrossRef Full Text | Google Scholar

Candi, E., Schmidt, R., and Melino, G. (2005). The cornified envelope: a model of cell death in the skin. Nat. Rev. Mol. Cell Biol. 6, 328–340. doi: 10.1038/nrm1619

PubMed Abstract | CrossRef Full Text | Google Scholar

Carrer, D. C., Vermehren, C., and Bagatolli, L. A. (2008). Pig skin structure and transdermal delivery of liposomes: a two photon microscopy study. J. Control. Release 132, 12–20. doi: 10.1016/j.jconrel.2008.08.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Carugo, D., Aron, M., Sezgin, E., Bernardino de la Serna, J., Kuimova, M. K., Eggeling, C., et al. (2017). Modulation of the molecular arrangement in artificial and biological membranes by phospholipid-shelled microbubbles. Biomaterials 113, 105–117. doi: 10.1016/j.biomaterials.2016.10.034

PubMed Abstract | CrossRef Full Text | Google Scholar

Charest, J. L., Jennings, J. M., King, W. P., Kowalczyk, A. P., and García, A. J. (2009). Cadherin-mediated cell–cell contact regulates keratinocyte differentiation. J. Invest. Dermatol. 129, 564–572. doi: 10.1038/jid.2008.265

PubMed Abstract | CrossRef Full Text | Google Scholar

Chon, S. H., Tannahill, R., Yao, X., Southall, M. D., and Pappas, A. (2015). Keratinocyte differentiation and upregulation of ceramide synthesis induced by an oat lipid extract via the activation of PPAR pathways. Exp. Dermatol. 24, 290–295. doi: 10.1111/exd.12658

PubMed Abstract | CrossRef Full Text | Google Scholar

Dale, B. A., Haugen Scofield, J. A., Hennings, H., Stanley, J. R., and Yuspa, S. H. (1983). Identification of filaggrin in cultured mouse keratinocytes and its regulation by calcium. J. Invest. Dermatol. 81 (1 Suppl), S90–S95. doi: 10.1111/1523-1747.ep12540769

PubMed Abstract | CrossRef Full Text | Google Scholar

de la Serna, J. B., Orädd, G., Bagatolli, L. A., Simonsen, A. C., Marsh, D., Lindblom, G., et al. (2009). Segregated phases in pulmonary surfactant membranes do not show coexistence of lipid populations with differentiated dynamic properties. Biophys. J. 97, 1381–1389. doi: 10.1016/j.bpj.2009.06.040

PubMed Abstract | CrossRef Full Text | Google Scholar

Eckhart, L., Lippens, S., Tschachler, E., and Declercq, W. (2013). Cell death by cornification. Biochim. Biophys. Acta 1833, 3471–3480. doi: 10.1016/j.bbamcr.2013.06.010

PubMed Abstract | CrossRef Full Text | Google Scholar

Elias, P. M., Brown, B. E., Crumrine, D., Feingold, K. R., and Ahn, S. K. (2002). Origin of the epidermal calcium gradient: regulation by barrier status and role of active vs passive mechanisms. J. Invest. Dermatol. 119, 1269–1274. doi: 10.1046/j.1523-1747.2002.19622.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Elias, P. M., Gruber, R., Crumrine, D., Menon, G., Williams, M. L., Wakefield, J. S., et al. (2014). Formation and functions of the corneocyte lipid envelope (CLE). Biochim. Biophys. Acta 1841, 314–318. doi: 10.1016/j.bbalip.2013.09.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Elias, P. M., Nau, P., Hanley, K., Cullander, C., Crumrine, D., Bench, G., et al. (1998). Formation of the epidermal calcium gradient coincides with key milestones of barrier ontogenesis in the rodent. J. Invest. Dermatol. 110, 399–404. doi: 10.1046/j.1523-1747.1998.00151.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Eyerich, S., Eyerich, K., Traidl-Hoffmann, C., and Biedermann, T. (2018). Cutaneous barriers and skin immunity: differentiating a connected network. Trends Immunol. 39, 315–327. doi: 10.1016/j.it.2018.02.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Flohr, C., Perkin, M., Logan, K., Marrs, T., Radulovic, S., Campbell, L. E., et al. (2014). Atopic dermatitis and disease severity are the main risk factors for food sensitization in exclusively breastfed infants. J. Invest. Dermatol. 134, 345–350. doi: 10.1038/jid.2013.298

PubMed Abstract | CrossRef Full Text | Google Scholar

Gaus, K., Gratton, E., Kable, E. P. W., Jones, A. S., Gelissen, I., Kritharides, L., et al. (2003). Visualizing lipid structure and raft domains in living cells with two-photon microscopy. Proc. Natl. Acad. Sci. U.S.A. 100, 15554–15559. doi: 10.1073/pnas.2534386100

PubMed Abstract | CrossRef Full Text | Google Scholar

Gdula, M. R., Poterlowicz, K., Mardaryev, A. N., Sharov, A. A., Peng, Y., Fessing, M. Y., et al. (2013). Remodeling of three-dimensional organization of the nucleus during terminal keratinocyte differentiation in the epidermis. J. Invest. Dermatol. 133, 2191–2201. doi: 10.1038/jid.2013.66

PubMed Abstract | CrossRef Full Text | Google Scholar

Gutowska-Owsiak, D., de La Serna, J. B., Fritzsche, M., Naeem, A., Podobas, E. I., Leeming, M., et al. (2018). Orchestrated control of filaggrin–actin scaffolds underpins cornification. Cell Death Dis. 9:412. doi: 10.1038/s41419-018-0407-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Hennings, H., Michael, D., Cheng, C., Steinert, P., Holbrook, K., and Yuspa, S. H. (1980). Calcium regulation of growth and differentiation of mouse epidermal cells in culture. Cell 19, 245–254. doi: 10.1016/0092-8674(80)90406-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Horimukai, K., Morita, K., Narita, M., Kondo, M., Kitazawa, H., Nozaki, M., et al. (2014). Application of moisturizer to neonates prevents development of atopic dermatitis. J. Allergy Clin. Immunol. 134, 824–830.e826. doi: 10.1016/j.jaci.2014.07.060

PubMed Abstract | CrossRef Full Text | Google Scholar

Ipponjima, S., Umino, Y., Nagayama, M., and Denda, M. (2020). Live imaging of alterations in cellular morphology and organelles during cornification using an epidermal equivalent model. Sci. Rep. 10:5515. doi: 10.1038/s41598-020-62240-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Iwai, I., Han, H., den Hollander, L., Svensson, S., Ofverstedt, L. G., Anwar, J., et al. (2012). The human skin barrier is organized as stacked bilayers of fully extended ceramides with cholesterol molecules associated with the ceramide sphingoid moiety. J. Invest. Dermatol. 132, 2215–2225. doi: 10.1038/jid.2012.43

PubMed Abstract | CrossRef Full Text | Google Scholar

Janssens, M., van Smeden, J., Gooris, G. S., Bras, W., Portale, G., Caspers, P. J., et al. (2012). Increase in short-chain ceramides correlates with an altered lipid organization and decreased barrier function in atopic eczema patients. J. Lipid Res. 53, 2755–2766. doi: 10.1194/jlr.P030338

PubMed Abstract | CrossRef Full Text | Google Scholar

Kezic, S., Kammeyer, A., Calkoen, F., Fluhr, J. W., and Bos, J. D. (2009). Natural moisturizing factor components in the stratum corneum as biomarkers of filaggrin genotype: evaluation of minimally invasive methods. Br. J. Dermatol. 161, 1098–1104. doi: 10.1111/j.1365-2133.2009.09342.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Kubo, A., Nagao, K., and Amagai, M. (2012). Epidermal barrier dysfunction and cutaneous sensitization in atopic diseases. J. Clin. Invest. 122, 440–447. doi: 10.1172/JCI57416

PubMed Abstract | CrossRef Full Text | Google Scholar

Lakowicz, J. R. (1999). Principles of Fluorescence Spectroscopy. New York, NY: Plenum Publishers. Available online at: https://www.springer.com/gp/book/9780387312781#aboutBook

Google Scholar

López, O., Cócera, M., Wertz, P. W., López-Iglesias, C., and de la Maza, A. (2007). New arrangement of proteins and lipids in the stratum corneum cornified envelope. Biochim. Biophys. Acta. 1768, 521–529. doi: 10.1016/j.bbamem.2006.11.023

PubMed Abstract | CrossRef Full Text | Google Scholar

Mao-Qiang, M., Fowler, A. J., Schmuth, M., Lau, P., Chang, S., Brown, B. E., et al. (2004). Peroxisome-proliferator-activated receptor (PPAR)-gamma activation stimulates keratinocyte differentiation. J. Invest. Dermatol. 123, 305–312. doi: 10.1111/j.0022-202X.2004.23235.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Marshall, D., Hardman, M. J., Nield, K. M., and Byrne, C. (2001). Differentially expressed late constituents of the epidermal cornified envelope. Proc. Natl. Acad. Sci. U.S.A. 98, 13031–13036. doi: 10.1073/pnas.231489198

PubMed Abstract | CrossRef Full Text | Google Scholar

Menon, G. K., Elias, P. M., Lee, S. H., and Feingold, K. R. (1992). Localization of calcium in murine epidermis following disruption and repair of the permeability barrier. Cell Tissue Res. 270, 503–512. doi: 10.1007/BF00645052

PubMed Abstract | CrossRef Full Text | Google Scholar

Menon, G. K., Grayson, S., and Elias, P. M. (1985). Ionic calcium reservoirs in mammalian epidermis: ultrastructural localization by ion-capture cytochemistry. J. Invest. Dermatol. 84, 508–512. doi: 10.1111/1523-1747.ep12273485

PubMed Abstract | CrossRef Full Text | Google Scholar

Narangifard, A., den Hollander, L., Wennberg, C. L., Lundborg, M., Lindahl, E., Iwai, I., et al. (2018). Human skin barrier formation takes place via a cubic to lamellar lipid phase transition as analyzed by cryo-electron microscopy and EM-simulation. Exp. Cell Res. 366, 139–151. doi: 10.1016/j.yexcr.2018.03.010

PubMed Abstract | CrossRef Full Text | Google Scholar

Owen, D. M., Rentero, C., Magenau, A., Abu-Siniyeh, A., and Gaus, K. (2011). Quantitative imaging of membrane lipid order in cells and organisms. Nat. Protoc. 7, 24–35. doi: 10.1038/nprot.2011.419

PubMed Abstract | CrossRef Full Text | Google Scholar

Owens, D. W., Brunton, V. G., Parkinson, E. K., and Frame, M. C. (2000). E-cadherin at the cell periphery is a determinant of keratinocyte differentiation in vitro. Biochem. Biophys. Res. Commun. 269, 369–376. doi: 10.1006/bbrc.2000.2292

PubMed Abstract | CrossRef Full Text | Google Scholar

Parasassi, T., De Stasio, G., d'Ubaldo, A., and Gratton, E. (1990). Phase fluctuation in phospholipid membranes revealed by Laurdan fluorescence. Biophys. J. 57, 1179–1186. doi: 10.1016/S0006-3495(90)82637-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Parasassi, T., Di Stefano, M., Loiero, M., Ravagnan, G., and Gratton, E. (1994). Cholesterol modifies water concentration and dynamics in phospholipid bilayers: a fluorescence study using laurdan probe. Biophys. J. 66 (3 Pt. 1), 763–768. doi: 10.1016/S0006-3495(94)80852-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Parasassi, T., Gratton, E., Yu, W. M., Wilson, P., and Levi, M. (1997). Two-photon fluorescence microscopy of laurdan generalized polarization domains in model and natural membranes. Biophys. J. 72, 2413–2429. doi: 10.1016/S0006-3495(97)78887-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Pillai, S., Bikle, D. D., Mancianti, M. L., Cline, P., and Hincenbergs, M. (1990). Calcium regulation of growth and differentiation of normal human keratinocytes: modulation of differentiation competence by stages of growth and extracellular calcium. J. Cell. Physiol. 143, 294–302. doi: 10.1002/jcp.1041430213

PubMed Abstract | CrossRef Full Text | Google Scholar

Plasencia, I., Norlén, L., and Bagatolli, L. A. (2007). Direct visualization of lipid domains in human skin stratum corneum's lipid membranes: effect of pH and temperature. Biophys. J. 93, 3142–3155. doi: 10.1529/biophysj.106.096164

PubMed Abstract | CrossRef Full Text | Google Scholar

Ponec, M., Kempenaar, J., Weerheim, A., and Boonstra, J. (1987). Differentiation of human keratinocytes: changes in lipid synthesis, plasma membrane lipid composition, and 125I-EGF binding upon administration of 25-hydroxycholesterol and mevinolin. J. Cell. Physiol. 133, 358–364. doi: 10.1002/jcp.1041330221

PubMed Abstract | CrossRef Full Text | Google Scholar

Ponec, M., Weerheim, A., Kempenaar, J., Mommaas, A. M., and Nugteren, D. H. (1988). Lipid composition of cultured human keratinocytes in relation to their differentiation. J. Lipid Res. 29, 949–961.

PubMed Abstract | Google Scholar

Reynier, M., Sari, H., d'Anglebermes, M., Kye, E. A., and Pasero, L. (1991). Differences in lipid characteristics of undifferentiated and enterocytic-differentiated HT29 human colonic cells. Cancer Res. 51, 1270–1277.

PubMed Abstract | Google Scholar

Rogerson, C., Bergamaschi, D., and O'Shaughnessy, R. F. L. (2018). Uncovering mechanisms of nuclear degradation in keratinocytes: a paradigm for nuclear degradation in other tissues. Nucleus 9, 56–64. doi: 10.1080/19491034.2017.1412027

PubMed Abstract | CrossRef Full Text | Google Scholar

Sandilands, A., Sutherland, C., Irvine, A. D., and McLean, W. H. (2009). Filaggrin in the frontline: role in skin barrier function and disease. J. Cell Sci. 122 (Pt. 9), 1285–1294. doi: 10.1242/jcs.033969

PubMed Abstract | CrossRef Full Text | Google Scholar

Schmidt, R., Parish, E. J., Dionisius, V., Cathelineau, C., Michel, S., Shroot, B., et al. (1991). Modulation of cellular cholesterol and its effect on cornified envelope formation in cultured human epidermal keratinocytes. J. Invest. Dermatol. 97, 771–775. doi: 10.1111/1523-1747.ep12486720

PubMed Abstract | CrossRef Full Text | Google Scholar

Sezgin, E., Waithe, D., Bernardino de la Serna, J., and Eggeling, C. (2015). Spectral imaging to measure heterogeneity in membrane lipid packing. Chemphyschem 16, 1387–1394. doi: 10.1002/cphc.201402794

PubMed Abstract | CrossRef Full Text | Google Scholar

Venkataraman, D., Soto-Ramírez, N., Kurukulaaratchy, R. J., Holloway, J. W., Karmaus, W., Ewart, S. L., et al. (2014). Filaggrin loss-of-function mutations are associated with food allergy in childhood and adolescence. J. Allergy Clin. Immunol. 134, 876–882.e874. doi: 10.1016/j.jaci.2014.07.033

PubMed Abstract | CrossRef Full Text | Google Scholar

Weber, G., and Farris, F. J. (1979). Synthesis and spectral properties of a hydrophobic fluorescent probe: 6-propionyl-2-(dimethylamino)naphthalene. Biochemistry. 14, 3075–3078. doi: 10.1021/bi00581a025

PubMed Abstract | CrossRef Full Text

Wennberg, C. L., Narangifard, A., Lundborg, M., Norlén, L., and Lindahl, E. (2018). Structural transitions in ceramide cubic phases during formation of the human skin barrier. Biophys. J. 114, 1116–1127. doi: 10.1016/j.bpj.2017.12.039

PubMed Abstract | CrossRef Full Text | Google Scholar

Wong, M., Xu, G., Park, D., Barboza, M., and Lebrilla, C. B. (2018). Intact glycosphingolipidomic analysis of the cell membrane during differentiation yields extensive glycan and lipid changes. Sci. Rep. 8:10993. doi: 10.1038/s41598-018-29324-7

PubMed Abstract | CrossRef Full Text | Google Scholar