Involucrin Modulates Vitamin D Receptor Activity in the Epidermis.

Terminally differentiated keratinocytes are critical for epidermal function and are surrounded by involucrin (IVL). Increased IVL expression is associated with a near-selective sweep in European populations compared with those in Africa. This positive selection for increased IVL in the epidermis identifies human adaptation outside of Africa. The functional significance is unclear. We hypothesize that IVL modulates the environmentally sensitive vitamin D receptor (VDR) in the epidermis. We investigated VDR activity in Ivl‒/‒ and wild-type mice using vitamin D agonist (MC903) treatment and comprehensively determined the inflammatory response using single-cell RNA sequencing and associated skin microbiome changes using 16S bacterial phylotyping. VDR activity and target gene expression were reduced in Ivl‒/‒ mouse skin, with decreased MC903-mediated skin inflammation and significant reductions in CD4+ T cells, basophils, macrophages, monocytes, and type II basal keratinocytes and an increase in suprabasal keratinocytes. Coinciding with the dampened MC903-mediated inflammation, the skin microbiota of Ivl‒/‒ mice was more stable than that of the wild-type mice, which exhibited an MC903-responsive increase in Bacteroidetes and a decrease in Firmicutes. Together, our studies in Ivl‒/‒ mice identify a functional role for IVL to positively impact VDR activity and suggest an emerging IVL/VDR paradigm for adaptation in the human epidermis.

Selective sweep for an enhancer involucrin allele identifies skin barrier adaptation out of Africa.

The genetic modules that contribute to human evolution are poorly understood. Here we investigate positive selection in the Epidermal Differentiation Complex locus for skin barrier adaptation in diverse HapMap human populations (CEU, JPT/CHB, and YRI). Using Composite of Multiple Signals and iSAFE, we identify selective sweeps for LCE1A-SMCP and involucrin (IVL) haplotypes associated with human migration out-of-Africa, reaching near fixation in European populations. CEU-IVL is associated with increased IVL expression and a known epidermis-specific enhancer. CRISPR/Cas9 deletion of the orthologous mouse enhancer in vivo reveals a functional requirement for the enhancer to regulate Ivl expression in cis. Reporter assays confirm increased regulatory and additive enhancer effects of CEU-specific polymorphisms identified at predicted IRF1 and NFIC binding sites in the IVL enhancer (rs4845327) and its promoter (rs1854779). Together, our results identify a selective sweep for a cis regulatory module for CEU-IVL, highlighting human skin barrier evolution for increased IVL expression out-of-Africa.

Current understanding of epigenetics in atopic dermatitis.

Atopic dermatitis (AD) is an inflammatory skin disorder affecting up to 20% of the paediatric population worldwide. AD patients commonly exhibit dry skin and pruritus and are at a higher risk for developing asthma as well as allergic rhinitis. Filaggrin loss-of-function variants are the most widely replicated genetic risk factor among >40 genes associated with AD susceptibility. The prevalence of AD has tripled in the past 30 years in industrial countries around the world. This urgent public health issue has prompted the field to more thoroughly investigate the mechanisms that underlie AD pathogenesis amidst environmental exposures. Epigenetics is the study of heritable, yet reversible, modifications to the genome that affect gene expression. The past decade has seen an emergence of exciting studies identifying a role for epigenetic regulation associated with AD and at the interface of environmental factors. Such epigenetic studies have been empowered by sequencing technologies and human genome variation and epigenome maps. miRNAs that post-transcriptionally modify gene expression and circRNAs have also been discovered to be associated with AD. Here, we review our current understanding of epigenetics associated with atopic dermatitis. We discuss studies identifying distinct DNA methylation changes in keratinocytes and T cells, eQTLs as DNA methylation switches that impact gene expression, and histone modification changes associated with AD-related microbial dysbiosis. We further highlight the need for integrative and collaborative analyses to elucidate the impact of these epigenetic findings as potential drivers for AD pathogenesis and the translation of this new knowledge to develop newer targeted treatments.

