Next Article in Journal
Effects of a Nutrition Education Programme in Stage IV Cardiac Rehabilitation Patients: A 3-Arm Randomised Controlled Trial
Next Article in Special Issue
Mechanism, Diagnostic and Therapeutic Novelties in Dermatology: Where Are We Going?
Previous Article in Journal
In Vitro Effects of Cypermethrin and Glyphosate on LPS-Induced Immune Cell Activation
Previous Article in Special Issue
Laser Therapy for Vulvar Lichen Sclerosus, a Systematic Review
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Life
Volume 14
Issue 1
10.3390/life14010060
life-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Genetic Association and Differential RNA Expression of Histone (De)Acetylation-Related Genes in Pemphigus Foliaceus—A Possible Epigenetic Effect in the Autoimmune Response

by
Maiara Sulzbach Denardin
1,
Valéria Bumiller-Bini Hoch
1,2,
Amanda Salviano-Silva
1,2,3,
Sara Cristina Lobo-Alves
1,4,
Gabriel Adelman Cipolla
1,
Danielle Malheiros
1,2,
Danillo G. Augusto
1,5,
Michael Wittig
6,
Andre Franke
6,
Claudia Pföhler
7,
Margitta Worm
8,
Nina van Beek
9,
Matthias Goebeler
10,
Miklós Sárdy
11,12,
Saleh Ibrahim
13,14,
Hauke Busch
14,
Enno Schmidt
9,14,
Jennifer Elisabeth Hundt
14,
Maria Luiza Petzl-Erler
1,2 and
Angelica Beate Winter Boldt
1,2,*
1
Laboratory of Human Molecular Genetics, Department of Genetics, Federal University of Paraná (UFPR), Curitiba 81531-980, Brazil
2
Postgraduate Program in Genetics, Department of Genetics, Federal University of Paraná (UFPR), Curitiba 81531-980, Brazil
3
Department of Neurosurgery, University Medical Center Hamburg-Eppendorf, 20251 Hamburg, Germany
4
Research Institut Pelé Pequeno Príncipe, Curitiba 80250-060, Brazil
5
Department of Biological Sciences, The University of North Carolina at Charlotte, Charlotte, NC 28223, USA
6
Institute of Clinical Molecular Biology (IKMB), Christian-Albrechts-University of Kiel, 24105 Kiel, Germany
7
Department of Dermatology, Saarland University Medical Center, 66421 Homburg, Germany
8
Division of Allergy and Immunology, Department of Dermatology, Venerology and Allergy, Charité-Universitätsmedizin Berlin, 10117 Berlin, Germany
9
Department of Dermatology, University of Lübeck, 23562 Lübeck, Germany
10
Department of Dermatology, Venereology and Allergology, University Hospital Würzburg, 97080 Würzburg, Germany
11
Department of Dermatology and Allergy, University Hospital, LMU Munich, 80539 Munich, Germany
12
Department of Dermatology, Venereology and Dermatooncology, Semmelweis University, 1085 Budapest, Hungary
13
College of Medicine and Health Sciences, Khalifa University, Abu Dhabi 127788, United Arab Emirates
14
Lübeck Institute of Experimental Dermatology (LIED), University of Lübeck, 23562 Lübeck, Germany
 
add Show full affiliation list
 
remove Hide full affiliation list
*
Author to whom correspondence should be addressed.
Life 2024, 14(1), 60; https://doi.org/10.3390/life14010060
Submission received: 28 November 2023 / Revised: 23 December 2023 / Accepted: 25 December 2023 / Published: 29 December 2023
(This article belongs to the Special Issue Mechanism, Diagnostic and Therapeutic Novelties in Dermatology)
Browse Figures
Versions Notes

Abstract

:
Pemphigus foliaceus (PF) is an autoimmune skin blistering disease characterized by antidesmoglein-1 IgG production, with an endemic form (EPF) in Brazil. Genetic and epigenetic factors have been associated with EPF, but its etiology is still not fully understood. To evaluate the genetic association of histone (de)acetylation-related genes with EPF susceptibility, we evaluated 785 polymorphisms from 144 genes, for 227 EPF patients and 194 controls. Carriers of HDAC4_rs4852054*A were more susceptible (OR = 1.79, p = 0.0038), whereas those with GSE1_rs13339618*A (OR = 0.57, p = 0.0011) and homozygotes for PHF21A_rs4756055*A (OR = 0.39, p = 0.0006) were less susceptible to EPF. These variants were not associated with sporadic PF (SPF) in German samples of 75 SPF patients and 150 controls, possibly reflecting differences in SPF and EPF pathophysiology. We further evaluated the expression of histone (de)acetylation-related genes in CD4+ T lymphocytes, using RNAseq. In these cells, we found a higher expression of KAT2B, PHF20, and ZEB2 and lower expression of KAT14 and JAD1 in patients with active EPF without treatment compared to controls from endemic regions. The encoded proteins cause epigenetic modifications related to immune cell differentiation and cell death, possibly affecting the immune response in patients with PF.

1. Introduction

Pemphigus foliaceus (PF) is an autoimmune blistering skin disease endemic in Brazil, where it is also known as fogo selvagem (“wild fire”) [1]. PF pathophysiology is distinguished by IgG1 and IgG4 pathogenic autoantibodies produced mainly against desmoglein 1 (DSG1), a transmembrane cell adhesion protein found in keratinocyte cell junction (desmosomes) [1,2]. The skin blisters in PF are caused by keratinocyte separation, mostly in the epidermal granular layer, a process known as acantholysis [2]. Studies trying to unveil the breakdown of the immunological tolerance in PF etiology indicate associations with environmental, genetic, and epigenetic factors [3,4,5]. Continued exposure to mosquito bites is among the candidate environmental factors for EPF [3]. Haplotypes of major histocompatibility complex (MHC) class II genes stand out as one of the most important genetic components: DR1_DQw1 and DR4_DQw3 haplotypes have been associated with susceptibility to EPF, while DR7_DQw2 and DR3_DQw2 have been associated with resistance against EPF [6]. Hence, certain alleles of single nucleotide polymorphisms (SNPs) in immune system genes have been reported to alter PF susceptibility [4]. Long noncoding RNAs, several of which are involved in epigenetic regulation of gene expression, were also associated with PF susceptibility [7,8,9]. Furthermore, allelic variants of genes involved in epigenetic regulation, whose products are histone lysine demethylase and methyltransferases (KDM4C, SETD7, MECOM, and PRDM16), have been associated with PF susceptibility [5].
Epigenetic modifications can alter chromatin structure, triggering different patterns of gene expression. Chromatin comprises DNA and histone (H1, H2A, H2B, H3, H4, as well as histone variants) and nonhistone proteins. Covalent post-translational modifications (PTMs) on histone tails, such as (de)methylation, (de)acetylation, and (de)phosphorylation, are established and removed by chromatin-modifying enzymes that integrate multiunit complexes [10]. Among those enzymes, the activity of well-characterized histone acetyltransferases (HATs) and histone deacetylases (HDACs) is correlated with gene activation and repression, respectively [11].
Dynamic epigenetic modifications control immune cells’ development, activation, and differentiation, which affect the development of autoimmune diseases [12]. Abnormal histone acetylation (altered H3/H4 acetylation and/or HDAC expression patterns) has been implicated in autoimmune diseases affecting the skin, such as systemic lupus erythematosus [13,14], alopecia areata [15], and pemphigus vulgaris [16]. In this work, we investigated whether allelic variants of genes encoding HATs, HDACs, and members of complexes that interact with these enzymes influence PF susceptibility and whether these genes are differentially expressed in CD4+ T cells. Genetic association and differential RNA expression of histone (de)acetylation-related genes in EPF may indicate a possible epigenetic effect in its autoimmune response, especially related to immune cell differentiation and cell death.

2. Materials and Methods

2.1. Association Analysis Samples

Two samples were enrolled in this study—one of EPF and another of sporadic PF (SPF) (Figure 1). DNA from 227 EPF patients and 194 controls was extracted from peripheral blood mononuclear cells (PBMCs) by the phenol–chloroform–isoamyl alcohol method. They lived in Brazilian Midwest, Southeast, and South regions and were recruited in the Hospital Adventista do Pênfigo (Campo Grande, Mato Grosso do Sul), Lar da Caridade—Hospital do Fogo Selvagem (Uberaba, Minas Gerais), Hospital das Clínicas—Universidade de São Paulo (Ribeirão Preto, São Paulo), Hospital de Clínicas—Universidade Federal do Paraná, Hospital de Dermatologia Sanitária São Roque, and Hospital Santa Casa de Misericórdia (Curitiba, Paraná). Patients and controls were unrelated, of predominantly European ancestry [17], and had a similar sex ratio and age distribution: in both groups, 52% were women and the mean age of patients was 40.9 (6–83) years and of controls, 44.8 (11–86) years. No patient had a history of other autoimmune diseases. None of the controls presented autoimmune diseases.
DNA from 75 SPF patients and 150 controls was isolated from whole blood using the QIAamp DNA Maxi Blood Kit (Qiagen, Hilden, Germany). Patients and controls were all of European ancestry, predominantly of German origin, and recruited in German hospitals by the German Autoimmune Bullous Diseases Genetic Study Group. They had a similar sex ratio and age distribution: 46% of patients and 51% of controls were women and the mean age of patients was 60.1 (25–88) years, and of controls, 60.3 (21–77) years.

2.2. RNA Expression Analysis Samples

Total RNA was extracted from CD4+ T lymphocytes isolated from peripheral blood samples from four patients (three females and one male), before the start of oral corticosteroid treatment, who constantly presented new skin lesions and whose mean age was 36 (15–54) years, as well as from five controls (four females and one male) from EPF-endemic regions, whose mean age was 39.6 (22–58) years. The selection process of CD4+ T cells has been previously described by Bumiller-Bini Hoch et al. [18].

2.3. SNP Selection and Genotyping

Genotype data were available from a previous study using microarray hybridization (CoreExome-24, version 1.1, Illumina, San Diego, CA, USA) [19]. The genomic positions of 144 genes, encoding HATs (12) and HDACs (18) and members of histone acetylation (59) and deacetylation (55) complexes, were identified, considering 1000 base pairs downstream and upstream from the sequence for the longest transcript, according to the human genome version GRC37/hg19 (NCBI gene) [20] (Table S1). A total of 2486 SNPs located in these gene sequences were extracted from DNA microarray data of EPF samples.
SNPs associated with EPF were genotyped in the SPF samples using the iPLEX platform of the MassARRAY system (Agena Bioscience, San Diego, CA, USA). The primer sequences are available in Table S2. MassARRAY Typer software (v4.0) (Agena Bioscience, San Diego, CA, USA) with standard settings was used to call the genotypes. The genotype distribution followed the Hardy–Weinberg equilibrium in patients, controls, and the set of participants (data not shown), except for rs13339618 (pcontrols = 0.049; pall participants = 0.005).

