Next Article in Journal
Novel Regulatory Roles of Hydrogen Sulfide in Health and Disease
Previous Article in Journal
Deletion of Growth Hormone Secretagogue Receptor in Kisspeptin Neurons in Female Mice Blocks Diet-Induced Obesity
Previous Article in Special Issue
Anti-Inflammatory microRNAs for Treating Inflammatory Skin Diseases
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Transcriptome Meta-Analysis Confirms the Hidradenitis Suppurativa Pathogenic Triad: Upregulated Inflammation, Altered Epithelial Organization, and Dysregulated Metabolic Signaling

by
Ana Sofia Lima Estevao de Oliveira
1,*,†,
Giovanna Bloise
2,3,†,
Chiara Moltrasio
4,5,
Antonio Coelho
3,
Almerinda Agrelli
6,
Ronald Moura
7,
Paola Maura Tricarico
7,
Stéphane Jamain
8,
Angelo Valerio Marzano
4,9,
Sergio Crovella
10,† and
Lucas André Cavalcanti Brandão
2,†
1
Laboratory of Immunopathology Keizo Asami-LIKA, Federal University of Pernambuco, Recife 50670-901, Brazil
2
Department of Pathology, Federal University of Pernambuco, Recife 50670-901, Brazil
3
Hospital Israelita Albert Einstein, São Paulo 05652-000, Brazil
4
Dermatology Unit, Fondazione IRCCS Ca’ Granda Ospedale Maggiore Policlinico, 20122 Milan, Italy
5
Department of Medical Surgical and Health Sciences, University of Trieste, 34137 Trieste, Italy
6
Laboratory of Nanostructured Materials (LMNANO), Center for Strategic Technologies Northeastern (CETENE), Av. Prof. Luís Freire, 1-Cidade Universitária, Recife 50740-545, Brazil
7
Department of Advanced Diagnostics, Institute for Maternal and Child Health-IRCCS “Burlo Garofolo”, 34137 Trieste, Italy
8
Translational Neuropsychiatry, Univ. Paris Est Créteil, Inserm, IMRB, 94010 Créteil, France
9
Department of Pathophysiology and Transplantation, Università degli Studi di Milano, 20122 Milan, Italy
10
Biological Science Program, Department of Biological and Environmental Sciences, College of Arts and Sciences, University of Qatar, Doha 2713, Qatar
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Biomolecules 2022, 12(10), 1371; https://doi.org/10.3390/biom12101371
Submission received: 15 August 2022 / Revised: 15 September 2022 / Accepted: 20 September 2022 / Published: 25 September 2022

Abstract

:
Hidradenitis suppurativa (HS) is an inflammatory skin condition clinically characterized by recurrent painful deep-seated nodules, abscesses, and sinus tracks in areas bearing apocrine glands, such as axillae, breasts, groins, and buttocks. Despite many recent advances, the pathophysiological landscape of HS still demands further clarification. To elucidate HS pathogenesis, we performed a meta-analysis, set analysis, and a variant calling on selected RNA-Sequencing (RNA-Seq) studies on HS skin. Our findings corroborate the HS triad composed of upregulated inflammation, altered epithelial differentiation, and dysregulated metabolism signaling. Upregulation of specific genes, such as KRT6, KRT16, serpin-family genes, and SPRR3 confirms the early involvement of hair follicles and the impairment of barrier function in HS lesioned skin. In addition, our results suggest that adipokines could be regarded as biomarkers of HS and metabolic-related disorders. Finally, the RNA-Seq variant calling identified several mutations in HS patients, suggesting potential new HS-related genes associated with the sporadic form of this disease. Overall, this study provides insights into the molecular pathways involved in HS and identifies potential HS-related biomarkers.

1. Introduction

Hidradenitis suppurativa (HS), also known as acne inversa, is a chronic inflammatory skin disease involving the pilosebaceous unit, with a worldwide prevalence varying from 0.03% to 4% [1,2,3,4,5,6,7]. HS lesioned skin is typically associated with apocrine gland-bearing body areas, such as axillae, breasts, groins, and buttocks [1,2]. HS is clinically characterized by painful and recurrent deep-seated inflammatory nodules and abscesses [3] that, with disease progression, may evolve into dermal tunnels with malodorous discharge and disfiguring scarring, severely affecting patients’ quality of life [4,5].
HS pathogenesis is still poorly understood, although hereditary factors have demonstrated the ability to increase risks associated with the development of the skin condition. [6]. It has been proposed that the primum movens of the disease are represented by keratin follicular plugging progressing to follicle rupture, with subsequent upregulation of the immune system response that leads to an overexpression of pro-inflammatory cytokines/chemokines, as well as several other pro-inflammatory mediators [7,8,9]. Several genetic changes have been involved in HS pathogenesis, including mutations in the γ-secretase complex genes, both in familial and syndromic forms of HS [10,11,12]. On the other hand, the γ-secretase gene mutations only occur in roughly 6% of sporadic cases of HS [13]. The γ-secretase complex has the notch signaling pathway as a substrate, which has been shown to alter the keratinocyte differentiation program leading to uncontrolled cell proliferation when downregulated [14]. Moreover, both γ-secretase genes and notch signaling pathway play a crucial role in the development of epidermal cysts and comedones, two phenotypic traits of HS [15]. With several new gene pathogenic variants identified, HS can be considered a multifactorial, polygenic, autoinflammatory condition [12].
In the -omic era, the next-generation sequencing (NGS) approach made major contributions to the study of HS [12,16]. Herein, we performed a meta-analysis of selected RNA sequencing (RNA-Seq) studies, comparing HS lesioned and healthy skin. With this integrated meta-analysis approach, it [17] was possible to combine results from independent studies to increase statistical power and obtain a more precise estimate of differentially expressed genes (DEGs). In addition, we carried out pathway analysis to characterize the DEGs in HS pathogenesis. Lastly, a variant calling was performed allowing the identification of expression quantitative trait loci (eQTL) associated with DEGs. Through these analyses performed we explored HS pathogenesis and created a link between candidate genes and their transcripts, genotypes, and clinical phenotypes.

2. Materials and Methods

2.1. Selection of Publicly Available Studies

We searched Sequence Read Archive (SRA) [18] and Gene Expression Omnibus (GEO) [19] datasets on 11 July 2022, to find studies involving RNA-Seq of HS tissue samples. The search keywords were: ((Hidradenitis Suppurativa) OR (Acne Inversa)) AND ((RNA-Seq) OR (Transcriptome)) [20]. After removing duplicates, we included studies that met the criteria of: experiments carried out in HS patients’ skin (Homo sapiens organism filter); had matching healthy skin controls; and had raw data (fastq files) available for each sample. Biopsies collected from lesioned skin of HS patients were considered cases, whereas biopsies collected from the skin of patients without HS, or any other skin condition, were accepted as controls.

2.2. RNA-Seq Data Collection, Processing, and Analysis

Figure 1 summarizes the methodology. Data collection, processing, and analysis were conducted as previously described by [21,22]. Briefly, SR Adb package [23] for R software version 4.1.0 [24] was used to download all raw. fastq files. Only the sequencing reads that met the study criteria were downloaded. Then, we reprocessed all reads using the same standard workflow to avoid bias due to heterogeneous pipelines of the original studies. The workflow used Trimmomatic v0.39 [25] to trim Illumina adapters, and to exclude low-quality reads (Q < 30) and reads counting shorter than 25 bases (length < 25). Then, the remaining reads were mapped on the National Center for Biotechnology (NCBI) human GRCh38/hg38 reference genome and sorted by coordinates using STAR aligner [26]. Aligned reads (BAM files) were imported into R software and processed with an annotation file from the reference genome with Rsubread package [27]. As a result, a gene counts table was created for each sample. These tables were converted into DESeq2 package objects, where we screened for DEGs [28]. Genes with log2(fold change) > 1 and false discovery rate (FDR)-adjusted p-values < 0.05 were acknowledged as statistically significant DEGs.

2.3. Meta-Analysis

We integrated the results produced by the independent groups with a meta-analysis using the RankProd package for R software. The package performs a nonparametric approach based on ranks of fold changes (FC) that detect differentially consistently expressed genes (mDEGs) from independent and replicated experiments. Overall, the meta-analysis approach was based on Lee et al. methodology (2019) [21].

2.4. Pathway Analysis

We performed a pathway enrichment analysis with the mDEGs. We searched the REACTOME database, using ReactomePA package [24] for R software, and PANTHER v17.0 [29,30]. As in the previous analyses, results with FDR-adjusted enrichment test p-value < 0.05 were considered significant.

2.5. RNA-Seq Variant Calling

We performed a variant calling sourcing from the aligned RNA-Seq data. Post-processing was applied using Genome Analysis Toolkit v4.1.2.0 [31] to the aligned reads: duplicates were marked and removed with Mark Duplicates and Split NCigar Reads, and base recalibration was completed with Base Reca librator and Apply BQSR. Variant Calling was performed using Strelka v2.9.10 [32], and single nucleotide polymorphisms (SNPs) and insertions and deletions (INDELs) were separated in two variant calling format (VCF) files using select variants. The resulting variants were removed using Variant Filtration if SNPs had QD < 2.0, FS > 60.0, MQ < 40.0, MQ Rank Sum < −12.5, Read Pos Rank Sum < −8.0 and GQ < 20.0, and if INDELs had QD < 2.0, FS > 200.0, Read Pos Rank Sum < −20.0 and GQ < 20.0. Annotation was performed using Annovar v2019Oct24 [33]. Finally, variants that were present in HS patients were only kept if they had depth (DP > 10) and quality score (QUAL > 30).
To infer the connection between genetic variants and HS pathogenesis, we conducted several approaches. Initially we searched for variants in genes previously associated with HS. Then, we sought variants in new potential HS genes. To unravel the new ones, we selected variants based on some filters detailed as follows. Estimates of HS prevalence vary between 0.03% and 4% among different populations [34,35]. Therefore, it has not been well defined if HS is a rare or a neglected disease [36]. Since rare variants are commonly responsible for disease appearance, herein, for analysis purposes only, we considered HS a rare skin condition (AF ≤ 0.01). Thus, potentially deleterious variants found in genes previously associated with HS or with a CADD (phred) score > 15, possessing clinical significance different from benign or likely benign, based on the Clin Var database, and having statistically significant (p-value < 0.05) differences in the genotype were selected as potential HS genetic markers. The median value of the RNA expression level was calculated for each genotype; large differences between the genotypes were considered statistically significant according to the Wilcoxon test.
Figure 1. Analysis flowchart: 1: Publicly available studies were selected from the Sequence Read Archive (SRA) and Gene Expression Omnibus (GEO); 2: Three studies passed our selection criteria and were downloaded. SRAdb package for R software version 4.1.0 was used to download all raw. fastq files. Trimmomatic v0.39 was used to trim Illumina adapters and to exclude reads. Resulting gene counts tables were converted into DESeq2 package objects [37,38,39]; 3: Lists of DEGs found through the meta-analysis (mDEGs) or the set analysis (saDEGs) (See supplementary data) were created; 4: List of common DEGs (cDEGs) found in both analyses was made (See supplementary data); 5: RNA-Seq variant calling analysis was performed by Strelka v2.9.10; 6: A list of potential HS-related variants originated based on deleterious capacity, allele frequency, and significant genotypic variation based on Wilcoxon Test (p value < 0.5).
Figure 1. Analysis flowchart: 1: Publicly available studies were selected from the Sequence Read Archive (SRA) and Gene Expression Omnibus (GEO); 2: Three studies passed our selection criteria and were downloaded. SRAdb package for R software version 4.1.0 was used to download all raw. fastq files. Trimmomatic v0.39 was used to trim Illumina adapters and to exclude reads. Resulting gene counts tables were converted into DESeq2 package objects [37,38,39]; 3: Lists of DEGs found through the meta-analysis (mDEGs) or the set analysis (saDEGs) (See supplementary data) were created; 4: List of common DEGs (cDEGs) found in both analyses was made (See supplementary data); 5: RNA-Seq variant calling analysis was performed by Strelka v2.9.10; 6: A list of potential HS-related variants originated based on deleterious capacity, allele frequency, and significant genotypic variation based on Wilcoxon Test (p value < 0.5).
Biomolecules 12 01371 g001

2.6. Statistics

