Figures
Abstract
Purpose
Storage of cultured human oral keratinocytes (HOK) allows for transportation of cultured transplants to eye clinics worldwide. In a previous study, one-week storage of cultured HOK was found to be superior with regard to viability and morphology at 12°C compared to 4°C and 37°C. To understand more of how storage temperature affects cell phenotype, gene expression of HOK before and after storage at 4°C, 12°C, and 37°C was assessed.
Materials and Methods
Cultured HOK were stored in HEPES- and sodium bicarbonate-buffered Minimum Essential Medium at 4°C, 12°C, and 37°C for one week. Total RNA was isolated and the gene expression profile was determined using DNA microarrays and analyzed with Partek Genomics Suite software and Ingenuity Pathway Analysis. Differentially expressed genes (fold change > 1.5 and P < 0.05) were identified by one-way ANOVA. Key genes were validated using qPCR.
Results
Gene expression of cultures stored at 4°C and 12°C clustered close to the unstored control cultures. Cultures stored at 37°C displayed substantial change in gene expression compared to the other groups. In comparison with 12°C, 2,981 genes were differentially expressed at 37°C. In contrast, only 67 genes were differentially expressed between the unstored control and the cells stored at 12°C. The 12°C and 37°C culture groups differed most significantly with regard to the expression of differentiation markers. The Hedgehog signaling pathway was significantly downregulated at 37°C compared to 12°C.
Conclusion
HOK cultures stored at 37°C showed considerably larger changes in gene expression compared to unstored cells than cultured HOK stored at 4°C and 12°C. The changes observed at 37°C consisted of differentiation of the cells towards a squamous epithelium-specific phenotype. Storing cultured ocular surface transplants at 37°C is therefore not recommended. This is particularly interesting as 37°C is the standard incubation temperature used for cell culture.
Citation: Utheim TP, Islam R, Fostad IG, Eidet JR, Sehic A, Olstad OK, et al. (2016) Storage Temperature Alters the Expression of Differentiation-Related Genes in Cultured Oral Keratinocytes. PLoS ONE 11(3): e0152526. https://doi.org/10.1371/journal.pone.0152526
Editor: Andrzej T. Slominski, University of Alabama at Birmingham, UNITED STATES
Received: December 19, 2015; Accepted: March 15, 2016; Published: March 29, 2016
Copyright: © 2016 Utheim et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Data Availability: All relevant data are within the paper. Supporting Information files are not provided.
Funding: Funding was received from Department of Oral Biology, Faculty of Dentistry, University of Oslo and Department of Medical Biochemistry, Oslo University Hospital, Oslo, Norway.
Competing interests: The authors have declared that no competing interests exist.
Introduction
The stem cells of the cornea are located in the periphery, in a region known as the limbus. Limbal stem cells can be destroyed by a multitude of diseases, including certain autoimmune diseases and genetic conditions [1]. These cells can also be damaged by external factors, such as chemical or thermal burns, ultraviolet radiation, and infections (e.g. trachoma). Contingent upon the extent of damage to limbal stem cells, various clinical presentations of limbal stem cell deficiency (LSCD) may develop. In the most serious cases, patients may become blind and experience substantial pain.
In 1997, LSCD was for the first time successfully treated by transplantation of ex vivo cultured limbal stem cells [2]. In unilateral LSCD, autologous limbal stem cells can be harvested from the contralateral healthy cornea, but this is generally not feasible in bilateral LSCD, which is by far the most common form. If allogeneic limbal stem cells are applied, immunosuppression, which can have severe adverse effects [3], is required at least for a certain period of time [1]. This has urged researchers to the search for alternative autologous cell sources.
In 2004, oral keratinocytes were shown to be effective for treating LSCD in humans [4, 5]. Since then, there have been 20 clinical reports confirming their potential to treat LSCD [6]. Except for conjunctival cells, oral keratinocytes are the only non-limbal cell type that has been used clinically. Accumulating evidence of the rationale for transplanting cultured oral keratinocytes in LSCD substantiates the need to make this regenerative medicine technology available worldwide. Currently, the treatment is restricted to a few centers of expertise [6]. Increasingly stricter regulations for cell therapy will likely lead to the centralization of culture units [7]. Centralization requires effective transportation strategies [8], which calls for a practical method for storage of cultured cells outside the incubator (Fig 1).
An oral mucosa biopsy is removed from the mouth (A) and sent to a laboratory (B, C, D). Oral keratinocytes are then cultured in an incubator for six days before the generated cell sheet is transferred to a sealed storage container where it can be preserved for up to one week. This allows the cultured tissue to be returned to the patient (E) for transplantation onto the diseased eye (F). Courtesy of Amer Sehic, Department of Oral Biology, University of Oslo.
Storage in a small sealed container for some days offers a number of advantages. These include: 1) sufficient time for phenotypic assessment of the cultured transplants prior to surgery, which has become increasingly important with recent knowledge about the critical role of the phenotype of transplants for good clinical outcome [9]; 2) microbiological assessment after aspiration of a storage medium sample from the septum of the hermetically sealed storage container; 3) increased flexibility for the surgeon in scheduling surgeries, which may be convenient if unforeseen factors related to the patient or cultured cells should occur [10, 11], and importantly; 4) transportation of transplants to reach eye clinics worldwide.
In a previous study, one-week storage of cultured human oral keratinocytes (HOK) at 12°C was superior with regard to viability and morphology compared to storage at 4°C and 37°C [12]. In the present study, we have used genome-wide analysis of gene expression to: 1) investigate whether differences in temperature following one-week storage result in phenotypic changes, and 2) explore potential mechanisms behind these differences. In the previous study, phenotype was analysed by immunocytochemistry and found to be preserved at both 4°C and 12°C, but not at 37°C [12]. Based on this finding, we hypothesize that cells stored at 37°C show substantial differences in gene expression compared to cells stored at 12°C.
Materials and Methods
First passage normal HOK, oral keratinocyte medium (OKM), oral keratinocyte growth supplement (OKGS), and penicillin/streptomycin solution (P/S) were purchased from ScienCell Research Laboratories (San Diego, Carlsbad, CA, USA). Nunclon Δ-Cell culture flasks and pipettes were purchased from VWR International (West Chester, PA, USA). The Minimum Essential Medium (MEM) was obtained from Invitrogen (Carlsbad, CA, USA). Phosphate-buffered saline (PBS), 4-(2-hydroxyethyl)-1-piperazineethanesulfonic acid (HEPES), and sodium bicarbonate were all purchased from Sigma-Aldrich (St. Louis, MO). RNeasy Plus Mini Kit and the QIAzol Lysis Reagent were provided by Qiagen (Hilden, Germany).
Culture and Storage of Human Oral Keratinocytes
Human oral keratinocytes were grown to confluence in T25 cell culture flasks in complete OKM (made by adding 5 mL OKGS and 5 mL P/S to 500 mL OKM), in a 37°C humidified incubator with 5% CO2 supply. The HOK were cultured in the dark, and the culture medium was changed every other day. All cells were cultured for six days. Control cells were immediately processed for analysis, while the rest were randomized to storage at either 4°C, 12°C or 37°C. These were stored for one week before being analyzed.
On day six of culture, when confluent cultures were obtained, the OKM was removed and the cultures were rinsed with PBS before adding the storage medium. The storage medium consisted of 70 mL MEM, 25 mM HEPES, 600 mg/L sodium bicarbonate, and 50 μg/mL gentamycin (hereafter named MEM). The screw caps of the T-25 flasks were tightened to reduce air exchange and evaporation. The cultures were randomized for storage at three temperatures (4°C, 12°C, and 37°C) for one week. Cells cultured for six days, but not subjected to storage, served as controls in all experiments. Custom-built storage cabinets with a very small standard deviation (±0.4°C) for the set temperatures were used for regulating temperature during storage [13]. The temperature inside each storage container was monitored throughout all experiments. The storage cabinets were kept in a cold room maintaining an ambient temperature of 4°C. Each cabinet was equipped with a light bulb functioning as a heater, which increased the temperature inside the box from the ambient room temperature (4°C) to the desired storage temperature. The light bulbs were continuously regulated by a highly sensitive thermometer, and the storage containers were equipped with a small fan that ensured a homogeneous temperature inside the box by circulating the air. The light bulbs were separated from the cells by dark walls, which ensured that the cells were not directly exposed to light, and minimized indirect light exposure. However, we cannot exclude the possibility that cells stored at higher temperatures (12°C and 37°C) were to some extent exposed to the light.
Isolation of RNA
Cultured HOK stored for one week at 4°C, 12°C, and 37°C, and control cultures that had not been subjected to storage, were rinsed with PBS and directly lysed with QIAzol Lysis Reagent. According to the manufacturer’s protocol, the fractions of total RNA were isolated using miRNeasy Mini Kit (Qiagen). The concentrations of purified RNA were assayed using a Nanodrop ND-1000 Spectrophotometer (Thermo Fisher Scientific, Wilmington, DE, USA). This yielded RNA fractions exhibiting absorbance ratios—A260/280 and A260/230 –of at least 1.8 and 2.0, respectively. The quality of RNA in solutions was assessed using the Agilent-Bioanalyzer 2100 System and RNA 6000 Nano Assay (Agilent Technologies, Santa Clara, CA, USA). All solutions used had RNA integrity number (RIN) values of > 8.5.
Microarray Analysis
The Affymetrix GeneChip Human Gene 1.0 ST Microarrays (Affymetrix, Santa Clara, CA, USA) used in this study included approximately 28,000 gene transcripts. Microarray analysis was carried out in triplicate using cultured HOK stored for one week at 4°C, 12°C, and 37°C. Unstored control cultures served as control. Preparation of complementary DNA (cDNA) was carried out using GeneChip HT One-Cycle cDNA Synthesis Kit (Affymetrix). Each of three microarrays was hybridized with cDNA prepared from 150 ng of total RNA from each resulting solution. Biotinylated and fragmented single stranded cDNAs were hybridized to the GeneChips. The arrays were washed and stained using FS-450 fluidics station (Affymetrix).
Signal intensities were detected by Hewlett Packard Gene Array Scanner 3000 7G (Hewlett Packard, Palo Alto, CA, USA). The scanned images were processed using the AGCC (Affymetrix GeneChip Command Console) software and the CEL files were imported into Partek Genomics Suite software (Partek, Inc. MO, USA). The Robust Multichip Analysis (RMA) algorithm was applied to generate signal values and normalization. Gene transcripts with maximal signal values of less than 32 across all arrays were removed to filter for low and non-expressed genes, reducing the number of gene transcripts to 17,684. For expression comparisons of different groups, profiles were compared using a one-way ANOVA model. The results were expressed as fold changes (FC) with corresponding P-values.
Bioinformatic Analysis
Bioinformatic analysis using Ingenuity Pathways Analysis (IPA) (Ingenuity Inc, IL) was carried out to find molecular and cellular functions and canonical pathways that were significantly associated with differentially expressed genes. Briefly, the data set containing gene identifiers and corresponding fold changes and P-values was uploaded onto the web-delivered application and each gene identifier was mapped to its corresponding gene object in the Ingenuity Pathways Knowledge Base (IPKB). Functional analysis identified the biological functions and/or diseases that were significantly associated with the data sets. Fisher’s exact test was performed to calculate a P-value determining the probability that each biological function and/or disease assigned to the data set was due to chance alone. The data sets were mined for significant pathways with the IPA library of canonical pathways, using IPA generated networks as graphical representations of the molecular relationships between genes and gene products.
Validation of Microarray Results by Quantitative Real-Time PCR
The differential gene expression data were validated for selected transcripts using TaqMan® Gene Expression Assays and the Applied Biosystems®ViiA™ 7 Real-Time PCR system (Applied Biosystems, Life Technologies, Carlsbad, CA, USA). The genes encoding heat shock 22kDa protein 8 (HSPB8), tumor protein p63 (TP63), and keratin 10 (KRT10) were selected for validation. Briefly, 200 ng of total RNA was reverse transcribed using qScript™ cDNA Super Mix (Quanta Biosciences Gaithersburg, MD) following the manufacturer’s instructions. After completion of cDNA synthesis, 1/10th of the first strand reaction was used for PCR amplification. A total amount of 9 μl of diluted cDNA (diluted in H2O), 1 μl of selected primer/ probes TaqMan® Gene Expression Assays (Life Technologies), and 10 μl TaqMan® Universal Master Mix (Life Technologies) were used, as per the manufacturer’s instructions. Transuducin-like enhancer of split 1 (TLE1) was used as an endogenous control due to the low coefficient of variation (CV) (0.444) in the Affymetrix study. Each gene was run in duplicates. TaqMan® Gene Expression Assays (Life Technologies) used assays detecting HSPB8 (HSPB8-Hs00205056_m1), TP63 (TP63-Hs00978343_m1), KRT10 (KRT10-Hs01043114_g1), and TLE1 (Hs00270768_m1).
P-values were calculated using Student's t-test in Microsoft Excel (Redmond, WA, USA) using delta Ct-values. Normalized target gene expression levels were calculated using the formula: 2^(–ΔΔCt).
Results
Global Perspective of Microarray Results
Gene expression of cultures stored at 4°C and 12°C were clustered close to those of fresh cultures that had not been subjected to storage (control group) (Fig 2). Cultures stored at 37°C displayed substantial change in gene expression compared to the other groups (Fig 3). In comparison with 12°C, 2,981 genes were differentially expressed at 37°C (Table 1). In contrast, only 67 and 117 genes were differentially expressed when comparing the 12°C group to the control and the 4°C group, respectively. While only 67 genes were differentially regulated at 12°C compared to control cells, almost twice as many (124) were differentially regulated at 4°C compared to the control. Given the relatively small differences between the control, 4°C and 12°C but large differences between 12°C and 37°C, we have chosen to direct our focus primarily on the differential gene expression between 12°C and 37°C.
In contrast, gene expression of cultures stored at 37°C (green) showed a distant clustering compared to the other experimental groups.
The Most Differentially Regulated Genes
Repetin (RPTN) was the most differentially regulated gene, with a 136.9-fold upregulation at 37°C compared to 12°C. Repetin is a constituent of the epidermal differentiation complex and functions in the cornified cell envelope formation [14]. Desmoglein (DSG1) was the second most upregulated gene at 37°C compared to 12°C, with a 94.6-fold upregulation at this temperature. It is a constituent of the desmosome, providing cell-cell adhesion [15]. Keratinocyte differentiation-associated protein (KDAP) was upregulated 54.7-fold at 37°C compared to 12°C. It is a regulatory protein of keratinocyte differentiation and influences the stratification of epithelia [16, 17]. Keratin 10 (KRT10), involved in the differentiation of human oral keratinocytes [18], was upregulated 45.6-fold at 37°C compared to 12°C. Lipase K (LIPK), a gene which is highly specific for the last step of keratinocyte differentiation [19], was upregulated 43.9-fold at 37°C compared to 12°C. Cornulin (CRNN), another marker of late epidermal differentiation [20], was upregulated 43.2-fold at 37°C compared to 12°C (Table 2). Hence, the most differentially regulated proteins at 37°C compared to 12°C are directly associated with differentiation of epithelia.
Melatonin receptor 1A (MTNR1A) was significantly upregulated in the 12°C storage group compared to all other groups: 5.4-fold compared to control cultures, 4.8-fold compared to 4°C and 5.9-fold compared to 37°C (Table 3). It was the single most upregulated gene when comparing the 12°C group to the control and the 4°C group. Expression of MTNR1A in the skin is modified by several factors, including UVB exposure [21]. A significant association between MTNR1A polymorphisms and oral carcinogenesis has been demonstrated [22], and MTNR1A has been designated a putative tumor suppressor [23].
Cysteine-serine-rich nuclear protein 1 (CSRNP1) was upregulated 2.4-fold after storage at 12°C compared to control cultures (Table 3). This result is in line with recent findings from our research group, indicating that CSRNP1 is the second most upregulated gene (12.7-fold increase) in retinal pigment epithelial cells, when stored at 16°C compared to unstored cells (unpublished data).
Late cornified envelope 3D (LCE3D) was downregulated 5.9-fold at 4°C and 6.9-fold at 12°C compared to the control. Similarly, late cornified envelope 3E (LCE3E) and keratinocyte proline rich protein (KPRP) were downregulated 3.4-fold and 2.7-fold at 4°C compared to the control, respectively. Genes that code for late cornified envelope proteins are enriched and clustered within the epidermal differentiation complex [24, 25]. The relatively low differences in gene expression at 4°C and 12°C compared to the control stand in sharp contrast to the much greater changes observed when comparing 12°C cultures to 37°C.
Expression of Differentiation Markers
In addition to RPTN, KDAP, KRT10, LIPK, and CRNN, as presented in the previous section, the following genes associated with differentiation were upregulated: First, the oral mucosal differentiation marker keratin 4 (KRT4) [26] was upregulated 10.2-fold at 37°C storage compared to 12°C (Table 4 and Fig 4), indicative of a lower degree of differentiation in 12°C cultures. Second, keratin 6B, a specific marker of oral mucosal cells [26], was upregulated 2.5-fold at 37°C compared to 12°C cultures. Third, expression of keratin 19 (KRT19), a marker of undifferentiated cells [27], was 1.5-fold higher in cells stored at 12°C compared to those stored at 37°C (Table 4). Taken together, a total of 17 keratins were differentially regulated at 37°C compared to 12°C; 14 of these were upregulated.
The color scale illustrates the relative expression level of mRNAs: red color represents a high expression level; blue color represents a low expression level.
Other structural markers of keratinocyte differentiation had also changed during storage at 37°C. Filaggrin (FLG), which aggregates keratin intermediate filaments in mammalian epidermis [28], was upregulated 18.9-fold at 37°C compared to 12°C. Involucrin (IVL), a marker of differentiated keratinocytes [17, 18, 27], was upregulated 8.8-fold in cells stored at 37°C compared to those stored at 12°C, indicating increased differentiation of cells stored at 37°C. The lipase M (LIPM) gene, closely related to lipase K and exclusively expressed in the epidermis [19], was upregulated 3.7-fold at 37°C compared to 12°C (Table 4).
The cornified cell envelope is an insoluble protein layer that provides barrier function to stratified squamous epithelial cells [29]. Small proline-rich proteins (SPRRs) are constituents of this structure, and their expression is restricted to terminally differentiating squamous cells [18, 30]. Eight SPRR genes were upregulated between 2- and 11-fold at 37°C compared to 12°C, further indicating a more differentiated phenotype of cells stored at this temperature (Table 4). Apart from a 1.7-fold downregulation of TP63, a marker of undifferentiated cells, at 37°C compared to 12°C, few stem cell markers seemed to be affected by storage temperature. Neither OCT-4, FGF2, nor Nanog were differentially expressed at 12°C compared to 37°C, suggesting no significant impact of storage temperature on these stem cell related genes.
Regulation of Cell-Cell Contact
Identified claudins (CLDN) 1, 4, 7, 9, and 16 were upregulated between 1.5 and 3.6-fold at 37°C compared to 12°C. Genes encoding tight junction proteins 1 (TJP1) and 3 (TJP3) were both upregulated 1.8 and 1.6-fold at 37°C compared to 12°C, respectively (Table 4 and Fig 4). These changes indicate an increased synthesis of tight junctions in cells stored at 37°C.
Regulatory changes of the constituents of the desmosomal adherens junction were noted (Table 4 and Fig 4). Desmosomes are intercellular junctions that link the intermediate filaments of the cytoskeleton of neighboring epithelial cells and consist of desmocollins, desmogleins, and cadherins [15]. Desmocollins (DSC) 1, 2, and 3 were upregulated 5.5, 4.5, and 2.2-fold at 37°C compared to 12°C, respectively. Loss of function of these genes is associated with skin barrier defects [31] and metastasis of cancer cells [32]. Desmocollin 3 has been described as a tumor suppressor of several types of cancer [33–36]. Desmogleins 1 and 3 were upregulated 94.6-fold and 2.8-fold, respectively. Of the many identified cadherins in our material, only cadherins 2 and 4 were differentially regulated, with a 3.9-fold upregulation and a 1.7-fold downregulation at 37°C, respectively (Table 4). Our findings point in the direction of increased adherence between cells stored at 37°C compared to other temperatures.
Regulation of Cellular Stress Responses
Very few genetic markers of the oxidative stress response were significantly altered when comparing cells stored at 12°C and 37°C, indicating little difference in oxidative insult between these temperatures (Fig 5). The heat shock protein family members comprise an important cellular defense pathway [37]. The heat shock protein encoding gene HSPB8 was upregulated 27.6-fold at 37°C compared to 12°C, which might indicate cell stress (Table 4). Cornulin, also known as squamous epithelial heat shock protein 53, was upregulated 43.2-fold at 37°C (Table 2). This protein may play a role in the mucosal/epithelial immune response in addition to its role in epidermal differentiation [38]. An additional three heat shock proteins were upregulated between 1.6 and 2.3-fold at 37°C, while two were downregulated 1.8 and 1.9-fold at the same temperature (Table 4).
The color scale illustrates the relative expression level of mRNAs: red color represents a high expression level; blue color represents a low expression level.
Various environmental stressors have been shown to induce glucocorticoid production in epidermal keratinocytes [39]. Several genes coding for enzymes, receptors and transport proteins involved in the production of corticosteroids were investigated in the present material. Their expression was either not detected (CYP11B1, CYP11A1, 3βHSD, CYP17, CYP21A2, MC2, StAR) or not significantly altered (NR3C1, encoding the glucocorticoid receptor) when comparing any of the culture groups.
Collectively, these findings suggest that storage at 37°C, compared to 12°C, induces a heat shock response, but does not trigger oxidative stress.
Regulation of Signaling Pathways
Expression levels of the Wnt, BMP, Hedgehog, JAK/STAT, Notch, and TGF-β signaling pathways were analyzed. Cultures stored at 12°C showed no changes in expression levels of either pathway compared to control cells, and cultures stored at 4°C did not differ from the 12°C cultures. However, cells stored at 37°C expressed slight changes in regulation of various elements of all these pathways compared to the 12°C group. These changes included both up- and downregulation of different elements in all signaling pathways except the Hedgehog pathway. In this pathway, regulation was exclusively negative at 37°C and was comprised of pathway elements cyclin B1 (CCNB1, downregulated 4.9-fold), patched 2 (PTCH2, downregulated 8.6-fold), arrestin, β2 (ARRB2, downregulated 1.7-fold), and serine/threonine kinase 36 (STK36, downregulated 1.8-fold) (Table 4, Figs 4 and 6).
Pathway elements marked in red are significantly downregulated at 37°C. There is no upregulation of pathway elements.
Regulation of Proliferation, Apoptosis, and Cell Death
The effect of storage temperature on the expression of proliferation and cell death markers was also studied. The proliferation marker Ki-67 [18] was upregulated 7.9-fold at 12°C compared to 37°C cultures, indicative of increased proliferative potential (not shown). Expression of proliferation markers ABCG2 and PCNA was not significantly altered between culture groups (not shown). Slight regulatory changes in the expression of cell death markers were noted between the culture groups, but few genes were markedly changed. Thirteen cell death-related genes were upregulated between 1.6 and 13.7-fold at 37°C compared to 12°C, while 15 genes were downregulated between 1.5 and 6.1-fold (Table 4, Fig 5). Caspase 4 (CASP4) was the only caspase-encoding gene to be differentially expressed, with a 1.7-fold upregulation at 37°C compared to 12°C. Lipase H (LIPH), which is selectively upregulated in lung cancer and associated with increased survival in lung cancer patients [40], was upregulated 13.7-fold at 37°C compared to 12°C.
Of note, of the 114 genes important for cellular function that were significantly regulated at 37°C compared to 12°C (Table 4), only one was differentially regulated when comparing 4°C storage to 12°C storage (Growth arrest and DNA-damage-inducible beta (GADD45B); downregulated 2.81-fold), control cultures to 12°C storage (GADD45B; downregulated 3.6-fold), and 4°C storage to the control (member RAS oncogene family (RAB3B); downregulated 1.53-fold). GADD45 is induced by environmental stress or DNA damage [41], while RAB3B localizes to tight junctions where it has been suggested to contribute to the polarization of epithelia [42, 43]. Hence, 4°C and 12°C storage do not induce notable changes in the regulation of the important genes analyzed herein, and the two storage groups seem to offer equivalent results.
Quantitative Real-Time PCR Validation of Microarray Data
HSPB8, TP63 and KRT10 were selected for validation by qPCR (Table 5). The expression of HSPB8 was substantially upregulated at 37°C compared to 12°C; a 123.6-fold upregulation by PCR compared to a 27.6-fold upregulation by microarray. The expression of TP63 was significantly downregulated 3.6-fold at 37°C compared to 12°C, which is in line with the microarray results demonstrating a 1.7-fold downregulation at 37°C. Keratin 10 expression was significantly upregulated 6.7-fold at 37°C compared to 12°C. Upregulation of this gene was higher in the microarray analysis (45.6-fold). Expression levels of HSPB8, TP63, and KRT10 at the remaining temperatures were consistent between qPCR and microarray results, showing no significant differential regulation between temperatures (Table 5 and Fig 7).
Black bars indicate microarray expression values; grey bars represent PCR verification values. *P < 0.05.
Discussion
The current study investigated the effects of storage temperature on gene expression in cultured HOK using microarray analysis. The temperatures selected included 4°C, the standard temperature of a refrigerator; 12°C, a temperature which previously gave the most optimal results with regard to both morphology and viability of stored HOK; and 37°C, the temperature of a standard cell culture incubator.
Five of the six most differentially regulated proteins at 37°C compared to 12°C are directly associated with epithelial differentiation. An epidermal differentiation profile of these HOK cells is regarded as a disadvantage when used for treating LSCD [9]. The RPTN gene is very active during the final steps of epidermal keratinocyte differentiation, since the repetin protein is associated with the keratin intermediate filaments that are present in mature epidermal cells [14]. Keratinocyte differentiation-associated protein localizes to the stratum corneum of normal skin, but is expressed in suprabasal keratinocytes in psoriatic lesions [17]. Expression of the KDAP gene is markedly upregulated during keratinocyte differentiation in vitro [17]. Filaggrin contributes to the hydration and pH homeostasis of the stratum corneum [28], and mutations of the filaggrin gene are associated with ichtyosis vulgaris [44] and eczema [45]. Cornulin is a squamous cell-specific polypeptide [37] which is downregulated in eczema [46] and is a component of the epithelial innate immune response [38]. This protein is a constituent of the heat shock response of esophageal squamous epithelial tissue, where it is known as the squamous epithelial heat shock protein 53 (SEP53) [37]. Cornulin’s expression increases markedly as a consequence of heat shock, which might indicate activation of this cellular defense pathway in cell cultures stored at 37°C compared to those stored at 12°C. However, the activation of cornulin as a component of the differentiation process, and not primarily as a heat shock modulating protein, is an alternative interpretation.
The oral mucosal marker keratin 6b [26] was upregulated 2.5-fold in cultures stored at 37°C compared to 12°C. Its relative downregulation in 12°C cultures might indicate a dedifferentiation of cells at these temperatures. The upregulation of keratin 6b has also been demonstrated in conjunctival epithelium of patients with Sjögren’s syndrome [47]. Keratin 6b is also used as a marker of activated keratinocytes, and its expression can be induced both by proliferative signals and the proinflammatory cytokine TNF-α [48]. TNF was upregulated 1.4-fold at 37°C compared to 12°C (not shown), and several TNF-related proteins are differentially regulated at 37°C (Table 4). Given the stable expression of proliferative markers at 37°C compared to 12°C, it is more likely that the keratin 6b upregulation might be a consequence of inflammatory responses rather than proliferative signals.
CSRNP1 was upregulated 2.4-fold after storage at 12°C compared to control cultures. CSRNP1 has been described as a tumor suppressor gene, its expression level decreased in several types of cancers [49]. Overexpression has been reported to halt cell cycle progression at mitosis [50]. Its function is essential for normal development of the brain [51], and it is a known negative regulator of the Wnt pathway [52].
Both adhesion and tight junction-related genes were significantly upregulated at 37°C compared to 12°C and control cultures. The upregulation of DSC1 in cultures stored at 37°C is in line with our findings in retinal pigment epithelial cells [53]. These findings indicate a functional change towards a more tightly woven, squamous-like epithelium after storage at 37°C.
The expression of stress-related genes was also evaluated. Except for a few heat shock proteins, few genetic markers of the oxidative stress response were significantly altered when comparing cells stored at 12°C and 37°C. Recent studies have described that epidermal keratinocytes subjected to various environmental stressors can respond by upregulating their glucocorticoid production [39, 54]. Expression of genes encoding enzymes, receptors and transport proteins necessary for the production and effect of corticosteroids were either not detected or were unaltered in our material.
Cultures stored at 37°C show both up- and downregulation of several important molecules of the Wnt, BMP, JAK/STAT, Notch, and TGF-β signaling pathways, rendering the effects of these regulational changes inconclusive. However, the Hedgehog pathway was exclusively downregulated in the 37°C cell group compared to the 12°C group. The Hedgehog signaling pathway is instrumental for vertebrate embryogenesis and has been demonstrated to regulate cell fate, proliferation, and survival in multiple cell types, especially those of neuroectodermal origin such as cells of the retina and optic nerve [55]. The pathway also regulates adult stem cells in several self-renewing organs including the subventricular zone of the brain [56, 57]. The deactivating effect of 37°C on this pathway may therefore perturb cellular function. The sonic hedgehog protein, Shh, binds Patched, a transmembrane receptor of the target cell [58]. Patched functions as a tumor suppressor in the hedgehog signaling pathway [58] and mutations of the gene have been detected in basal cell carcinomas and medulloblastomas, among other cancers [55, 59]. The downregulation of Patched in the 37°C culture group compared to 12°C might destabilize cellular quiescence, an undesirable event in cell preservation.
The patched protein interacts with cyclin B1 and participates in determining its cellular localization [60]. Cyclin B1 is a regulatory protein involved in the promotion of mitosis, transitioning the cell from the G2 to M phase. The effect of Patched on Cyclin 1 is inhibition of cellular proliferation [60]. Arrestin β2, which was downregulated at 37°C, is a ubiquitously distributed protein with a critical role in the regulation of several important signaling pathways, including Hedgehog [61]. The serine/threonine kinase 36 (STK36), also downregulated at 37°C, plays a key role in the Hedgehog signaling pathway. Hence, several crucial constituents of the pathway are downregulated at 37°C, contributing to a reduced activity of the pathway as a whole. The finding of perturbed signaling pathway elements of several major pathways in cultures stored at 37°C is in concordance with findings of similar changes in retinal pigment epithelial cells stored at 37°C (unpublished data).
The analysis of cell proliferation and death markers after storage at different temperatures indicates a disturbance of the regulation of several cell death-related genes at 37°C. These changes may be partly responsible for the reduced viability found in cultures stored at this temperature.
Expression of the melatonin receptor MTNR1A was significantly upregulated at 12°C compared to the other groups. Melatonin exerts numerous effects on a multitude of organ systems, including the skin. Its effects are mediated through both receptor dependent and independent mechanisms [62]. Importantly, the skin cannot be regarded a passive target for the effects of melatonin, but a vibrant site of its synthesis and metabolism [63, 64]. There is accumulating evidence that the stringently regulated effects of melatonin in the skin are organized through an interlaced local neuroendocrine system exploiting both auto- and paracrine mechanisms [63, 64]. Melatonin and its metabolites exert broad antioxidant effects, and melatonin is able to activate cytoprotective molecules and enzymes, including glutathione [63, 65, 66]. It may also protect DNA from oxidative damage, thereby providing anti-apoptotic and anti-carcinogenic effects [63, 64, 67]. Specifically, both the initiation and promotion of skin carcinogenesis can be decreased by melatonin [68]. There is also evidence that its oncostatic effects are dependent on the MTNR1A receptor [69–71], and that the tumor suppressive effect of the MTNR1A gene is silenced in oral squamous cell carcinoma [23]. A significant association between MTNR1A polymorphisms and oral carcinogenesis has been demonstrated [22], in which environmental factors (betel quid chewing and cigarette smoking) are required to increase the susceptibility to oral cancer in individuals with MTNR1A gene polymorphisms.
It has also been demonstrated that melatonin and its metabolites protect keratinocytes from UV radiation [64, 72]. Following UVB exposure, MTNR1A expression has been shown to be upregulated in normal neonatal epidermal melanocytes and downregulated in melanoma lines [73]. In the current study, HOK were cultured in the dark except when being handled, but cultures stored at 12°C and 37°C were probably exposed to small amounts of light during heating of the storage containers. This light exposure might have contributed in inducing MTNR1A expression in the 12°C group, but it does not explain why a similar upregulation in the 37°C storage group could not be detected.
While there are clear indications to the upregulation of differentiation-related genes at 37°C compared to 12°C, the gene regulation changes in the 4°C and 12°C groups compared to the control are not as clear. However, the downregulation of some constituents of the epidermal differentiation complex at 4°C might indicate some degree of de-differentiation in these cultures, contrary to the effect of 37°C storage. For cells stored at 12°C, the evidence toward a less differentiated phenotype was not as strong.
In conclusion, HOK cultures stored at 37°C demonstrated considerably larger changes in both the amount of genes affected as well as their differential regulation levels compared to unstored cells than cultured HOK stored at 4°C and 12°C. Temporary storage of cell cultures in sealed containers at 37°C, rather than at 12°C, appears to promote differentiation similarly to conventional cell culture, which employs a humidified incubator at 37°C with CO2 supply. Storage at 37°C may reduce the stemness, and thereby the therapeutic potential, of cultured cells intended for the treatment of LSCD [9]. Storage at 12°C also maintains the regulation of genes closer to control levels than storage at 4°C. However, storage at 4°C might steer cells toward a less differentiated phenotype. Nevertheless, we conclude that storage at 4°C and 12°C are more suitable than storage at 37°C for preserving cultured HOK for transplantation. Our findings are in line with a recent study [12], demonstrating superior viability of HOK when stored at 12°C compared to 4°C and 37°C. Thus, collectively, 12°C seems to be the most ideal storage temperature among those investigated.
Author Contributions
Conceived and designed the experiments: TPU IF JRE DAD EM MG LP. Performed the experiments: TPU RI OKO. Analyzed the data: TPU AS OKO LP. Contributed reagents/materials/analysis tools: TPU RI OKO EM. Wrote the paper: TPU IF JRE AS OKO DAD EM MG LP.
References
- 1. Utheim TP. Limbal epithelial cell therapy: past, present, and future. Methods in molecular biology (Clifton, NJ). 2013;1014:3–43. Epub 2013/05/22. pmid:23690002.
- 2. Pellegrini G, Traverso CE, Franzi AT, Zingirian M, Cancedda R, De Luca M. Long-term restoration of damaged corneal surfaces with autologous cultivated corneal epithelium. Lancet (London, England). 1997;349(9057):990–3. Epub 1997/04/05. pmid:9100626.
- 3. Niethammer D, Kummerle-Deschner J, Dannecker GE. Side-effects of long-term immunosuppression versus morbidity in autologous stem cell rescue: striking the balance. Rheumatology (Oxford). 1999;38(8):747–50. Epub 1999/09/29. pmid:10501425.
- 4. Nakamura T, Inatomi T, Sotozono C, Amemiya T, Kanamura N, Kinoshita S. Transplantation of cultivated autologous oral mucosal epithelial cells in patients with severe ocular surface disorders. The British journal of ophthalmology. 2004;88(10):1280–4. Epub 2004/09/21. pmid:15377551; PubMed Central PMCID: PMCPMC1772364.
- 5. Nishida K, Yamato M, Hayashida Y, Watanabe K, Yamamoto K, Adachi E, et al. Corneal reconstruction with tissue-engineered cell sheets composed of autologous oral mucosal epithelium. The New England journal of medicine. 2004;351(12):1187–96. Epub 2004/09/17. pmid:15371576.
- 6. Utheim TP. Concise review: transplantation of cultured oral mucosal epithelial cells for treating limbal stem cell deficiency-current status and future perspectives. Stem Cells. 2015;33(6):1685–95. pmid:25786664.
- 7. Daniels JT, Secker GA, Shortt AJ, Tuft SJ, Seetharaman S. Stem cell therapy delivery: treading the regulatory tightrope. Regen Med. 2006;1(5):715–19. pmid:17465740
- 8. Ahmad S, Osei-Bempong C, Dana R, Jurkunas U. The culture and transplantation of human limbal stem cells. Journal of Cellular Physiology. 2010;225(1):15–9. pmid:20506173
- 9. Rama P, Matuska S, Paganoni G, Spinelli A, De Luca M, Pellegrini G. Limbal stem-cell therapy and long-term corneal regeneration. The New England journal of medicine. 2010;363(2):147–55. Epub 2010/06/25. pmid:20573916.
- 10. O'Callaghan AR, Daniels JT. Limbal epithelial stem cell therapy: controversies and challenges. Stem Cells. 2011;29(12):1923–32. pmid:21997829
- 11. Utheim TP, Raeder S, Olstad OK, Utheim OA, de la Paz MF, Cheng R, et al. Comparison of the histology, gene expression profile, and phenotype of cultured human limbal epithelial cells from different limbal regions. InvestOphthalmolVisSci. 2009;50(11):5165–72. iovs.08-2884 [pii];
- 12. Islam R, Jackson C, Eidet JR, Messelt EB, Corraya RM, Lyberg T, et al. Effect of Storage Temperature on Structure and Function of Cultured Human Oral Keratinocytes. PLOS ONE. 2015;10(6):e0128306. pmid:26052937; PubMed Central PMCID: PMCPMC4459984.
- 13. Pasovic L, Utheim TP, Maria R, Lyberg T, Messelt EB, Aabel P, et al. Optimization of Storage Temperature for Cultured ARPE-19 Cells. Journal of Ophthalmology. 2013;2013:11.
- 14. Huber M, Siegenthaler G, Mirancea N, Marenholz I, Nizetic D, Breitkreutz D, et al. Isolation and Characterization of Human Repetin, a Member of the Fused Gene Family of the Epidermal Differentiation Complex. J Invest Dermatol. 2005;124:998–1007. pmid:15854042
- 15. Garrod DR, Merritt AJ, Nie Z. Desmosomal cadherins. Current Opinion in Cell Biology. 2002;14:537–45. pmid:12231347
- 16. Oomizu S, Sahuc F, Asahina K, Inamatsu M, Matsuzaki T, Sasaki M, et al. Kdap, a novel gene associated with the stratification of the epithelium. Gene. 2000;256:19–27. pmid:11054531
- 17. Tsuchida S, Bonkobara M, McMillan JR, Akiyama M, Yudate T, Aragane Y, et al. Characterization of Kdap, a Protein Secreted by Keratinocytes. J Invest Dermatol. 2004;122:1225–34. pmid:15140226
- 18. Gibbs D, Ponec M. Intrinsic regulation of differentiation markers in human epidermis, hard palate and buccal mucosa. Arch Oral Biochem. 2000;45:149–58.
- 19. Toulza E, Mattiuzzo NR, Galliano MF, Jonca N, Dossat C, Jacob D, et al. Large-scale identification of human genes implicated in epidermal barrier function. Genome biology. 2007;8(6):R107. pmid:17562024; PubMed Central PMCID: PMC2394760.
- 20. Contzler R, Favre B, Huber M, Hohl D. Cornulin, a New Member of the ‘‘Fused Gene” Family, Is Expressed During Epidermal Differentiation. J Invest Dermatol. 2005;124:990–97. pmid:15854041
- 21. Slominski A, Tobin DJ, Zmijewski MA, Wortsman J, Paus R. Melatonin in the skin: synthesis, metabolism and functions. Trends Endocrinol Metab. 2008;19(1):17–24. pmid:18155917.
- 22. Lin FY, Lin CW, Yang SF, Lee WJ, Lin YW, Lee LM, et al. Interactions between environmental factors and melatonin receptor type 1A polymorphism in relation to oral cancer susceptibility and clinicopathologic development. PLOS ONE. 2015;10(3):e0121677. pmid:25806809; PubMed Central PMCID: PMCPMC4373723.
- 23. Nakamura E, Kozaki K, Tsuda H, Suzuki E, Pimkhaokham A, Yamamoto G, et al. Frequent silencing of a putative tumor suppressor gene melatonin receptor 1 A (MTNR1A) in oral squamous-cell carcinoma. Cancer Sci. 2008;99(7):1390–400. pmid:18452558.
- 24. Backendorf C, Hohl D. A common origin for cornified envelope proteins? Nature genetics. 1992;2:91. pmid:1303269
- 25. Jackson B, Tilli CM, Hardman MJ, Avilion AA, MacLeod MC, Ashcroft GS, et al. Late cornified envelope family in differentiating epithelia—response to calcium and ultraviolet irradiation. J Invest Dermatol. 2005;124(5):1062–70. pmid:15854049.
- 26. Sugiyama H, Yamato M, Nishida K, Okano T. Evidence of the survival of ectopically transplanted oral mucosal epithelial stem cells after repeated wounding of cornea. Molecular therapy: the journal of the American Society of Gene Therapy. 2014;22(8):1544–55. pmid:24769908.
- 27. Calenic B, Ishkitiev N, Yaegaki K, Imai T, Costache M, Tovaru M, et al. Characterization of oral keratinocyte stem cells and prospects of its differentiation to oral epithelial equivalents. Rom J Morphol Embryol. 2010;51(4):641–45. pmid:21103620
- 28. Thyssen JP, Kezic S. Causes of epidermal filaggrin reduction and their role in the pathogenesis of atopic dermatitis. The Journal of allergy and clinical immunology. 2014;134(4):792–9. pmid:25065719.
- 29. Steinert PM, Candi E, Kartasova T, Marekov L. Small proline-rich proteins are cross-bridging proteins in the cornified cell envelopes of stratified squamous epithelia. Journal of structural biology. 1998;122(1–2):76–85. Epub 1998/09/02. pmid:9724607.
- 30. Hohl D, de Viragh PA, Amiguet-Barras F, Gibbs S, Backendorf C, Huber M. The small proline-rich proteins constitute a multigene family of differentially regulated cornified cell envelope precursor proteins. J Invest Dermatol. 1995;104(6):902–9. pmid:7769256
- 31. Chidgey M, Brakebusch C, Gustafsson E, Cruchley A, Hail C, Kirk S, et al. Mice lacking desmocollin 1 show epidermal fragility accompanied by barrier defects and abnormal differentiation. J Cell Biol. 2001;155(5):821–32. Epub 2001/11/21. pmid:11714727; PubMed Central PMCID: PMC2150874.
- 32. Shinohara M, Hiraki A, Ikebe T, Nakamura S, Kurahara S-I, Shirasuna K, et al. Immunohistochemical study of desmosomes in oral squamous cell carcinoma: correlation with cytokeratin and E-cadherin staining, and with tumour behaviour. J Pathol. 1998;184:369–81. pmid:9664902
- 33. Cui T, Chen Y, Yang L, Knosel T, Huber O, Pacyna-Gengelbach M, et al. The p53 target gene desmocollin 3 acts as a novel tumor suppressor through inhibiting EGFR/ERK pathway in human lung cancer. Carcinogenesis. 2012;33(12):2326–33. Epub 2012/09/04. pmid:22941060.
- 34. Oshiro MM, Kim CJ, Wozniak RJ, Junk DJ, Munoz-Rodriguez JL, Burr JA, et al. Epigenetic silencing of DSC3 is a common event in human breast cancer. Breast Cancer Res. 2005;7(5):R669–80. Epub 2005/09/20. pmid:16168112; PubMed Central PMCID: PMC1242132.
- 35. Wang L, Liu T, Wang Y, Cao L, Nishioka M, Aguirre RL, et al. Altered expression of desmocollin 3, desmoglein 3, and beta-catenin in oral squamous cell carcinoma: correlation with lymph node metastasis and cell proliferation. Virchows Arch. 2007;451(5):959–66. Epub 2007/09/12. pmid:17846785.
- 36. Cui T, Chen Y, Yang L, Knosel T, Zoller K, Huber O, et al. DSC3 expression is regulated by p53, and methylation of DSC3 DNA is a prognostic marker in human colorectal cancer. Br J Cancer. 2011;104(6):1013–9. Epub 2011/03/03. pmid:21364582; PubMed Central PMCID: PMC3065270.
- 37. Yagui-Beltran A, Craig AL, Lawrie L, Thompson D, Posipilova S, Johnston D, et al. The human oesophageal squamous epithelium exhibits a novel type of heat shock protein response. Eur J Biochem. 2001;268:5343–55. pmid:11606197
- 38. Little TJ, Nelson L, Hupp T. Adaptive evolution of a stress response protein. PLOS ONE. 2007;2(10):e1003. pmid:17925851; PubMed Central PMCID: PMC1994589.
- 39. Zhu G, Janjetovic Z, Slominski A. On the role of environmental humidity on cortisol production by epidermal keratinocytes. Exp Dermatol. 2014;23(1):15–7. pmid:24372648; PubMed Central PMCID: PMCPMC3923494.
- 40. Seki Y, Yoshida Y, Ishimine H, Shinozaki-Ushiku A, Ito Y, Sumitomo K, et al. Lipase member H is a novel secreted protein selectively upregulated in human lung adenocarcinomas and bronchioloalveolar carcinomas. Biochem Biophys Res Commun. 2014;443(4):1141–7. pmid:24380866.
- 41. Takekawa M, Saito H. A Family of Stress-Inducible GADD45-like Proteins Mediate Activation of the Stress-Responsive MTK1/MEKK4 MAPKKK. Cell. 1998;95(4):521–30. pmid:9827804
- 42. Weber E, Berta G, Tousson A, John PS, Green MW, Gopalokrishnan U, et al. Expression and Polarized Targeting of a Rab3 Isoform in Epithelial Cells. J Cell Biol. 1994;125(3):583–94. pmid:8175882
- 43. Yamamoto Y, Nishimura N, Morimoto S, Kitamura H, Manabe S, Kanayama H-O, et al. Distinct roles of Rab3B and Rab13 in the polarized transport of apical, basolateral, and tight junctional membrane proteins to the plasma membrane. Biochemical and Biophysical Research Communications. 2003;308(2):270–5. pmid:12901864
- 44. Smith FJ, Irvine AD, Terron-Kwiatkowski A, Sandilands A, Campbell LE, Zhao Y, et al. Loss-of-function mutations in the gene encoding filaggrin cause ichthyosis vulgaris. Nature genetics. 2006;38(3):337–42. pmid:16444271.
- 45. Palmer CN, Irvine AD, Terron-Kwiatkowski A, Zhao Y, Liao H, Lee SP, et al. Common loss-of-function variants of the epidermal barrier protein filaggrin are a major predisposing factor for atopic dermatitis. Nature genetics. 2006;38(4):441–6. pmid:16550169.
- 46. Lieden A, Ekelund E, Kuo IC, Kockum I, Huang CH, Mallbris L, et al. Cornulin, a marker of late epidermal differentiation, is down-regulated in eczema. Allergy. 2009;64(2):304–11. pmid:19133922.
- 47. Kawasaki S, Kawamoto S, Yokoi N, Connon C, Minesaki Y, Kinoshita S, et al. Up-regulated gene expression in the conjunctival epithelium of patients with Sjögren's syndrome. Exp Eye Res. 2003;77(1):17–26. pmid:12823984
- 48. Komine M, Rao LS, Kaneko T, Tomic-Canic M, Tamaki K, Freedberg IM, et al. Inflammatory versus proliferative processes in epidermis. Tumor necrosis factor alpha induces K6b keratin synthesis through a transcriptional complex containing NFkappa B and C/EBPbeta. J Biol Chem. 2000;275(41):32077–88. pmid:10887174.
- 49. Ishiguro H, Tsunoda T, Tanaka T, Fujii Y, Nakamura Y, Furukawa Y. Identification of AXUD1, a novel human gene induced by AXIN1 and its reduced expression in human carcinomas of the lung, liver, colon and kidney. Oncogene. 2001;20:5062–66. pmid:11526492
- 50. Glavic A, Molnar C, Cotoras D, de Celis JF. Drosophila Axud1 is involved in the control of proliferation and displays pro-apoptotic activity. Mech Dev. 2009;126(3–4):184–97. Epub 2008/12/17. pmid:19084594.
- 51. Feijoo CG, Sarrazin AF, Allende ML, Glavic A. Cystein-serine-rich nuclear protein 1, Axud1/Csrnp1, is essential for cephalic neural progenitor proliferation and survival in zebrafish. Dev Dyn. 2009;238(8):2034–43. Epub 2009/06/23. pmid:19544579.
- 52. Wang K, Ling T, Wu H, Zhang J. Screening of candidate tumor-suppressor genes in 3p21.3 and investigation of the methylation of gene promoters in oral squamous cell carcinoma. Oncol Rep. 2013;29(3):1175–82. Epub 2013/01/08. pmid:23292452.
- 53. Pasovic L, Eidet JR, Brusletto BS, Lyberg T, Utheim TP. Effect of Storage Temperature on Key Functions of Cultured Retinal Pigment Epithelial Cells. Journal of ophthalmology. 2015;2015:263756. Epub 2015/10/09. pmid:26448872; PubMed Central PMCID: PMCPmc4584032.
- 54. Takei K, Denda S, Kumamoto J, Denda M. Low environmental humidity induces synthesis and release of cortisol in an epidermal organotypic culture system. Exp Dermatol. 2013;22(10):662–4. pmid:24079737.
- 55. Varjosalo M, Taipale J. Hedgehog: functions and mechanisms. Genes Dev. 2008;22(18):2454–72. pmid:18794343.
- 56. Ihrie RA, Shah JK, Harwell CC, Levine JH, Guinto CD, Lezameta M, et al. Persistent sonic hedgehog signaling in adult brain determines neural stem cell positional identity. Neuron. 2011;71(2):250–62. pmid:21791285; PubMed Central PMCID: PMC3346180.
- 57. Palma V, Lim DA, Dahmane N, Sanchez P, Brionne TC, Herzberg CD, et al. Sonic hedgehog controls stem cell behavior in the postnatal and adult brain. Development. 2005;132(2):335–44. pmid:15604099; PubMed Central PMCID: PMC1431583.
- 58. Stone DM, Hynes M, Armanini M, Swanson TA, Gu Q, Johnson RL, et al. The tumour-suppressor gene patched encodes a candidate receptor for Sonic hedgehog. Nature. 1996;384:129–33. pmid:8906787
- 59. Villavicencio EH, Walterhouse DO, Iannaccone PM. The Sonic Hedgehog–Patched–Gli Pathway in Human Development and Disease. Am J Hum Genet. 2000;67:1047–54. pmid:11001584
- 60. Barnes E, Kong M, Ollendorff V, Donoghue DJ. Patched1 interacts with cyclin B1 to regulate cell cycle progression. Embo J. 2001;20(9):2214–23. pmid:11331587
- 61. Kovacs JJ, Hara MR, Davenport CL, Kim J, Lefkowitz RJ. Arrestin development: emerging roles for beta-arrestins in developmental signaling pathways. Developmental cell. 2009;17(4):443–58. pmid:19853559; PubMed Central PMCID: PMC3221601.
- 62. Slominski RM, Reiter RJ, Schlabritz-Loutsevitch N, Ostrom RS, Slominski AT. Melatonin membrane receptors in peripheral tissues: distribution and functions. Mol Cell Endocrinol. 2012;351(2):152–66. pmid:22245784; PubMed Central PMCID: PMCPMC3288509.
- 63. Slominski A, Wortsman J, Tobin DJ. The cutaneous serotoninergic/melatoninergic system: securing a place under the sun. FASEB J. 2005;19(2):176–94. pmid:15677341.
- 64. Slominski AT, Kleszczynski K, Semak I, Janjetovic Z, Zmijewski MA, Kim TK, et al. Local melatoninergic system as the protector of skin integrity. Int J Mol Sci. 2014;15(10):17705–32. pmid:25272227; PubMed Central PMCID: PMCPMC4227185.
- 65. Rodriguez C, Mayo JC, Sainz RM, Antolín I, Herrera F, Martín V, et al. Regulation of antioxidant enzymes: a significant role for melatonin. J Pineal Res. 2004;36(1):1–9. pmid:14675124
- 66. Martín M, Macías M, Escames G, León J, Acuña-Castroviejo D. Melatonin but not vitamins C and E maintains glutathione homeostasis in t-butyl hydroperoxide-induced mitochondrial oxidative stress. FASEB J. 2000;14(12):1677–9. pmid:10973915
- 67. Karbownik M, Lewinski A, Reiter RJ. Anticarcinogenic actions of melatonin which involve antioxidative processes: comparison with other antioxidants. Int J Biochem Cell Biol. 2001;33(8):735–53. pmid:11404179
- 68. Kumar CA, Das UN. Effect of melatonin on two stage skin carcinogenesis in Swiss mice. Med Sci Monit. 2000;6(3):471–75. pmid:11208355
- 69. Jones MP, Melan MA, Witt-Enderby PA. Melatonin decreases cell proliferation and transformation in a melatonin receptor-dependent manner. Cancer Lett. 2000;151(2):133–43. pmid:10738107
- 70. Ram PT, Dai J, Yuan L, Dong C, Kiefer TL, Lai L, et al. Involvement of the mt1 melatonin receptor in human breast cancer. Cancer Lett. 2002;178(2):141–50.
- 71. Hill SM, Belancio VP, Dauchy RT, Xiang S, Brimer S, Mao L, et al. Melatonin: an inhibitor of breast cancer. Endocr Relat Cancer. 2015;22(3):R183–204. pmid:25876649; PubMed Central PMCID: PMCPMC4457700.
- 72. Janjetovic Z, Nahmias ZP, Hanna S, Jarrett SG, Kim TK, Reiter RJ, et al. Melatonin and its metabolites ameliorate ultraviolet B-induced damage in human epidermal keratinocytes. J Pineal Res. 2014;57(1):90–102. pmid:24867336; PubMed Central PMCID: PMCPMC4106994.
- 73. Slominski A, Fischer TW, Zmijewski MA, Wortsman J, Semak I, Zbytek B, et al. On the Role of Melatonin in Skin Physiology and Pathology. Endocrine. 2005;27(2):137–47. pmid:16217127