2.4. Association Analysis

Rare variants (minor allele frequency < 0.01) deviating from the Hardy–Weinberg equilibrium for controls (p < 0.05) or in high linkage disequilibrium (LD) with another SNP from the dataset (r2 ≥ 0.80) were excluded. After filtering, 785 SNPs remained for subsequent logistic regression association analysis with allele frequencies (additive model), frequencies of homozygotes for the minor allele (recessive model), and the sum of the frequencies of heterozygotes and homozygotes for the minor allele (dominant model) (PLINK version 1.1.9), correcting possible associations resulting from population stratification for two principal components of analysis and for sex. p < 0.005 was established as significant [21].

2.5. In Silico Analysis

Functional annotation available in reference public databases was used to explore the potential effects of the SNPs associated with PF. Gene and SNP annotations were performed with UCSC [22] and Ensembl [23] genome browsers. Gene expression was entrained with Protein Atlas [24], which, in addition to presenting its own database, also exhibits the Monaco dataset [25] and Schmiedel dataset [26]. Allelic variants’ epigenetic annotations were entrained with the HaploReg v.4.1 database [27] and SNPNexus tool [28], which filters many databases, including Roadmap and Encode [29,30]. LD between alleles was evaluated using LDLink [31], with European as a reference. The expression and splicing quantitative trait loci (eQTL and sQTL) effects were observed in Genehopper Qtlizer [32], which compiles information from several databases, such as Blood eQTL Browser [33], GRASP 2 Catalog [34], GTEx v8 [35], The Cardiogenics Project [36,37], and Zeller et al. [38]. Clinical significance was evaluated with ClinVar/NCBI [39].

2.6. RNAseq

Gene expression of CD4+ T lymphocytes was evaluated by RNA sequencing performed with the Illumina Hi-seq platform (San Diego, CA, USA) using a paired-end 150 bp protocol. The quality of RNAseq reads was investigated using FASTQC version 0.11.5. Reads were pseudoaligned to the human transcriptome (Ensembl version 93, GRCh38), using coding and non-coding genes. Expression was quantified using Salmon version 0.11.0 with 30 bootstrap cycles. Quantification files were imported to SLEUTH version 0.29.0 and pairwise comparisons between controls and EPF patients were performed. The Log2 fold change was obtained with the Wald test. All statistical analyses were performed using R version 3.4.4. p < 0.005 was established as significant [21].

3. Results

3.1. Genetic Association of Histone (De)Acetylation-Related Genes in EPF

Three intronic SNPs in different genes were associated with EPF. One of these genes encodes an HDAC (HDAC4) and the other two, HDAC complex members—GSE1 (Gse1 coiled-coil protein) and PHF21A (PHD zinc-finger protein 21A). Carriers of the HDAC4 rs4852054*A intronic allele (OR = 1.79 [95%CI = 1.21–2.67], p = 0.004) presented increased susceptibility to EPF. The GSE1 rs13339618*A regulatory allele (OR = 0.57 [95%CI = 0.41–0.80], p = 0.001) and the PHF21A rs4756055 A/A intronic genotype (OR = 0.391 [95%CI = 0.228–0.672], p = 0.0006) were associated with decreased EPF susceptibility. These associations were not replicated in SPF samples, suggesting different etiologies for EPF and SPF forms (Table 1). This may indicate distinctions in the pathogenesis of EPF and SPF.

3.2. Predicted Consequences of SNPs in Histone (De)Acetylation-Related Genes Associated with EPF

Concerning functional annotation, HDAC4_rs4852054*A and GSE1_rs13339618*A present regulatory features, being located on a promoter flanking region and an enhancer, respectively [23].
HDAC4 presents two similar transcripts (differing by five codons), expressed at similar levels in most tissues. In the blood, HDAC4 is mostly expressed in granulocytes (neutrophils—47.0 normalized transcript per million, nTPM, and basophils—32.0 nTPM on the HPA dataset, 26.0 nTPM in the Monaco dataset), and much less so in the myeloid and lymphocytic lineages (very low in HPA, between 5 and 12 nTPM in monocytes and 1 and 10 nTPM in different kinds of T cells of the Monaco dataset) [24,25]. Its expression in the skin is rather low, at 11.1 nTPM in 1734 suprabasal keratinocyte cells [24]. Even so, the region containing HDAC4_rs4852054*A is enriched for enhancer- (H3K4me1 and H3K27ac) and promoter-associated histone marks (H3K9ac) and/or open-chromatin (H3K36, DNAse) in immune (common myeloid progenitor CD34+ cells, monocytes, and T cells) and skin cells (fibroblasts, normal human keratinocytes, and melanocyte primary cells) (Table S3, with data from Haploreg v.4.1, Roadmap and ENCODE [27,29,30]). The region with this allelic variant is occupied by the NFIC regulatory protein in the myeloid K562 cell lineage [30]. Furthermore, HDAC4_rs4852054*A is in high LD (r2 > 0.90) with six other HDAC4 polymorphisms [31]. An HDAC4 haplotype (composed of six SNPs in almost absolute LD) is associated with higher HDAC4-AS1 (HDAC4 antisense RNA 1) gene expression in skin and/or blood cells [35] (Table 2).
GSE1 encodes five transcripts in the skin, whose expression is not significantly affected by sun exposure [35]. In the skin, GSE1 is expressed in suprabasal keratinocytes (29.1 nTPM in 1734 cells), B cells (28.8 nTPM in 268 skin B cells), and macrophages (31.2 nTPM in 961 skin macrophages, 24.5 nTPM in 1776 Langerhans cells). In lymph nodes, its expression in two B cell populations doubles (50.7–53.6 nTPM, 2283–2025 B cells) [24]. The region containing GSE1_rs13339618*A is enriched with enhancer- (H3K4me1 and H3K27ac) and/or promoter-associated histone marks (H3K4m3) in skin cells (fibroblasts and melanocyte primary cells) (Table S3). It is also occupied by the ZBTB33 and REST regulatory proteins in the K562 myeloid cell lineage [30]. GSE1_rs13339618*A presents high LD with seven other GSE1 alleles (Table 2) [31]. These alleles, except GSE1_rs8060638*A and GSE1_rs28671512*T, are associated with five times higher GSE1 gene expression in blood cells [33]. GSE1_rs13329722*C and GSE1_rs7498141*A are also associated with TRIM35 (tripartite motif containing 35) gene expression in monocytes [38] (Table 2).
PH21A presents a moderate expression of five transcript variants in the skin, regardless of sun exposure. Skin-resident Langerhans cells, T and B cells, as well as macrophages, all present similar expression of this gene (varying from 12.6 to 16.5 nTPM in Langerhans cells and T cells, respectively). Its expression is relatively higher in lymph node B cells (27.4–48.1 nTPM) and T cells (20.0–34.2 nTPM) [24]. The PHF21A_rs4756055*A region is enriched with enhancers (H3K4me1 and H3K27ac), promoters (H3K4me3 and H3K9ac), and/or open-chromatin (DNase) in immune (neutrophils, monocytes, and mononuclear, natural killer, T and B cells) and skin cells (fibroblasts, keratinocytes, and melanocytes) [27,29] (Table S3). This region is also occupied by myocyte enhancer factor 2B (MEF2B), paired box 8 (PAX8), and nuclear factor of activated T cells 3 (NFATC3) regulatory proteins in the GM12878 cell line [30]. PHF21A_rs4756055*A is in high LD with seven other PHF21 alleles (Table 2) [31]. These SNPs, excepting PHF21A_rs11374563*T, are associated with slightly higher HNRNPH1 (heterogeneous nuclear ribonucleoprotein H1) gene expression in macrophages [37]. Three SNPs in strong LD, including the one associated with lower susceptibility to EPF, are also associated with three times higher CRY2 (cryptochrome circadian regulator 2) expression in blood cells [33]. Three other SNPs displaying a somewhat lower LD (0.8 < r2 < 0.9) are associated with slightly lower CRY2 expression in monocytes [36]. PHF21A_rs923530*T, PHF21A_rs7107550*C, and PHF21A_rs74366855*A are associated with lower SLC37A1 (solute carrier family 37 member 1) gene expression in macrophages [37]. PHF21A_rs4756055*A is also associated with lower PEX16 (peroxisomal biogenesis factor 16) gene expression [33] and altered PHF21A gene expression [34] in blood cells (Table 2).
To date, no clinical significance has been reported for these variants [39].

3.3. Differential RNA Expression of Histone (De)Acetylation-Related Genes in EPF

Of the 144 genes encoding HATs and HDACs, and members of histone (de)acetylation complexes, 139 were analyzed by RNAseq. Five genes were found to be differently expressed in CD4+ T cells from patients with EPF without oral corticoid treatment and presenting active disease, compared to controls from PF endemic regions: ZEB2 (zinc finger E-box binding homeobox 2), KAT2B (lysine acetyltransferase 2B), and PHF20 (PHD finger protein 20) were overexpressed in patients, while KAT14 (lysine acetyltransferase 14) and JADE1 (jade family PHD finger 1) were underexpressed (Table 3). Interestingly, ZEB2 expression is elevated in T cells from bone marrow (229.6 nTPM), as well as from peripheral blood (248.6 nTPM). In turn, KAT2B is highly expressed in memory CD4+ T helper (Th) 1/Th17 cells, in contrast to PHF20 (303.5 vs. 47 TPM, respectively, in the Schmiedel dataset [26]). Both may reach 80–100 nTPM in T cells of vascular tissue. JADE1 also is highly expressed in memory CD4+ Th1/Th17 cells (100.2 nTPM in the Monaco dataset [25], although its expression in memory CD4+ T cells in the HPA dataset was lower—18.4 nTPM). KAT14, however, is less expressed in T cells (in these cells, its higher expression occurs in naïve CD4+ cells—37.3 nTPM) [24].
None of them coincided with genes presenting EPF-associated variants. This is not surprising, since the associated polymorphisms do not have an eQTL effect in lymphocytes. However, HDAC4 and KAT2B are known to co-occur (may physically interact). Furthermore, HDAC4, KAT2B, and ZEB2 are known to be coexpressed, as well as HNRNPH1, PHF21A, and GSE1 (Figure 2).

4. Discussion