The Shapiro-Wilk test was used to assess if the counts were normally distributed. To uncover if genotypes were associated with differential RNA expression, we used an appropriate test for normally distributed counts based on the number of observed genotypes (normal and two genotypes: Sztudent’s t-test for independent samples; not normal and two genotypes: Wilcoxon-Mann-Whitney test; normal and three genotypes: ANOVA; not normal and three genotypes: Kruskal-Wallis test. In the case of three-genotypes tests, we performed a post-hoc paired Wilcoxon test with FDR correction to identify which genotype was responsible for the overall difference).

3. Results

3.1. Dataset Selection

The search strategy conducted by three authors independently (A.S.L.E.O, G.B., and A.A.) resulted in 25 datasets. Among these, seven studies used microarrays; eight used RNA-Seq; six used quantitative real time polymerase chain reaction (PCR); three used single cell (sc) RNA-Seq; and one used target capture. Since we were searching for RNA-Seq studies only, 17 studies were excluded. Out of the eight RNA-Seq studies, six contained analysis of skin biopsies. Of these six studies, two of them did not have matching healthy skin control, and another analyzed familial HS cases only. Therefore, only three studies were included in our meta-analysis and set analysis. Herein, a total of 51 unique HS lesioned samples, and 30 unique healthy skin samples were considered (Table 1).
In the first study, Lowe et al. (SRAdb id: SRP273719) [37] aimed to uncover genes related to HS immunogenesis and their influence on patients’ responsiveness to anti-tumor necrosis factor alpha (TNF-α) therapy. This study was conducted through transcriptome sequencing of lesioned and non-lesioned skin samples before and after anti-TNF-α therapy. Herein, we focused our attention on the RNA-Seq from samples collected before anti-TNF-α therapy. In summary, 26 HS lesioned skin samples and 16 healthy skin controls were included in our analyses.
In the second study, Gudjonsson et al. (SRAdb id: SRP272659) [38] focused on outlining the major dysregulated cell types and inflammatory pathways in HS. In order to evaluate abnormal inflammatory pathways, they performed skin and whole blood RNA-Seq, and scRNA-Seq. In addition, they compared their findings with RNA-Seq data from 28 psoriasis patients and 32 atopic dermatitis patients. In summary, they found 4797 DEGs within the HS skin RNA-Seq compared to healthy controls, and 332 DEGs within HS patients’ whole blood. From this study, we incorporated into our analysis 32 skin RNA-Seq samples; 22 from HS lesioned skin and 10 from healthy skin controls.
In the third study, Witte-Händel et al. (SRAdb id: SRP168930) [39] investigated HS skin lesions’ cytokine milieu and found IL-1β to be highly expressed. First, they quantified the expression levels of approximately 30 mediators in HS lesioned skin and compared the results with the expression patterns from these same mediators in healthy skin control and in psoriasis lesioned skin. After finding IL-1β to be upregulated in HS, they compared the results of HS lesioned skin’s RNA-Seq with the transcriptome of various cell types isolated from healthy patients’ biopsies that were exposed to IL-1β and observed transcriptomic overlapping profiles. RT-qPCR and ELISA analyses were also conducted. Our analysis included samples from three HS patients’ lesioned skin biopsy and four healthy individuals skin biopsy not previously exposed with IL-1β.

3.2. Meta-Analisis

We performed a meta-analysis, which resulted in 6445 annotated genes, and 3628 mDEGs with an FDR < 0.05 and |log2(fold change)| > 1. From these mDEGs, 1699 were upregulated and 1929 were downregulated (Table S1). Briefly, upregulated mDEGs that play a role in the maintenance of skin inflammation and structure, namely DEFB4A/B [40,41], TCN1 [9], MMP1 [42], S100A7A [43], and MMP3 [44] appeared at the top of the overall meta-analysis rank product list. Alongside these mDEGs, other immune-related genes such as TNF-α, IFN-γ, IL-1β, IL-36, IL-17, MZB1, CD19, CD79A, CXCL13, RETN [45], and RARRES2 [46], and skin structure-related genes particularly PI3 [47], SERPINB3/4 [48], SERPINA1 [49], KRT6/16 [50], and SPRR1/2/3 [51] were observed. Conversely, the highest positions on the overall rank product list for the down-regulated mDEGs included THRSP, AWAT2, DCD, UGT3A2, and WIF1, genes associated with energetic metabolism and skin’s barrier function [52,53,54,55,56]. Down-regulation of these mDEGs as well as other genes, especially AQP2/4/5/6/7/8, FOXA1, and ADIPOQ, may impair these pathways.
Pathway enrichment analysis of the mDEGs by reactome resulted in 76 statistically significant enriched pathways (Table S2). We then selected the 20 pathways with the lowest FDR values and summarized them in Table 2. Pathways such as keratinization (FDR = 0.0049, gene ratio = 54/1663), extracellular matrix (ECM) organization (FDR = 2.42 × 1014, gene ratio 109/1663), neutrophil degranulation (FDR = 5.25 × 109, gene ratio = 132/1663), PD-1 signaling (FDR = 1.5532 × 109, gene ratio = 18/1663), interferon gamma signaling (FDR = 0.00003, gene ratio = 58/1663), signaling by interleukins (FDR = 2.3270 × 106, gene ratio = 118/1663), complement cascade (p value = 0.00003, gene ratio = 25/1663), antimicrobial peptides (p value = 3.476 × 106, gene ratio = 30/1663), cell surface interactions at the vascular wall (p value = 8.1106 × 106, gene ratio = 47/1663), and g-protein coupled receptor (GPCR) ligand binding (p value = 5.2463 × 109, gene ratio = 126/1663) appeared enriched.
In addition, a pathway enrichment analysis by PANTHER database was conducted. Summarily, pathways with the highest percentage of upregulated mDEGs enclosed inflammation mediated by chemokine and cytokine, integrin and interleukin signaling, T and B cell activation, angiogenesis, and apoptosis (Figure S1). Where, pathways containing the highest percentage of down-regulated mDEGs included wnt signaling, angiogenesis, G-protein signaling, gonadotropin-releasing hormone receptor, and cadherin signaling (Figure S2).
We then performed a set analysis to extract DEGs commonly found in all three selected studies (saDEGs), as well as in both analyses (cDEGs). These analyses can be located in the supplementary data.

3.3. RNA-Seq Variant Calling

We performed the variant calling in the HS lesioned samples and controls. The pipeline identified 323,827 unique variants across all samples. Among them, 190,845 variants were annotated as known SNPs and Indels.
After filtering, a total of 73,134 genetic variants present in 12,665 genes were found only in patients, and thus selected as candidate variants for HS pathogenesis. Up to 39,709 novel variants not present on dbSNP and Gnomad were found and marked by an asterisk at the end of the information regarding amino acids change and reference gene at the column “AAChange.refGene”. These results are summarized in Table S3. Briefly, 63,235 single nucleotide variants (SNVs) and 9899 insertion/deletion variants (Indels) were detected across the patients’ transcriptomes. They were distributed as follows: 39,577 in 3′-UTR, 3577 in 5′-UTR, 2325 in splicing sites, 25,577 exonic, 22 non-coding RNA splicing, and 2056 non-coding RNA exonic. Among the exonic variants present only in patients, 11,890 were synonymous, 12,266 were nonsynonymous, 35 were start-loss, 140 were stop gain, 203 were non-frameshift deletions, 205 were non-frameshift insertions, 120 were frameshift deletions, and 651 were frameshift insertions.
The variant calling analysis identified 34 variants in genes of the γ-secretase complex. This complex, composed of 4 genes (presenilins (PSEN), nicastrin (NCSTN), anterior pharynx defective 1 (APH1), and presenilin enhancer 2 (PSENEN)), modulates NOTCH pathway by cleaving and releasing NOTCH’s intracellular domain [14], thus, being important for skin and immune response’s homeostasis and HS pathogenesis [57,58]. The 34 mutations were distributed as follows: 12 in the NCSTN gene, 3 in the APH1A, 6 in the APH1B, 8 in the PSEN1, 3 in the PSEN2, and 2 in the PSENEN. One variant was already reported in HS patients: a stop gain mutation NM_172341.4: c.168T>G (p.Tyr56X) of the PSENEN gene previously associated with concomitant HS and Dowling Degos Disease (DDD), a rare genodermatosis, classically characterized by acquired reticular hyperpigmentation in flexural sites [59]. Another non-synonymous variant NM_000243.2:c.2177T>C (p.Val726Ala) in the MEFV gene, previously associated with HS, has also been spotted [60]. In Table 3 we summarized these two mutations, along with other variants found in coding regions of genes already linked with HS [60,61,62].
Apart from variants in genes already considered related to the disease, a second table containing potential HS new genetic markers was created (Table 4). Overall, variants were observed in genes associated with cellular metabolism, immune response, and skin homeostasis. For instance, ACSF3 plays a crucial role in fatty acid synthesis [63], and its dysregulation may affect the ability of the body to properly process certain proteins and lipids. In conjunction, we also observed mutations in genes, such as BTN2A1, FNIP2, AKR1C3, ALDH6A1, YTHDF1, and GALNT7 that are similarly involved in the metabolism of lipids, proteins, fatty-acids, and sterols. Variants associated with immune response were spotted in genes such as GSDMD, SIT1, WDR92, and GPANK1. GSDMD, a gene that encodes protein Gasdermin D, is a positive regulator of IL-1β involved in pyroptosis cell death [64] and was found upregulated by our meta-analysis. SIT1 is involved in the regulation of T cell activation [65], while WDR92 is involved in apoptosis via activation of caspase-3 [66]. Variants that impair apoptosis of keratinocytes may be associated with the hyperkeratosis responsible for hair follicle occlusion [67]. In the context of skin regulation, variants were detected in genes such as KLF4, FLOT2, PHACTR4, and CORO1B. KLF4 is required for the normal development of the barrier function of skin. Moreover, KLF4 is also part of the NOTCH super pathway. FLOT2 may act as a scaffolding protein and may be involved in epidermal adhesion, structure, and function, while PHACTR4 and CORO1B behave as actin-binding proteins. Thus, variants in these genes could explain dysfunction of the skin barrier and structure. Finally, MYOF plays a role in the cell membrane repair mechanism of endothelial cells that permits rapid resealing of membranes disrupted by mechanical stress [68]. A decreased number of circulating progenitor endothelial cells and endothelial dysfunction in HS has been previously described [69,70].
The variants found by the variant calling corroborated abnormal pathways found by the meta-analysis. Therefore, loss or gain of function of these genes should be closely investigated to better understand their influence in the HS pathogenesis.

4. Discussion

Here, we operated an RNA-Seq meta-analysis to explore HS lesioned skin transcription profiles among three independent studies. These studies exploited HS immunopathogenesis, yet the re-evaluation of these transcriptomes resulted in more than immunological regulation, they highlighted other important DEGs in the context of skin homeostasis and energetic metabolism. The main findings are summarized in Figure 2, which briefly represents the three main outcomes of the meta-analysis at the molecular level: immune dysregulation, skin homeostasis perturbation, and abnormal metabolic signaling. In addition, we performed a variant calling from RNA-Seq data that may potentially contribute to novel insights for HS pathogenesis.

4.1. Immune Dysregulation

