Signaling sphingolipids are biomarkers for atopic dermatitis prone to disseminated viral infections.

BACKGROUND: Life-threatening viral diseases such as eczema herpeticum (EH) and eczema vaccinatum (EV) occur in <5% of individuals with atopic dermatitis (AD). The diagnosis of AD, however, excludes all individuals with AD from smallpox vaccination. OBJECTIVES: We sought to identify circulatory and skin lipid biomarkers associated with EH and EV. METHODS: Stratum corneum and plasma samples from 15 subjects with AD and a history of EH, 13 age- and gender-matched subjects with AD and without EH history, and 13 healthy nonatopic (NA) controls were analyzed by liquid chromatography tandem mass spectrometry for sphingolipid content. Sphingosine-1-phosphate (S1P) and ceramide levels were validated in plasma samples from the Atopic Dermatitis Vaccinia Network/Atopic Dermatitis Research Network repository (12 NA, 12 AD, 23 EH) and plasma from 7 subjects with EV and 7 matched subjects with AD. S1P lyase was downregulated in human primary keratinocytes to evaluate its effect on herpes simplex virus 1 (HSV-1) replication in vitro. RESULTS: The stratum corneum of patients with EH demonstrated significantly higher levels of free sphingoid bases than those in patients who were NA, indicating enhanced sphingolipid turnover in keratinocytes (P < .05). Plasma from 2 independent cohorts of patients with EH had a significantly increased S1P/ceramide ratio in subjects with EH versus those with AD and or who were NA (P < .01). The S1P level in plasma from subjects with EV was twice the level in plasma from subjects with AD (mean = 1,533 vs 732 pmol/mL; P < .001). Downregulation of S1P lyase expression with silencing RNA led to an increased S1P level and doubled HSV-1 titer in keratinocytes. CONCLUSIONS: Our data point to long-term abnormalities in the S1P signaling system as a biomarker for previous disseminated viral diseases and a potential treatment target in recurring infections.

Dupilumab significantly improves skin barrier function in patients with moderate-to-severe atopic dermatitis.

BACKGROUND: Atopic dermatitis (AD) is characterized by abnormal skin lipids that are largely driven by hyperactivated type 2 immune responses. The antibody to the α-subunit of interleukin (IL)-4 receptor, dupilumab, was recently approved to treat AD and demonstrated strong efficacy. However, the role of dupilumab therapy in the regulation of skin barrier structure and function has not been fully explored. METHODS: We have evaluated the content of lipids and transepidermal water loss (TEWL) in lesional and non-lesional skin of adults and adolescents with moderate-to-severe AD over the course of 16-week treatment with dupilumab and compared those values with that of matched healthy volunteers. RESULTS: Dupilumab treatment provided a significant decrease in TEWL in AD lesions, lowering it almost to the levels seen in the skin of healthy subjects. Blocking IL-4/IL-13 signaling with dupilumab normalized lipid composition (decreased levels of ceramides with non-hydroxy fatty acids and C18-sphingosine and increased the level of esterified omega-hydroxy fatty acid-containing ceramides) and increased ceramide chain length in lesional as well as non-lesional stratum corneum of AD patients. Partial changes for these parameters were already observed after 2 weeks, with a full response achieved after 8 weeks of dupilumab treatment. CONCLUSIONS: Inhibition of IL-4/IL-13 signaling by dupilumab allows restoration of skin lipid composition and barrier function in patients with moderate-to-severe AD.

Whole genome sequencing identifies novel genetic mutations in patients with eczema herpeticum.

BACKGROUND: Eczema herpeticum (EH) is a rare complication of atopic dermatitis (AD) caused by disseminated herpes simplex virus (HSV) infection. The role of rare and/or deleterious genetic variants in disease etiology is largely unknown. This study aimed to identify genes that harbor damaging genetic variants associated with HSV infection in AD with a history of recurrent eczema herpeticum (ADEH+). METHODS: Whole genome sequencing (WGS) was performed on 49 recurrent ADEH+ (≥3 EH episodes), 491 AD without a history of eczema herpeticum (ADEH-) and 237 non-atopic control (NA) subjects. Variants were annotated, and a gene-based approach (SKAT-O) was used to identify genes harboring damaging genetic variants associated with ADEH+. Genes identified through WGS were studied for effects on HSV responses and keratinocyte differentiation. RESULTS: Eight genes were identified in the comparison of recurrent ADEH+to ADEH-and NA subjects: SIDT2, CLEC7A, GSTZ1, TPSG1, SP110, RBBP8NL, TRIM15, and FRMD3. Silencing SIDT2 and RBBP8NL in normal human primary keratinocytes (NHPKs) led to significantly increased HSV-1 replication. SIDT2-silenced NHPKs had decreased gene expression of IFNk and IL1b in response to HSV-1 infection. RBBP8NL-silenced NHPKs had decreased gene expression of IFNk, but increased IL1b. Additionally, silencing SIDT2 and RBBP8NL also inhibited gene expression of keratinocyte differentiation markers keratin 10 (KRT10) and loricrin (LOR). CONCLUSION: SIDT2 and RBBP8NL participate in keratinocyte's response to HSV-1 infection. SIDT2 and RBBP8NL also regulate expression of keratinocyte differentiation genes of KRT10 and LOR.