Dysregulated homeostasis in autoimmune diseases is closely associated with epigenetic modifications [13]. Altered DNA methylation and histone PTMs are known to break immune tolerance [40], but their role is still poorly understood in PF. In this work, we identified HATs, HDACs, and interacting proteins either genetically associated with EPF or differentially expressed in CD4+ T cells (Figure 3). Three variants within three genes related to histone deacetylase were associated with altered EPF susceptibility. The lack of association in SPF may be related to different pathophysiological pathways (reviewed in [41]).
One of the herein associated genes encodes HDAC4, a member of the human class II histone deacetylases, known to silence the transcription of genes related to the immune response [42,43,44]. It represses myocyte enhancer factors 2 A (MEF2A), C (MEF2C), and D (MEF2D) [43,44,45]. Moreover, the caspase-cleaved amino-terminal fragment of HDAC4 (independently from the HDAC4 deacetylase domain) has been shown to repress MEF2C, triggering cell death [46]. It modulates thymocyte apoptosis and T and B cell differentiation and trafficking, by regulating the expression of the steroid receptor Nur77 [45,47] and KLF2 (Krüppel-like factor 2) [48,49]. These HDAC4 interactions indicate a possible role for cell death in EPF physiopathology, reinforced by our results and by others [50,51]. HDAC4 is also associated with B-cell lymphoma 6 (BCL6), a repressor protein essential for germinal center B cells and T follicular helper cells (Tfh), known to be increased in the blood of bullous pemphigoid patients [52,53,54]. By the way, DNA binding sites for MEF2A and BCL6 were significantly enriched in systemic lupus erythematosus susceptibility loci [55]. HDAC4 is also associated with regulatory factor X-associated ankyrin-containing protein and class II MHC transactivator (CIITA) protein, and represses CIITA-mediated MHC class II gene expression [56]. Moreover, CIITA_rs3087456*G carriers present an increased susceptibility to EPF [57]. This variant is also an eQTL for lower CIITA expression in sun-exposed skin [35]. Finally, HDAC4 forms a complex with GATA binding protein 3 and yin-yang 1 protein and downregulates IL5 (interleukin 5) expression [58], thus influencing Th2 cell differentiation and function. In the initial phase of the EPF immune response, Th1 and Th2 cytokines have been proposed to drive IgG1 and IgG4 anti-DSG1 production, while during the second phase, Th2 cytokines seem to drive IgG4 anti-DSG1 production [59,60].
HDAC4_rs4852054*A carriers were approximately 1.8 times more susceptible to EPF. This variant occurs within a five SNP haplotype consistently associated with higher expression of the HDAC4-AS1 lncRNA in blood and skin cells [35]. On hypoxia, this antisense inhibits HDAC4 expression in the human retinal pigment epithelial ARPE-19 cell line [61]. If this regulation also occurs in immune cells, higher HDAC4-AS1 expression in these individuals may ultimately decrease HDAC4 expression and activity, maintaining aberrant gene activation in the autoimmune response.
The other two polymorphisms associated with decreased EPF susceptibility are also associated with altered expression of histone deacetylase complex members, namely of two BRAF35-HDAC chromatin-remodeling complex components in blood: GSE1 and PHF21A [62,63]. GSE1_rs13339618*A is associated with greatly enhanced expression of GSE1 in blood cells, probably of the myeloid lineage. This higher gene expression may increase the efficiency of acetyl group removal from histones in inflammatory genes, reducing the odds of autoimmunity. GSE1_rs13339618*A is also associated with altered TRIM35 expression [38]. Interestingly, TRIM35 negatively regulates type I interferon response by interacting with IRF7, promoting its ubiquitination and subsequent degradation [64]. Cells and pathways with a pathogenic role in autoimmunity can be activated by type I interferon [65].
PHF21A regulates IL-1-induced gene expression [66] and is a core protein of the CoREST complex [67], whose inhibition impairs the suppressive function of regulatory T (Treg) cells [68]. These findings support the notion that dysfunctional Treg cells and cytokine production contribute to EPF pathophysiology. Furthermore, PHF21A_rs4756055*A and seven SNPs at LD are also associated with higher expression of the circadian CRY2 and HNRNPH1 genes in blood cells and macrophages, respectively, and with lower expression of PEX16 and SLC37A1 genes in the same types of cells. Thus, this region has a regulatory influence on several genes in immune cells, which may also be implicated in EPF. CRY2 expression is also reduced in celiac disease [69], and Cry double knockout mice have an overactivated B cell receptor signaling pathway, contributing to high serum IgG concentrations and antinuclear antibodies, with severe lung and kidney involvement [70]. Thus, higher CRY2 expression in PHF21A_rs4756055*A carriers may improve resistance to autoimmunity. Moreover, HNRNPH1 has been shown to reduce IL8 expression [71]. Interestingly, higher levels of IL-8 were found in the serum and blister fluids of patients with pemphigus [72]. The possible involvement of the other genes remains elusive.
Five genes were differentially expressed in CD4+ T cells from patients with active EPF without oral corticoid treatment. Three of these were overexpressed in patients—ZEB2, KAT2B, and PHF20, while the other two were underexpressed—KAT14 and JADE1.
ZEB2, also known as Smad-interacting protein 1, is a member of two-handed zinc finger/homeodomain transcription factors [73], which binds E box sequences and interacts with several co-repressors containing HDACs [74]. ZEB2 is essential for normal hematopoietic stem cell and hematopoietic progenitor cell differentiation and mobilization [75]. It interacts with T-Box Transcription Factor 21 and represses IL7R and IL2 expression, promoting terminal differentiation of lymphocytes, especially of CD8+ T cells [76,77]. It also regulates dendritic cell differentiation [78]. The consequences of its overexpression in CD4+ T cells from EPF patients warrant further investigation.
Consistent with our findings, the HDAC KAT2B gene is also upregulated in PBMCs from rheumatoid arthritis patients [79], and its overexpression increases histone H4 acetylation on the enhancer of MHC class I promoters [80]. It is nevertheless also important in Treg cell development [68,81,82]. KAT2B interacts with KLF10 and induces FOXP3 transcription [68]. Indeed, decreased levels of CD4+ FOXP3+ T cells occur in Kat2b-deficient mice [82,83]. Moreover, KAT2B acetylates the IL2 promoter, protecting Treg cells from undergoing apoptosis upon T cell receptor stimulation [81]. It even acetylates FOXP3 itself [81] and DEK protein [84]. When overexpressed, KAT2B acetylates DEK, reducing its affinity for DNA and impairing its transcription repression function [84]. Interestingly, anti-DEK autoantibodies have been identified in juvenile rheumatoid arthritis [85] and other inflammatory diseases [86]. Thus, KAT2B overexpression in T lymphocytes may maintain a deregulated immune response in EPF patients.
In contrast to ZEB2 and KAT2B, which are related to histone deacetylation, PHF20 is a component of the NSL (non-specific lethal) complex, which has HAT activity [87]. Phf20 knockout impairs thymocyte differentiation and maturation in mice [88]. Moreover, the second Tudor domain and the PHD finger of PHF20 recognize dimethylated lysines on p53, promoting nuclear factor kappa B (NF-κB) transcriptional activity [88,89,90], by stabilizing p53 [91] and avoiding its acetylation by the NSL complex [90]. PHF20 further regulates autophagy genes [92]. Its higher expression possibly implicates this protein in EPF autoimmune activation.
Among the two genes with lower expression in CD4+ T lymphocytes of patients, KAT14 is a subunit of the ADA2A-contaning acetyltransferase complex, essential for mammalian development, exerting a role in cell cycle progression and apoptosis [93]. Its impact on EPF autoimmune activation should be further investigated.
Jade family PHD finger 1 (JADE1, also known as PHF17) belongs to a subfamily of PHD zinc finger proteins, which contain two PHD zinc finger domains. JADE1 is subjected to post-transcriptional regulation and polyadenylation, resulting in multiple transcripts and two protein isoforms: the major isoform JADE1L and the truncated isoform JADE1S [94]. JADE1 is a transcriptional co-activator associated with HAT activity specific for histone H4 [95]. It interacts with lysine acetyltransferase 5 and 7 (KAT5 and KAT7, known as TIP60 and HBO1, respectively), which are close homologs, but is only able to induce transcription with KAT7 [95,96]. Essential roles of JADE1 in cell cycle regulation [97], and also as a consequence on apoptosis [98], have been described.

5. Conclusions

This is the first study suggesting that disruption of histone (and other proteins) (de)acetylation is involved in EPF pathophysiology, indicating a possible epigenetic effect in the autoimmune response.
We suggest that by turning chromatin into a close state, deacetylation-related proteins associated with increased EPF susceptibility (HDAC4) or overexpressed in patients’ CD4+ T cells (ZEB2) repress the expression of protective genes. Conversely, those associated with decreased EPF susceptibility (GSE1 and PHF21A) may repress the expression of genes conferring susceptibility to the disease. In the opposite direction, by turning chromatin into an open state, acetylation-related proteins overexpressed in patients’ CD4+ T cells (KAT2B and PHF20) may induce the expression of genes conferring susceptibility to the disease. Furthermore, those underexpressed (JADE1 and KAT14) induce the expression of protective genes. Some of these genes, whose expression is supposedly regulated by histone (de)acetylase-related proteins, seem to be essential genes in immune cell differentiation and cell death processes (Figure 3).
In conclusion, our findings highlight the importance of histone (de)acetylation dynamic processes, and its crosstalk with other histone PTMs, such as (de)methylation, in the autoimmune pathophysiology of EPF. Experiments evaluating the chromatin state may increase our understanding of epigenetic alterations in PF.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/life14010060/s1, Table S1: Genomic position of HATs, HDACs and histone (de)acetylation complexes members; Table S2: Primer sequences used for investigating genetic association in sporadic pemphigus foliaceus samples; Table S3: Epigenomic annotation of allelic variants associated with PF in immune and skin cells.

Author Contributions

Conceptualization, A.B.W.B.; methodology, M.S.D., V.B.-B.H., A.S.-S., S.C.L.-A., G.A.C., D.G.A., M.W. (Michael Wittig), A.F. and A.B.W.B.; resources, D.M., D.G.A., G.A.C., M.L.P.-E., C.P., M.W. (Margitta Worm), N.v.B., M.G., M.S., S.I. and E.S.; validation and visualization, M.S.D., V.B.-B.H. and A.B.W.B.; formal analysis and investigation, M.S.D., V.B.-B.H., H.B., S.C.L.-A. and A.B.W.B.; writing—original draft preparation, M.S.D., V.B.-B.H. and A.B.W.B.; writing—review and editing, M.S.D., V.B.-B.H., A.B.W.B., M.L.P.-E., A.S.-S., S.C.L.-A., D.G.A., G.A.C., D.M., C.P., M.W. (Margitta Worm), N.v.B., M.G. and H.B.; supervision and project administration, M.L.P.-E. and A.B.W.B.; funding acquisition, J.E.H., H.B., M.L.P.-E., E.S., S.I. and A.B.W.B. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by grants of the following funding agencies: Programa de Apoio a Núcleos de Excelência (PRONEX)—Fundação Araucária (FA) and Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq), Convênio 116/2018—Protocolo 50530, which provided support to M.L.P.-E.; Projeto 48.025—Chamada de Projetos 01/2016 (Fundação Araucária-PR/SESA-PR/CNPq/MS-Decit) for A.B.W.B.; the Swiss National Science Foundation CRSII5_202301/1; the University of Lübeck; the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) under Germany‘s Excellence Strategy—EXC 22167-390884018 provided financial support to H.B. and E.S. V.B.B.H. received a post-doc scholarship from Fundação Araucária/CNPq/UFPR (JDT2022271000049); G.A.C. received a postdoctoral scholarship (CAPES protocol 88882.306040/2018-01); M.L.P.-E. and A.B.W.B. received CNPq research scholarships (protocols 308783/2019-0 for M.L.P.-E. and 313741/2021 for A.B.W.B., respectively). The funding agencies had no role in study design, sample collection, data analysis and interpretation, or manuscript drafting and submission.