The importance of immunological pathways in HS is well accepted and their dysregulation may predispose the development of a strong initial inflammatory reaction commonly seen in HS patients [9]. Here, we provide evidence that many immune-associated DEGs are present in HS lesioned skin and absent in healthy control skin. These DEGs are crucial for commonly HS-associated enriched pathways, namely neutrophil degranulation, signaling by interleukins, complement cascade, immunoregulatory interactions between a lymphoid and a non-lymphoid cell, antimicrobial peptides, and interferon signaling.
Overall, our results were in agreement with the three studies used in our analyses. Lowe et al. [37], when comparing lesioned skin to healthy control, identified IFN-γ and IL-1β as highly upregulated DEGs. These genes were found to be more than 4 times overexpressed in the meta-analysis. Our results suggested that upregulation of both proinflammatory molecules was essential to local inflammation in the skin condition. In HS lesions, IFN-γ may function as the primary activator of macrophages, which were the most numerous immune cells in the inflammatory infiltrates [71]. In normal conditions, the recruitment of macrophages was crucial not only for the immune response but also for the removal of cellular debris, the promotion of healing, and reorganization of areas within inflammation [72]. However, dysregulated activation and proliferation of macrophages may contribute to an elevated secretion of pro-inflammatory cytokines, such as IL-1β [71]. Upregulation of IL-1β exacerbates the inflammation contributing to the recruitment of other immune cells to the lesion, and to the pus formation seen in HS abscess [73]. Indeed, interferon alpha/beta/gamma signaling, as well as neutrophil degranulation pathway appeared enriched by the reactome. Witte-Händel et al. [39], who also found IL-1β highly expressed in HS lesions, postulates that its pathways-related are of paramount importance for HS phenotype. Corroborating this theory, a potential gain-of-function mutation in a positive regulator of IL-1β (GSDMD), was found by the variant calling, and upregulated expression of important molecules for the IL-1 signaling, as well as genes that stimulate or are stimulated by the expression of this cytokine were found in the mDEGs.
Gudjonsson et al. [38] obtained similar results for IFN-γ, but not for IL-1β. Instead, they found upregulation of IL36A and IL36G genes. Still, these genes have been recently reported as agonistic cytokines in the interleukin IL-1 superfamily through the activation of nuclear factor-κB (NF-κB) and mitogen-activated protein kinase (MAPK) [74]. IL-36A and IL-36G were also found upregulated in the meta-analysis. Hessam S. et al. [74] provide evidence for a distinctive IL-36 pro-inflammatory role in the development of an inflammatory loop commonly seen in HS phenotypes. Besides, IL-36 was also found to induce HS like acanthosis and hyperkeratosis in transgenic mice [75].
In our analyses, we found the upregulation of IL-17A and IL-17F. In HS, the release of follicular debris in the dermis results in the activation of an immune response mediated by Th17 [76]. IL-17 induces proinflammatory cytokine expression by keratinocytes, contributing to immune cell infiltration by neutrophils, dendritic cells, and memory T cells to the lesioned area [77]. However, Kelly et al. proposed that Th17 cells were present in HS skin prior to the formation of an active lesion, indicating that they could participate in lesion development. Indeed, genes that are influenced by IL-17, such as DEFB4A/B, MMP1, and MMP9 were also found upregulated, and had the highest-ranking product positions in the meta-analysis. Other skin genes, namely psoriasin (S100A7) and calprotectin (S100A8/9), whose expression augmented in the presence of IL-17, were also found upregulated in the meta-analysis. Since keratinocytes are believed to have an inflammatory profile in HS skin [78], these findings suggested that vicious interactions between keratinocytes and TH17/IL-17 contributed more to chronic inflammation than to enhancing a tissue protective response.
Enhanced B cell signatures were found by Lowe et al. [37] and Gudjonsson et al. [38] when they compared HS lesioned skin with the healthy control and other skin conditions, such as psoriasis. Although Witte-Händel et al. [39] did not explore B cell signaling, they mentioned finding upregulation of CXCL13, a B cell chemoattractant [79]. In our analyses, several genes associated with B cells, such as MZB1, CD19, CD79A, and CXCL13 were found to be expressed 10 times more in HS lesions when compared to controls. This may suggest that B cell alteration was a common event for HS development. B cells are important for HS pathogenesis since they contribute to fibrosis, stromal remodeling, and therapeutic responses [80]. The latter still remaining controversial; reports using rituximab, a monoclonal antibody that depletes B cells from the circulation targeting CD20 molecules, has shown to ameliorate HS phenotype [81] in some HS patients, but to induce HS-like lesions appearance in others [82]. When we analyzed the upregulated mDEGs in the PANTHER gene analysis tool, the B cell activation (P00010) pathway was enriched with 18 genes, suggesting an important role of B cells for HS lesions maintenance.

4.2. Skin Homeostasis

The main enriched pathways associated with skin regulation were ECM organization, keratinization, formation of the cornified envelope, epidermal cell differentiation, and collagen formation and degradation. Moreover, several variants associated with disruption of these pathways were identified, suggesting potential HS-related genes. Zouboulis C.C. et al. [9] demonstrated that the inflammatory process in HS was linked with abnormalities in these signaling routes. Our results confirmed these connections. For instance, elafin (PI3), an elastase-specific inhibitor that acted as an antimicrobial peptide and was expressed by epithelial and certain immune cells, appeared eight times more expressed. Besides acting against gram-positive and gram-negative bacteria, fungal pathogens, and being involved in NF-κB pathway modulation, cytokine secretion and cell recruitment [83], the gene encoding elafin has also been correlated to abnormal epithelial differentiation in a context of hyperproliferation [47,84]. Additionally, we found dysregulation in some genes of the serine protease inhibitor superfamily, serpin. Among them, SERPINB4, SERPINB3, and SERPINA1 stood out in our ranking product, corroborating the involvement of the hair follicle unit [85]. Serpins overexpression has been previously spotted at the internal and epithelial root sheath of hair follicles in HS patients [9]. In general, serpins are associated with epidermal barrier homeostasis, and chronic skin inflammation [48]. These roles have been corroborated in atopic dermatitis and psoriasis [49,86], two inflammatory skin diseases that overlap molecular pathways with HS [87,88]. Thus, upregulation of PI3 and SERPINs may be associated with keratinocyte hyperproliferation responsible for hair follicle plugin, a hallmark of HS.
Additionally, several cytokeratin appeared dysregulated, which may affect keratinization and formation of the cornified envelope. KRT6 and KRT16 are particularly interesting since they are expressed together in the outer root sheath of the hair follicles under stressful conditions [50]. In contrast, under homeostatic conditions, KRT6-KRT16 are co-expressed to respond to barrier breach, stimulating hyperproliferation of interfollicular keratinocytes and modulating the inflammatory response to wounds [89]. Dysregulation of these KRTs have been previously assigned as psoriasis biomarkers [90,91]. In parallel, dysregulation of other mDEGs, such as SPRR3, and several subtypes of SPRR2 and SPRR1 also provided evidence of impairment of barrier functions in HS. SPRRs proteins play a structural role in the cornified envelope, and SPRR3 is generally not detected in normal skin [92]. It had been suggested that an altered expression of SPRR3 could impact barrier function through altered production of cornified envelope scaffold that impaired the supramolecular organization of lamellar body-derived lipids into normal bilayer structures [93,94]. This type of dysregulation may be associated with thinning of the cornified envelope [95]. Thus, upregulation of these genes may be linked to the impaired structure and fragility of the cornified envelope and follicular infundibulum that easily break and start an immunological reaction.
Equally important to hair follicle fragility is the family of matrix metalloproteinases (MMPs), a zinc-dependent extracellular protease that break down and remodel the ECM [96]. MMPs play an important role in HS pathogenesis because it is involved in inflammation, sinus formation, and in the thinning of the basement membrane surrounding the hair follicle unit [97]. Interestingly, MMPs are also known to play a critical role in the conversion of IL-1β into its active form [98], which respond by enhancing MMP expression [99], potentially creating a long-lasting inflammatory response in HS lesioned skin. Here, we found up to 12 types of MMPs upregulated, of which MMP1 and MMP3 composed some of the highest positions in the meta-analysis overall ranking, suggesting that skin intrinsic components may represent a critical step in disease progression.
Finally, sweat glands are believed to regulate epidermal homeostasis and wound repair [100,101]. Genes that are relevant for a proper sweat-gland function namely WIF1, AQP5, FOXA1, and DCD were found highly decreased in HS skin by our meta-analysis. Respectively, the genes appeared approximately 10, 2, 6, and 24 times less expressed when compared to healthy skin. In normal conditions, sweat glands’ multipotent progenitor cells contribute to repair and skin re-epithelialization [100,101]. Thus, down-regulation of these genes may affect the capacity of skin repair, contributing to the clinical HS phenotype of non-healing wound-like environment. Similar results were found by Coates et al. (2019) [102], corroborating the possibility of impaired sweat gland function contributing to HS pathogenesis. Moreover, alongside AQP5, several other aquaporin-related genes such as AQP7, AQP6, AQP8, AQP2, AQP4, AQP7P3, AQP7P1, and AQP4-AS1 were also found down-regulated in the mDEGs. Aquaporins are membrane channel proteins that through the bidirectional transport of water, glycerol, and small solutes across the membrane, serve as critical players of the skin barrier [103]. Besides their essential role for correct skin barrier activity, AQPs, mainly AQP7, also operated in the modulation of skin’s inflammatory responses [104]. In this context, the downregulation of several aquaporins may be associated with the impaired skin barrier function evidenced by our analyses and commonly seen in inflammatory skin diseases, as well as atypical skin immune surveillance.

4.3. Energy Metabolism

Metabolic dysregulation in HS was detected through several downregulated DEGs and pathways including the metabolism of lipids, triglyceride and fatty acid, and glucocorticoid biosynthesis. Additionally, variants found by the RNA variant calling may be associated with several metabolic-related disorders, including metabolic syndrome (MetS). MetS concomitant with HS had already been described [105,106,107]. Nevertheless, it is not clear if inflammation induced by MetS leads to the initiation of HS or if the systemic inflammation in HS leads to manifestations of MetS [106]. MetS is characterized by a combination of clinical conditions including central obesity, hyperglycemia, dyslipidemia, and/or hypertension [108,109]. Although not necessarily present, central obesity plays a critical role in MetS’ development [109], and phenotype severity in HS [110,111]. In fact, at the molecular level, metabolic disorders and obesity may be considered a primary risk factor in HS [112].
Body mass index (BMI) data was available for two studies. The data confirmed that most of the patients had a BMI > 25. Obesity is supposed to favor HS skin alteration in a few manners [113]. From the molecular point of view, obesity induces low levels of systemic inflammation and metabolic changes [97]. Inflammatory immune cells including M1-type macrophages and T cells infiltrate the hypertrophic and damaged adipose tissue, producing inflammatory cytokines and inducing a dysregulated pattern of soluble mediators called adipokines [114]. Adipokines seem to drive metabolic alterations at the same time it represents a mechanistic link in the interaction between skin and metabolic comorbidities [115]. They are regularly expressed by several skin cells, such as keratinocytes, melanocytes, sebocytes, and fibroblasts, and their abnormal expression has been associated with various inflammatory skin conditions [116]. A few types of adipokines have been described and can have good and bad properties regarding the skin and body’s homeostasis [116]. The “bad” adipokines include resistin (RETN), chemerin (RARRES2), and classical pro-inflammatory cytokines, such as TNF-α, and IL-1β [115]. These metabolic villains were found upregulated in our meta-analysis. Besides their role in skin inflammation with immune cell tissue infiltration and cell dysfunction, they also drive insulin resistance, disturbance of glucose and lipid metabolism, and vascular dysfunction [117,118,119,120]. The dysregulation of these pathways was corroborated by our findings, suggesting concomitant HS and metabolic-related disorders in the patients. On the other hand, “good” adipokines such as adiponectin (ADIPOQ) that have properties, such as control of fat metabolism and insulin sensitivity, with direct anti-diabetic, anti-atherogenic, and anti-inflammatory activities [117,121,122,123], were found ten times less expressed in HS lesions. This may represent an association between the downregulation of ADIPOQ and susceptibility of HS and metabolic comorbidities. The exploration of adipokines’ role in skin and cellular metabolism is still in an early phase.

5. Limitations

Not all the studies provided clinical information of the patients, such as weight, smoking history, and race, nor general information about lifestyle, comorbidities, and different HS phenotypes. Therefore, clinical information was not taken into consideration in our analysis. Still, independent studies provided similar results, confirming the importance of the DEGs found here for HS pathogenesis. Moreover, whole exome or genome sequencing data from the same samples were not available. Thus, we could not confirm the results found by our variant calling, nor could we confirm them by PCR or sanger sequencing. However, known HS-related variants were identified substantiating our findings. Ultimately, our study design involves re-analyzing previously published transcriptomic data and aggregating evidence of transcriptional changes, thus we were not able to confirm our results on the translational level. Notwithstanding, we hope our results and conclusions stimulate new functional studies that focus on confirming our findings at the protein/tissue level. In general, our results are in compliance with what is reported in the literature.

6. Conclusions