High-throughput Identification of Gene Regulatory Sequences Using Next-generation Sequencing of Circular Chromosome Conformation Capture (4C-seq).

The identification of regulatory elements for a given target gene poses a significant technical challenge owing to the variability in the positioning and effect sizes of regulatory elements to a target gene. Some progress has been made with the bioinformatic prediction of the existence and function of proximal epigenetic modifications associated with activated gene expression using conserved transcription factor binding sites. Chromatin conformation capture studies have revolutionized our ability to discover physical chromatin contacts between sequences and even within an entire genome. Circular chromatin conformation capture coupled with next-generation sequencing (4C-seq), in particular, is designed to discover all possible physical chromatin interactions for a given sequence of interest (viewpoint), such as a target gene or a regulatory enhancer. Current 4C-seq strategies directly sequence from within the viewpoint but require numerous and diverse viewpoints to be simultaneously sequenced to avoid the technical challenges of uniform base calling (imaging) with next generation sequencing platforms. This volume of experiments may not be practical for many laboratories. Here, we report a modified approach to the 4C-seq protocol that incorporates both an additional restriction enzyme digest and qPCR-based amplification steps that are designed to facilitate a greater capture of diverse sequence reads and mitigate the potential for PCR bias, respectively. Our modified 4C method is amenable to the standard molecular biology lab for assessing chromatin architecture.

Tiled array-based sequencing identifies enrichment of loss-of-function variants in the highly homologous filaggrin gene in African-American children with severe atopic dermatitis.

Filaggrin (FLG) loss-of-function (LOF) variants are a major risk factor for the common inflammatory skin disease, atopic dermatitis (AD) and are often population-specific. African-American (AA) children are disproportionately affected with AD, often later developing asthma and/or allergic rhinitis and comprise an atopy health disparity group for which the role of FLG LOF is not well known. Discovery of FLG LOF using exome sequencing is challenging given the known difficulties for accurate short-read alignment to FLG's high homology repeat variation. Here, we employed an array-based sequencing approach to tile across each FLG repeat and discover FLG LOF in a well-characterized cohort of AA children with moderate-to-severe AD. Five FLG LOF were identified in 23% of our cohort. Two novel FLG LOF singletons, c.488delG and p.S3101*, were discovered as well as p.R501*, p.R826* and p.S3316* previously reported for AD. p.S3316* (rs149484917) is likely an African ancestral FLG LOF, reported in African individuals in ExAC (Exome Aggregation Consortium), Exome Variant Server (ESP), and 4 African 1000G population databases (ESN, MSL, ASW, and ACB). The proportion of FLG LOF (11.5%) among the total FLG alleles in our cohort was significantly higher in comparisons with FLG LOF reported for African individuals in ExAC (2.5%; P = 4.3 × 10-4 ) and ESP (1.7%; P = 3.5 × 10-5 ) suggesting a disease-enrichment effect for FLG LOF. Our results demonstrate the utility of array-based sequencing in discovering FLG LOF, including novel and population-specific, which are of higher prevalence in our AA severe AD group than previously reported.

Recent evolution of the human skin barrier.

The skin is the first line of defense against the environment, with the epidermis as the outermost tissue providing much of the barrier function. Given its direct exposure to and encounters with the environment, the epidermis must evolve to provide an optimal barrier for the survival of an organism. Recent advances in genomics have identified a number of genes for the human skin barrier that have undergone evolutionary changes since humans diverged from chimpanzees. Here, we highlight a selection of key and innovative genetic findings for skin barrier evolution in our divergence from our primate ancestors and among modern human populations.

Electrophilic properties of itaconate and derivatives regulate the IκBζ-ATF3 inflammatory axis.

