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.

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.

Translational profiling of hypocretin neurons identifies candidate molecules for sleep regulation.

Hypocretin (orexin; Hcrt)-containing neurons of the hypothalamus are essential for the normal regulation of sleep and wake behaviors and have been implicated in feeding, anxiety, depression, and reward. The absence of these neurons causes narcolepsy in humans and model organisms. However, little is known about the molecular phenotype of these cells; previous attempts at comprehensive profiling had only limited sensitivity or were inaccurate. We generated a Hcrt translating ribosome affinity purification (bacTRAP) line for comprehensive translational profiling of all ribosome-bound transcripts in these neurons in vivo. From this profile, we identified >6000 transcripts detectably expressed above background and 188 transcripts that are highly enriched in these neurons, including all known markers of the cells. Blinded analysis of in situ hybridization databases suggests that ~60% of these are expressed in a Hcrt marker-like pattern. Fifteen of these were confirmed with double labeling and microscopy, including the transcription factor Lhx9. Ablation of this gene results in a >30% loss specifically of Hcrt neurons, without a general disruption of hypothalamic development. Polysomnography and activity monitoring revealed a profound hypersomnolence in these mice. These data provide an in-depth and accurate profile of Hcrt neuron gene expression and suggest that Lhx9 may be important for specification or survival of a subset of these cells.

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.

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.

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.

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.

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.

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.

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.

Progressive activation of T(H)2/T(H)22 cytokines and selective epidermal proteins characterizes acute and chronic atopic dermatitis.

BACKGROUND: Atopic dermatitis (AD) is a common disease with an increasing prevalence. The primary pathogenesis of the disease is still elusive, resulting in the lack of specific treatments. AD is currently considered a biphasic disease, with T(H)2 predominating in acute disease and a switch to T(H)1 characterizing chronic disease. Elucidation of the molecular factors that participate in the onset of new lesions and maintenance of chronic disease is critical for the development of targeted therapeutics. OBJECTIVES: We sought to characterize the mechanisms underlying the onset and maintenance of AD. METHODS: We investigated intrapersonal sets of transcriptomes from nonlesional skin and acute and chronic lesions of 10 patients with AD through genomic, molecular, and cellular profiling. RESULTS: Our study associated the onset of acute lesions with a striking increase in a subset of terminal differentiation proteins, specifically the cytokine-modulated S100A7, S100A8, and S100A9. Acute disease was also associated with significant increases in gene expression levels of major T(H)22 and T(H)2 cytokines and smaller increases in IL-17 levels. A lesser induction of T(H)1-associated genes was detected in acute disease, although some were significantly upregulated in chronic disease. Further significant intensification of major T(H)22 and T(H)2 cytokines was observed between acute and chronic lesions. CONCLUSIONS: Our data identified increased S100A7, S100A8, and S100A9 gene expression with AD initiation and concomitant activation of T(H)2 and T(H)22 cytokines. Our findings support a model of progressive activation of T(H)2 and T(H)22 immune axes from the acute to chronic phases, expanding the prevailing view of pathogenesis with important therapeutic implications.

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.

A milieu of regulatory elements in the epidermal differentiation complex syntenic block: implications for atopic dermatitis and psoriasis.

Two common inflammatory skin disorders with impaired barrier, atopic dermatitis (AD) and psoriasis, share distinct genetic linkage to the Epidermal Differentiation Complex (EDC) locus on 1q21. The EDC is comprised of tandemly arrayed gene families encoding proteins involved in skin cell differentiation. Discovery of semi-dominant mutations in filaggrin (FLG) associated with AD and a copy number variation within the LCE genes associated with psoriasis provide compelling evidence for the role of EDC genes in the pathogenesis of these diseases. To date, little is known about the potentially complex regulatory landscape within the EDC. Here, we report a computational approach to identify conserved non-coding elements (CNEs) in the EDC queried for regulatory function. Coordinate expression of EDC genes during mouse embryonic skin development and a striking degree of synteny and linearity in the EDC locus across a wide range of mammalian (placental and marsupial) genomes suggests an evolutionary conserved regulatory milieu in the EDC. CNEs identified by comparative genomics exhibit dynamic regulatory activity (enhancer or repressor) in differentiating or proliferating conditions. We further demonstrate epidermal-specific, developmental in vivo enhancer activities (DNaseI and transgenic mouse assays) in CNEs, including one within the psoriasis-associated deletion, LCE3C_LCE3B-del. Together, our multidisciplinary study features a network of regulatory elements coordinating developmental EDC gene expression as an unexplored resource for genetic variants in skin diseases.

Lipid defect underlies selective skin barrier impairment of an epidermal-specific deletion of Gata-3.

Skin lies at the interface between the complex physiology of the body and the external environment. This essential epidermal barrier, composed of cornified proteins encased in lipids, prevents both water loss and entry of infectious or toxic substances. We uncover that the transcription factor GATA-3 is required to establish the epidermal barrier and survive in the ex utero environment. Analysis of Gata-3 mutant transcriptional profiles at three critical developmental stages identifies a specific defect in lipid biosynthesis and a delay in differentiation. Genomic analysis identifies highly conserved GATA-3 binding sites bound in vivo by GATA-3 in the first intron of the lipid acyltransferase gene AGPAT5. Skin from both Gata-3-/- and previously characterized barrier-deficient Kruppel-like factor 4-/- newborns up-regulate antimicrobial peptides, effectors of innate immunity. Comparison of these animal models illustrates how impairment of the skin barrier by two genetically distinct mechanisms leads to innate immune responses, as observed in the common human skin disorders psoriasis and atopic dermatitis.

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.

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.

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.