Herein, we confirm the HS pathogenic triad composed by upregulated inflammation, altered epithelial differentiation, and dysregulated metabolism signaling. The variant calling and meta-analysis verify this highly interconnected network. The upregulation of PI3, TNF-α, IL-1β, and IFN-γ suggest that the inflammatory process in HS may be linked with keratinocyte hyperproliferation. Moreover, the upregulation of KRT6, KRT16, serpin-family genes, and SPRR3 alongside downregulation of aquaporin-family genes confirm the involvement of hair follicles and the impairment of barrier function in early phases of HS pathogenesis. Finally, our results demonstrate that HS and metabolic-related syndromes such as MetS and obesity may share transcription profiles and suggest that adipokines may be potential biomarkers for this interaction; down-regulation of ADIPOQ may be associated with diseases’ co-occurrence susceptibility. Several mutations associated with the HS pathogenic triad were found, highlighting potentially new HS-related genes associated with the sporadic form of this disease.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/biom12101371/s1, Figure S1: Percent of upregulated gene hits found in the mDEGs against total pathway hits, p-value < 0.05; Figure S2: Percent of down-regulated gene hits found in the mDEGs against total pathway hits, p-value < 0.05; Figure S3: Venn Diagram analysis representing saDEGs found in the individual studies; (a) saDEGs of the included studies; (b) upregulated saDEGs; (c) downregulated saDEGs. Figure S4: Percent of upregulated gene hits against total pathway hits; Table S1: Total mDEGs; Table S2: Reactome enriched pathways analysis associated with mDEGs (FDR < 0.05); Table S3: Variants found in patients only by the RNA-Seq variant calling; Table S4: Enriched pathways associated with saDEGs, and therefore, cDEGs, filtered by FDR < 0.05 found in the set analysis. Among the upregulated cDEGs, 91 genes had major roles in immunological pathways. For instance, 34 genes, such as CYBB (p-value = 1.758 × 10−20, logfc = 3.1303), CD3E (p value = 2.421 × 10−15, logfc = 2.6583), CD8A (p value = 3.734, logfc = 2.2508), and CD19 (p value = 2.485 × 10−30, logfc = 5.5543) were associated with adaptive immune response and their cellular signaling networks, while 44 genes, including Matrix Metallopeptidases (MMPs) family genes, DEFB4 (p value = 2.278 × 10−58, logfc = 7.6490), TLR8 (p value = 3.003 × 10−19, logfc = 2.8510), and BIN2 (p value = 2.563 × 10−13, logfc = 2.4414), were linked with innate immune system. In other important pathways, such as extracellular matrix organization, 23 genes such as ADAM12 (p value = 1.492 × 10−22, logfc = 4.1200), ADAMTS2 (p value = 2.6091, logfc = 1.111 × 10−14), and COL3A1 (p value = 1.433 × 10−17, logfc = 2.7178), were found. On the other hand, down-regulated cDEGs, such as FA2H (p value = 5.394 × 10−20, logfc = −3.7949), FADS2 (p value =1.073 × 10−24, logfc = −4.61), GUCY2EP (p value = 5.673, logfc = −1.9362), HSD3B1 (p value = 4.861 × 10−30, logfc = −5.2634), IYD (p value = 1.337 × 10−20, logfc = −2.4972), MOGAT2 (p value = 1.59 × 10−23, logfc = −3.7538), MPPED1 (p value = 1.519 × 10−25, logfc = −3.1882), PNPLA5 (p value = 3.314 × 10−29, logfc = −4.2126), SOX9-AS1 (p value = 8.702 × 10−17, logfc = −2.3383), and THRSP (p value = 4.101 × 10−48, logfc = −5.4498) play a crucial role in metabolic pathways. In particular, FA2H, FADS2, MOGAT2, HSD3B1, THRSP, and PNPLA5 are important for the metabolism of lipids while IYD participates in thyroxine biosynthesis and metabolism of amine-derived hormones.

Author Contributions