Metabolic regulation has been recognized as a powerful principle guiding immune responses. Inflammatory macrophages undergo extensive metabolic rewiring 1 marked by the production of substantial amounts of itaconate, which has recently been described as an immunoregulatory metabolite 2 . Itaconate and its membrane-permeable derivative dimethyl itaconate (DI) selectively inhibit a subset of cytokines 2 , including IL-6 and IL-12 but not TNF. The major effects of itaconate on cellular metabolism during macrophage activation have been attributed to the inhibition of succinate dehydrogenase2,3, yet this inhibition alone is not sufficient to account for the pronounced immunoregulatory effects observed in the case of DI. Furthermore, the regulatory pathway responsible for such selective effects of itaconate and DI on the inflammatory program has not been defined. Here we show that itaconate and DI induce electrophilic stress, react with glutathione and subsequently induce both Nrf2 (also known as NFE2L2)-dependent and -independent responses. We find that electrophilic stress can selectively regulate secondary, but not primary, transcriptional responses to toll-like receptor stimulation via inhibition of IκBζ protein induction. The regulation of IκBζ is independent of Nrf2, and we identify ATF3 as its key mediator. The inhibitory effect is conserved across species and cell types, and the in vivo administration of DI can ameliorate IL-17-IκBζ-driven skin pathology in a mouse model of psoriasis, highlighting the therapeutic potential of this regulatory pathway. Our results demonstrate that targeting the DI-IκBζ regulatory axis could be an important new strategy for the treatment of IL-17-IκBζ-mediated autoimmune diseases.

Novel oxytocin receptor variants in laboring women requiring high doses of oxytocin.

BACKGROUND: Although oxytocin commonly is used to augment or induce labor, it is difficult to predict its effectiveness because oxytocin dose requirements vary significantly among women. One possibility is that women requiring high or low doses of oxytocin have variations in the oxytocin receptor gene. OBJECTIVES: To identify oxytocin receptor gene variants in laboring women with low and high oxytocin dosage requirements. STUDY DESIGN: Term, nulliparous women requiring oxytocin doses of ≤4 mU/min (low-dose-requiring, n = 83) or ≥20 mU/min (high-dose-requiring, n = 104) for labor augmentation or induction provided consent to a postpartum blood draw as a source of genomic DNA. Targeted-amplicon sequencing (coverage >30×) with MiSeq (Illumina) was performed to discover variants in the coding exons of the oxytocin receptor gene. Baseline relevant clinical history, outcomes, demographics, and oxytocin receptor gene sequence variants and their allele frequencies were compared between low-dose-requiring and high-dose-requiring women. The Scale-Invariant Feature Transform algorithm was used to predict the effect of variants on oxytocin receptor function. The Fisher exact or χ2 tests were used for categorical variables, and Student t tests or Wilcoxon rank sum tests were used for continuous variables. A P value < .05 was considered statistically significant. RESULTS: The high-dose-requiring women had greater rates of obesity and diabetes and were more likely to have undergone labor induction and required prostaglandins. High-dose-requiring women were more likely to undergo cesarean delivery for first-stage arrest and less likely to undergo cesarean delivery for nonreassuring fetal status. Targeted sequencing of the oxytocin receptor gene in the total cohort (n = 187) revealed 30 distinct coding variants: 17 nonsynonymous, 11 synonymous, and 2 small structural variants. One novel variant (A243T) was found in both the low- and high-dose-requiring groups. Three novel variants (Y106H, A240_A249del, and P197delfs*206) resulting in an amino acid substitution, loss of 9 amino acids, and a frameshift stop mutation, respectively, were identified only in low-dose-requiring women. Nine nonsynonymous variants were unique to the high-dose-requiring group. These included 3 known variants (R151C, G221S, and W228C) and 6 novel variants (M133V, R150L, H173R, A248V, G253R, and I266V). Of these, R150L, R151C, and H173R were predicted by Scale-Invariant Feature Transform algorithm to damage oxytocin receptor function. There was no statistically significant association between the numbers of synonymous and nonsynonymous substitutions in the patient groups. CONCLUSION: Obesity, diabetes, and labor induction were associated with the requirement for high doses of oxytocin. We did not identify significant differences in the prevalence of oxytocin receptor variants between low-dose-requiring and high-dose-requiring women, but novel oxytocin receptor variants were enriched in the high-dose-requiring women. We also found 3 oxytocin receptor variants (2 novel, 1 known) that were predicted to damage oxytocin receptor function and would likely increase an individual's risk for requiring a high oxytocin dose. Further investigation of oxytocin receptor variants and their effects on protein function will inform precision medicine in pregnant women.