Institutional Review Board Statement

The study was conducted in accordance with the Declaration of Helsinki and approved by the Brazilian National Ethics Committee (CONEP) under the protocol number CAAE 02727412.4.0000.0096 (approval number 505988 from 2 January 2014, amended on 5 May 2017, number 2045739) and by the Ethics Committee of the University of Lübeck (approval number 08-156 from 19 November 2008, amended on 24 April 2009, and 12 December 2012).

Informed Consent Statement

Written informed consent has been obtained from all subjects involved in the study.

Data Availability Statement

The data presented in this study are available on request from the corresponding author. The data are not publicly available due to ethical restrictions.

Acknowledgments

We gratefully acknowledge the patients and all other volunteers for their participation in this study. We also thank the staff of the “Laboratório de Genética Molecular Humana/UFPR” for their assistance and helpful discussions, as well as the investigators of the German Autoimmune Bullous Diseases Study Group and funding sources.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Kasperkiewicz, M.; Ellebrecht, C.T.; Takahashi, H.; Yamagami, J.; Zillikens, D.; Payne, A.S.; Amagai, M. Pemphigus. Nat. Rev. Dis. Primers 2017, 3, 17026. [Google Scholar] [CrossRef] [PubMed]
  2. Schmidt, E.; Kasperkiewicz, M.; Joly, P. Pemphigus. Lancet 2019, 394, 882–894. [Google Scholar] [CrossRef] [PubMed]
  3. Qian, Y.; Jeong, J.S.; Maldonado, M.; Valenzuela, J.G.; Gomes, R.; Teixeira, C.; Evangelista, F.; Qaqish, B.; Aoki, V.; Hans, G.; et al. Cutting Edge: Brazilian Pemphigus Foliaceus Anti-Desmoglein 1 Autoantibodies Cross-React with Sand Fly Salivary LJM11 Antigen. J. Immunol. 2012, 189, 1535–1539. [Google Scholar] [CrossRef] [PubMed]
  4. Petzl-Erler, M.L. Beyond the HLA Polymorphism: A Complex Pattern of Genetic Susceptibility to Pemphigus. Genet. Mol. Biol. 2020, 43, e20190369. [Google Scholar] [CrossRef] [PubMed]
  5. Spadoni, M.B.; Bumiller-Bini, V.; Petzl-Erler, M.L.; Augusto, D.G.; Boldt, A.B.W. First Glimpse of Epigenetic Effects on Pemphigus Foliaceus. J. Investig. Dermatol. 2020, 140, 488–491.e1. [Google Scholar] [CrossRef] [PubMed]
  6. Petzl-Erler, M.L.; Santamaria, J. Are HLA Class II Genes Controlling Susceptibility and Resistance to Brazilian Pemphigus Foliaceus (Fogo Selvagem)? Tissue Antigens 1988, 33, 408–414. [Google Scholar] [CrossRef] [PubMed]
  7. Lobo-Alves, S.C.; Augusto, D.G.; Magalhães, W.C.S.; Tarazona-Santos, E.; Lima-Costa, M.F.; Barreto, M.L.; Horta, B.L.; de Almeida, R.C.; Petzl-Erler, M.L. Long Noncoding RNA Polymorphisms Influence Susceptibility to Endemic Pemphigus Foliaceus. Br. J. Dermatol. 2019, 181, 324–331. [Google Scholar] [CrossRef]
  8. Salviano-Silva, A.; Farias, T.D.J.; Bumiller-Bini, V.; Castro, M.D.S.; Lobo-Alves, S.C.; Busch, H.; Pföhler, C.; Worm, M.; Goebeler, M.; van Beek, N.; et al. Genetic Variability of Immune-Related LncRNAs: Polymorphisms in LINC-PINT and LY86-AS1 Are Associated with Pemphigus Foliaceus Susceptibility. Exp. Dermatol. 2021, 30, 831–840. [Google Scholar] [CrossRef]
  9. Salviano-Silva, A.; Becker, M.; Augusto, D.G.; Busch, H.; Adelman Cipolla, G.; Farias, T.D.J.; Bumiller-Bini, V.; Calonga-Solís, V.; Munz, M.; Franke, A.; et al. Genetic Association and Differential Expression of HLA Complex Group LncRNAs in Pemphigus. J. Autoimmun. 2021, 123, 102705. [Google Scholar] [CrossRef]
  10. Allis, C.D.; Caparros, M.-L.; Jenuwein, T.; Lachner, M.; Reinberg, D. Overview and Concepts. In Epigenetics; Allis, C.D., Caparros, M.-L., Jenuwein, T., Reinberg, D., Lachner, M., Eds.; Cold Spring Harbor Laboratory Press: New York, NY, USA, 2015. [Google Scholar]
  11. Marmorstein, R.; Zhou, M.M. Writers and Readers of Histone Acetylation: Structure, Mechanism, and Inhibition. Cold Spring Harb. Perspect. Biol. 2014, 6, a018762. [Google Scholar] [CrossRef]
  12. Zhao, M.; Wang, Z.; Yung, S.; Lu, Q. Epigenetic Dynamics in Immunity and Autoimmunity. Int. J. Biochem. Cell Biol. 2015, 67, 65–74. [Google Scholar] [CrossRef] [PubMed]
  13. Zhou, Y.; Su, Y.; Lu, Q.; Hu, N.; Qiu, X.; Luo, Y.; Yuan, J.; Li, Y.; Lei, W.; Zhang, G. Abnormal Histone Modification Patterns in Lupus CD4+ T Cells. J. Rheumatol. 2008, 35, 804–810. [Google Scholar]
  14. Zhang, Z.; Song, L.; Maurer, K.; Petri, M.A.; Sullivan, K.E. Global H4 Acetylation Analysis by ChIP-Chip in Systemic Lupus Erythematosus Monocytes. Genes. Immun. 2010, 11, 124–133. [Google Scholar] [CrossRef] [PubMed]
  15. Zhao, M.; Liang, G.; Wu, X.; Wang, S.; Zhang, P.; Su, Y.; Yin, H.; Tan, Y.; Zhang, J.; Lu, Q. Abnormal Epigenetic Modifications in Peripheral Blood Mononuclear Cells from Patients with Alopecia Areata. Br. J. Dermatol. 2012, 166, 266–273. [Google Scholar] [CrossRef] [PubMed]
  16. Zhao, M.; Huang, W.; Zhang, Q.; Gao, F.; Wang, L.; Zhang, G.; Su, Y.; Xiao, R.; Zhang, J.; Tang, M.; et al. Aberrant Epigenetic Modifications in Peripheral Blood Mononuclear Cells from Patients with Pemphigus Vulgaris. Br. J. Dermatol. 2012, 167, 523–531. [Google Scholar] [CrossRef] [PubMed]
  17. Calonga-Solís, V.; Amorim, L.M.; Farias, T.D.J.; Petzl-Erler, M.L.; Malheiros, D.; Augusto, D.G. Variation in Genes Implicated in B-Cell Development and Antibody Production Affects Susceptibility to Pemphigus. Immunology 2021, 162, 58–67. [Google Scholar] [CrossRef] [PubMed]
  18. Hoch, V.B.B.; Kohler, A.F.; Augusto, D.G.; Lobo-Alves, S.C.; Malheiros, D.; Cipolla, G.A.; Boldt, A.B.W.; Braun-Prado, K.; Wittig, M.; Franke, A.; et al. Genetic Associations and Differential MRNA Expression Levels of Host Genes Suggest a Viral Trigger for Endemic Pemphigus Foliaceus. Viruses 2022, 14, 879. [Google Scholar] [CrossRef] [PubMed]
  19. Augusto, D.G.; de Almeida, R.C.; Farias, T.D.J.; Magalhães, W.C.S.; Malheiros, D.; Lima-Costa, M.F.; Barreto, M.L.; Horta, B.L.; Kumar, V.; Wittig, M.; et al. Unsuspected Associations of Variants within the Genes NOTCH4 and STEAP2-AS1 Uncovered by a GWAS in Endemic Pemphigus Foliaceus. J. Investig. Dermatol. 2021, 141, 2741–2744. [Google Scholar] [CrossRef]
  20. Maglott, D.; Ostell, J.; Pruitt, K.D.; Tatusova, T. Entrez Gene: Gene-centered information at NCBI. Nucleic Acids Res. 2005, 33, D54–D58. [Google Scholar] [CrossRef]
  21. Benjamin, D.J.; Berger, J.O.; Johannesson, M.; Nosek, B.A.; Wagenmakers, E.J.; Berk, R.; Bollen, K.A.; Brembs, B.; Brown, L.; Camerer, C.; et al. Redefine Statistical Significance. Nat. Hum. Behav. 2018, 2, 6–10. [Google Scholar] [CrossRef]
  22. Kent, W.J.; Sugnet, C.W.; Furey, T.S.; Roskin, K.M.; Pringle, T.H.; Zahler, A.M.; Haussler, A.D. The Human Genome Browser at UCSC. Genome Res. 2002, 12, 996–1006. [Google Scholar] [CrossRef] [PubMed]
  23. Yates, A.D.; Achuthan, P.; Akanni, W.; Allen, J.; Allen, J.; Alvarez-Jarreta, J.; Amode, M.R.; Armean, I.M.; Azov, A.G.; Bennett, R.; et al. Ensembl 2020. Nucleic Acids Res. 2020, 48, D682–D688. [Google Scholar] [CrossRef] [PubMed]
  24. Uhlén, M.; Fagerberg, L.; Hallström, B.M.; Lindskog, C.; Oksvold, P.; Mardinoglu, A.; Sivertsson, Å.; Kampf, C.; Sjöstedt, E.; Asplund, A.; et al. Tissue-Based Map of the Human Proteome. Science 2015, 347, 1260419. [Google Scholar] [CrossRef] [PubMed]
  25. Monaco, G.; Lee, B.; Xu, W.; Mustafah, S.; Hwang, Y.Y.; Carré, C.; Burdin, N.; Visan, L.; Ceccarelli, M.; Poidinger, M.; et al. RNA-Seq Signatures Normalized by MRNA Abundance Allow Absolute Deconvolution of Human Immune Cell Types. Cell Rep. 2019, 26, 1627–1640. [Google Scholar] [CrossRef] [PubMed]
  26. Schmiedel, B.J.; Singh, D.; Madrigal, A.; Valdovino-Gonzalez, A.G.; White, B.M.; Zapardiel-Gonzalo, J.; Ha, B.; Altay, G.; Greenbaum, J.A.; McVicker, G.; et al. Impact of Genetic Polymorphisms on Human Immune Cell Gene Expression. Cell 2018, 175, 1701–1715. [Google Scholar] [CrossRef] [PubMed]
  27. Ward, L.D.; Kellis, M. HaploReg: A Resource for Exploring Chromatin States, Conservation, and Regulatory Motif Alterations within Sets of Genetically Linked Variants. Nucleic Acids Res. 2012, 40, D930–D934. [Google Scholar] [CrossRef] [PubMed]
  28. Oscanoa, J.; Sivapalan, L.; Gadaleta, E.; Dayem Ullah, A.Z.; Lemoine, N.R.; Chelala, C. SNPnexus: A Web Server for Functional Annotation of Human Genome Sequence Variation (2020 Update). Nucleic Acids Res. 2020, 48, W185–W192. [Google Scholar] [CrossRef] [PubMed]
  29. Roadmap Epigenomics Consortium; Kundaje, A.; Meuleman, W.; Ernst, J.; Bilenky, M.; Yen, A.; Heravi-Moussavi, A.; Kheradpour, P.; Zhang, Z.; Wang, J.; et al. Integrative Analysis of 111 Reference Human Epigenomes. Nature 2015, 518, 317–329. [Google Scholar] [CrossRef]
  30. Dunham, I.; Kundaje, A.; Aldred, S.F.; Collins, P.J.; Davis, C.A.; Doyle, F.; Epstein, C.B.; Frietze, S.; Harrow, J.; Kaul, R.; et al. An Integrated Encyclopedia of DNA Elements in the Human Genome. Nature 2012, 489, 57–74. [Google Scholar] [CrossRef]
  31. Myers, T.A.; Chanock, S.J.; Machiela, M.J. LDlinkR: An R Package for Rapidly Calculating Linkage Disequilibrium Statistics in Diverse Populations. Front. Genet. 2020, 11, 157. [Google Scholar] [CrossRef]
  32. Munz, M.; Tönnies, S.; Balke, W.T.; Simon, E. Multidimensional Gene Search with Genehopper. Nucleic Acids Res. 2015, 43, W98–W103. [Google Scholar] [CrossRef] [PubMed]
  33. Westra, H.J.; Peters, M.J.; Esko, T.; Yaghootkar, H.; Schurmann, C.; Kettunen, J.; Christiansen, M.W.; Fairfax, B.P.; Schramm, K.; Powell, J.E.; et al. Systematic Identification of Trans EQTLs as Putative Drivers of Known Disease Associations. Nat. Genet. 2013, 45, 1238–1243. [Google Scholar] [CrossRef] [PubMed]
  34. Fehrmann, R.S.N.; Jansen, R.C.; Veldink, J.H.; Westra, H.J.; Arends, D.; Bonder, M.J.; Fu, J.; Deelen, P.; Groen, H.J.M.; Smolonska, A.; et al. Trans-Eqtls Reveal That Independent Genetic Variants Associated with a Complex Phenotype Converge on Intermediate Genes, with a Major Role for the Hla. PLoS Genet. 2011, 7, e1002197. [Google Scholar] [CrossRef] [PubMed]
  35. Lonsdale, J.; Thomas, J.; Salvatore, M.; Phillips, R.; Lo, E.; Shad, S.; Hasz, R.; Walters, G.; Garcia, F.; Young, N.; et al. The Genotype-Tissue Expression (GTEx) Project. Nat. Genet. 2013, 45, 580–585. [Google Scholar] [CrossRef] [PubMed]
  36. Garnier, S.; Truong, V.; Brocheton, J.; Zeller, T.; Rovital, M.; Wild, P.S.; Ziegler, A.; Munzel, T.; Tiret, L.; Blankenberg, S.; et al. Genome-Wide Haplotype Analysis of Cis Expression Quantitative Trait Loci in Monocytes. PLoS Genet. 2013, 9, e1003240. [Google Scholar] [CrossRef] [PubMed]
  37. Codoni, V.; Blum, Y.; Civelek, M.; Proust, C.; Franzén, O.; Consortium, C.; Björkegren, J.L.M.; Le Goff, W.; Cambien, F.; Lusis, A.J.; et al. Preservation Analysis of Macrophage Gene Coexpression between Human and Mouse Identifies PARK2 as a Genetically Controlled Master Regulator of Oxidative Phosphorylation in Humans. G3 Genes. Genomes Genet. 2016, 6, 3361–3371. [Google Scholar] [CrossRef]
  38. Zeller, T.; Wild, P.; Szymczak, S.; Rotival, M.; Schillert, A.; Castagne, R.; Maouche, S.; Germain, M.; Lackner, K.; Rossmann, H.; et al. Genetics and beyond—the Transcriptome of Human Monocytes and Disease Susceptibility. PLoS ONE 2010, 5, e10693. [Google Scholar] [CrossRef]
  39. Landrum, M.J.; Chitipiralla, S.; Brown, G.R.; Chen, C.; Gu, B.; Hart, J.; Hoffman, D.; Jang, W.; Kaur, K.; Liu, C.; et al. ClinVar: Improvements to Accessing Data. Nucleic Acids Res. 2020, 48, D835–D844. [Google Scholar] [CrossRef]
  40. Meda, F.; Folci, M.; Baccarelli, A.; Selmi, C. The Epigenetics of Autoimmunity. Cell Mol. Immunol. 2011, 8, 226–236. [Google Scholar] [CrossRef]
  41. Bumiller-Bini, V.; Cipolla, G.A.; de Almeida, R.C.; Petzl-Erler, M.L.; Augusto, D.G.; Boldt, A.B.W. Sparking Fire under the Skin? Answers from the Association of Complement Genes with Pemphigus Foliaceus. Front. Immunol. 2018, 9, 695. [Google Scholar] [CrossRef]
  42. Grozinger, C.M.; Hassig, C.A.; Schreiber, S.L. Three Proteins Define a Class of Human Histone Deacetylases Related to Yeast Hda1p. Proc. Natl. Acad. Sci. USA 1999, 96, 4868–4873. [Google Scholar] [CrossRef] [PubMed]
  43. Miska, E.A.; Karlsson, C.; Langley, E.; Nielsen, S.J.; Pines, J.; Kouzarides, T. HDAC4 Deacetylase Associates with and Represses the MEF2 Transcription Factor. EMBO J. 1999, 18, 5099–5107. [Google Scholar] [CrossRef] [PubMed]
  44. Wang, A.H.; Bertos, N.R.; Vezmar, M.; Pelletier, N.; Crosato, M.; Heng, H.H.; Th’ng, J.; Han, J.; Yang, X.-J. HDAC4, a Human Histone Deacetylase Related to Yeast HDA1, Is a Transcriptional Corepressor. Mol. Cell Biol. 1999, 19, 7816–7827. [Google Scholar] [CrossRef] [PubMed]
  45. Youn, H.D.; Grozinger, C.M.; Liu, J.O. Calcium Regulates Transcriptional Repression of Myocyte Enhancer Factor 2 by Histone Deacetylase 4. J. Biol. Chem. 2000, 275, 22563–22567. [Google Scholar] [CrossRef] [PubMed]
  46. Paroni, G.; Mizzau, M.; Henderson, C.; Giannino; Sal, D.; Schneider, C.; Brancolini, C. Caspase-Dependent Regulation of Histone Deacetylase 4 Nuclear-Cytoplasmic Shuttling Promotes Apoptosis. Mol. Biol. Cell 2004, 15, 2804–2818. [Google Scholar] [CrossRef]
  47. Liu, Z.-G.; Smith, S.W.; McLaughlin, K.A.; Schwartz, L.M.; Osborne, B.A. Apoptotic Signals Delivered through the T-Cell Receptor of a T-Cell Hybrid Require the Immediate–Early Gene Nur77. Nature 1994, 367, 281–284. [Google Scholar] [CrossRef] [PubMed]
  48. Kumar, A.; Lin, Z.; SenBanerjee, S.; Jain, M.K. Tumor Necrosis Factor Alpha-Mediated Reduction of KLF2 Is Due to Inhibition of MEF2 by NF-ΚB and Histone Deacetylases. Mol. Cell Biol. 2005, 25, 5893–5903. [Google Scholar] [CrossRef]
  49. Jha, P.; Das, H. KLF2 in Regulation of NF-ΚB-Mediated Immune Cell Function and Inflammation. Int. J. Mol. Sci. 2017, 18, 2383. [Google Scholar] [CrossRef]
  50. Bumiller-Bini, V.; Cipolla, G.A.; Spadoni, M.B.; Augusto, D.G.; Petzl-Erler, M.L.; Beltrame, M.H.; Boldt, A.B.W. Condemned or Not to Die? Gene Polymorphisms Associated with Cell Death in Pemphigus Foliaceus. Front. Immunol. 2019, 10, 2416. [Google Scholar] [CrossRef]
  51. Hoch, V.B.B.; Schneider, L.; Pumpe, A.E.; Lüders, E.; Hundt, J.E.; Boldt, A.B.W. Marked to Die-Cell Death Mechanisms for Keratinocyte Acantholysis in Pemphigus Diseases. Life 2022, 12, 329. [Google Scholar] [CrossRef]
  52. Hatzi, K.; Philip Nance, J.; Kroenke, M.A.; Bothwell, M.; Haddad, E.K.; Melnick, A.; Crotty, S. BCL6 Orchestrates Tfh Cell Differentiation via Multiple Distinct Mechanisms. J. Exp. Med. 2015, 212, 539–553. [Google Scholar] [CrossRef] [PubMed]
  53. Lemercier, C.; Brocard, M.P.; Puvion-Dutilleul, F.; Kao, H.Y.; Albagli, O.; Khochbin, S. Class II Histone Deacetylases Are Directly Recruited by BCL6 Transcriptional Repressor. J. Biol. Chem. 2002, 277, 22045–22052. [Google Scholar] [CrossRef] [PubMed]
  54. Li, Q.; Liu, Z.; Dang, E.; Jin, L.; He, Z.; Yang, L.; Shi, X.; Wang, G. Follicular Helper T Cells (Tfh) and IL-21 Involvement in the Pathogenesis of Bullous Pemphigoid. PLoS ONE 2013, 8, e68145. [Google Scholar] [CrossRef] [PubMed]
  55. Molineros, J.E.; Yang, W.; Zhou, X.J.; Sun, C.; Okada, Y.; Zhang, H.; Chua, K.H.; Lau, Y.L.; Kochi, Y.; Suzuki, A.; et al. Confirmation of Five Novel Susceptibility Loci for Systemic Lupus Erythematosus (SLE) and Integrated Network Analysis of 82 SLE Susceptibility Loci. Hum. Mol. Genet. 2017, 26, 1205–1216. [Google Scholar] [CrossRef] [PubMed]
  56. McKinsey, T.A.; Kuwahara, K.; Bezprozvannaya, S.; Olson, E.N. Class II Histone Deacetylases Confer Signal Responsiveness to the Ankyrin-Repeat Proteins ANKRA2 and RFXANK. Mol. Biol. Cell 2006, 17, 438–447. [Google Scholar] [CrossRef] [PubMed]
  57. Piovezan, B.Z.; Petzl-Erler, M.L. Both Qualitative and Quantitative Genetic Variation of MHC Class II Molecules May Influence Susceptibility to Autoimmune Diseases: The Case of Endemic Pemphigus Foliaceus. Hum. Immunol. 2013, 74, 1134–1140. [Google Scholar] [CrossRef] [PubMed]
  58. Han, S.; Lu, J.; Zhang, Y.; Cheng, C.; Han, L.; Wang, X.; Li, N.; Liu, C.; Huang, B. Recruitment of Histone Deacetylase 4 by Transcription Factors Represses Interleukin-5 Transcription. Biochem. J. 2006, 400, 439–448. [Google Scholar] [CrossRef] [PubMed]
  59. Warren, S.J.P.; Arteaga, L.A.; Rivitti, E.A.; Aoki, V.; Hans-Filho, G.; Qaqish, B.F.; Lin, M.S.; Giudice, G.J.; Diaz, L.A. The Role of Subclass Switching in the Pathogenesis of Endemic Pemphigus Foliaceus. J. Invest. Dermatol. 2003, 120, 1–5. [Google Scholar] [CrossRef]
  60. Aoki, V.; Millikan, R.C.; Rivitti, E.A.; Hans-Filho, G.; Eaton, D.P.; Warren, S.J.P.; Li, N.; Hilario-Vargas, J.; Hoffmann, R.G.; Diaz, L.A. Environmental Risk Factors in Endemic Pemphigus Foliaceus (Fogo Selvagem). JID Symp. Proc. 2004, 9, 34–40. [Google Scholar] [CrossRef]
  61. Pan, J.; Zhao, L. Long Non-Coding RNA Histone Deacetylase 4 Antisense RNA 1 (HDAC4-AS1) Inhibits HDAC4 Expression in Human ARPE-19 Cells with Hypoxic Stress. Bioengineered 2021, 12, 2228–2237. [Google Scholar] [CrossRef]
  62. Hakimi, M.A.; Dong, Y.; Lane, W.S.; Speicher, D.W.; Shiekhattar, R. A Candidate X-Linked Mental Retardation Gene Is a Component of a New Family of Histone Deacetylase-Containing Complexes. J. Biol. Chem. 2003, 278, 7234–7239. [Google Scholar] [CrossRef]
  63. Hakimi, M.-A.; Bochar, D.A.; Chenoweth, J.; Lane, W.S.; Mandel, G.; Shiekhattar, R.; Rosenfeld, M.G. A Core-BRAF35 Complex Containing Histone Deacetylase Mediates Repression of Neuronal-Specific Genes. PNAS 2002, 99, 7420–7425. [Google Scholar] [CrossRef] [PubMed]
  64. Wang, Y.; Yan, S.; Yang, B.; Wang, Y.; Zhou, H.; Lian, Q.; Sun, B. TRIM35 Negatively Regulates TLR7- and TLR9-Mediated Type I Interferon Production by Targeting IRF7. FEBS Lett. 2015, 589, 1322–1330. [Google Scholar] [CrossRef] [PubMed]
  65. Hall, J.C.; Rosen, A. Type I Interferons: Crucial Participants in Disease Amplification in Autoimmunity. Nat. Rev. Rheumatol. 2010, 6, 40–49. [Google Scholar] [CrossRef] [PubMed]
  66. Meier-Soelch, J.; Jurida, L.; Weber, A.; Newel, D.; Kim, J.; Braun, T.; Lienhard Schmitz, M.; Kracht, M. RNAi-Based Identification of Gene-Specific Nuclear Cofactor Networks Regulating Interleukin-1 Target Genes. Front. Immunol. 2018, 9, 775. [Google Scholar] [CrossRef] [PubMed]
  67. Meier, K.; Brehm, A. Chromatin Regulation: How Complex Does It Get? Epigenetics 2014, 9, 1485–1495. [Google Scholar] [CrossRef] [PubMed]
  68. Xiong, Y.; Wang, L.; Giorgio, E.D.; Akimova, T.; Beier, U.H.; Han, R.; Trevisanut, M.; Kalin, J.H.; Cole, P.A.; Hancock, W.W. Inhibiting the Coregulator CoREST Impairs Foxp3+ Treg Function and Promotes Antitumor Immunity. J. Clin. Investig. 2020, 130, 1830–1842. [Google Scholar] [CrossRef]
  69. Weintraub, Y.; Cohen, S.; Yerushalmy-Feler, A.; Chapnik, N.; Tsameret, S.; Anafy, A.; Damari, E.; Ben-Tov, A.; Shamir, R.; Froy, O. Circadian Clock Gene Disruption in White Blood Cells of Patients with Celiac Disease. Biochimie 2023, 219, 51–54. [Google Scholar] [CrossRef]
  70. Cao, Q.; Zhao, X.; Bai, J.; Gery, S.; Sun, H.; Lin, D.C.; Chen, Q.; Chen, Z.; Mack, L.; Yang, H.; et al. Circadian Clock Cryptochrome Proteins Regulate Autoimmunity. Proc. Natl. Acad. Sci. USA 2017, 114, 12548–12553. [Google Scholar] [CrossRef]
  71. Tanu, T.; Taniue, K.; Imamura, K.; Onoguchi-Mizutani, R.; Han, H.; Jensen, T.H.; Akimitsu, N. HnRNPH1-MTR4 Complex-Mediated Regulation of NEAT1v2 Stability Is Critical for IL8 Expression. RNA Biol. 2021, 18, 537–547. [Google Scholar] [CrossRef]
  72. Kowalski, E.H.; Kneibner, D.; Kridin, K.; Amber, K.T. Serum and Blister Fluid Levels of Cytokines and Chemokines in Pemphigus and Bullous Pemphigoid. Autoimmun. Rev. 2019, 18, 526–534. [Google Scholar] [CrossRef] [PubMed]
  73. Verschueren, K.; Remacle, J.E.; Collart, C.; Kraft, H.; Baker, B.S.; Tylzanowski, P.; Nelles, L.; Wuytens, G.; Su, M.T.; Bodmer, R.; et al. SIP1, a Novel Zinc Finger/Homeodomain Repressor, Interacts with Smad Proteins and Binds to 5′-CACCT Sequences in Candidate Target Genes. J. Biol. Chem. 1999, 274, 20489–20498. [Google Scholar] [CrossRef] [PubMed]
  74. Postigo, A.A.; Dean, D.C.; Kipnis, D.M. Differential Expression and Function of Members of the Zfh-1 Family of Zinc Fingerhomeodomain Repressors. Proc. Natl. Acad. Sci. USA 2000, 97, 6391–6396. [Google Scholar] [CrossRef] [PubMed]
  75. Goossens, S.; Janzen, V.; Bartunkova, S.; Yokomizo, T.; Drogat, B.; Crisan, M.; Haigh, K.; Seuntjens, E.; Umans, L.; Riedt, T.; et al. The EMT Regulator Zeb2/Sip1 Is Essential for Murine Embryonic Hematopoietic Stem/Progenitor Cell Differentiation and Mobilization. Blood J. Am. Soc. Hematol. 2011, 117, 5620–5630. [Google Scholar] [CrossRef] [PubMed]
  76. Omilusik, K.D.; Adam Best, J.; Yu, B.; Goossens, S.; Weidemann, A.; Nguyen, J.V.; Seuntjens, E.; Stryjewska, A.; Zweier, C.; Roychoudhuri, R.; et al. Transcriptional Repressor ZEB2 Promotes Terminal Differentiation of CD8+ Effector and Memory T Cell Populations during Infection. J. Exp. Med. 2015, 212, 2027–2039. [Google Scholar] [CrossRef]
  77. Dominguez, C.X.; Amezquita, R.A.; Guan, T.; Marshall, H.D.; Joshi, N.S.; Kleinstein, S.H.; Kaech, S.M. The Transcription Factors ZEB2 and T-Bet Cooperate to Program Cytotoxic T Cell Terminal Differentiation in Response to LCMV Viral Infection. J. Exp. Med. 2015, 212, 2041–2056. [Google Scholar] [CrossRef]
  78. Scott, C.L.; Soen, B.; Martens, L.; Skrypek, N.; Saelens, W.; Taminau, J.; Blancke, G.; Van Isterdael, G.; Huylebroeck, D.; Haigh, J.; et al. The Transcription Factor Zeb2 Regulates Development of Conventional and Plasmacytoid DCs by Repressing Id2. J. Exp. Med. 2016, 213, 897–911. [Google Scholar] [CrossRef]
  79. Ghasemi, A.; Farazmand, A.; Hassanzadeh, V.; Poursani, S.; Soltani, S.; Akhtari, M.; Akhlaghi, M.; Farhadi, E.; Jamshidi, A.; Mahmoudi, M. Upregulation of KAT2B and ESCO2 Gene Expression Level in Patients with Rheumatoid Arthritis. Clin. Rheumatol. 2023, 42, 253–259. [Google Scholar] [CrossRef]
  80. Brockmann, D.; Lehmkühler, O.; Schmücker, U.; Esche, H. The Histone Acetyltransferase Activity of PCAF Cooperates with the Brahma/SWI2-Related Protein BRG-1 in the Activation of the Enhancer A of the MHC Class I Promoter. Gene 2001, 277, 111–120. [Google Scholar] [CrossRef]
  81. Liu, Y.; Bao, C.; Wang, L.; Han, R.; Beier, U.H.; Akimova, T.; Cole, P.A.; Dent, S.Y.R.; Hancock, W.W. Complementary Roles of Gcn5 and Pcaf in Foxp3+ T-Regulatory Cells. Cancers 2019, 11, 554. [Google Scholar] [CrossRef]
  82. de Jong, A.; de Jong, R.C.M.; Peters, E.A.; Arens, R.; Jukema, J.W.; de Vries, M.R.; Quax, P.H.A. P300/CBP Associated Factor (PCAF) Deficiency Enhances Diet-Induced Atherosclerosis in ApoE3*Leiden Mice via Systemic Inhibition of Regulatory T Cells. Front. Cardiovasc. Med. 2021, 7, 604821. [Google Scholar] [CrossRef] [PubMed]
  83. Bastiaansen, A.J.N.M.; Ewing, M.M.; De Boer, H.C.; Van Der Pouw Kraan, T.C.; De Vries, M.R.; Peters, E.A.B.; Welten, S.M.J.; Arens, R.; Moore, S.M.; Faber, J.E.; et al. Lysine Acetyltransferase PCAF Is a Key Regulator of Arteriogenesis. Arter. Thromb. Vasc. Biol. 2013, 33, 1902–1910. [Google Scholar] [CrossRef] [PubMed]
  84. Cleary, J.; Sitwala, K.V.; Khodadoust, M.S.; Kwok, R.P.S.; Mor-Vaknin, N.; Cebrat, M.; Cole, P.A.; Markovitz, D.M. P300/CBP-Associated Factor Drives DEK into Interchromatin Granule Clusters. J. Biol. Chem. 2005, 280, 31760–31767. [Google Scholar] [CrossRef] [PubMed]
  85. Sierakowska, H.; Williams, K.R.; Szer, I.S.; Szer, W. The Putative Oncoprotein DEK, Part Ofa Chimera Protein Associated with Acute Myeloid Leukaemia, Is an Autoantigen in Juvenile Rheumatoid Arthritis. Clin. Exp. Immunol. 1993, 94, 435–439. [Google Scholar] [CrossRef] [PubMed]
  86. Dong, X.; Wang, J.; Kabir, F.N.; Shaw, M.; Reed, A.M.; Stein, L.; Andrade, L.E.C.; Trevisani, V.F.M.; Miller, M.L.; Fujii, T.; et al. Autoantibodies to Dek Oncoprotein in Human Inflammatory Disease. Arthritis Rheum. Off. J. Am. Coll. Rheumatol. 2000, 43, 85–93. [Google Scholar] [CrossRef]
  87. Cai, Y.; Jin, J.; Swanson, S.K.; Cole, M.D.; Choi, S.H.; Florens, L.; Washburn, M.P.; Conaway, J.W.; Conaway, R.C. Subunit Composition and Substrate Specificity of a MOF-Containing Histone Acetyltransferase Distinct from the Male-Specific Lethal (MSL) Complex. J. Biol. Chem. 2010, 285, 4268–4272. [Google Scholar] [CrossRef] [PubMed]
  88. Badeaux, A.I.; Yang, Y.; Cardenas, K.; Vemulapalli, V.; Chen, K.; Kusewitt, D.; Richie, E.; Li, W.; Bedford, M.T. Loss of the Methyl Lysine Effector Protein PHF20 Impacts the Expression of Genes Regulated by the Lysine Acetyltransferase MOF. J. Biol. Chem. 2012, 287, 429–437. [Google Scholar] [CrossRef]
  89. Adams-Cioaba, M.A.; Li, Z.; Tempel, W.; Guo, Y.; Bian, C.; Li, Y.; Lam, R.; Min, J. Crystal Structures of the Tudor Domains of Human PHF20 Reveal Novel Structural Variations on the Royal Family of Proteins. FEBS Lett. 2012, 586, 859–865. [Google Scholar] [CrossRef]
  90. Klein, B.J.; Wang, X.; Cui, G.; Yuan, C.; Botuyan, M.V.; Lin, K.; Lu, Y.; Wang, X.; Zhao, Y.; Bruns, C.J.; et al. PHF20 Readers Link Methylation of Histone H3K4 and P53 with H4K16 Acetylation. Cell Rep. 2016, 17, 1158–1170. [Google Scholar] [CrossRef]
  91. Cui, G.; Park, S.; Badeaux, A.I.; Kim, D.; Lee, J.; Thompson, J.R.; Yan, F.; Kaneko, S.; Yuan, Z.; Botuyan, M.V.; et al. PHF20 Is an Effector Protein of P53 Double Lysine Methylation That Stabilizes and Activates P53. Nat. Struct. Mol. Biol. 2012, 19, 916–924. [Google Scholar] [CrossRef]
  92. Park, S.W.; Kim, J.; Oh, S.; Lee, J.; Cha, J.; Lee, H.S.; Kim, K.I.; Park, D.; Baek, S.H. PHF20 Is Crucial for Epigenetic Control of Starvation-Induced Autophagy through Enhancer Activation. Nucleic Acids Res. 2022, 50, 7856–7872. [Google Scholar] [CrossRef] [PubMed]
  93. Guelman, S.; Kozuka, K.; Mao, Y.; Pham, V.; Solloway, M.J.; Wang, J.; Wu, J.; Lill, J.R.; Zha, J. The Double-Histone-Acetyltransferase Complex ATAC Is Essential for Mammalian Development. Mol. Cell Biol. 2009, 29, 1176–1188. [Google Scholar] [CrossRef] [PubMed]
  94. Tzouanacou, E.; Tweedie, S.; Wilson, V. Identification of Jade1, a Gene Encoding a PHD Zinc Finger Protein, in a Gene Trap Mutagenesis Screen for Genes Involved in Anteroposterior Axis Development. Mol. Cell Biol. 2003, 23, 8553–8562. [Google Scholar] [CrossRef] [PubMed]
  95. Panchenko, M.V.; Zhou, M.I.; Cohen, H.T. Von Hippel-Lindau Partner Jade-1 Is a Transcriptional Co-Activator Associated with Histone Acetyltransferase Activity. J. Biol. Chem. 2004, 279, 56032–56041. [Google Scholar] [CrossRef]
  96. Foy, R.L.; Ihn, Y.S.; Chitalia, V.C.; Cohen, H.T.; Saksouk, N.; Cayrou, C.; Vaziri, C.; Côté, J.; Panchenko, M.V. Role of Jade-1 in the Histone Acetyltransferase (HAT) HBO1 Complex. J. Biol. Chem. 2008, 283, 28817–28826. [Google Scholar] [CrossRef]
  97. Avvakumov, N.; Lalonde, M.-E.; Saksouk, N.; Paquet, E.; Glass, K.C.; Landry, A.-J.; Doyon, Y.; Cayrou, C.; Robitaille, G.A.; Richard, D.E.; et al. Conserved Molecular Interactions within the HBO1 Acetyltransferase Complexes Regulate Cell Proliferation. Mol. Cell Biol. 2012, 32, 689–703. [Google Scholar] [CrossRef]
  98. Zhou, M.I.; Foy, R.L.; Chitalia, V.C.; Zhao, J.; Panchenko, M.V.; Wang, H.; Cohen, H.T. Jade-1, a Candidate Renal Tumor Suppressor That Promotes Apoptosis. Proc. Natl. Acad. Sci. USA 2005, 102, 11035–11040. [Google Scholar] [CrossRef]