Conceptualization: A.S.L.E.d.O., A.V.M., P.M.T., A.A., S.J., C.M., L.A.C.B. and S.C.; methodology: G.B., R.M. and A.C.; software: R.M. and L.A.C.B.; validation: R.M. and L.A.C.B.; formal analysis: R.M., G.B. and A.C.; investigation, A.S.L.E.d.O.; resources, L.A.C.B. and R.M.; data curation: G.B. and A.C.; writing—original draft preparation: A.S.L.E.d.O.; writing—review and editing: A.V.M., P.M.T., A.A., S.J., L.A.C.B. and C.M.; visualization: A.S.L.E.d.O. and L.A.C.B.; supervision: S.C. and L.A.C.B.; project administration: S.C. and L.A.C.B.; funding acquisition: S.C., L.A.C.B. and S.J. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by “Coordenação de Aperfeiçoamento de Pessoal de Nível Superior—Brasil (CAPES)”—Finance Code 001, Fondation René Touraine, “Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)” (311415/2020-2 and 430353/2018-9), and EraPerMed 2018-17 European Community funds. L.A.C.B. is supported by CNPq (311415/2020-2). This work was also supported by the Italian Ministry of Health, through the contribution given to the Institute for Maternal and Child Health IRCCS Burlo Garofolo, Trieste, Italy for the Starting Grant (SG-2019-12369421) and for RC03/2020.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Transcriptomic data sets can be accessed at the Sequence Read Archive (SRA) database and Gene Expression Omnibus (GEO) database (SRP273719/GSE155176, SRP272659/GSE154773, and SRP168930/GSE122592), both from the National Center for Biotechnology Information.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Katoulis, A.C.; Liakou, A.I.; Rotsiamis, N.; Bonovas, S.; Bozi, E.; Rallis, E.; Christodoulou, C.; Rigopoulos, D. Descriptive Epidemiology of Hidradenitis Suppurativa in Greece: A Study of 152 Cases. Ski. Appendage Disord. 2017, 3, 197–201. [Google Scholar] [CrossRef] [PubMed]
  2. Ralf Paus, L.; Kurzen, H.; Kurokawa, I.; Jemec, G.B.E.; Emtestam, L.; Sellheyer, K.; Giamarellos-Bourboulis, E.J.; Nagy, I.; Bechara, F.G.; Sartorius, K.; et al. What Causes Hidradenitis Suppurativa? Exp. Dermatol. 2008, 17, 455–456. [Google Scholar] [CrossRef]
  3. Dunstan, R.W.; Salte, K.M.; Todorović, V.; Lowe, M.; Wetter, J.B.; Harms, P.W.; Burney, R.E.; Scott, V.E.; Smith, K.M.; Rosenblum, M.D.; et al. Histologic Progression of Acne Inversa/Hidradenitis Suppurativa: Implications for Future Investigations and Therapeutic Intervention. Exp. Dermatol. 2021, 30, 820–830. [Google Scholar] [CrossRef]
  4. Jemec, G.B.E. Hidradenitis Suppurativa. N. Engl. J. Med. 2012, 366, 158–164. [Google Scholar] [CrossRef] [PubMed]
  5. Matusiak, Ł. Profound Consequences of Hidradenitis Suppurativa: A Review. Br. J. Dermatol. 2020, 183, e171–e177. [Google Scholar] [CrossRef] [PubMed]
  6. Dufour, D.N.; Emtestam, L.; Jemec, G.B. Hidradenitis Suppurativa: A Common and Burdensome, yet under-Recognised, Inflammatory Skin Disease. Postgrad. Med. J. 2014, 90, 216–221. [Google Scholar] [CrossRef] [PubMed]
  7. von Laffert, M.; Helmbold, P.; Wohlrab, J.; Fiedler, E.; Stadie, V.; Marsch, W.C. Hidradenitis Suppurativa (Acne Inversa): Early Inflammatory Events at Terminal Follicles and at Interfollicular Epidermis*. Exp. Dermatol. 2009, 19, 533–537. [Google Scholar] [CrossRef]
  8. von Laffert, M.; Stadie, V.; Wohlrab, J.; Marsch, W.C. Hidradenitis Suppurativa/Acne Inversa: Bilocated Epithelial Hyperplasia with Very Different Sequelae. Br. J. Dermatol. 2011, 164, 367–371. [Google Scholar] [CrossRef]
  9. Zouboulis, C.C.; Nogueira da Costa, A.; Makrantonaki, E.; Hou, X.X.; Almansouri, D.; Dudley, J.T.; Edwards, H.; Readhead, B.; Balthasar, O.; Jemec, G.B.E.; et al. Alterations in Innate Immunity and Epithelial Cell Differentiation Are the Molecular Pillars of Hidradenitis Suppurativa. J. Eur. Acad. Dermatol. Venereol. 2020, 34, 846–861. [Google Scholar] [CrossRef]
  10. Duchatelet, S.; Miskinyte, S.; Delage, M.; Ungeheuer, M.-N.; Lam, T.; Benhadou, F.; del Marmol, V.; Vossen, A.R.J.V.; Prens, E.P.; Cogrel, O.; et al. Low Prevalence of GSC Gene Mutations in a Large Cohort of Predominantly Caucasian Patients with Hidradenitis Suppurativa. J. Investig. Dermatol. 2020, 140, 2085–2088.e14. [Google Scholar] [CrossRef]
  11. Pink, A.E.; Simpson, M.A.; Desai, N.; Dafou, D.; Hills, A.; Mortimer, P.; Smith, C.H.; Trembath, R.C.; Barker, J.N.W. Mutations in the γ-Secretase Genes NCSTN, PSENEN, and PSEN1 Underlie Rare Forms of Hidradenitis Suppurativa (Acne Inversa). J. Investig. Dermatol. 2012, 132, 2459–2461. [Google Scholar] [CrossRef] [PubMed]
  12. Marzano, A.V.; Genovese, G.; Moltrasio, C.; Tricarico, P.M.; Gratton, R.; Piaserico, S.; Garcovich, S.; Boniotto, M.; Brandão, L.; Moura, R.; et al. Whole-Exome Sequencing in 10 Unrelated Patients with Syndromic Hidradenitis Suppurativa: A Preliminary Step for a Genotype-Phenotype Correlation. Dermatology 2022, 238, 860–869. [Google Scholar] [CrossRef] [PubMed]
  13. Zouboulis, C.C.; Benhadou, F.; Byrd, A.S.; Chandran, N.S.; Giamarellos-Bourboulis, E.J.; Fabbrocini, G.; Frew, J.W.; Fujita, H.; González-López, M.A.; Guillem, P.; et al. What Causes Hidradenitis Suppurativa?—15 Years After. Exp. Dermatol. 2020, 29, 1154–1170. [Google Scholar] [CrossRef] [PubMed]
  14. Pan, Y.; Lin, M.-H.; Tian, X.; Cheng, H.-T.; Gridley, T.; Shen, J.; Kopan, R. γ-Secretase Functions through Notch Signaling to Maintain Skin Appendages but Is Not Required for Their Patterning or Initial Morphogenesis. Dev. Cell 2004, 7, 731–743. [Google Scholar] [CrossRef] [PubMed]
  15. Melnik, B.C.; Plewig, G. Impaired Notch-MKP-1 Signalling in Hidradenitis Suppurativa: An Approach to Pathogenesis by Evidence from Translational Biology. Exp. Dermatol. 2013, 22, 172–177. [Google Scholar] [CrossRef]
  16. Puente, X.S.; Pinyol, M.; Quesada, V.; Conde, L.; Ordóñez, G.R.; Villamor, N.; Escaramis, G.; Jares, P.; Beà, S.; González-Díaz, M.; et al. Whole-Genome Sequencing Identifies Recurrent Mutations in Chronic Lymphocytic Leukaemia. Nature 2011, 475, 101–105. [Google Scholar] [CrossRef]
  17. Ramasamy, A.; Mondry, A.; Holmes, C.C.; Altman, D.G. Key Issues in Conducting a Meta-Analysis of Gene Expression Microarray Datasets. PLoS Med. 2008, 5, e184. [Google Scholar] [CrossRef]
  18. Leinonen, R.; Sugawara, H.; Shumway, M. The Sequence Read Archive. Nucleic Acids Res. 2011, 39, D19–D21. [Google Scholar] [CrossRef]
  19. Clough, E.; Barrett, T. The Gene Expression Omnibus Database. Methods Mol. Biol. 2016, 1418, 93–110. [Google Scholar] [CrossRef]
  20. Brandão, L.A.C.; Tricarico, P.M.; Gratton, R.; Agrelli, A.; Zupin, L.; Abou-Saleh, H.; Moura, R.; Crovella, S. Multiomics Integration in Skin Diseases with Alterations in Notch Signaling Pathway: PlatOMICs Phase 1 Deployment. Int. J. Mol. Sci. 2021, 22, 1523. [Google Scholar] [CrossRef]
  21. Lee, S.Y.; Park, Y.K.; Yoon, C.-H.; Kim, K.; Kim, K.-C. Meta-Analysis of Gene Expression Profiles in Long-Term Non-Progressors Infected with HIV-1. BMC Med. Genom. 2019, 12, 3. [Google Scholar] [CrossRef] [PubMed]
  22. Coelho, A.V.C.; Gratton, R.; de Melo, J.P.B.; Andrade-Santos, J.L.; Guimarães, R.L.; Crovella, S.; Tricarico, P.M.; Brandão, L.A.C. HIV-1 Infection Transcriptomics: Meta-Analysis of CD4+ T Cells Gene Expression Profiles. Viruses 2021, 13, 244. [Google Scholar] [CrossRef]
  23. Zhu, Y.; Stephens, R.M.; Meltzer, P.S.; Davis, S.R. SRAdb: Query and Use Public next-Generation Sequencing Data from within R. BMC Bioinform. 2013, 14, 19. [Google Scholar] [CrossRef]
  24. The R Project for Statistical Computing. Available online: https://www.r-project.org/ (accessed on 9 June 2022).
  25. Bolger, A.M.; Lohse, M.; Usadel, B. Trimmomatic: A Flexible Trimmer for Illumina Sequence Data. Bioinformatics 2014, 30, 2114–2120. [Google Scholar] [CrossRef] [PubMed]
  26. Dobin, A.; Davis, C.A.; Schlesinger, F.; Drenkow, J.; Zaleski, C.; Jha, S.; Batut, P.; Chaisson, M.; Gingeras, T.R. STAR: Ultrafast Universal RNA-Seq Aligner. Bioinformatics 2013, 29, 15–21. [Google Scholar] [CrossRef]
  27. Liao, Y.; Smyth, G.K.; Shi, W. The R Package Rsubread Is Easier, Faster, Cheaper and Better for Alignment and Quantification of RNA Sequencing Reads. Nucleic Acids Res. 2019, 47, e47. [Google Scholar] [CrossRef]
  28. Love, M.I.; Huber, W.; Anders, S. Moderated Estimation of Fold Change and Dispersion for RNA-Seq Data with DESeq2. Genome Biol. 2014, 15, 550. [Google Scholar] [CrossRef] [PubMed]
  29. Mi, H.; Muruganujan, A.; Huang, X.; Ebert, D.; Mills, C.; Guo, X.; Thomas, P.D. Protocol Update for Large-Scale Genome and Gene Function Analysis with the PANTHER Classification System (v.14.0). Nat. Protoc. 2019, 14, 703–721. [Google Scholar] [CrossRef]
  30. Panther Classification System. Available online: http://www.pantherdb.org/ (accessed on 19 July 2022).
  31. McKenna, A.; Hanna, M.; Banks, E.; Sivachenko, A.; Cibulskis, K.; Kernytsky, A.; Garimella, K.; Altshuler, D.; Gabriel, S.; Daly, M.; et al. The Genome Analysis Toolkit: A MapReduce Framework for Analyzing next-Generation DNA Sequencing Data. Genome Res. 2010, 20, 1297–1303. [Google Scholar] [CrossRef]
  32. Kim, S.; Scheffler, K.; Halpern, A.L.; Bekritsky, M.A.; Noh, E.; Källberg, M.; Chen, X.; Kim, Y.; Beyter, D.; Krusche, P.; et al. Strelka2: Fast and Accurate Calling of Germline and Somatic Variants. Nat. Methods 2018, 15, 591–594. [Google Scholar] [CrossRef]
  33. Wang, K.; Li, M.; Hakonarson, H. ANNOVAR: Functional Annotation of Genetic Variants from High-Throughput Sequencing Data. Nucleic Acids Res. 2010, 38, e164. [Google Scholar] [CrossRef] [PubMed]
  34. Garg, A.; Wertenteil, S.; Baltz, R.; Strunk, A.; Finelt, N. Prevalence Estimates for Hidradenitis Suppurativa among Children and Adolescents in the United States: A Gender- and Age-Adjusted Population Analysis. J. Investig. Dermatol. 2018, 138, 2152–2156. [Google Scholar] [CrossRef] [PubMed]
  35. Jemec, G.B.E.; Heidenheim, M.; Nielsen, N.H. The Prevalence of Hidradenitis Suppurativa and Its Potential Precursor Lesions. J. Am. Acad. Dermatol. 1996, 35, 191–194. [Google Scholar] [CrossRef]
  36. Kirby, J.S.; Zaenglein, A.L. Recognizing the Effects and Disparities of Pediatric Hidradenitis Suppurativa. JAMA Dermatol. 2021, 157, 379. [Google Scholar] [CrossRef]
  37. Lowe, M.M.; Naik, H.B.; Clancy, S.; Pauli, M.; Smith, K.M.; Bi, Y.; Dunstan, R.; Gudjonsson, J.E.; Paul, M.; Harris, H.; et al. Immunopathogenesis of Hidradenitis Suppurativa and Response to Anti–TNF-α Therapy. JCI Insight 2020, 5, e139932. [Google Scholar] [CrossRef]
  38. Gudjonsson, J.E.; Tsoi, L.C.; Ma, F.; Billi, A.C.; van Straalen, K.R.; Vossen, A.R.J.V.; van der Zee, H.H.; Harms, P.W.; Wasikowski, R.; Yee, C.M.; et al. Contribution of Plasma Cells and B Cells to Hidradenitis Suppurativa Pathogenesis. JCI Insight 2020, 5, e139930. [Google Scholar] [CrossRef]
  39. Witte-Händel, E.; Wolk, K.; Tsaousi, A.; Irmer, M.L.; Mößner, R.; Shomroni, O.; Lingner, T.; Witte, K.; Kunkel, D.; Salinas, G.; et al. The IL-1 Pathway Is Hyperactive in Hidradenitis Suppurativa and Contributes to Skin Infiltration and Destruction. J. Investig. Dermatol. 2019, 139, 1294–1305. [Google Scholar] [CrossRef]
  40. Kolbinger, F.; Loesche, C.; Valentin, M.-A.; Jiang, X.; Cheng, Y.; Jarvis, P.; Peters, T.; Calonder, C.; Bruin, G.; Polus, F.; et al. β-Defensin 2 Is a Responsive Biomarker of IL-17A–Driven Skin Pathology in Patients with Psoriasis. J. Allergy Clin. Immunol. 2017, 139, 923–932. [Google Scholar] [CrossRef]
  41. Clausen, M.-L.; Jungersted, J.M.; Andersen, P.S.; Slotved, H.-C.; Krogfelt, K.A.; Agner, T. Human β-Defensin-2 as a Marker for Disease Severity and Skin Barrier Properties in Atopic Dermatitis. Br. J. Dermatol. 2013, 169, 587–593. [Google Scholar] [CrossRef]
  42. Cao, Y.; Harvey, B.P.; Hong, F.; Ruzek, M.; Wang, J.; Murphy, E.R.; Kaymakcalan, Z. Adalimumab Induces a Wound Healing Profile in Patients with Hidradenitis Suppurativa by Regulating Macrophage Differentiation and Matrix Metalloproteinase Expression. J. Investig. Dermatol. 2021, 141, 2730–2740.e9. [Google Scholar] [CrossRef]
  43. Alowami, S.; Qing, G.; Emberley, E.; Snell, L.; Watson, P.H. Psoriasin (S100A7) Expression Is Altered during Skin Tumorigenesis. BMC Dermatol. 2003, 3, 1. [Google Scholar] [CrossRef]
  44. Mirastschijski, U.; Lupše, B.; Maedler, K.; Sarma, B.; Radtke, A.; Belge, G.; Dorsch, M.; Wedekind, D.; McCawley, L.J.; Boehm, G.; et al. Matrix Metalloproteinase-3 Is Key Effector of TNF-α-Induced Collagen Degradation in Skin. Int. J. Mol. Sci. 2019, 20, 5234. [Google Scholar] [CrossRef] [PubMed]
  45. Pine, G.M.; Batugedara, H.M.; Nair, M.G. Here, There and Everywhere: Resistin-like Molecules in Infection, Inflammation, and Metabolic Disorders. Cytokine 2018, 110, 442–451. [Google Scholar] [CrossRef] [PubMed]
  46. Er, L.-K.; Wu, S.; Hsu, L.-A.; Teng, M.-S.; Sun, Y.-C.; Ko, Y.-L. Pleiotropic Associations of RARRES2 Gene Variants and Circulating Chemerin Levels: Potential Roles of Chemerin Involved in the Metabolic and Inflammation-Related Diseases. Mediat. Inflamm. 2018, 2018, 4670521. [Google Scholar] [CrossRef] [PubMed]
  47. Pol, A.; Pfundt, R.; Zeeuwen, P.; Molhuizen, H.; Schalkwijk, J. Transcriptional Regulation of the Elafin Gene in Human Keratinocytes. J. Investig. Dermatol. 2003, 120, 301–307. [Google Scholar] [CrossRef]
  48. Sivaprasad, U.; Kinker, K.G.; Ericksen, M.B.; Lindsey, M.; Gibson, A.M.; Bass, S.A.; Hershey, N.S.; Deng, J.; Medvedovic, M.; Khurana Hershey, G.K. SERPINB3/B4 Contributes to Early Inflammation and Barrier Dysfunction in an Experimental Murine Model of Atopic Dermatitis. J. Investig. Dermatol. 2015, 135, 160–169. [Google Scholar] [CrossRef] [PubMed]
  49. Kantaputra, P.; Chaowattanapanit, S.; Kiratikanon, S.; Chaiwarith, R.; Choonhakarn, C.; Intachai, W.; Quarto, N.; Tongsima, S.; Ketudat Cairns, J.R.; Ngamphiw, C.; et al. SERPINA1, Generalized Pustular Psoriasis, and Adult-onset Immunodeficiency. J. Dermatol. 2021, 48, 1597–1601. [Google Scholar] [CrossRef]
  50. Lessard, J.C.; Piña-Paz, S.; Rotty, J.D.; Hickerson, R.P.; Kaspar, R.L.; Balmain, A.; Coulombe, P.A. Keratin 16 Regulates Innate Immunity in Response to Epidermal Barrier Breach. Proc. Natl. Acad. Sci. USA 2013, 110, 19537–19542. [Google Scholar] [CrossRef]
  51. Kainu, K.; Kivinen, K.; Zucchelli, M.; Suomela, S.; Kere, J.; Inerot, A.; Baker, B.S.; Powles, A.V.; Fry, L.; Samuelsson, L.; et al. Association of Psoriasis to PGLYRP and SPRR Genes at PSORS4 Locus on 1q Shows Heterogeneity between Finnish, Swedish and Irish Families. Exp. Dermatol. 2009, 18, 109–115. [Google Scholar] [CrossRef]
  52. Keil, K.P.; Mehta, V.; Branam, A.M.; Abler, L.L.; Buresh-Stiemke, R.A.; Joshi, P.S.; Schmitz, C.T.; Marker, P.C.; Vezina, C.M. Wnt Inhibitory Factor 1 (Wif1) Is Regulated by Androgens and Enhances Androgen-Dependent Prostate Development. Endocrinology 2012, 153, 6091–6103. [Google Scholar] [CrossRef] [Green Version]
  53. Yen, C.-L.E.; Brown, C.H.; Monetti, M.; Farese, R.V. A Human Skin Multifunctional O-Acyltransferase That Catalyzes the Synthesis of Acylglycerols, Waxes, and Retinyl Esters. J. Lipid Res. 2005, 46, 2388–2397. [Google Scholar] [CrossRef] [PubMed]
  54. Kim, K.-A.; Ka, S.-O.; Moon, W.S.; Yi, H.-K.; Lee, Y.-H.; Kwon, K.-B.; Park, J.-W.; Park, B.-H. Effect of Dermcidin, an Antimicrobial Peptide, on Body Fat Mobilization in Normal Mice. J. Endocrinol. 2008, 198, 111–118. [Google Scholar] [CrossRef] [PubMed]
  55. Ahonen, M.A.; Höring, M.; Nguyen, V.D.; Qadri, S.; Taskinen, J.H.; Nagaraj, M.; Wabitsch, M.; Fischer-Posovszky, P.; Zhou, Y.; Liebisch, G.; et al. Insulin-Inducible THRSP Maintains Mitochondrial Function and Regulates Sphingolipid Metabolism in Human Adipocytes. Mol. Med. 2022, 28, 68. [Google Scholar] [CrossRef] [PubMed]
  56. MacKenzie, P.I.; Rogers, A.; Elliot, D.J.; Chau, N.; Hulin, J.-A.; Miners, J.O.; Meech, R. The Novel UDP Glycosyltransferase 3A2: Cloning, Catalytic Properties, and Tissue Distribution. Mol. Pharmacol. 2011, 79, 472–478. [Google Scholar] [CrossRef] [PubMed]
  57. Massi, D.; Panelos, J. Notch Signaling and the Developing Skin Epidermis. In Notch Signaling in Embryology and Cancer; Springer: Berlin/Heidelberg, Germany, 2012; pp. 131–141. [Google Scholar]
  58. Radtke, F.; Fasnacht, N.; MacDonald, H.R. Notch Signaling in the Immune System. Immunity 2010, 32, 14–27. [Google Scholar] [CrossRef] [PubMed]
  59. Pavlovsky, M.; Sarig, O.; Eskin-Schwartz, M.; Malchin, N.; Bochner, R.; Mohamad, J.; Gat, A.; Peled, A.; Hafner, A.; Sprecher, E. A Phenotype Combining Hidradenitis Suppurativa with Dowling-Degos Disease Caused by a Founder Mutation in PSENEN. Br. J. Dermatol. 2018, 178, 502–508. [Google Scholar] [CrossRef]
  60. Jfri, A.; Litvinov, I.V.; Netchiporouk, E.; O’Brien, E. Novel Variants of MEFV and NOD2 Genes in Familial Hidradenitis Suppurativa: A Case Report. SAGE Open Med. Case Rep. 2020, 8, 2050313X2095311. [Google Scholar] [CrossRef]
  61. Higgins, R.; Pink, A.; Hunger, R.; Yawalkar, N.; Navarini, A.A. Generalized Comedones, Acne, and Hidradenitis Suppurativa in a Patient with an FGFR2 Missense Mutation. Front. Med. 2017, 4, 16. [Google Scholar] [CrossRef]
  62. Saito, N.; Minami-Hori, M.; Nagahata, H.; Nozaki, H.; Iinuma, S.; Igawa, S.; Kanno, K.; Kishibe, M.; Kanazawa, N.; Ishida-Yamamoto, A. Novel PSTPIP1 Gene Mutation in Pyoderma Gangrenosum, Acne and Suppurative Hidradenitis Syndrome. J. Dermatol. 2018, 45, e213–e214. [Google Scholar] [CrossRef]
  63. Witkowski, A.; Thweatt, J.; Smith, S. Mammalian ACSF3 Protein Is a Malonyl-CoA Synthetase That Supplies the Chain Extender Units for Mitochondrial Fatty Acid Synthesis. J. Biol. Chem. 2011, 286, 33729–33736. [Google Scholar] [CrossRef] [Green Version]
  64. Fenini, G.; Contassot, E.; French, L.E. Potential of IL-1, IL-18 and Inflammasome Inhibition for the Treatment of Inflammatory Skin Diseases. Front. Pharmacol. 2017, 8, 278. [Google Scholar] [CrossRef] [PubMed]
  65. Arndt, B.; Krieger, T.; Kalinski, T.; Thielitz, A.; Reinhold, D.; Roessner, A.; Schraven, B.; Simeoni, L. The Transmembrane Adaptor Protein SIT Inhibits TCR-Mediated Signaling. PLoS ONE 2011, 6, e23761. [Google Scholar] [CrossRef]
  66. Saeki, M.; Irie, Y.; Ni, L.; Yoshida, M.; Itsuki, Y.; Kamisaki, Y. Monad, a WD40 Repeat Protein, Promotes Apoptosis Induced by TNF-α. Biochem. Biophys. Res. Commun. 2006, 342, 568–572. [Google Scholar] [CrossRef] [PubMed]
  67. Raj, D.; Brash, D.E.; Grossman, D. Keratinocyte Apoptosis in Epidermal Development and Disease. J. Investig. Dermatol. 2006, 126, 243–257. [Google Scholar] [CrossRef]
  68. Dong, Y.; Kang, H.; Liu, H.; Wang, J.; Guo, Q.; Song, C.; Sun, Y.; Zhang, Y.; Zhang, H.; Zhang, Z.; et al. Myoferlin, a Membrane Protein with Emerging Oncogenic Roles. BioMed Res. Int. 2019, 2019, 7365913. [Google Scholar] [CrossRef] [PubMed]
  69. González-López, M.A.; Ocejo-Viñals, J.G.; López-Sundh, A.E.; Guiral, S.; Ruiz-Solana, M.; Mata, C.; Portilla, V.; Corrales, A.; Blanco, R.; Hernández, J.L. Biomarkers of Endothelial Dysfunction and Atherosclerosis in Hidradenitis Suppurativa. J. Dermatol. 2022. [Google Scholar] [CrossRef] [PubMed]
  70. Matusiak, Ł.; Batycka-Baran, A.; Bieniek, A.; Szepietowski, J.C. Decreased Number of Circulating Endothelial Progenitor Cells in Hidradenitis Suppurativa Patients. Dermatology 2015, 230, 228–233. [Google Scholar] [CrossRef]
  71. Shah, A.; Alhusayen, R.; Amini-Nik, S. The Critical Role of Macrophages in the Pathogenesis of Hidradenitis Suppurativa. Inflamm. Res. 2017, 66, 931–945. [Google Scholar] [CrossRef]
  72. Krzyszczyk, P.; Schloss, R.; Palmer, A.; Berthiaume, F. The Role of Macrophages in Acute and Chronic Wound Healing and Interventions to Promote Pro-Wound Healing Phenotypes. Front. Physiol. 2018, 9, 419. [Google Scholar] [CrossRef]
  73. Banerjee, A.; McNish, S.; Shanmugam, V.K. Interferon-Gamma (IFN-γ) Is Elevated in Wound Exudate from Hidradenitis Suppurativa. Immunol. Investig. 2017, 46, 149–158. [Google Scholar] [CrossRef] [Green Version]
  74. Hessam, S.; Sand, M.; Gambichler, T.; Skrygan, M.; Rüddel, I.; Bechara, F.G. Interleukin-36 in Hidradenitis Suppurativa: Evidence for a Distinctive Proinflammatory Role and a Key Factor in the Development of an Inflammatory Loop. Br. J. Dermatol. 2018, 178, 761–767. [Google Scholar] [CrossRef] [PubMed]
  75. Tortola, L.; Rosenwald, E.; Abel, B.; Blumberg, H.; Schäfer, M.; Coyle, A.J.; Renauld, J.-C.; Werner, S.; Kisielow, J.; Kopf, M. Psoriasiform Dermatitis Is Driven by IL-36–Mediated DC-Keratinocyte Crosstalk. J. Clin. Investig. 2012, 122, 3965–3976. [Google Scholar] [CrossRef]
  76. Vossen, A.R.J.V.; van der Zee, H.H.; Prens, E.P. Hidradenitis Suppurativa: A Systematic Review Integrating Inflammatory Pathways Into a Cohesive Pathogenic Model. Front. Immunol. 2018, 9, 2965. [Google Scholar] [CrossRef] [PubMed]
  77. Keijsers, R.R.M.C.; Joosten, I.; van Erp, P.E.J.; Koenen, H.J.P.M.; van de Kerkhof, P.C.M. Cellular Sources of IL-17 in Psoriasis: A Paradigm Shift? Exp. Dermatol. 2014, 23, 799–803. [Google Scholar] [CrossRef] [PubMed]
  78. Hotz, C.; Boniotto, M.; Guguin, A.; Surenaud, M.; Jean-Louis, F.; Tisserand, P.; Ortonne, N.; Hersant, B.; Bosc, R.; Poli, F.; et al. Intrinsic Defect in Keratinocyte Function Leads to Inflammation in Hidradenitis Suppurativa. J. Investig. Dermatol. 2016, 136, 1768–1780. [Google Scholar] [CrossRef] [PubMed]
  79. Kowarik, M.C.; Cepok, S.; Sellner, J.; Grummel, V.; Weber, M.S.; Korn, T.; Berthele, A.; Hemmer, B. CXCL13 Is the Major Determinant for B Cell Recruitment to the CSF during Neuroinflammation. J. Neuroinflammation 2012, 9, 93. [Google Scholar] [CrossRef]
  80. Frew, J.W.; Grand, D.; Navrazhina, K.; Krueger, J.G. Beyond Antibodies: B Cells in Hidradenitis Suppurativa: Bystanders, Contributors or Therapeutic Targets? Exp. Dermatol. 2020, 29, 509–515. [Google Scholar] [CrossRef]
  81. Takahashi, K.; Yanagi, T.; Kitamura, S.; Hata, H.; Imafuku, K.; Iwami, D.; Hotta, K.; Morita, K.; Shinohara, N.; Shimizu, H. Successful Treatment of Hidradenitis Suppurativa with Rituximab for a Patient with Idiopathic Carpotarsal Osteolysis and Chronic Active Antibody-Mediated Rejection. J. Dermatol. 2018, 45, e116–e117. [Google Scholar] [CrossRef]
  82. Faivre, C.; Villani, A.P.; Aubin, F.; Lipsker, D.; Bottaro, M.; Cohen, J.-D.; Durupt, F.; Jeudy, G.; Sbidian, E.; Toussirot, E.; et al. Hidradenitis Suppurativa (HS): An Unrecognized Paradoxical Effect of Biologic Agents (BA) Used in Chronic Inflammatory Diseases. J. Am. Acad. Dermatol. 2016, 74, 1153–1159. [Google Scholar] [CrossRef]
  83. Verrier, T.; Solhonne, B.; Sallenave, J.-M.; Garcia-Verdugo, I. The WAP Protein Trappin-2/Elafin: A Handyman in the Regulation of Inflammatory and Immune Responses. Int. J. Biochem. Cell Biol. 2012, 44, 1377–1380. [Google Scholar] [CrossRef]
  84. Wiedow, O.; Schröder, J.M.; Gregory, H.; Young, J.A.; Christophers, E. Elafin: An Elastase-Specific Inhibitor of Human Skin. Purification, Characterization, and Complete Amino Acid Sequence. J. Biol. Chem. 1990, 265, 14791–14795. [Google Scholar] [CrossRef]
  85. Jensen, P.J.; Yang, T.; Yu, D.-W.; Baker, M.S.; Risse, B.; Sun, T.-T.; Lavker, R.M. Serpins in the Human Hair Follicle. J. Investig. Dermatol. 2000, 114, 917–922. [Google Scholar] [CrossRef] [PubMed]
  86. Tian, S.; Krueger, J.G.; Li, K.; Jabbari, A.; Brodmerkel, C.; Lowes, M.A.; Suárez-Fariñas, M. Meta-Analysis Derived (MAD) Transcriptome of Psoriasis Defines the “Core” Pathogenesis of Disease. PLoS ONE 2012, 7, e44274. [Google Scholar] [CrossRef]
  87. Kaakati, R.N.; Tanaka, J.; Liu, B.; Ward, R.; Macleod, A.S.; Green, C.L.; Jaleel, T. Atopic Dermatitis Is Associated with Hidradenitis Suppurativa Diagnosis: A Single Institution Retrospective Cohort Study. JAAD Int. 2021, 4, 18–24. [Google Scholar] [CrossRef]
  88. Fletcher, J.M.; Moran, B.; Petrasca, A.; Smith, C.M. IL-17 in Inflammatory Skin Diseases Psoriasis and Hidradenitis Suppurativa. Clin. Exp. Immunol. 2020, 201, 121–134. [Google Scholar] [CrossRef] [PubMed]
  89. Zhang, X.; Yin, M.; Zhang, L. Keratin 6, 16 and 17—Critical Barrier Alarmin Molecules in Skin Wounds and Psoriasis. Cells 2019, 8, 807. [Google Scholar] [CrossRef] [PubMed]
  90. Jiang, S.; Hinchliffe, T.E.; Wu, T. Biomarkers of An Autoimmune Skin Disease—Psoriasis. Genom. Proteom. Bioinform. 2015, 13, 224–233. [Google Scholar] [CrossRef]
  91. Mommers, J.M.; van Rossum, M.M.; van Erp, P.E.J.; van de Kerkhof, P.C.M. Changes in Keratin 6 and Keratin 10 (Co-)Expression in Lesional and Symptomless Skin of Spreading Psoriasis. Dermatology 2000, 201, 15–20. [Google Scholar] [CrossRef]
  92. Tesfaigzi, J.; Carlson, D.M. Expression, Regulation, and Function of the SPR Family of Proteins. Cell Biochem. Biophys. 1999, 30, 243–265.e1. [Google Scholar] [CrossRef]
  93. Elias, P.M.; Wakefield, J.S. Mechanisms of Abnormal Lamellar Body Secretion and the Dysfunctional Skin Barrier in Patients with Atopic Dermatitis. J. Allergy Clin. Immunol. 2014, 134, 781–791. [Google Scholar] [CrossRef] [Green Version]
  94. Carregaro, F.; Stefanini, A.C.B.; Henrique, T.; Tajara, E.H. Study of Small Proline-Rich Proteins (SPRRs) in Health and Disease: A Review of the Literature. Arch. Dermatol. Res. 2013, 305, 857–866. [Google Scholar] [CrossRef] [PubMed]
  95. Guttman-Yassky, E.; Suárez-Fariñas, M.; Chiricozzi, A.; Nograles, K.E.; Shemer, A.; Fuentes-Duculan, J.; Cardinale, I.; Lin, P.; Bergman, R.; Bowcock, A.M.; et al. Broad Defects in Epidermal Cornification in Atopic Dermatitis Identified through Genomic Analysis. J. Allergy Clin. Immunol. 2009, 124, 1235–1244.e58. [Google Scholar] [CrossRef] [PubMed]
  96. Lukes, A.; Mun-Bryce, S.; Lukes, M.; Rosenberg, G.A. Extracellular Matrix Degradation by Metalloproteinases and Central Nervous System Diseases. Mol. Neurobiol. 1999, 19, 267–284. [Google Scholar] [CrossRef] [PubMed]
  97. Wolk, K.; Join-Lambert, O.; Sabat, R. Aetiology and Pathogenesis of Hidradenitis Suppurativa. Br. J. Dermatol. 2020, 183, 999–1010. [Google Scholar] [CrossRef]
  98. Opdenakker, G.; van den Steen, P.E.; Dubois, B.; Nelissen, I.; van Coillie, E.; Masure, S.; Proost, P.; van Damme, J. Gelatinase B Functions as Regulator and Effector in Leukocyte Biology. J. Leukoc. Biol. 2001, 69, 851–859. [Google Scholar] [PubMed]
  99. Sanchez, J.; le Jan, S.; Muller, C.; François, C.; Renard, Y.; Durlach, A.; Bernard, P.; Reguiai, Z.; Antonicelli, F. Matrix Remodelling and MMP Expression/Activation Are Associated with Hidradenitis Suppurativa Skin Inflammation. Exp. Dermatol. 2019, 28, 593–600. [Google Scholar] [CrossRef]
  100. Rittié, L.; Sachs, D.L.; Orringer, J.S.; Voorhees, J.J.; Fisher, G.J. Eccrine Sweat Glands Are Major Contributors to Reepithelialization of Human Wounds. Am. J. Pathol. 2013, 182, 163–171. [Google Scholar] [CrossRef]
  101. Lu, C.P.; Polak, L.; Rocha, A.S.; Pasolli, H.A.; Chen, S.-C.; Sharma, N.; Blanpain, C.; Fuchs, E. Identification of Stem Cell Populations in Sweat Glands and Ducts Reveals Roles in Homeostasis and Wound Repair. Cell 2012, 150, 136–150. [Google Scholar] [CrossRef]
  102. Coates, M.; Mariottoni, P.; Corcoran, D.L.; Kirshner, H.F.; Jaleel, T.; Brown, D.A.; Brooks, S.R.; Murray, J.; Morasso, M.I.; MacLeod, A.S. The Skin Transcriptome in Hidradenitis Suppurativa Uncovers an Antimicrobial and Sweat Gland Gene Signature Which Has Distinct Overlap with Wounded Skin. PLoS ONE 2019, 14, e0216249. [Google Scholar] [CrossRef]
  103. Guo, L.; Chen, H.; Li, Y.; Zhou, Q.; Sui, Y. An Aquaporin 3-Notch1 Axis in Keratinocyte Differentiation and Inflammation. PLoS ONE 2013, 8, e80179. [Google Scholar] [CrossRef]
  104. Tricarico, P.M.; Mentino, D.; de Marco, A.; del Vecchio, C.; Garra, S.; Cazzato, G.; Foti, C.; Crovella, S.; Calamita, G. Aquaporins Are One of the Critical Factors in the Disruption of the Skin Barrier in Inflammatory Skin Diseases. Int. J. Mol. Sci. 2022, 23, 4020. [Google Scholar] [CrossRef] [PubMed]
  105. Sabat, R.; Chanwangpong, A.; Schneider-Burrus, S.; Metternich, D.; Kokolakis, G.; Kurek, A.; Philipp, S.; Uribe, D.; Wolk, K.; Sterry, W. Increased Prevalence of Metabolic Syndrome in Patients with Acne Inversa. PLoS ONE 2012, 7, e31810. [Google Scholar] [CrossRef] [PubMed]
  106. Mintoff, D.; Benhadou, F.; Pace, N.P.; Frew, J.W. Metabolic Syndrome and Hidradenitis Suppurativa: Epidemiological, Molecular, and Therapeutic Aspects. Int. J. Dermatol. 2021, 61, 1175–1186. [Google Scholar] [CrossRef] [PubMed]
  107. Shalom, G.; Freud, T.; Harman-Boehm, I.; Polishchuk, I.; Cohen, A.D. Hidradenitis Suppurativa and Metabolic Syndrome: A Comparative Cross-sectional Study of 3207 Patients. Br. J. Dermatol. 2015, 173, 464–470. [Google Scholar] [CrossRef]
  108. Alberti, K.G.M.; Zimmet, P.; Shaw, J. The Metabolic Syndrome—A New Worldwide Definition. Lancet 2005, 366, 1059–1062. [Google Scholar] [CrossRef]
  109. Cameron, A.J.; Boyko, E.J.; Sicree, R.A.; Zimmet, P.Z.; Söderberg, S.; Alberti, K.G.M.M.; Tuomilehto, J.; Chitson, P.; Shaw, J.E. Central Obesity as a Precursor to the Metabolic Syndrome in the AusDiab Study and Mauritius. Obesity 2008, 16, 2707–2716. [Google Scholar] [CrossRef]
  110. Kromann, C.; Ibler, K.; Kristiansen, V.; Jemec, G. The Influence of Body Weight on the Prevalence and Severity of Hidradenitis Suppurativa. Acta Derm. Venereol. 2014, 94, 553–557. [Google Scholar] [CrossRef]
  111. Balgobind, A.; Finelt, N.; Strunk, A.; Garg, A. Association between Obesity and Hidradenitis Suppurativa among Children and Adolescents: A Population-Based Analysis in the United States. J. Am. Acad. Dermatol. 2020, 82, 502–504. [Google Scholar] [CrossRef]
  112. Kaleta, K.P.; Nikolakis, G.; Hossini, A.M.; Balthasar, O.; Almansouri, D.; Vaiopoulos, A.; Knolle, J.; Boguslawska, A.; Wojas-Pelc, A.; Zouboulis, C.C. Metabolic Disorders/Obesity Is a Primary Risk Factor in Hidradenitis Suppurativa: An Immunohistochemical Real-World Approach. Dermatology 2022, 238, 251–259. [Google Scholar] [CrossRef]
  113. Zouboulis, V.A.; Zouboulis, K.C.; Zouboulis, C.C. Hidradenitis Suppurativa and Comorbid Disorder Biomarkers, Druggable Genes, New Drugs and Drug Repurposing—A Molecular Meta-Analysis. Pharmaceutics 2021, 14, 44. [Google Scholar] [CrossRef]
  114. Curat, C.A.; Wegner, V.; Sengenès, C.; Miranville, A.; Tonus, C.; Busse, R.; Bouloumié, A. Macrophages in Human Visceral Adipose Tissue: Increased Accumulation in Obesity and a Source of Resistin and Visfatin. Diabetologia 2006, 49, 744–747. [Google Scholar] [CrossRef]
  115. Wolk, K.; Sabat, R. Adipokines in Psoriasis: An Important Link between Skin Inflammation and Metabolic Alterations. Rev. Endocr. Metab. Disord. 2016, 17, 305–317. [Google Scholar] [CrossRef] [PubMed]
  116. Kovács, D.; Fazekas, F.; Oláh, A.; Törőcsik, D. Adipokines in the Skin and in Dermatological Diseases. Int. J. Mol. Sci. 2020, 21, 9048. [Google Scholar] [CrossRef] [PubMed]
  117. Yamauchi, T.; Kamon, J.; Waki, H.; Terauchi, Y.; Kubota, N.; Hara, K.; Mori, Y.; Ide, T.; Murakami, K.; Tsuboyama-Kasaoka, N.; et al. The Fat-Derived Hormone Adiponectin Reverses Insulin Resistance Associated with Both Lipoatrophy and Obesity. Nat. Med. 2001, 7, 941–946. [Google Scholar] [CrossRef] [PubMed]
  118. Gustafson, B.; Smith, U. Cytokines Promote Wnt Signaling and Inflammation and Impair the Normal Differentiation and Lipid Accumulation in 3T3-L1 Preadipocytes. J. Biol. Chem. 2006, 281, 9507–9516. [Google Scholar] [CrossRef]
  119. Gustafson, B.; Hammarstedt, A.; Andersson, C.X.; Smith, U. Inflamed Adipose Tissue. Arterioscler. Thromb. Vasc. Biol. 2007, 27, 2276–2283. [Google Scholar] [CrossRef]
  120. Andersson, C.X.; Gustafson, B.; Hammarstedt, A.; Hedjazifar, S.; Smith, U. Inflamed Adipose Tissue, Insulin Resistance and Vascular Injury. Diabetes. Metab. Res. Rev. 2008, 24, 595–603. [Google Scholar] [CrossRef] [PubMed]
  121. Ouchi, N.; Kihara, S.; Arita, Y.; Nishida, M.; Matsuyama, A.; Okamoto, Y.; Ishigami, M.; Kuriyama, H.; Kishida, K.; Nishizawa, H.; et al. Adipocyte-Derived Plasma Protein, Adiponectin, Suppresses Lipid Accumulation and Class A Scavenger Receptor Expression in Human Monocyte-Derived Macrophages. Circulation 2001, 103, 1057–1063. [Google Scholar] [CrossRef]
  122. Yamaguchi, N.; Argueta, J.G.M.; Masuhiro, Y.; Kagishita, M.; Nonaka, K.; Saito, T.; Hanazawa, S.; Yamashita, Y. Adiponectin Inhibits Toll-like Receptor Family-Induced Signaling. FEBS Lett. 2005, 579, 6821–6826. [Google Scholar] [CrossRef]
  123. Wolf, A.M.; Wolf, D.; Rumpold, H.; Enrich, B.; Tilg, H. Adiponectin Induces the Anti-Inflammatory Cytokines IL-10 and IL-1RA in Human Leukocytes. Biochem. Biophys. Res. Commun. 2004, 323, 630–635. [Google Scholar] [CrossRef]