The Molecular Revolution in Cutaneous Biology: EDC and Locus Control.

The epidermal differentiation complex (EDC) locus consists of a cluster of genes important for the terminal differentiation of the epidermis. While early studies identified the functional importance of individual EDC genes, the recognition of the EDC genes as a cluster with its shared biology, homology, and physical linkage was pivotal to later studies that investigated the transcriptional regulation of the locus. Evolutionary conservation of the EDC and the transcriptional activation during epidermal differentiation suggested a cis-regulatory mechanism via conserved noncoding elements or enhancers. This line of pursuit led to the identification of CNE 923, an epidermal-specific enhancer that was found to mediate chromatin remodeling of the EDC in an AP-1 dependent manner. These genomic studies, as well as the advent of high-throughput sequencing and genome engineering techniques, have paved the way for future investigation into enhancer-mediated regulatory networks in cutaneous biology.

Recent Positive Selection in Genes of the Mammalian Epidermal Differentiation Complex Locus.

The epidermal differentiation complex (EDC) is the most rapidly evolving locus in the human genome compared to that of the chimpanzee. Yet the EDC genes that are undergoing positive selection across mammals and in humans are not known. We sought to identify the positively selected genetic variants and determine the evolutionary events of the EDC using mammalian-wide and clade-specific branch- and branch-site likelihood ratio tests and a genetic algorithm (GA) branch test. Significant non-synonymous substitutions were found in filaggrin, SPRR4, LELP1, and S100A2 genes across 14 mammals. By contrast, we identified recent positive selection in SPRR4 in primates. Additionally, the GA branch test discovered lineage-specific evolution for distinct EDC genes occurring in each of the nodes in the 14-mammal phylogenetic tree. Multiple instances of positive selection for FLG, TCHHL1, SPRR4, LELP1, and S100A2 were noted among the primate branch nodes. Branch-site likelihood ratio tests further revealed positive selection in specific sites in SPRR4, LELP1, filaggrin, and repetin across 14 mammals. However, in addition to continuous evolution of SPRR4, site-specific positive selection was also found in S100A11, KPRP, SPRR1A, S100A7L2, and S100A3 in primates and filaggrin, filaggrin2, and S100A8 in great apes. Very recent human positive selection was identified in the filaggrin2 L41 site that was present in Neanderthal. Together, our results identifying recent positive selection in distinct EDC genes reveal an underappreciated evolution of epidermal skin barrier function in primates and humans.

The Asian atopic dermatitis phenotype combines features of atopic dermatitis and psoriasis with increased TH17 polarization.