Life 14 00060 g001
Figure 1. Flowchart of the study samples. The samples, composed of EPF and SPF patients from Brazil and Germany, respectively, and of controls were enrolled for experiments of genetic association with single nucleotide polymorphisms of (de)acetylation-related genes. EPF samples were also enrolled for RNA expression analyses of (de)acetylation-related genes.
Figure 1. Flowchart of the study samples. The samples, composed of EPF and SPF patients from Brazil and Germany, respectively, and of controls were enrolled for experiments of genetic association with single nucleotide polymorphisms of (de)acetylation-related genes. EPF samples were also enrolled for RNA expression analyses of (de)acetylation-related genes.
Life 14 00060 g001
Life 14 00060 g002
Figure 2. Interactions among proteins encoded by genes whose polymorphisms were associated with EPF susceptibility (HDAC4, GSE1, and PHF21A), have an eQTL effect (HNRNPH1), or which were differently expressed in EPF patients’ CD4+ T cells (KAT14, ZEB2, KAT2B). The interactions map was drawn using STRING (stringdb.org; accessed on 28 November 2023). The edges present a protein–protein association. Blue and purple edges indicate interactions with known co-occurrence and experimental evidence, respectively. Black and yellow edges indicate co-expression and textmining evidence, respectively. All the 3D protein structures are known or predicted.
Figure 2. Interactions among proteins encoded by genes whose polymorphisms were associated with EPF susceptibility (HDAC4, GSE1, and PHF21A), have an eQTL effect (HNRNPH1), or which were differently expressed in EPF patients’ CD4+ T cells (KAT14, ZEB2, KAT2B). The interactions map was drawn using STRING (stringdb.org; accessed on 28 November 2023). The edges present a protein–protein association. Blue and purple edges indicate interactions with known co-occurrence and experimental evidence, respectively. Black and yellow edges indicate co-expression and textmining evidence, respectively. All the 3D protein structures are known or predicted.
Life 14 00060 g002
Life 14 00060 g003
Figure 3. HDAC4_rs4852054*A was associated with increased EPF susceptibility, whereas GSE1_rs13339618*A and PHF21A_ rs4756055 A/A were associated with decreased EPF susceptibility. In CD4+ T cells from EPF patients, ZEB2, KAT2B, and PHF20 were overexpressed, and instead, KAT14 and JADE1 were underexpressed. By turning chromatin into a close state, deacetylation-related proteins whose polymorphisms were associated with increased EPF susceptibility (HDAC4) or which were overexpressed in patients’ CD4+ T cells (ZEB2) are thought to repress the expression of genes conferring protection to EPF in patients, whereas those associated with decreased EPF susceptibility (GSE1 and PHF21A) are thought to repress the expression of genes conferring susceptibility to the disease in controls. Conversely, in turning chromatin into an open state, acetylation-related proteins encoded by genes overexpressed in patients’ CD4+ T cells (KAT2B and PHF20) are supposed to induce expression of genes conferring susceptibility to the disease in EPF patients, while those underexpressed (JADE1 and KAT14) are supposed to induce expression of genes conferring protection to EPF in controls.
Figure 3. HDAC4_rs4852054*A was associated with increased EPF susceptibility, whereas GSE1_rs13339618*A and PHF21A_ rs4756055 A/A were associated with decreased EPF susceptibility. In CD4+ T cells from EPF patients, ZEB2, KAT2B, and PHF20 were overexpressed, and instead, KAT14 and JADE1 were underexpressed. By turning chromatin into a close state, deacetylation-related proteins whose polymorphisms were associated with increased EPF susceptibility (HDAC4) or which were overexpressed in patients’ CD4+ T cells (ZEB2) are thought to repress the expression of genes conferring protection to EPF in patients, whereas those associated with decreased EPF susceptibility (GSE1 and PHF21A) are thought to repress the expression of genes conferring susceptibility to the disease in controls. Conversely, in turning chromatin into an open state, acetylation-related proteins encoded by genes overexpressed in patients’ CD4+ T cells (KAT2B and PHF20) are supposed to induce expression of genes conferring susceptibility to the disease in EPF patients, while those underexpressed (JADE1 and KAT14) are supposed to induce expression of genes conferring protection to EPF in controls.
Life 14 00060 g003
Table 1. Genetic association analysis with endemic and sporadic pemphigus foliaceus.
Table 1. Genetic association analysis with endemic and sporadic pemphigus foliaceus.
GeneSNPEndemic Pemphigus Foliaceus Sporadic Pemphigus Foliaceus
MAF (%)ModelControlPatientsOR95% CIpMAF (%)ModelControlPatientsOR95% CIp
ControlsPatientsControlsPatients
HDAC4rs485205430.6739.82ADD119/269180/2721.525[1.14–2.05]0.005118.7123.03ADD52/22635/1171.299[0.80–2.11]0.2891
2q37.3G > a REC19/17532/1941.536[0.84–2.81]0.1647 REC05/13404/721.489[0.39–5.72]0.5621
intron 1 DOM78/14894/1001.794[1.21–2.67]0.0038 DOM47/9231/451.348[0.76–2.40]0.3097
GSE1rs1333961829.6419.42ADD115/27387/3610.571[0.41–0.80]0.001119.0623.65ADD53/22535/1331.257[0.81–1.96]0.3108
16q24.1G > a REC16/1785/2190.261[0.09–0.73]0.0102 REC09/13007/671.509[0.54–4.23]0.4339
regulatory region downstream DOM99/9582/1420.578[0.39–0.86]0.0064 DOM44/9528/461.314[0.73–2.37]0.3643
PHF21Ars475605538.0542.82ADD186/198172/2800.671[0.50–0.90]0.007945.3237.50ADD126/15257/950.740[0.50–1.10]0.1322
11p11.2g > A REC50/14226/2000.391[0.23–0.67]0.0007 REC31/10812/640.653[0.38–1.36]0.2559
intron 1 DOM136/56146/800.792[0.52–1.21]0.2831 DOM95/4445/310.672[0.38–1.20]0.1801
In bold, significant associations (p < 0.005). The minor alleles in our sample are given in lowercase; they are the reference for the associations. Abbreviations: ADD, additive model; CI, confidence interval; DOM, dominant model; GSE1, Gse1 coiled-coil protein; HDAC4, histone deacetylase 4; MAF, minor allele frequency; OR, odds ratio; PHF21A, PHD finger protein 21A; REC, recessive model; SNP, single nucleotide polymorphism.
Table 2. Alleles in high LD with HDAC4, GSE1, and PHF21A alleles associated with EPF and their eQTL effect in relevant tissues and cells for PF pathophysiology.
Table 2. Alleles in high LD with HDAC4, GSE1, and PHF21A alleles associated with EPF and their eQTL effect in relevant tissues and cells for PF pathophysiology.
SNPsr2eQTLTissues and Cells
HDAC4_rs4852054*A (Susceptibility)-HDAC4-AS1Blood (ES = 0.15; p = 1.5 × 10−4)
Skin not sun-exposed (ES = 0.17; p = 3.4 × 10−3)
HDAC4_rs4852053*C1.0HDAC4-AS1Blood (ES = 0.17, p = 4.7 × 10−5)
Skin sun-exposed (ES = 0.16; p = 5.0 × 10−3)
HDAC4_rs55932933*T0.99HDAC4-AS1Blood (ES = 0.17; p = 6.3 × 10−5)
Skin sun-exposed (ES = 0.16; p = 5.0 × 10−3)
HDAC4_rs56064197*T0.99HDAC4-AS1Blood (ES = 0.17; p = 5.8 × 10−5)
Skin sun-exposed (ES = 0.16; p = 5.0 × 10−3)
HDCA4_rs1476321*T0.99HDAC4-AS1Blood (ES = 0.15; p = 1.9 × 10−4)
Skin not sun-exposed (ES = 0.17; p = 3.4 × 10−3)
HDCA4_rs10182344*A0.99HDAC4-AS1Blood (ES = 0.14; p = 3.5 × 10−4)
HDCA4_rs58332998*A0.93--
GSE1_rs13339618*A (Protection) GSE1Blood (ES = 5.37; p = 7.9 × 10−8)
GSE1_rs13339626*A1.0GSE1Blood (ES = 5.37; p = 7.9 × 10−8)
GSE1_rs13329722*C1.0GSE1Blood (ES = 5.40; p = 6.9 × 10−8)
TRIM35Monocytes (EF = -; p = 9.2 × 10−6)
GSE1_rs7498141*A1.0GSE1Blood (EF = 5.40; p = 6.9 × 10−8)
TRIM35Monocytes (ES = -; p = 9.0 × 10−6)
GSE1_rs4843505*G1.0GSE1Blood (ES = 5.40; p = 6.9 × 10−8)
GSE1_rs13332576*A0.99GSE1Blood (ES = 5.31; p = 1.1 × 10−7)
GSE1_rs8060638*A0.89--
GSE1_rs28671512*T0.99--
PHF21A_rs4756055*A/A (Protection)-CRY2Blood (ES = 3.64; p = 2.7 × 10−4)
PEX16Blood (ES = −3.16; p = 1.0 × 10−3)
PHF21ABlood (ES = -; p = 2.8 × 10−4)
HNRNPH1Macrophages (ES = 0.08; p = 5.8 × 10−6)
PHF21A_rs1976182*A0.98HNRNPH1Macrophages (ES = 0.08; p = 5.4 × 10−6)
PHF21A_rs923530*T0.87CRY2Monocytes (ES = −0.04, p = 3.8 × 10−4)
HNRNPH1Macrophages (ES = 0.08, p = 7.6 × 10−6)
SLC37A1Macrophages (EF = −0.05, p = 2.3 × 10−6)
PHF21A_rs11374563*T0.87--
PHF21A_rs7109480*T0.93CRY2Blood (ES = 3.46, p = 5.4 × 10−4)
HNRNPH1Macrophages (EF = 0.08; p = 6.8 × 10−6)
PHF21A_rs7107550*C0.86CRY2Monocytes (ES = −0.04, p = 3.2 × 10−4)
HNRNPH1Macrophages (ES = 0.08, p = 9.3 × 10−6)
SLC37A1Macrophages (EF = −0.05, p = 4.3 × 10−6)
PHF21A_rs950105*A1.0CRY2Blood (ES = 3.41; p = 6.4 × 10−4)
HNRNPH1Macrophages (EF = 0.08; p = 9.1 × 10−6)
PHF21A_rs74366855*A0.87CRY2Monocytes (ES = −0.04, p = 5.8 × 10−4)
HNRNPH1Macrophages (ES = 0.08, p = 8.0 × 10−6)
SLC37A1Macrophages (EF = −0.05, p = 5.7 × 10−6)
High LD (r2 > 0.80) between SNPs in the European population was evaluated with LDLink [31]. eQTL effect (p ≤ 5.0 × 10−3) was observed in Genehopper Qtlizer [32], which compiles information from Blood eQTL [33], GRASP 2 Catalog [34], GTEx v8 [35], The Cardiogenics Project [36,37], and Zeller et al. [38]. “-” indicates missing/unfound data. Abbreviations: CRY2, cryptochrome circadian regulator 2; ES, effect size; EPF, endemic pemphigus foliaceus; eQTL, expression quantitative trait loci; GSE1, Gse1 coiled-coil protein; HDAC4, histone deacetylase 4; HDAC4-AS1, HDAC4 antisense RNA 1; HNRNPH1, heterogeneous nuclear ribonucleoprotein H1; LD, linkage disequilibrium; PEX16, peroxisomal biogenesis factor 16; PF, pemphigus foliaceus; PHF21A, PHD finger protein 21A; SLC37A1, solute carrier family 37 member 1; SNP, single nucleotide polymorphism; TRIM35, tripartite motif containing 35.
Table 3. (De)acetylation-related genes differently expressed in CD4+ T cells of endemic pemphigus foliaceus patients compared to controls.
Table 3. (De)acetylation-related genes differently expressed in CD4+ T cells of endemic pemphigus foliaceus patients compared to controls.
GeneFull NameFold Changep-Value
ZEB2zinc finger E-box binding homeobox 21.086.2 × 10−6
KAT2Blysine acetyltransferase 2B0.552.8 × 10−5
PHF20PHD finger protein 200.418.6 × 10−5
KAT14lysine acetyltransferase 14−0.592.0 × 10−4
JADE1jade family PHD finger 1−0.312.6 × 10−3
PF patients presented new skin lesions coming up constantly, before starting oral corticosteroid treatment. Controls were from endemic areas for PF.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Sulzbach Denardin, M.; Bumiller-Bini Hoch, V.; Salviano-Silva, A.; Lobo-Alves, S.C.; Adelman Cipolla, G.; Malheiros, D.; Augusto, D.G.; Wittig, M.; Franke, A.; Pföhler, C.; et al. Genetic Association and Differential RNA Expression of Histone (De)Acetylation-Related Genes in Pemphigus Foliaceus—A Possible Epigenetic Effect in the Autoimmune Response. Life 2024, 14, 60. https://doi.org/10.3390/life14010060