Figure 2. Representation of the HS pathogenesis based on the three main outcomes found by the meta-analysis: upregulated inflammation, altered epithelial organization, and dysregulated metabolic signaling; 1: It is believed that keratinocytes may have an abnormal inflammatory profile in HS lesions. These keratinocytes expressing cytokines such as TNF-α, IL-1β, IFN-γ, alongside pro-inflammatory adipokines, may induce infundibular hyperkeratosis and subsequent perifollicular immune cell infiltration; 2: Hyperkeratosis occludes the hair follicle, resulting in follicular hyperplasia and nodule formation demonstrated by upregulation of KRT6-KRT16 and SPRR3; 3: Non-stop follicular dilatation leads to follicle rupture and impairment of the skin barrier function hinted at due to the upregulation of serpins-family genes, exacerbating inflammation. Recurrent injured nodules may evolve into dermal tunnels with impaired wound healing capacity, a consequence of abnormal eccrine gland function suggested by downregulation of WIF1, AQP5, FOXA1, and DCD.
Figure 2. Representation of the HS pathogenesis based on the three main outcomes found by the meta-analysis: upregulated inflammation, altered epithelial organization, and dysregulated metabolic signaling; 1: It is believed that keratinocytes may have an abnormal inflammatory profile in HS lesions. These keratinocytes expressing cytokines such as TNF-α, IL-1β, IFN-γ, alongside pro-inflammatory adipokines, may induce infundibular hyperkeratosis and subsequent perifollicular immune cell infiltration; 2: Hyperkeratosis occludes the hair follicle, resulting in follicular hyperplasia and nodule formation demonstrated by upregulation of KRT6-KRT16 and SPRR3; 3: Non-stop follicular dilatation leads to follicle rupture and impairment of the skin barrier function hinted at due to the upregulation of serpins-family genes, exacerbating inflammation. Recurrent injured nodules may evolve into dermal tunnels with impaired wound healing capacity, a consequence of abnormal eccrine gland function suggested by downregulation of WIF1, AQP5, FOXA1, and DCD.
Biomolecules 12 01371 g002
Table 1. Detailed information regarding the three selected studies matching our study criteria.
Table 1. Detailed information regarding the three selected studies matching our study criteria.
SRATitleSamples Included in Our StudyMain Findings
SRP273719Immunopathogenesis of hidradenitis suppurativa and response to anti–TNF-α therapy42 samples (HS skin lesion pre-TNF = 19: HS skin lesion mild-moderate HS = 7: healthy skin control = 16)Highly enriched pathways in HS lesioned skin are immune related. Signatures of complement activation, B cell signaling, and pathways involving phagocytosis were found to be unique to HS. TNF-α–regulated genes, IFN-γ, and IL-1β were selected as the major drivers of the inflammatory pathways in HS skin lesions. Nonetheless, IL-1 receptor antagonist, IL-1RN, and IL-10RA, 2 potent immunoregulatory molecules, were relatively reduced in HS skin. Alongside, α-catenin and sirtuin 1, both important for regulation of cell proliferation and survival, were reduced in HS skin.
SRP272659Contribution of plasma cells and B cells to hidradenitis suppurativa pathogenesis32 samples (HS skin lesion = 22: healthy skin control = 10)Several upregulated genes in the skin were associated with B cell responses, including immunoglobulin genes such as IGLV3-27, CD19, and CD79a. Other important genes found were the antimicrobial gene DEFB4A; CXCL13, a B cell chemoattractant, and the neutrophil chemokine CXCL1. In summary, they found B cells, and in particular plasma cells, as a potential therapeutic target in HS.
SRP168930The IL-1 pathway is hyperactive in hidradenitis suppurativa and contributes to skin infiltration and destruction7 samples (HS skin lesion = 3: healthy skin control = 4)IL-1β is highly active in HS, contributing to local and systemic inflammation. IL-1β induces expression of many molecules involved in extracellular matrix destruction including MMPs, ADAM12, serpinA1, COL3A1, and COL10A1, and immune cell infiltration such as CXCL1, CXCL6, CCL7, CXCL10, CXCL16, CXCL13, CCL24, CCL2, CCL8, and CCL20. IL-1β and, therefore, MMP1, MMP3, MMP9, MMP10, CCL2, CXCL1, IL-6, and IL-32 were upregulated when compared with healthy control skin and psoriasis lesions.
Table 2. Reactome analysis highlighting the 20 pathways with the lowest FDR, which are associated with upregulated mDEGs.
Table 2. Reactome analysis highlighting the 20 pathways with the lowest FDR, which are associated with upregulated mDEGs.
Pathway IdentifierPathway NameEntities FoundFDR
R-HSA-198933Immunoregulatory interactions between a lymphoid and a non-lymphoid cell817.2839 × 10−31
R-HSA-1474244Extracellular matrix organization1092.4259 × 10−14
R-HSA-380108Chemokine receptors bind chemokines351.4862 × 10−11
R-HSA-373076Class A/1 (Rhodopsin-like receptors)1062.2223 × 10−11
R-HSA-202430Translocation of ZAP-70 to immunological synapse161.2195 × 10−9
R-HSA-389948PD-1 signaling181.55 × 10−9
R-HSA-6798695Neutrophil degranulation1325.2463 × 10−9
R-HSA-500792GPCR ligand binding1265.25 × 10−9
R-HSA-909733Interferon alpha/beta signaling335.2463 × 10−9
R-HSA-202433Generation of second messenger molecules225.2463 × 10−9
R-HSA-202427Phosphorylation of CD3 and TCR zeta chains175.2463 × 10−9
R-HSA-1474228Degradation of the extracellular matrix551.2723 × 10−8
R-HSA-877300Interferon gamma signaling403.5328 × 0−8
R-HSA-1442490Collagen degradation326.6811 × 10−8
R-HSA-375276Peptide ligand-binding receptors677.395 × 10−8
R-HSA-6809371Formation of the cornified envelope508.368 × 10−8
R-HSA-6785807Interleukin-4 and Interleukin-13 signaling441.1813 × 10−7
R-HSA-6783783Interleukin-10 signaling252.22 × 10−7
R-HSA-1474290Collagen formation377.0322 × 10−7
R-HSA-418594G alpha (i) signaling events1081.4245 × 10−6
Table 3. Distribution of potential pathogenic variants with AF ≤ 0.01 among genes that have been already associated with HS.
Table 3. Distribution of potential pathogenic variants with AF ≤ 0.01 among genes that have been already associated with HS.
GeneDEGsSNP IDRefAltDistribution of Genotypes among HS PatientsWilcoxon TestExonic FunctionHGVSAF
NCSTN-rs35603924GCGG (50)/GC (1)0.1087nonsynonymous SNVc.G231C:p.E77D0.00432
APH1A-rs996158631ATAA (50)/AT (1)0.5630nonsynonymous SNVc.T123A:p.D41E0.00010
APH1B-rs142676640CTCC (50)/CT (1)0.2229nonsynonymous SNVc.C640T:p.R2140.000699
PSEN2-rs143912759CACC (50)/CA (1)0.4288nonsynonymous SNVc.C1139A:p.T380K0.00026
PSEN1-rs1174374799-T−(50)/−T (1)0.4542frameshift insertionc.526dupT:p.S178Ffs * 100.000006573
PSENEN-rs751542345TGTT (50)/TG (1)0.1519stopgainc.T168G:p.Y56X0.000008
FGFR2-rs56226109GAGG (49)/GA (2)0.2413nonsynonymous SNVc.C170T:p.S57L0.003722
MEFVuprs28940579AGAA (50)/AG (1)0.0215nonsynonymous SNVc.T2177C:p.V726A0.001440
MEFVuprs104895094TCTT (50)/TC (1)0.0996nonsynonymous SNVc.A2084G:p.K695R0.005245
NOD2-rs104895452CACC (50)/CA (1)0.9148nonsynonymous SNVc.C2672A:p.A891D0.000707
NOD2-rs5743279GAGG (50)/GA (1)0.1775nonsynonymous SNVc.G2288A:p.R763Q0.001217
NOD2-rs5743272AGAA (49)/AG (2)0.0699nonsynonymous SNVc.A974G:p.H325R0.000392
NOD2-rs35285618GAGG (50)/GA (1)0.2845nonsynonymous SNVc.G2042A:p.R681H0.00198
NOD2-rs2066847-C−(50)/−C (1)0.6377frameshift insertionc.2936dupC:p.L980Pfs * 20.015002
NOD2-rs34684955GAGG (50)/GA (1)0.2845nonsynonymous SNVc.G337A:p.A113T0.00251
NOD2-rs5743278CGCC (50)/CG (1)0.1038nonsynonymous SNVc.C2093G:p.A698G0.00371
NOD2-rs576658764CTCC (50)/CT (1)0.2845nonsynonymous SNVc.C1540T:p.R514W0.00007
PSTPIP1uprs34240327GCGG (49)/GC (2)0.4935nonsynonymous SNVc.G773C:p.G258A0.00461
Ref = reference; Alt = altered; HGVS = Human Genome Variation Society; AF = allelic frequency.
Table 4. Potential Clin Var pathogenic non-synonymous variants found by the RNA-Seq variant calling analysis in HS patients with an allele frequency ≤ 0.01, CADD score > 15, and genotypic differences that were statistically significant (p value < 0.05).
Table 4. Potential Clin Var pathogenic non-synonymous variants found by the RNA-Seq variant calling analysis in HS patients with an allele frequency ≤ 0.01, CADD score > 15, and genotypic differences that were statistically significant (p value < 0.05).
GenesDEGsSNP IDRefAltDistribution of Genotypes among HS PatientsWilcoxon TestHGVSAFCADD Score
ACSF3-rs144681140GAGG (49)/GA (2)0.0295c.G1406A:p.R469Q0.003222.9
KLF4-rs139237114GAGG (49)/GA (2)0.0385c.C859T:p.H287Y0.001623.8
DUSP23-rs11544443ATAA (49)/AT (2)0.0475c.A371T:p.E124V0.002225.1
BTN2A1-rs143104579GAGG (48)/GA (3)0.045c.G188A:p.R63H0.009618.01
FLOT2-rs3736238CTCC (49)/CT (2)0.0284c.G982A:p.A328T0.011918.44
GPANK1-rs35265780GAGG (49)/GA (2)0.04c.C335T:p.A112V0.009632
FNIP2-rs62001914CACC (49)/CA (2)0.0476c.C1653A:p.S551R0.009224.8
CORO1B-rs145707942CGCC (49)/CG (2)0.0462c.G367C:p.E123Q0.000224.2
ADCY4-rs61745073TATT (49)/TA (2)0.0357c.A1358T:p.E453V0.002222.7
AKR1C3downrs34186955CTCC (49)/CT (2)0.0295c.C538T:p.P180S0.008623.3
ALDH6A1downrs139579994GAGG (49)/GA (2)0.0395c.C716T:p.P239L0.001827.9
GSDMDuprs62000416CACC (49)/CA (2)0.0315c.C556A:p.L186M0.005623.4
YTHDF1-rs141487890GAGG (49)/GA (2)0.0344c.C437T:p.A146V0.000824
MYOF-rs61861290GAGG (48)/GA (3)0.0207c.C4576T:p.P1526S0.006226.2
SIT1uprs138786883CACC (49)/CA (2)0.0496c.G520T:p.A174S0.003217.54
RBMXL1-rs139713926TCTT (49)/TC (2)0.0242c.A701G:p.Y234C0.002225.5
GALNT7-rs144873913CACC (48)/CA (2)0.0496c.C1585A:p.P529T0.001430
PHACTR4-rs72661785GCGG (49)/GC (2)0.043c.G1609C:p.A537P0.003626.3
WDR92-rs138784630CTCC (49)CT (2)0.04c.G841A:p.A281T0.007223.4
Ref = reference; Alt = altered; HGVS = Human Genome Variation Society sequence variant nomenclature; AF = allelic frequency; CADD score = mutation impact prediction.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