BACKGROUND: Atopic dermatitis (AD) shows very high prevalence in Asia, with a large unmet need for effective therapeutics. Direct comparisons between European American (EA) and Asian patients with AD are unavailable, but earlier blood studies detected increased IL-17(+)-producing cell counts in Asian patients with AD. OBJECTIVE: We sought to characterize the Asian AD skin phenotype and compare it with the EA AD skin phenotype. METHODS: We performed genomic profiling (real-time PCR) and immunohistochemistry on lesional and nonlesional biopsy specimens from 52 patients with AD (25 EAs and 27 Asians), 10 patients with psoriasis (all EAs), and 27 healthy subjects (12 EAs and 15 Asians). RESULTS: Although disease severity/SCORAD scores were similar between the AD groups (58.0 vs 56.7, P = .77), greater acanthosis, higher Ki67 counts, and frequent parakeratosis were characteristics of lesional epidermis from Asian patients with AD (P < .05). Most (24/27) Asian patients had high IgE levels. A principal component analysis using real-time PCR data clustered the Asian AD phenotype between the EA AD and psoriasis phenotypes. TH2 skewing characterized both Asian and EA patients with AD but not patients with psoriasis. Significantly higher TH17 and TH22 (IL17A, IL19, and S100A12 in lesional and IL-22 in nonlesional skin; P < .05) and lower TH1/interferon (CXCL9, CXCL10, MX1, and IFNG in nonlesional skin; P < .05) gene induction typified AD skin in Asian patients. CONCLUSION: The Asian AD phenotype presents (even in the presence of increased IgE levels) a blended phenotype between that of EA patients with AD and those with psoriasis, including increased hyperplasia, parakeratosis, higher TH17 activation, and a strong TH2 component. The relative pathogenic contributions of the TH17 and TH2 axes in creating the Asian AD phenotype need to be tested in future clinical trials with appropriate targeted therapeutics.

Epigenetic coordination of embryonic heart transcription by dynamically regulated long noncoding RNAs.

The vast majority of mammalian DNA does not encode for proteins but instead is transcribed into noncoding (nc)RNAs having diverse regulatory functions. The poorly characterized subclass of long ncRNAs (lncRNAs) can epigenetically regulate protein-coding genes by interacting locally in cis or distally in trans. A few reports have implicated specific lncRNAs in cardiac development or failure, but precise details of lncRNAs expressed in hearts and how their expression may be altered during embryonic heart development or by adult heart disease is unknown. Using comprehensive quantitative RNA sequencing data from mouse hearts, livers, and skin cells, we identified 321 lncRNAs present in the heart, 117 of which exhibit a cardiac-enriched pattern of expression. By comparing lncRNA profiles of normal embryonic (∼E14), normal adult, and hypertrophied adult hearts, we defined a distinct fetal lncRNA abundance signature that includes 157 lncRNAs differentially expressed compared with adults (fold-change ≥ 50%, false discovery rate = 0.02) and that was only poorly recapitulated in hypertrophied hearts (17 differentially expressed lncRNAs; 13 of these observed in embryonic hearts). Analysis of protein-coding mRNAs from the same samples identified 22 concordantly and 11 reciprocally regulated mRNAs within 10 kb of dynamically expressed lncRNAs, and reciprocal relationships of lncRNA and mRNA levels were validated for the Mccc1 and Relb genes using in vitro lncRNA knockdown in C2C12 cells. Network analysis suggested a central role for lncRNAs in modulating NFκB- and CREB1-regulated genes during embryonic heart growth and identified multiple mRNAs within these pathways that are also regulated, but independently of lncRNAs.

Cyclosporine in patients with atopic dermatitis modulates activated inflammatory pathways and reverses epidermal pathology.

BACKGROUND: Atopic dermatitis (AD) is the most common inflammatory disease. Evolving disease models link changes in epidermal growth and differentiation to T(H)2/T(H)22 cytokine activation. However, these models have not been tested by in vivo suppression of T-cell cytokines. Cyclosporine (CsA) is an immunosuppressant that is highly effective for severe disease, but its mechanism in AD skin lesions has not been studied. OBJECTIVE: We sought to establish the ability of a systemic immunosuppressant to modulate immune and epidermal alterations that form the pathogenic disease phenotype and to correlate changes with clinical improvement. METHODS: CsA's effects on AD skin pathology were evaluated by using gene expression and immunohistochemistry studies in baseline, week 2, and week 12 lesional and nonlesional biopsy specimens from 19 patients treated with 5 mg/kg/d CsA for 12 weeks. RESULTS: After 2 and 12 weeks of treatment, we observed significant reductions of 51% and 72%, respectively, in SCORAD scores. Clinical improvements were associated with significant gene expression changes in lesional but also nonlesional skin, particularly reductions in levels of T(H)2-, T(H)22-, and some T(H)17-related molecules (ie, IL-13, IL-22, CCL17, S100As, and elafin/peptidase inhibitor 3), and modulation of epidermal hyperplasia and differentiation measures. CONCLUSIONS: This is the first study that establishes a relationship between cytokine activation and molecular epidermal alterations, as well as correlations between disease biomarkers in the skin and clinical improvement. The reversal of the molecular phenotype with CsA and the associated biomarkers can serve as a reference for the successful modulation of tissue inflammation with specific immune antagonists in future studies, contributing to the understanding of the specific cytokines involved in epidermal pathology.