Ankyrin repeat domain 1 regulates innate immune responses against herpes simplex virus 1: A potential role in eczema herpeticum.

BACKGROUND: Atopic dermatitis (AD) is a common inflammatory skin disease. A subset of patients with AD are susceptible to disseminated herpes simplex virus (HSV) infection, a complication termed eczema herpeticum (ADEH+). The immune mechanisms causing ADEH+ remain elusive. Using RNA sequencing, we recently found that ankyrin repeat domain 1 (ANKRD1) was significantly induced in human PBMCs upon HSV-1 stimulation, and its induction in patients with ADEH+ was significantly reduced compared with that seen in AD patients without a history of eczema herpeticum (ADEH-). OBJECTIVE: We sought to validate ANKRD1 gene expression in nonatopic (NA) subjects, patients with ADEH-, and patients with ADEH+ and to delineate the biological function of ANKRD1 and the signaling pathway or pathways involved. METHODS: Purification of human PBMCs, monocytes, B cells, dendritic cells, T cells, and natural killer cells; RNA extraction and quantitative RT-PCR; small interfering RNA technique; co-immunoprecipitation; and Western blot assays were used. RESULTS: ANKRD1 expression was significantly reduced in PBMCs from patients with ADEH+ after HSV-1 stimulation compared with PBMCs from patients with ADEH-. We found that the induction of ANKRD1 by HSV-1 and multiple pattern recognition receptor agonists are mediated by inflammatory cytokines. Silencing ANKRD1 gene expression in antigen-presenting cells led to increased viral load and reduced IFNB1 and IL29 production. Using co-immunoprecipitation methods, we demonstrated that ANKRD1 formed protein complexes with interferon regulatory factor (IRF) 3 and IRF7, which are important transcription factors regulating signaling transduction of pattern recognition receptors. Overexpression of ANKRD1 enhanced the IRF3-mediated signaling pathways. CONCLUSION: ANKRD1 is involved in IRF3-mediated antiviral innate immune signaling pathways. Its reduced expression in patients with ADEH+ might contribute to the pathogenesis of ADEH+.

Minimally invasive skin tape strip RNA sequencing identifies novel characteristics of the type 2-high atopic dermatitis disease endotype.

BACKGROUND: Expression profiling of skin biopsy specimens has established molecular features of the skin in patients with atopic dermatitis (AD). The invasiveness of biopsies has prevented their use in defining individual-level AD pathobiological mechanisms (endotypes) in large research studies. OBJECTIVE: We sought to determine whether minimally invasive skin tape strip transcriptome analysis identifies gene expression dysregulation in AD and molecular disease endotypes. METHODS: We sampled nonlesional and lesional skin tape strips and biopsy specimens from white adult patients with AD (18 male and 12 female patients; age [mean ± SE], 36.3 ± 2.2 years) and healthy control subjects (9 male and 16 female subjects; age [mean ± SE], 34.8 ± 2.2 years). AmpliSeq whole-transcriptome sequencing was performed on extracted RNA. Differential expression, clustering/pathway analyses, immunostaining of skin biopsy specimens, and clinical trait correlations were performed. RESULTS: Skin tape expression profiles were distinct from skin biopsy profiles and better sampled epidermal differentiation complex genes. Skin tape expression of 29 immune and epidermis-related genes (false discovery rate < 5%) separated patients with AD from healthy subjects. Agnostic gene set analyses and clustering revealed 50% of patients with AD exhibited a type 2 inflammatory signature (type 2-high endotype) characterized by differential expression of 656 genes, including overexpression of IL13, IL4R, CCL22, CCR4 (log2 fold change = 5.5, 2.0, 4.0, and 4.1, respectively) and at a pathway level by TH2/dendritic cell activation. Both expression and immunostaining of skin biopsy specimens indicated this type 2-high group was enriched for inflammatory, type 2-skewed dendritic cells expressing FcεRI. The type 2-high endotype group exhibited more severe disease by using both the Eczema Area and Severity Index score and body surface area covered by lesions. CONCLUSION: Minimally invasive expression profiling of nonlesional skin reveals stratification in AD molecular pathology by type 2 inflammation that correlates with disease severity.

Identification of novel gene signatures in patients with atopic dermatitis complicated by eczema herpeticum.

BACKGROUND: A subset of patients with atopic dermatitis (AD) is prone to disseminated herpes simplex virus (HSV) infection (ie, atopic dermatitis with a history of eczema herpeticum [ADEH+]). Biomarkers that identify ADEH+ are lacking. OBJECTIVE: We sought to search for novel ADEH+ gene signatures in PBMCs. METHODS: An RNA-sequencing approach was applied to evaluate global transcriptional changes by using PBMCs from patients with ADEH+ and patients with atopic dermatitis without a history of eczema herpeticum (ADEH-). Candidate genes were confirmed by means of quantitative PCR or ELISA. RESULTS: PBMCs from patients with ADEH+ had distinct changes to the transcriptome when compared with those from patients with ADEH- after HSV-1 stimulation: 792 genes were differentially expressed at a false discovery rate of less than 0.05 (ANOVA), and 15 type I and type III interferon genes were among the top 20 most downregulated genes in patients with ADEH+. We further validated that IFN-α and IL-29 mRNA and protein levels were significantly decreased in HSV-1-stimulated PBMCs from patients with ADEH+ compared with those from patients with ADEH- and healthy subjects. Ingenuity Pathway Analysis demonstrated that the upstream regulators of type I and type III interferons, interferon regulatory factor (IRF) 3 and IRF7, were significantly inhibited in patients with ADEH+ based on the downregulation of their target genes. Furthermore, we found that gene expression of IRF3 and IRF7 was significantly decreased in HSV-1-stimulated PBMCs from patients with ADEH+. CONCLUSIONS: PBMCs from patients with ADEH+ have a distinct immune response after HSV-1 exposure compared with those from patients with ADEH-. Inhibition of the IRF3 and IRF7 innate immune pathways in patients with ADEH+ might be an important mechanism for increased susceptibility to disseminated viral infection.

Vitamin D inhibits monocyte/macrophage proinflammatory cytokine production by targeting MAPK phosphatase-1.

It is estimated that 1 billion people around the world are vitamin D deficient. Vitamin D deficiency has been linked to various inflammatory diseases. However, the mechanism by which vitamin D reduces inflammation remains poorly understood. In this study, we investigated the inhibitory effects of physiologic levels of vitamin D on LPS-stimulated inflammatory response in human blood monocytes and explored potential mechanisms of vitamin D action. We observed that two forms of the vitamin D, 1,25(OH)(2)D(3), and 25(OH)D(3), dose dependently inhibited LPS-induced p38 phosphorylation at physiologic concentrations, IL-6 and TNF-α production by human monocytes. Upon vitamin D treatment, the expression of MAPK phosphatase-1 (MKP-1) was significantly upregulated in human monocytes and murine bone marrow-derived macrophages (BMM). Increased binding of the vitamin D receptor and increased histone H4 acetylation at the identified vitamin D response element of the murine and human MKP-1 promoters were demonstrated. Moreover, in BMM from MKP1(-/-) mice, the inhibition of LPS-induced p38 phosphorylation by vitamin D was completely abolished. Vitamin D inhibition of LPS-induced IL-6 and TNF-α production by BMM from MKP-1(-/-) mice was significantly reduced as compared with wild-type mice. In conclusion, this study identified the upregulation of MKP-1 by vitamin D as a novel pathway by which vitamin D inhibits LPS-induced p38 activation and cytokine production in monocytes/macrophages.

Steroid requirements and immune associations with vitamin D are stronger in children than adults with asthma.

BACKGROUND: The effects of serum vitamin D status on atopy, steroid requirement, and functional responsiveness to corticosteroids in children versus adults with asthma have not been studied systematically. OBJECTIVE: We sought to explore the age-specific effects of vitamin D in asthmatic patients. METHODS: Serum vitamin D levels were examined in a prospective study of adults and children (102 healthy control subjects and 103 asthmatic patients). PBMCs were cultured for 3 hours with or without 100 nmol/L dexamethasone, and the expression of corticosteroid-regulated genes was detected by using real-time PCR. Serum IgE levels were measured, and information about asthmatic patients' steroid requirements was collected. RESULTS: Deficient serum vitamin D levels (<20 ng/mL) were found in 47.6% of asthmatic patients and 56.8% of healthy control subjects, with means ± SDs of 20.7 ± 9.8 and 19.2 ± 7.7 ng/mL, respectively. In multivariate regression models a significant positive correlation between serum vitamin D levels and the expression of vitamin D-regulated targets, cytochrome P450, family 24, subfamily a (cyp24a) expression by PBMCs (P = .0084, pediatric asthma group only) and serum LL-37 levels (P = .0006 in the pediatric group but P = .0067 in the adult asthma group), was found. An inverse association between vitamin D and serum IgE levels was observed in the pediatric (P = .006) asthma group. Serum vitamin D level (P = .05), as well as PBMC cyp24a expression (P = .0312), demonstrated a significant inverse relationship with daily inhaled corticosteroid dose in the pediatric asthma group only. Cyp24a expression in PBMCs correlated positively with in vitro suppression of TNF-α by dexamethasone (P = .05) and IL-13 (P = .0094) in PBMCs in the pediatric asthma group only. CONCLUSIONS: This study demonstrated significant associations between serum vitamin D status and steroid requirement and in vitro responsiveness to corticosteroids in the pediatric but not the adult asthma group. Vitamin D was also related to IgE levels in children but not in adults.

Atopic dermatitis complicated by recurrent eczema herpeticum is characterized by multiple, concurrent epidermal inflammatory endotypes.

BACKGROUND: A subgroup of atopic dermatitis (AD) patients suffer from recurrent, disseminated herpes simplex virus (HSV) skin infections, termed eczema herpeticum (EH), which can be life-threatening and contribute to AD morbidity. The pathobiology underlying ADEH is unknown. OBJECTIVE: To determine transcriptional mechanisms of skin and immune system pathobiology that underlie ADEH disease. METHODS: We performed whole transcriptome RNA-sequencing of non-lesional skin samples (epidermis, dermis) of AD patients with (ADEH + , n=15) and without (ADEH - , n=13) recurrent EH history, and healthy controls (HC, n=15). We also performed RNA-sequencing on plasmacytoid dendritic cells (pDCs) collected from these participants and infected in vitro with HSV-1. Differential expression, gene set enrichment, and endotyping analyses were performed. RESULTS: ADEH + disease was characterized by dysregulation in skin gene expression, which was limited in dermis (differentially expressed genes [DEGs]=14) and widespread in epidermis (DEGs=129). ADEH + -upregulated epidermal DEGs were enriched in type 2 cytokine (T2) ( IL4R, CCL22, CRLF2, IL7R ), interferon ( CXCL10, ICAM1, IFI44 , and IRF7) , and IL-36γ ( IL36G ) inflammatory pathway genes. At a person-level, all ADEH + participants exhibited T2 and interferon endotypes and 87% were IL36G-high. In contrast, these endotypes were more variably expressed among ADEH - participants. ADEH + patient skin also exhibited dysregulation in epidermal differentiation complex (EDC) genes within the LCE, S100 , and SPRR families, which are involved in skin barrier function, inflammation, and antimicrobial activities. pDC transcriptional responses to HSV-1 infection were not altered by ADEH status. CONCLUSIONS: ADEH + pathobiology is characterized by a unique, multi-faceted epidermal inflammation that accompanies dysregulation in the expression of EDC genes. Key Messages: AD patients with a history of recurrent EH exhibit molecular skin pathobiology that is similar in form, but more severe in degree, than in AD patients without this complication. Non-lesional skin of ADEH + patients concurrently exhibits excessive type 2 cytokine, interferon, and IL-36γ-driven epidermal inflammation. Expression of these inflammatory skin endotypes among ADEH + patients is associated with dysregulation in expression of epidermal differentiation complex genes involved in barrier function, inflammation, and antimicrobial activity. Capsule Summary: AD patients with a history of recurrent disseminated HSV-1 skin infections form a unique molecular skin endotype group that concurrently exhibits type 2 cytokine, interferon, and IL-36γ-driven skin inflammation, accompanied by dysregulation in expression of epidermal differentiation complex genes involved in barrier function, inflammation, and antimicrobial activity.

Dupilumab Inhibits Vascular Leakage of Blood Proteins Into Atopic Dermatitis Skin.

BACKGROUND: Atopic dermatitis (AD) skin lesions are associated with oozing, bleeding, and erythema. This suggests that AD is associated with vascular changes. Dupilumab is an antibody to the alpha subunit of IL-4 receptor that demonstrates strong efficacy in the treatment of AD. IL-4 is known to reduce the permeability barrier function of vascular endothelium. OBJECTIVE: To examine the effects of dupilumab on vascular barrier function in AD skin. METHODS: Using proteomic analysis, we evaluated the plasma protein composition in skin tapes of lesional and nonlesional skin of adults and adolescents with moderate to severe AD over the course of a 16-week treatment with dupilumab and compared those with matched healthy subjects. RESULTS: At baseline, 115 plasma proteins were detected in AD skin and globally increased (1.5-fold or greater) compared with healthy skin. Functionally, these proteins included immunoglobulins, proteins involved in the coagulation process, enzymes, protease inhibitors, transport proteins, acute-phase proteins, complement proteins, and other pleiotropic proteins. Noteworthy, fibrinogens, fibronectin, and heme-binding proteins haptoglobin and hemopexin were among the top proteins originating from plasma and were increased in AD lesional versus healthy skin at baseline (P < .0001). Dupilumab treatment resulted in significantly reduced levels of plasma proteins in AD skin (P < .0001), with most dropping to levels seen in healthy skin or no longer detectable at week 16. CONCLUSIONS: Inhibition of IL-4/IL-13 action by dupilumab significantly reduces the efflux of plasma proteins into AD skin. Several of these proteins, such as fibrinogens and fibronectin, are known to enhance Staphylococcus aureus colonization and are associated with AD skin severity.

Interferon Kappa Is Important for Keratinocyte Host Defense against Herpes Simplex Virus-1.

Type I interferon kappa (IFNκ) is selectively expressed in human keratinocytes. Herpes simplex virus-1 (HSV-1) is a human pathogen that infects keratinocytes and causes lytic skin lesions. Whether IFNκ plays a role in keratinocyte host defense against HSV-1 has not been investigated. In this study, we found that IFNκ mRNA expression was induced by addition of recombinant IFNκ and poly (I:C); and its expression level was significantly greater than IFNa2, IFNb1, and IFNL1 in both undifferentiated and differentiated normal human epidermal keratinocytes (NHEKs) under resting and stimulation conditions. Although IFNe was expressed at a relatively higher level than other IFNs in resting undifferentiated NHEK, its expression level did not change after stimulation with recombinant IFNκ and poly (I:C). HSV-1 infection inhibited gene expression of IFNκ and IFNe in NHEK. Silencing IFNκ in NHEK led to significantly enhanced HSV-1 replication in both undifferentiated and differentiated NHEK compared to scrambled siRNA-transfected cells, while the addition of recombinant IFNκ significantly reduced HSV-1 replication in NHEK. In addition, we found that IFNκ did not regulate protein expression of NHEK differentiation markers. Our results demonstrate that IFNκ is the dominant type of IFNs in keratinocytes and it has an important function for keratinocytes to combat HSV-1 infection.

The nonlesional skin surface distinguishes atopic dermatitis with food allergy as a unique endotype.

Skin barrier dysfunction has been reported in both atopic dermatitis (AD) and food allergy (FA). However, only one-third of patients with AD have FA. The purpose of this study was to use a minimally invasive skin tape strip sampling method and a multiomics approach to determine whether children with AD and FA (AD FA+) have stratum corneum (SC) abnormalities that distinguish them from AD without FA (AD FA-) and nonatopic (NA) controls. Transepidermal water loss was found to be increased in AD FA+. Filaggrin and the proportion of ω-hydroxy fatty acid sphingosine ceramide content in nonlesional skin of children with AD FA+ were substantially lower than in AD FA- and NA skin. These abnormalities correlated with morphologic changes in epidermal lamellar bilayer architecture responsible for barrier homeostasis. Shotgun metagenomic studies revealed that the nonlesional skin of AD FA+ had increased abundance of Staphylococcus aureus compared to NA. Increased expression of keratins 5, 14, and 16 indicative of hyperproliferative keratinocytes was observed in the SC of AD FA+. The skin transcriptome of AD FA+ had increased gene expression for dendritic cells and type 2 immune pathways. A network analysis revealed keratins 5, 14, and 16 were positively correlated with AD FA+, whereas filaggrin breakdown products were negatively correlated with AD FA+. These data suggest that the most superficial compartment of nonlesional skin in AD FA+ has unique properties associated with an immature skin barrier and type 2 immune activation.

Correction: Forkhead Box C1 Regulates Human Primary Keratinocyte Terminal Differentiation.

[This corrects the article DOI: 10.1371/journal.pone.0167392.].

Skin Wound Healing Is Accelerated by a Lipid Mixture Representing Major Lipid Components of Chamaecyparis obtusa Plant Extract.

In chronic nonhealing wounds, the healing process is disrupted and wounds are often infected with bacteria. About 85% of lower extremity amputations in diabetes are attributed to deep infection of foot ulcers. Therefore, infection control is critical for wound care. In this study, we analyzed lipid composition of Chamaecyparis obtusa extract, and we describe the wound-healing properties of its combination of 10 major lipid components. A 10-lipid mixture up-regulated HBD-3 and LL-37 through the olfactory receptor 2AT4 and induced phosphorylation of extracellular signal-regulated kinases and p38 mitogen-activated protein kinases in primary human keratinocytes. In addition, the 10-lipid mixture had direct bactericidal effects against Staphylococcus aureus and Streptococcus pyogenes and protected against staphylococcal α-toxin-induced keratinocyte cell death. In an animal model, the 10-lipid mixture accelerated skin wound healing and was also effective in healing wounds superinfected with S. aureus. We suggest that the 10-lipid mixture, because of its wound-healing and antimicrobial properties, can be beneficial for wound treatment.

Lipid abnormalities in atopic skin are driven by type 2 cytokines.

Lipids in the stratum corneum of atopic dermatitis (AD) patients differ substantially in composition from healthy subjects. We hypothesized that hyperactivated type 2 immune response alters AD skin lipid metabolism. We have analyzed stratum corneum lipids from nonlesional and lesional skin of AD subjects and IL-13 skin-specific Tg mice. We also directly examined the effects of IL-4/IL-13 on human keratinocytes in vitro. Mass spectrometric analysis of lesional stratum corneum from AD subjects and IL-13 Tg mice revealed an increased proportion of short-chain (N-14:0 to N-24:0) NS ceramides, sphingomyelins, and 14:0-22:0 lysophosphatidylcholines (14:0-22:0 LPC) with a simultaneous decline in the proportion of corresponding long-chain species (N-26:0 to N-32:0 sphingolipids and 24:0-30:0 LPC) when compared with healthy controls. An increase in short-chain LPC species was also observed in nonlesional AD skin. Similar changes were observed in IL-4/IL-13-driven responses in Ca2+-differentiated human keratinocytes in vitro, all being blocked by STAT6 silencing with siRNA. RNA sequencing analysis performed on stratum corneum of AD as compared with healthy subjects identified decreased expression of fatty acid elongases ELOVL3 and ELOVL6 that contributed to observed changes in atopic skin lipids. IL-4/IL-13 also inhibited ELOVL3 and ELOVL6 expression in keratinocyte cultures in a STAT6-dependent manner. Downregulation of ELOVL3/ELOVL6 expression in keratinocytes by siRNA decreased the proportion of long-chain fatty acids globally and in sphingolipids. Thus, our data strongly support the pathogenic role of type 2 immune activation in AD skin lipid metabolism.

Forkhead Box C1 Regulates Human Primary Keratinocyte Terminal Differentiation.

The epidermis serves as a critical protective barrier between the internal and external environment of the human body. Its remarkable barrier function is established through the keratinocyte (KC) terminal differentiation program. The transcription factors specifically regulating terminal differentiation remain largely unknown. Using a RNA-sequencing (RNA-seq) profiling approach, we found that forkhead box c 1 (FOXC1) was significantly up-regulated in human normal primary KC during the course of differentiation. This observation was validated in human normal primary KC from several different donors and human skin biopsies. Silencing FOXC1 in human normal primary KC undergoing differentiation led to significant down-regulation of late terminal differentiation genes markers including epidermal differentiation complex genes, keratinization genes, sphingolipid/ceramide metabolic process genes and epidermal specific cell-cell adhesion genes. We further demonstrated that FOXC1 works down-stream of ZNF750 and KLF4, and upstream of GRHL3. Thus, this study defines FOXC1 as a regulator specific for KC terminal differentiation and establishes its potential position in the genetic regulatory network.