AMA Style

Sulzbach Denardin M, Bumiller-Bini Hoch V, Salviano-Silva A, Lobo-Alves SC, Adelman Cipolla G, Malheiros D, Augusto DG, Wittig M, Franke A, Pföhler C, et al. Genetic Association and Differential RNA Expression of Histone (De)Acetylation-Related Genes in Pemphigus Foliaceus—A Possible Epigenetic Effect in the Autoimmune Response. Life. 2024; 14(1):60. https://doi.org/10.3390/life14010060

Chicago/Turabian Style

Sulzbach Denardin, Maiara, Valéria Bumiller-Bini Hoch, Amanda Salviano-Silva, Sara Cristina Lobo-Alves, Gabriel Adelman Cipolla, Danielle Malheiros, Danillo G. Augusto, Michael Wittig, Andre Franke, Claudia Pföhler, and et al. 2024. "Genetic Association and Differential RNA Expression of Histone (De)Acetylation-Related Genes in Pemphigus Foliaceus—A Possible Epigenetic Effect in the Autoimmune Response" Life 14, no. 1: 60. https://doi.org/10.3390/life14010060

APA Style

Sulzbach Denardin, M., Bumiller-Bini Hoch, V., Salviano-Silva, A., Lobo-Alves, S. C., Adelman Cipolla, G., Malheiros, D., Augusto, D. G., Wittig, M., Franke, A., Pföhler, C., Worm, M., van Beek, N., Goebeler, M., Sárdy, M., Ibrahim, S., Busch, H., Schmidt, E., Hundt, J. E., Petzl-Erler, M. L., & Beate Winter Boldt, A. (2024). Genetic Association and Differential RNA Expression of Histone (De)Acetylation-Related Genes in Pemphigus Foliaceus—A Possible Epigenetic Effect in the Autoimmune Response. Life, 14(1), 60. https://doi.org/10.3390/life14010060

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop