Whole-transcriptome gene expression profiling in an epidermolysis bullosa simplex Dowling-Meara model keratinocyte cell line uncovered novel, potential therapeutic targets and affected pathways
BMC Research Notes volume 8, Article number: 785 (2015)
To be able to develop effective therapeutics for epidermolysis bullosa simplex (EBS), it is necessary to elucidate the molecular pathomechanisms that give rise to the disease’s characteristic severe skin-blistering phenotype.
Starting with a whole-transcriptome microarray analysis of an EBS Dowling-Meara model cell line (KEB7), we identified 207 genes showing differential expression relative to control keratinocytes. A complementary qRT-PCR study of 156 candidates confirmed 76.58 % of the selected genes to be significantly up-regulated or down-regulated (p-value <0.05) within biological replicates. Our hit list contains previously identified genes involved in epithelial cell proliferation, cell-substrate adhesion, and responses to diverse biological stimuli. In addition, we identified novel candidate genes and potential affected pathways not previously considered as relevant to EBS pathology.
Our results broaden our understanding of the molecular processes dysregulated in EBS.
KRT5 and KRT14 are the main stress-absorbing keratins in basal keratinocytes of the epidermis and related stratified epithelia. These rod-shaped proteins form heterodimeric units that interact to build up the cytoskeletal intermediate filament (IF) network, a resilient yet adaptable scaffold that maintains cellular structural stability and, in turn, normal skin integrity and function. Disruption of IF-organization as a consequence of keratin mutation is the basis of a number of inherited skin fragility syndromes [1, 2]. In the case of epidermolysis bullosa simplex (EBS), which exhibits several clinical variants , specific phenotypes can be largely correlated with the positions of missense mutations in structurally sensitive portions of either KRT5 or KRT14 . Aberrant IF organization results in fragile epidermal basal cells that readily lyse following mild mechanical trauma or minor traction, leading to intraepithelial fluid accumulation and recurrent blister formation [5–8]. At the molecular level, this cytoskeletal collapse manifests as aggregates of misfolded keratins, along with activation of stress-response cascades [9, 10]. Detailed elucidation of the underlying pathomechanisms in EBS is an important prerequisite for developing innovative therapeutics; however, relatively few studies have focused on expression-profiling of mutant epidermal cells (summarized in Additional file 1: Table S1). In 2007, Lu et al.  reported the expression profile of KRT5−/− EBS mouse epidermis, concentrating mainly on the regulation of inflammatory cytokines. Liovic et al.  investigated the response of the EBS keratinocyte cell lines KEB4 (mild EBS-loc phenotype, KRT14 mutation V270 M) and KEB7 (severe EBS Dowling-Meara (EBS-DM) phenotype, K14 mutation R125P) to hypo-osmotic stress in comparison to the wild-type cell line NEB1, and found dual-specificity phosphatases and their downstream targets ERK and p38 to be differentially regulated in EBS cells. A subsequent study by the same group identified differences in the expression profiles of cell-junction components in EBS versus wild-type cell lines . More recent profiling studies reported aberrant expression of genes involved in keratinization, cell growth, proliferation, the immune response, and fatty acid metabolism in EBS . Because there is limited overlap of the genes identified in those studies, further efforts seem necessary to fully elucidate EBS-relevant genes to better understand EBS pathology. We recently described the gene expression profile of KEB7 cells after applying suppression subtractive hybridization, and found dysregulated genes involved in keratinocyte differentiation, migration and wound healing . Here, we follow-up that analysis with a more expansive expression profiling study, combining whole-transcriptome microarray examination with bioinformatics-assisted functional clustering, complementary qRT-PCR validation, and western blotting analysis of selected hits. Using this approach we were able to verify candidate genes previously described as being differentially expressed in EBS-DM, and to discover other differentially regulated genes not previously implicated in EBS pathology. A graphical summary of our experimental design is shown in Fig. 1.
Cell lines and culture conditions
The immortalized keratinocyte cell lines KEB7 (EBS-DM severe phenotype caused by the R125P mutation in KRT14) and NEB1 (wild-type control cell line)  were cultured in standard DMEM (HyClone Laboratories, GE Healthcare, South Logan, UT, USA) supplemented with 25 % Ham’s F12 nutrient mixture, 10 % fetal bovine serum (FBS), 1.8 × 10−4 M adenine, 0.4 µg/ml hydrocortisone, 5 µg/ml transferrin, 2 × 10−11 M liothyronine, 5 µg/ml insulin and 10 µg/ml EGF. Unless otherwise stated, all chemicals and reagents were obtained from Life Technologies™ (Karlsruhe, Germany) or Sigma-Aldrich Inc. (Taufkirchen, Germany). Cell lines were cultivated in standard 25 cm2 flasks (Techno Plastic Products AG, Trasadingen, Switzerland) at 37 °C under a 5 % CO2 atmosphere without fibroblast feeder cells.
Total RNA isolation from cultivated KEB7 and NEB1 cells
Total RNA from cells grown to approximately 80 % confluence was isolated using an RNeasy Mini Kit (QiagenGmbH, Hilden, Germany) according to the manufacturer’s instructions. RNA was resuspended in nuclease-free water and quantified spectrophotometrically at 260 nm (DS-11 Spectrophotometer, DeNovix, Wilmington, DE, USA). RNA preparations were considered as suitable only if samples were of sufficient yield, exhibited intact bands on agarose gels, and displayed no spurious peaks or RNA degradation artifacts on the UV absorption spectrum. RNA samples were stored in aliquots at −80 °C.
Whole-transcriptome microarray analysis
Total RNA, isolated from a pool of different passages of KEB7 versus NEB1 cells (p25 and p28–37), was used for probing a high-precision, whole-transcript Human Gene 1.0 ST array (Affymetrix, Santa Clara, CA, USA). The entire experimental pipeline, from sense target labeling to hybridization, washing, array scanning and final raw data capture, was performed at the Center of Excellence for Fluorescent Bioanalysis KFB, University of Regensburg, Germany, using standard Affymetrix protocols, reagents and instrumentation. Four microarrays per cell line were processed to account for technical and/or biological variability. The raw data were checked for quality, background-adjusted, quantile-normalized (by imposing the same distribution of gene signal intensities for each array used under the condition that the expression of most genes is relatively unaltered across mutant and wild-type cells), and statistically analyzed by the linear model for microarray data (LIMMA) method  using the Integromics Biomarker Discovery® platform (IBD, Integromics, Granada, Spain). Assuming experimental deviations/errors of less than a log2-fold change of +1/−1, only candidates showing regulation greater than +2/−2 were subjected to further analysis. The respective statistical parameters were also used for the subsequent enrichment and qRT-PCR studies.
Gene ontology enrichment and functional clustering of differentially expressed genes
Affymetrix IDs of differentially regulated genes were subjected to comprehensive bioinformatical analysis embedded within IBD. Gene ontology (GO) classification was performed to confirm enrichment of microarray-derived genes and to assign them to specific biological themes and functions. To ascertain GO-term enrichment of genes statistically overrepresented in the candidate hit list (log2-fold change of ≥ +1 and ≤ −1, adjusted p-value of <0.05), the software tool calculated a one-sided hypergeometric p-value (identical to the one-tailed version of Fisher’s exact test), so that terms with values below 0.05 can be considered to be significantly enriched.
cDNA synthesis, primer design, and complementary qRT-PCR
Total RNA isolated from one passage of NEB1 and KEB7 cells (each cell line in technical triplicates) was treated with DNase I (Sigma-Aldrich, Taufkirchen, Germany) and then reverse transcribed into cDNA using an iScript™ cDNA Synthesis Kit (Bio-Rad Laboratories Inc., Hercules, CA USA) according to the manufacturer’s instructions. The resultant cDNAs served as template for complementary quantitative real-time PCR (qRT-PCR) in a 96-well plate format using GoTaq®qPCR Master Mix (Promega Corporation, Madison, WI, USA) and a CFX96™ instrument (Bio-Rad Laboratories). Using Batch Entrez (http://www.ncbi.nlm.nih.gov/sites/batchentrez) and WIBR UTR extractor (http://jura.wi.mit.edu/bioc/tools/utrs/), we extracted the coding sequences (highlighted within the downloaded mRNAs) of the 79 most promising up-regulated and the 79 most promising down-regulated genes uncovered in the microarray analysis to automatically design oligonucleotides with similar properties (uniform GC content, comparable melting temperature, optimal length and similar product size; BATCH Primer 3 program, http://probes.pw.usda.gov/cgi-bin/batchprimer3/batchprimer3.cgi). Forward and reverse primer stocks (100 pmol/µl, see Additional file 2: Table S2) were purchased in 96-well format from Life Technologies GmbH (Karlsruhe, Germany). Differential expression of selected candidate genes was investigated by performing three independent PCR-runs (biological and technical replicates) in an in-house developed and established 96-well qRT-PCR array format. Seven constitutively expressed housekeeping genes (ACTB, ANXA1, B2 M, GAPDH, HPRT1, RPL13A, and TUBB) served as internal controls for data normalization as well as determination of experimental variance. Fold inductions and statistical significance were assessed via RealTimeStatMiner® from Integromics® that calculates the relative quantity (RQ) and log2-fold change using ΔΔCt (cycle threshold) and efficiency correction . A parametric LIMMA moderate t-test  and a Benjamini–Hochberg false discovery rate (FDR) p-value correction , included as standard default parameters of the software (log2-fold change of +1/−1, p-value <0.05), were applied in the data analysis. The meaning of an FDR-adjusted p-value is that, if all genes with a p-value below a threshold of 0.05 are picked as differentially expressed, then the expected proportion of false discoveries in the selected group is controlled to be less than the threshold value, in this case 5 %.
KEB7/NEB1 total protein lysates and Western blotting
For preparation of KEB7 versus NEB1 protein lysates (using three 75-cm2 flasks per cell line), the culture medium was removed, the cells were rinsed twice with sterile PBS, and then lysed in RIPA buffer (Radio-Immunoprecipitation Assay buffer; 150 mM NaCl, 1 % IGEPAL® CA-630, 0.5 % deoxycholate, 0.1 % SDS, 50 mM Tris–HCl, pH 8.0, Sigma-Aldrich) supplemented with 1× Protease Inhibitor Cocktail (GE Healthcare, USA) for 10 min at RT (room temperature) and a further 20 min at −80 °C. Cells were detached with a cell scraper, and the resultant total protein lysates transferred to 1.5-ml tubes and centrifuged at maximum speed for 15 s to remove cellular debris. Protein in the supernatants was photometrically quantified by a standard Bradford assay (Bio-Rad, Hercules, CA, USA), and the samples were frozen at −80 °C for long-term storage. 10 µg of protein lysate per slot were subjected to SDS-PAGE (sodium dodecyl sulfate polyacrylamide gel electrophoresis) and then transferred onto Amersham™ Hybond™ ECL nitrocellulose membranes (GE Healthcare). After blocking of the membrane in 3 % nonfat milk in Tris-buffered saline (TBS, pH 8.0, Sigma-Aldrich) containing 0.05 % Tween-20 for 1 h at RT, target-specific primary antibodies diluted 1:1000 in TBS with 0.05 % Tween-20 were added overnight at 4 °C with gentle agitation (Wnt-5a antibody, sc-365370, mouse monoclonal, Santa Cruz Biotechnology, Inc.; Cytokeratin 19 antibody, 11-120-C020, mouse monoclonal, EXBIO Praha a.s.; SOX2 antibody, D1951-1NB-1-1R1M1M3/1C16_110120, mouse monoclonal, ABMART; GAPDH antibody, TA302944-100, goat polyclonal, OriGene; antibodies were pre-tested for adequate performance in different human cell lysates; data not shown). The next day the membranes were incubated in horseradish-peroxidase (HRP)-labeled secondary antibodies (1:10,000, polyclonal rabbit anti-mouse IgG/HRP, polyclonal rabbit anti-goat IgG/HRP, Dako Denmark A/S, Glostrup, Denmark) for 1 h at RT. Blots were developed using Immobilon™ Western Chemiluminescent HRP Substrate (Merck/Millipore, Darmstadt, Germany) and signals detected with a ChemiDoc XRS system (Bio-Rad Laboratories). Band intensities were densitometrically quantified by using the image-processing program ImageJ (http://rsb.info.nih.gov/ij/).
Global transcriptional profiling uncovers 207 differentially expressed genes in the EBS-DM cell line KEB7
Transcriptional profiling using the Affymetrix Human Gene 1.0 ST gene chip platform offers coverage of nearly the entire transcriptome, defined by 28,869 annotated open reading frames/genes. We applied this platform to study the expression pattern of genes influenced by the cytoskeletal collapse caused by the R125P mutation in KRT14 in KEB7 cells. Table 1 summarizes the results of our microarray profiling (four technical replicates each) using a pool of different passages (biological replicates) of KEB7 cells relative to NEB1 wild-type control cells. A LIMMA statistical test comparing mutant versus wild-type cells identified 106 up-regulated and 101 down-regulated genes with significantly altered expression levels (log2-fold change ≥+1 and ≤−1, adjusted p-value <0.05) out of 382 differentially expressed genes (failed at log2-fold change and p-value set).
Bioinformatic enrichment and clustering suggest previously undescribed pathways and functions are affected in EBS-DM
We next assigned cellular themes to the differentially expressed genes by using IBD software. Functional annotation clustering allowed us to group highly related annotation terms into enriched functional categories, as shown in Table 2. The enrichment p-values of each GO term are presented as the ordered p-value of a one-sided hypergeometric test. Significance of enrichment is indicated by a p-value of <0.05. Apart from already well-known clusters such as the regulation of epithelial cell proliferation , cell adhesion and response to stress [12, 13], response to retinoic acid and UV [21–23], immune response  or fatty acid metabolism , our analysis also identified previously undescribed functional classes potentially implicated in EBS-DM. The latter include genes involved in morphogenesis, and genes that regulate BMP- and canonical Wnt-signaling. In addition, two enrichment clusters associated with protein tyrosine phosphatases and symporter activities, and one further cellular component cluster linked to neurofilaments, were uncovered.
A complementary qRT-PCR array approach confirms differential regulation of the majority of selected candidates
Typically in microarray investigations, just a few randomly selected candidates are re-analyzed by qRT-PCR analysis for confirmation of the first hits, but this can lead to a significant number of undetected false positives. In the present study we used an array-based qRT-PCR analysis to investigate 158 of the 207 (76 %) microarray hits. To reduce the number of transcripts to be able to fit our in-house 96-well qRT-PCR array arrangement, we examined 79 randomly selected up-regulated and 79 down-regulated genes (plus controls) in three independent PCR runs (biological and technical replicates). Furthermore, we based our fold-change cut-off values on the average of seven housekeeping genes (rather than on two, as has been customary in the literature [24, 25], so at to improve the reliability of the data obtained). Out of the 158 genes analyzed, significant (p-value <0.05) differences in expression (fold change set ≥+2 for up-regulated and ≤−2 for down-regulated genes) of 121 genes (~76.58 %) was confirmed, as highlighted in Table 1. Eleven genes (~6.96 %) turned out not to be significantly modulated, and 26 genes (~16.45 %) completely failed in all independent PCR runs, probably due to poorly designed primers; this would potentially give an underestimate of the number of reproducible hits.
Western blot analysis confirms differential expression at the protein level of 3 selected candidates
Western blot analysis using crude lysates of KEB7 and NEB1 cells confirmed differential regulation of some of the selected candidates at the protein level after normalization with glyceraldehyde-3-phosphate-dehydrogenase (GAPDH), as shown in Fig. 2. Sex determining region Y-box 2 (SOX2, a log2-fold change of −1.71 by microarray analysis, no data by qRT-PCR) and cytokeratin 19 (KRT19, log2-fold change of −5.73 by microarray analysis, 1633-fold down-regulation by qRT-PCR) were both found to be significantly down-regulated in KEB7 cells, whereas wingless-type MMTV integration site family member 5a (WNT5a) was found to be significantly up-regulated (log2-fold change of 1.94 by microarray analysis, 6.9-fold up-regulation by qRT-PCR). Based on densitometric evaluation of the protein blots, SOX2 is 25-fold down-regulated, KRT19 61.7-fold down-regulated, and WNT5a 123-fold up-regulated in KEB7 cells.
Global transcriptional profiling uncovers novel differentially expressed genes as well as affected pathways in the EBS-DM cell line KEB7
The present study employed global high-throughput microarray analysis together with bioinformatics-assisted functional clustering and enrichment to identify 207 differentially expressed genes in the EBS-DM cell line KEB7. A subset of 158 genes was subjected to further investigation by qRT-PCR, of which three quarters (76.58 %) could be validated as being differentially regulated. Parenthetically, the used human gene 1.0 ST Array from Affymetrix also has the ability to detect alternative mRNA isoforms. The mRNA source of the array transcripts arise from the reference sequence database (RefSeq). Each of these transcripts is read out by 26 different oligonucleotides (25mere), which are spread about the whole length of the transcription unit. Hence, a averaged probe set consists of 26 independent measure points per transcript, which are summated in only one signal. Because we have not taken into account this feature of the array to detect these alternative mRNA isoforms or splice variants we designed for qPCR validation only primer pairs for one transcript variant. Therefore, the number of our reproducible hits could potentially be underestimated.
In addition to verifying genes from previous studies (and thereby supporting their potential role in EBS pathology), a novel contribution of the present work is its uncovering of genes involved in the BMP-signaling pathway, the canonical Wnt-receptor signaling pathway, and the response to retinoic acid, suggesting that dysregulation of these processes as well may play a role in disease pathology. It is well demonstrated that keratinocyte differentiation is closely linked to hormonal action, particularly to calcitriol- [26, 27] and retinoic acid-signaling [21, 23]. In the present study, we identified four differentially regulated genes in KEB7 cells having roles in retinoic acid signaling: three genes, HSD17B2 (hydroxysteroid (17-β) dehydrogenase 2), SOX2 and MEST (mesoderm specific transcript), are down-regulated, and DKK1 (dickkopf Wnt signaling pathway inhibitor 1) is up-regulated. These potential disturbances to retinoic acid signaling may at least partially account for the observations by Wagner et al.  of aberrant differentiation processes in KEB7 cells. Generally, retinoids are considered to play a role in normalization of keratinocyte differentiation by down-regulating desmosomal proteins, exerting anti-proliferative effects, and regulating lipid-synthesis, growth factors, and cytokines [22, 23, 28, 29]. An interesting point in this regard is the ability of retinoic acid to act as an antagonist of Jun N-terminal kinase (JNK)-signaling, mediated by activating protein 1 (AP1) , since JNK/mitogen activated protein kinase (MAPK) signal transduction has previously been shown to be dysregulated in EBS-DM model keratinocytes [9, 31, 32]. Furthermore, we recently advanced a positive feedback model, suggesting that mutant KRT14 leads to an increase in AP1-dependent expression of KRT14, which in turn amplifies the level of aberrant JNK stress-signaling . Concerning the canonical Wnt-receptor signaling pathway, enrichment clustering applied in our present study identified four significantly dysregulated genes: GPC3 (glypican3) and SOX2 are both down-regulated, and WNT5a and PRICKLE1 (prickle homolog 1) are up-regulated. Interestingly, WNT5a is one of the so-called non-canonical Wnt ligands. During normal development, WNT5a is secreted and directs the migration of target cells along concentration gradients. However, deregulated WNT5a signaling facilitates invasion by multiple tumor types into contiguous tissues. EB per se is often associated with carcinoma development, mostly in EB dystrophicans (DEB) and junctionalis (JEB) [34–36], but rare cases of carcinoma in EBS have also been described . To date, the roles of WNT5a in cutaneous squamous cell carcinoma (SCC) and basal cell carcinoma (BCC) as well as the effect of WNT5a on keratinocyte migration have not been fully investigated, although Pourreyron et al.  recently demonstrated up-regulation of WNT5a in SCC/BCC and its localization to the leading edge of tumors as well as in tumor-associated fibroblasts. Regarding BMP signaling, enrichment clustering showed positive regulation of this pathway in EBS by highlighting three genes, GPC3, MSX2 (Msh homeobox 2) and ILK (Integrin-linked kinase), as being significantly down-regulated. BMPs are secreted signaling polypeptides belonging to the transforming growth factor-β (TGF-β) superfamily that act as important multifactor players in the development of vertebrate skin (and appendages). BMPs thereby not only regulate hair follicle (HF) morphogenesis and keratinocyte differentiation in both the HF and epidermis, but are, along with other interacting growth factor family members, also involved in normal postnatal tissue remodeling and homeostasis [38, 39]. Furthermore, BMP signaling acts as a potent tumor suppressor in the skin, inhibiting mainly epidermal- and HF-derived tumor formation [38, 40]. Consequently, dysregulation of this pathway can lead to abnormal skin development and tumor formation. Progress in this research area could be interesting to get a better knowledge of the disease pathology.
KRT19—a novel biomarker in EBS-DM?
A further example from our list of dysregulated genes is the smallest known acidic keratin, KRT19 (significantly down-regulated in our microarray analysis, independently confirmed by both qRT-PCR and western blotting). In our previous study, in which we applied a similar whole-transcriptome microarray methodology to calcitriol-stimulated human primary keratinocytes, we  found a significantly increased expression level of KRT19 after treatment with the indicated secosteroid. Hence, irrespective of the exact role of KRT19 in the disease mechanism, it could constitute another useful biomarker in EBS.
In summary, our study employed global high-throughput microarray analysis plus bioinformatics-assisted functional clustering and enrichment to uncover 207 differentially expressed genes in the EBS-DM cell line. Thereof, 158 genes were subjected to further investigation by a complementary qRT-PCR analysis, of which 76.58 % could be validated as being differentially regulated. We here present a pool of novel candidate genes and potential affected pathways, which may play a role in the disease pathology.
activating protein 1
bone morphogenetic protein
epidermolysis bullosa simplex Dowling-Meara
false discovery rate
Jun N-terminal kinase
linear model for microarray data
JNK/mitogen-activated protein kinase
quantitative real-time polymerase chain reaction
cutaneous squamous cell carcinoma
sex-determining region Y-box 2
wingless-type MMTV integration site family5a
Msh homeobox 2
Moll R, Divo M, Langbein L. The human keratins: biology and pathology. Histochem Cell Biol. 2008;129(6):705–33.
Wong P, Coulombe PA. Loss of keratin 6 (K6) proteins reveals a function for intermediate filaments during wound repair. J Cell Biol. 2003;163(2):327–37.
Fine JD, Eady RA, Bauer EA, Bauer JW, Bruckner-Tuderman L, Heagerty A, Hintner H, Hovnanian A, Jonkman MF, Leigh I, et al. The classification of inherited epidermolysis bullosa (EB): report of the third international consensus meeting on diagnosis and classification of EB. J Am Acad Dermatol. 2008;58(6):931–50.
Letai A, Coulombe PA, McCormick MB, Yu QC, Hutton E, Fuchs E. Disease severity correlates with position of keratin point mutations in patients with epidermolysis bullosa simplex. Proc Natl Acad Sci USA. 1993;90(8):3197–201.
Bonifas JM, Rothman AL, Epstein EH Jr. Epidermolysis bullosa simplex: evidence in two families for keratin gene abnormalities. Science. 1991;254(5035):1202–5.
Coulombe PA, Hutton ME, Letai A, Hebert A, Paller AS, Fuchs E. Point mutations in human keratin 14 genes of epidermolysis bullosa simplex patients: genetic and functional analyses. Cell. 1991;66(6):1301–11.
Irvine AD, McLean WH. Human keratin diseases: the increasing spectrum of disease and subtlety of the phenotype-genotype correlation. Br J Dermatol. 1999;140(5):815–28.
Lane EB, Rugg EL, Navsaria H, Leigh IM, Heagerty AH, Ishida-Yamamoto A, Eady RA. A mutation in the conserved helix termination peptide of keratin 5 in hereditary skin blistering. Nature. 1992;356(6366):244–6.
D’Alessandro M, Russell D, Morley SM, Davies AM, Lane EB. Keratin mutations of epidermolysis bullosa simplex alter the kinetics of stress response to osmotic shock. J Cell Sci. 2002;115(Pt 22):4341–51.
Russell D, Ross H, Lane EB. ERK involvement in resistance to apoptosis in keratinocytes with mutant keratin. J Invest Dermatol. 2010;130(3):671–81.
Lu H, Chen J, Planko L, Zigrino P, Klein-Hitpass L, Magin TM. Induction of inflammatory cytokines by a keratin mutation and their repression by a small molecule in a mouse model for EBS. J Invest Dermatol. 2007;127(12):2781–9.
Liovic M, Lee B, Tomic-Canic M, D’Alessandro M, Bolshakov VN, Lane EB. Dual-specificity phosphatases in the hypo-osmotic stress response of keratin-defective epithelial cell lines. Exp Cell Res. 2008;314(10):2066–75.
Liovic M, D’Alessandro M, Tomic-Canic M, Bolshakov VN, Coats SE, Lane EB. Severe keratin 5 and 14 mutations induce down-regulation of junction proteins in keratinocytes. Exp Cell Res. 2009;315(17):2995–3003.
Bchetnia M, Tremblay ML, Leclerc G, Duperee A, Powell J, McCuaig C, Morin C, Legendre-Guillemin V, Laprise C. Expression signature of epidermolysis bullosa simplex. Hum Genet. 2012;131:393–406.
Wagner M, Hintner H, Bauer JW, Onder K. Gene expression analysis of an epidermolysis bullosa simplex Dowling-Meara cell line by subtractive hybridization: recapitulation of cellular differentiation, migration and wound healing. Exp Dermatol. 2012;21(2):111–7.
Morley SM, D’Alessandro M, Sexton C, Rugg EL, Navsaria H, Shemanko CS, Huber M, Hohl D, Heagerty AI, Leigh IM, et al. Generation and characterization of epidermolysis bullosa simplex cell lines: scratch assays show faster migration with disruptive keratin mutations. Br J Dermatol. 2003;149(1):46–58.
Smyth GK. Linear models and empirical bayes methods for assessing differential expression in microarray experiments. Stat Appl Genet Mol Biol. 2004;3:Article3.
Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) method. Methods. 2001;25(4):402–8.
Benjamini Y, Drai D, Elmer G, Kafkafi N, Golani I. Controlling the false discovery rate in behavior genetics research. Behav Brain Res. 2001;125(1–2):279–84.
Alam H, Sehgal L, Kundu ST, Dalal SN, Vaidya MM. Novel function of keratins 5 and 14 in proliferation and differentiation of stratified epithelial cells. Mol Biol Cell. 2011;22(21):4068–78.
Green H, Watt FM. Regulation by vitamin A of envelope cross-linking in cultured keratinocytes derived from different human epithelia. Mol Cell Biol. 1982;2(9):1115–7.
Gatto H, Richard MH, Viac J, Charveron M, Schmitt D. Effects of retinoic acid on interleukin-1 alpha and -1 beta expression by normal human keratinocytes cultured in defined medium. Skin Pharmacol. 1993;6(1):10–9.
Torma H. Regulation of keratin expression by retinoids. Dermatoendocrinology. 2011;3(3):136–40.
de Jonge HJ, Fehrmann RS, de Bont ES, Hofstra RM, Gerbens F, Kamps WA, de Vries EG, van der Zee AG, te Meerman GJ, ter Elst A. Evidence based selection of housekeeping genes. PLoS One. 2007;2(9):e898.
Gilsbach R, Kouta M, Bonisch H, Bruss M. Comparison of in vitro and in vivo reference genes for internal standardization of real-time PCR data. Biotechniques. 2006;40(2):173–7.
Bikle DD. Vitamin D regulated keratinocyte differentiation. J Cell Biochem. 2004;92(3):436–44.
Bikle DD, Ng D, Tu CL, Oda Y, Xie Z. Calcium- and vitamin D-regulated keratinocyte differentiation. Mol Cell Endocrinol. 2001;177(1–2):161–71.
Humphries JD, Parry EJ, Watson RE, Garrod DR, Griffiths CE. All-trans retinoic acid compromises desmosome expression in human epidermis. Br J Dermatol. 1998;139(4):577–84.
Xiao JH, Feng X, Di W, Peng ZH, Li LA, Chambon P, Voorhees JJ. Identification of heparin-binding EGF-like growth factor as a target in intercellular regulation of epidermal basal cell growth by suprabasal retinoic acid receptors. EMBO J. 1999;18(6):1539–48.
Schule R, Rangarajan P, Yang N, Kliewer S, Ransone LJ, Bolado J, Verma IM, Evans RM. Retinoic acid is a negative regulator of AP-1-responsive genes. Proc Natl Acad Sci USA. 1991;88(14):6092–6.
Chamcheu JC, Navsaria H, Pihl-Lundin I, Liovic M, Vahlquist A, Torma H. Chemical chaperones protect epidermolysis bullosa simplex keratinocytes from heat stress-induced keratin aggregation: involvement of heat shock proteins and MAP kinases. J Invest Dermatol. 2011;131:1684–91.
Morley SM, Dundas SR, James JL, Gupta T, Brown RA, Sexton CJ, Navsaria HA, Leigh IM, Lane EB. Temperature sensitivity of the keratin cytoskeleton and delayed spreading of keratinocyte lines derived from EBS patients. J Cell Sci. 1995;108(Pt 11):3463–71.
Wagner M, Trost A, Hintner H, Bauer JW, Onder K. Imbalance of intermediate filament component keratin 14 contributes to increased stress signalling in epidermolysis bullosa simplex. Exp Dermatol. 2013;22(4):292–4.
Mallipeddi R. Epidermolysis bullosa and cancer. Clin Exp Dermatol. 2002;27(8):616–23.
Venugopal SS, Murrell DF. Treatment of skin cancers in epidermolysis bullosa. Dermatol Clin. 2010;28(2):283–7 (ix–x).
Schumann H, Roth W, Has C, Volz A, Erfurt-Berge C, Magin TM, Bruckner-Tuderman L. Verrucous carcinoma in epidermolysis bullosa simplex is possibly associated with a novel mutation in the keratin 5 gene. Br J Dermatol. 2012;167(4):929–36.
Pourreyron C, Reilly L, Proby C, Panteleyev A, Fleming C, McLean K, South AP, Foerster J. Wnt5a is strongly expressed at the leading edge in non-melanoma skin cancer, forming active gradients, while canonical Wnt signalling is repressed. PLoS One. 2012;7(2):e31827.
Sharov AA, Mardaryev AN, Sharova TY, Grachtchouk M, Atoyan R, Byers HR, Seykora JT, Overbeek P, Dlugosz A, Botchkarev VA. Bone morphogenetic protein antagonist noggin promotes skin tumorigenesis via stimulation of the Wnt and Shh signaling pathways. Am J Pathol. 2009;175(3):1303–14.
Sharov AA, Sharova TY, Mardaryev AN, Tommasi di Vignano A, Atoyan R, Weiner L, Yang S, Brissette JL, Dotto GP, Botchkarev VA. Bone morphogenetic protein signaling regulates the size of hair follicles and modulates the expression of cell cycle-associated genes. Proc Natl Acad Sci USA. 2006;103(48):18166–71.
Rishikaysh P, Dev K, Diaz D, Qureshi WM, Filip S, Mokry J. Signaling involved in hair follicle morphogenesis and development. Int J Mol Sci. 2014;15(1):1647–70.
Rid R, Wagner M, Maier CJ, Hundsberger H, Hintner H, Bauer JW, Onder K. Deciphering the calcitriol-induced transcriptomic response in keratinocytes: presentation of novel target genes. J Mol Endocrinol. 2013;50(2):131–49.
JH performed the research, analyzed the data and wrote the manuscript, KO designed the research study, MW and RR contributed their expertise in data analysis and condition of the manuscript, and JB, HH, and AE critically reviewed the manuscript. All authors read and approved the final manuscript.
This work was supported by the Dystrophic Epidermolysis Bullosa Research Association (DEBRA), Austria and Procomcure Biotech GmbH, Austria.
The authors declare that they have no competing interests.
Additional file 1: Table S1. Literature survey (2007 to date) focusing on differential gene expression profiling in EB, especially EBS. Overview of already identified deregulated candidate genes which might play a role in the severe pathologic phenotype observed in patients with this disease.
Additional file 2: Table S2. Primer sequences of genes regulated in the EBS-DM model keratinocyte cell line KEB7 identified by microarray analysis (r = revers).
About this article
Cite this article
Herzog, J., Rid, R., Wagner, M. et al. Whole-transcriptome gene expression profiling in an epidermolysis bullosa simplex Dowling-Meara model keratinocyte cell line uncovered novel, potential therapeutic targets and affected pathways. BMC Res Notes 8, 785 (2015). https://doi.org/10.1186/s13104-015-1783-7