Regulation of the dynamic chromatin architecture of the epidermal differentiation complex is mediated by a c-Jun/AP-1-modulated enhancer.

The epidermal differentiation complex (EDC) locus comprises a syntenic and linear cluster of genes whose concomitant expression is a hallmark feature of differentiation in the developing skin epidermis. Many of the EDC proteins are cross-linked together to form the cornified envelope, an essential and discrete unit of the mammalian skin barrier. The mechanism underlying coordinate transcriptional activation of the EDC is unknown. Within the human EDC, we identified an epidermal-specific regulatory enhancer, 923, which responded to the developmental and spatiotemporal cues at the onset of epidermal differentiation in the mouse embryo. Comparative chromosomal conformation capture assays in proliferating and differentiated primary mouse keratinocytes revealed multiple physiologically sensitive chromatin interactions between the 923 enhancer and EDC gene promoters, thus depicting the dynamic chromatin topology of the EDC. We elucidate a mechanistic link between c-Jun/AP-1 and 923, whereby AP-1- and 923-mediated EDC chromatin remodeling are required for functional EDC gene activation. Thus, we identify a critical enhancer/transcription factor axis governing the dynamic regulation of the EDC chromatin architecture and gene expression and provide a framework for future studies toward understanding gene regulation in cutaneous diseases.

Intrinsic atopic dermatitis shows similar TH2 and higher TH17 immune activation compared with extrinsic atopic dermatitis.

BACKGROUND: Atopic dermatitis (AD) is classified as extrinsic and intrinsic, representing approximately 80% and 20% of patients with the disease, respectively. Although sharing a similar clinical phenotype, only extrinsic AD is characterized by high serum IgE levels. Because most patients with AD exhibit high IgE levels, an "allergic"/IgE-mediated disease pathogenesis was hypothesized. However, current models associate AD with T-cell activation, particularly TH2/TH22 polarization, and epidermal barrier defects. OBJECTIVE: We sought to define whether both variants share a common pathogenesis. METHODS: We stratified 51 patients with severe AD into extrinsic AD (n = 42) and intrinsic AD (n = 9) groups (with similar mean disease activity/SCORAD scores) and analyzed the molecular and cellular skin pathology of lesional and nonlesional intrinsic AD and extrinsic AD by using gene expression (real-time PCR) and immunohistochemistry. RESULTS: A significant correlation between IgE levels and SCORAD scores (r = 0.76, P < 10(-5)) was found only in patients with extrinsic AD. Marked infiltrates of T cells and dendritic cells and corresponding epidermal alterations (keratin 16, Mki67, and S100A7/A8/A9) defined lesional skin of patients with both variants. However, higher activation of all inflammatory axes (including TH2) was detected in patients with intrinsic AD, particularly TH17 and TH22 cytokines. Positive correlations between TH17-related molecules and SCORAD scores were only found in patients with intrinsic AD, whereas only patients with extrinsic AD showed positive correlations between SCORAD scores and TH2 cytokine (IL-4 and IL-5) levels and negative correlations with differentiation products (loricrin and periplakin). CONCLUSIONS: Although differences in TH17 and TH22 activation exist between patients with intrinsic AD and those with extrinsic AD, we identified common disease-defining features of T-cell activation, production of polarized cytokines, and keratinocyte responses to immune products. Our data indicate that a TH2 bias is not the sole cause of high IgE levels in patients with extrinsic AD, with important implications for similar therapeutic interventions.