de Oliveira, A.S.L.E.; Bloise, G.; Moltrasio, C.; Coelho, A.; Agrelli, A.; Moura, R.; Tricarico, P.M.; Jamain, S.; Marzano, A.V.; Crovella, S.; et al. Transcriptome Meta-Analysis Confirms the Hidradenitis Suppurativa Pathogenic Triad: Upregulated Inflammation, Altered Epithelial Organization, and Dysregulated Metabolic Signaling. Biomolecules 2022, 12, 1371. https://doi.org/10.3390/biom12101371

AMA Style

de Oliveira ASLE, Bloise G, Moltrasio C, Coelho A, Agrelli A, Moura R, Tricarico PM, Jamain S, Marzano AV, Crovella S, et al. Transcriptome Meta-Analysis Confirms the Hidradenitis Suppurativa Pathogenic Triad: Upregulated Inflammation, Altered Epithelial Organization, and Dysregulated Metabolic Signaling. Biomolecules. 2022; 12(10):1371. https://doi.org/10.3390/biom12101371

Chicago/Turabian Style

de Oliveira, Ana Sofia Lima Estevao, Giovanna Bloise, Chiara Moltrasio, Antonio Coelho, Almerinda Agrelli, Ronald Moura, Paola Maura Tricarico, Stéphane Jamain, Angelo Valerio Marzano, Sergio Crovella, and et al. 2022. "Transcriptome Meta-Analysis Confirms the Hidradenitis Suppurativa Pathogenic Triad: Upregulated Inflammation, Altered Epithelial Organization, and Dysregulated Metabolic Signaling" Biomolecules 12, no. 10: 1371. https://doi.org/10.3390/biom12101371

APA Style

de Oliveira, A. S. L. E., Bloise, G., Moltrasio, C., Coelho, A., Agrelli, A., Moura, R., Tricarico, P. M., Jamain, S., Marzano, A. V., Crovella, S., & Cavalcanti Brandão, L. A. (2022). Transcriptome Meta-Analysis Confirms the Hidradenitis Suppurativa Pathogenic Triad: Upregulated Inflammation, Altered Epithelial Organization, and Dysregulated Metabolic Signaling. Biomolecules, 12(10), 1371. https://doi.org/10.3390/biom12101371

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