Analysis of Potential Genes and Pathways Involved in the Pathogenesis of Acne by Bioinformatics

Integrative bioinformatics and experimental validation of hub genetic markers in acne vulgaris: Toward personalized diagnostic and therapeutic strategies

BACKGROUND

Acne vulgaris is a widespread chronic inflammatory dermatological condition. The precise molecular and genetic mechanisms of its pathogenesis remain incompletely understood. This research synthesizes existing databases, targeting a comprehensive exploration of core genetic markers.

METHODS

Gene expression datasets (GSE6475, GSE108110, and GSE53795) were retrieved from the GEO. Differentially expressed genes (DEGs) were identified using the limma package. Enrichment analyses were conducted using GSVA for pathway assessment and clusterProfiler for GO and KEGG analyses. PPI networks and immune cell infiltration were analyzed using the STRING database and ssGSEA, respectively. We investigated the correlation between hub gene biomarkers and immune cell infiltration using Spearman's rank analysis. ROC curve analysis validated the hub genes' diagnostic accuracy. miRNet, TarBase v8.0, and ChEA3 identified miRNA/transcription factor-gene interactions, while DrugBank delineated drug-gene interactions. Experiments utilized HaCaT cells stimulated with Propionibacterium acnes, treated with retinoic acid and methotrexate, and evaluated using RT-qPCR, ELISA, western blot, lentiviral transduction, CCK-8, wound-healing, and transwell assays.

RESULTS

There were 104 genes with consistent differences across the three datasets of paired acne and normal skin. Functional analyses emphasized the significant enrichment of these DEGs in immune-related pathways. PPI network analysis pinpointed hub genes PTPRC, CXCL8, ITGB2, and MMP9 as central players in acne pathogenesis. Elevated levels of specific immune cell infiltration in acne lesions corroborated the inflammatory nature of the disease. ROC curve analysis identified the acne diagnostic potential of four hub genes. Key miRNAs, particularly hsa-mir-124-3p, and central transcription factors like TFEC were noted as significant regulators. In vitro validation using HaCaT cells confirmed the upregulation of hub genes following Propionibacterium acnes exposure, while CXCL8 knockdown reduced pro-inflammatory cytokines, cell proliferation, and migration. DrugBank insights led to the exploration of retinoic acid and methotrexate, both of which mitigated gene expression upsurge and inflammatory mediator secretion.

CONCLUSION

This comprehensive study elucidated pivotal genes associated with acne pathogenesis, notably PTPRC, CXCL8, ITGB2, and MMP9. The findings underscore potential biomarkers, therapeutic targets, and the therapeutic potential of agents like retinoic acid and methotrexate. The congruence between bioinformatics and experimental validations suggests promising avenues for personalized acne treatments.

Identification of biomarkers of acne based on transcriptome analysis and combined with network pharmacology to explore the therapeutic mechanism of Jinhuang ointment

The incidence of acne is on the rise due to unhealthy diet and living habits. Jinhuang ointment (JHO) is a classic prescription composed of 10 kinds of commonly used Chinese herbal medicine, which has been widely used in clinical prevention and treatment of skin inflammatory diseases since ancient times. However, the pharmacological mechanism and target of JHO are not clear. The acne microarray dataset was downloaded from gene expression omnibus database to identify differentially expressed genes (DEG). Immune infiltration was analyzed by CiberSort algorithm. HUB gene was identified by protein-protein interaction network. The gene expression omnibus dataset validates the biomarkers of acne with high diagnostic value. The potential active components and targets of JHO were obtained through Traditional Chinese Medicine Systems Pharmacology database, and the therapeutic targets were obtained by crossing with disease targets. R-packet is used for enrichment analysis. Molecular docking using Auto Dock Tools. A total of 202 DEGs were identified from 12 skin samples in the GSE6475. Immune infiltration analysis showed that there were a large number of macrophages and mast cells in acne skin. Gene set enrichment analysis analysis showed that DEGS was mainly involved in bacterial reaction, inflammatory reaction and so on. Six central genes and gene cluster modules were identified by Cytoscape software. A total of 185 JHO active components and 220 targets were obtained, of which 10 targets were potential targets for JHO in the treatment of acne. Kyoto encyclopedia of genes and genomes enrichment analysis showed that JHO treatment of acne was mainly related to Toll-like receptors, IL-17 and other signal pathways. The results of molecular docking showed that 5 active compounds in JHO had strong binding activity to the core protein receptor. IL-1 β, CXCL8, toll-like receptor 2, CXCL2, LCN2, and secretory phosphoprotein 1 may be potential biomarkers for early diagnosis of acne. JHO active components may regulate skin cell metabolism and inflammatory response and improve cellular immune microenvironment by acting on core targets (CXCL8, ESR1, IL-1 β, MMP1, MMP3, secretory phosphoprotein 1), thus achieving the purpose of treating acne. This is the result of the joint action of multiple targets and multiple pathways. It provides an idea for the development of a new combination of drugs for the treatment of acne.

Keratinocytes use FPR2 to detect Staphylococcus aureus and initiate antimicrobial skin defense

Introduction

Keratinocytes form a multilayer barrier that protects the skin from invaders or injuries. The barrier function of keratinocytes is in part mediated by the production of inflammatory modulators that promote immune responses and wound healing. Skin commensals and pathogens such as Staphylococcus aureus secrete high amounts of phenol-soluble modulin (PSM) peptides, agonists of formyl-peptide receptor 2 (FPR2). FPR2 is crucial for the recruitment of neutrophils to the sites of infection, and it can influence inflammation. FPR1 and FPR2 are also expressed by keratinocytes but the consequences of FPR activation in skin cells have remained unknown.

Methods

Since an inflammatory environment influences S. aureus colonization, e. g. in patients with atopic dermatitis (AD), we hypothesized that interference with FPRs may alter keratinocyte-induced inflammation, proliferation, and bacterial colonization of the skin. To assess this hypothesis, we investigated the effects of FPR activation and inhibition in keratinocytes with respect to chemokine and cytokine release as well as proliferation and skin wound gap closure.

Results

We observed that FPR activation induces the release of IL-8, IL-1α and promotes keratinocyte proliferation in a FPR-dependent manner. To elucidate the consequence of FPR modulation on skin colonization, we used an AD-simulating S. aureus skin colonization mouse model using wild-type (WT) or Fpr2-/- mice and demonstrate that inflammation enhances the eradication of S. aureus from the skin in a FPR2-dependent way. Consistently, inhibition of FPR2 in the mouse model or in human keratinocytes as well as human skin explants promoted S. aureus colonization.

Discussion

Our data indicate that FPR2 ligands promote inflammation and keratinocyte proliferation in a FPR2-dependent manner, which is necessary for eliminating S. aureus during skin colonization.

Network Pharmacology-Based Investigation on the Mechanism of the JinGuanLan Formula in Treating Acne Vulgaris

Background

JinGuanLan (JGL) formula is a traditional Chinese medicine (TCM) developed by the Department of Pharmacology at the First Hospital of Lanzhou University. The network pharmacology approach was applied to determine the potential active compounds, therapeutic targets, and main pathways of the JGL formula to evaluate its application value in acne vulgaris.

Methods

Data on the active compounds and their related targets were obtained from the Traditional Chinese Medicine Systems Pharmacology Database and Analysis Platform (TCMSP). Acne vulgaris-related targets were searched from the Online Mendelian Inheritance in Man (OMIM) database, GeneCards Database, Comparative Toxicogenomics Database (CTD), Therapeutic Target Database (TTD), and DisGeNET Database. Targets intersecting between JGL- and acne vulgaris-related targets were chosen as potential therapeutic targets. The protein-protein interaction (PPI) network of potential therapeutic targets was visualized using Cytoscape software based on the PPI data collected from the STRING database. Three topological features, namely, "Degree," "MCC," and "EPC" of each node in the PPI network were calculated using the cytoHubba plugin of Cytoscape to excavate the core targets. R program was used for the Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis of the potential therapeutic targets. Finally, the compound-target-pathway network was constructed.

Result

Among the 148 active compounds that were identified, quercetin and kaempferol showed the highest degree of target interaction and thus may play essential roles in the pharmacological effect of the JGL formula for acne treatment. Among the 97 potential therapeutic targets that were screened out, the 6 core targets were TNF, JUN, IL6, STAT3, MAPK1, and MAPK3. A total of 2260 terms of GO enrichment analysis were obtained, including 2090 for biological processes (BP), 37 for cellular components (CC), and 133 for molecular function (MF). A total of 156 enriched KEGG pathways were identified, including TNF, IL-17, Th17 cell differentiation, MAPK, PI3K-Akt, T cell receptor, and Toll-like receptor signalling pathways.

Conclusion

This work showed that the JGL formula might reverse the pathological changes associated with acne vulgaris through its antiinflammatory effect and regulate the excessive lipogenesis in sebaceous glands via different signalling pathways. This new drug has application value and is worthy of further research and development.

Characterization of a Cutibacterium acnes Camp Factor 1-Related Peptide as a New TLR-2 Modulator in In Vitro and Ex Vivo Models of Inflammation

Cutibacterium acnes (C. acnes) has been implicated in inflammatory acne where highly mutated Christie-Atkins-Munch-Petersen factor (CAMP)1 displays strong toll like receptor (TLR)-2 binding activity. Using specific antibodies, we showed that CAMP1 production was independent of C. acnes phylotype and involved in the induction of inflammation. We confirmed that TLR-2 bound both mutated and non-mutated recombinant CAMP1, and peptide array analysis showed that seven peptides (A14, A15, B1, B2, B3, C1 and C3) were involved in TLR-2 binding, located on the same side of the three-dimensional structure of CAMP1. Both mutated and non-mutated recombinant CAMP1 proteins induced the production of C-X-C motif chemokine ligand interleukin (CXCL)8/(IL)-8 in vitro in keratinocytes and that of granulocyte macrophage-colony stimulating factor (GM-CSF), tumor necrosis factor (TNF)-α, IL-1β and IL-10 in ex vivo human skin explants. Only A14, B1 and B2 inhibited the production of CXCL8/IL-8 by keratinocytes and that of (GM-CSF), TNF-α, IL-1β and IL-10 in human skin explants stimulated with rCAMP1 and C. acnes. Following pretreatment with B2, RNA sequencing on skin explants identified the 10 genes displaying the strongest differential expression as IL6, TNF, CXCL1, CXCL2, CXCL3, CXCL8, IL-1β, chemokine ligand (CCL)2, CCL4 and colony stimulating factor (CSF)2. We, thus, identified a new CAMP1-derived peptide as a TLR-2 modulator likely to be a good candidate for clinical evaluation.

Gene variants associated with acne vulgaris presentation and severity: a systematic review and meta-analysis

BACKGROUND

Multiple factors have been attributed to acne vulgaris predisposition and individual variations in the severity of skin symptoms, and genetics stood out as one of the major factors.

METHODS

We performed a systematic review on the genes and their variants that have been investigated for association with acne presentation and severity. A random-effect meta-analysis using the allele model (minor allele vs. major allele) was also conducted to provide an overall estimation of risk effects of frequently reported gene variants. This included a subset data of 982 acne cases and 846 controls extracted from our existing GWAS database on various allergic and skin diseases among Singapore Chinese.

RESULTS

Systematic review of 51 articles covering Asians and Caucasians found 60 genes/loci and their 100 variants implicated in acne; majority of them were in the intron, coding region/missense, and promoter regions. The commonly studied candidate genes/gene families include tumor necrosis factor (TNF), and the interleukin (IL) and cytochrome P450 (CYP) gene families. Our meta-analysis showed that most of the analyzed gene variants exhibited insignificant pooled odds ratio (pOR) and significant heterogeneity between studies. Nevertheless, we found that TNF rs1800629 A allele carriers and CYP17A1 rs743572 T allele carriers had significantly reduced mild acne risk [pOR: 0.60; 95% Confidence Interval (CI): 0.33-0.86] and severe acne risk (pOR: 0.59; 95% CI: 0.40-0.79), respectively, across populations. Overall, FST (follistatin) rs629725 A allele poses a significantly modest increased risk for acne presentation (pOR: 1.19, 95% CI: 1.14, 1.23), but neither TIMP2 (TIMP metallopeptidase inhibitor 2) rs8179090 nor CYP1A1 rs4646903 (pOR: 0.96, 95% CI: 0.80-1.12; pOR: 0.95, 95% CI: 0.83, 1.08), respectively. We discovered 15 novel SNPs in the 3' UTR region of the Toll-like Receptor 4 gene (TLR4) associated with acne presentation.

CONCLUSIONS

This systematic review and meta-analysis suggest that genes influencing inflammatory responses, specifically TNF, and genes influencing the function and activity of sebaceous glands, specifically CYP17A1 and FST, have potential risk variants for acne presentation and severity across populations. Understanding the genetic susceptibility factors and biological pathways involved in the pathogenesis of acne will help us to gain insights into developing effective acne treatments.

Altered Gene Expression in Acne Vulgaris Patients Treated by Oral Isotretinoin: A Preliminary Study

Background/Objective

The role of gene expression changes in acne patients treated by oral isotretinoin (ISO) and in influencing the ISO therapeutic effects is still unclear. In this study, we investigated the gene profiles of patients with severe acne who responded variously to ISO therapy.

Methods

The peripheral blood of 113 acne vulgaris patients (Pillsbury IV grade) was collected before treatment. After 8 weeks of oral ISO, nine acne patients were selected and divided into the following groups. A: effectively treated by ISO, group B: ineffectively treated by ISO, group C: ISO-induced acne flare-up, and 3 healthy subjects were included as control group D. The peripheral blood of patients pre- and post-treatment was subjected to high-throughput RNA sequencing technology and bioinformatics analysis of the separate groups (n = 3). The candidate genes were validated by qRT-PCR.

Results

Comparing pre- and post-oral ISO treatment, gene expression was changed as 39 genes in ISO-effective group, 345 genes in ISO-ineffective group, and 57 genes in ISO-induced acne flare-up group. Comparing the ISO-induced acne flare-up group with healthy control subjects revealed 34 upregulated genes and 23 downregulated genes, while comparing the ISO-induced acne flare-up group with ISO-ineffective patients identified 1835 changed genes. Expression of GATA2 (2.73 fold, P=0.024512), C4BPA (35.87 folds, P=0.038073), and CCR5 (2.48 folds, P=0.004681) increased in the ISO-induced acne flare-up patients. Meanwhile, the expression of DEFA3 (0.18 fold, P=0.041934), ELANE (0.14 fold, P=0.030767), MMP9 (0.41 fold, P=0.013383), and RPS4Y1 (0.00018 fold, P=0.000986) decreased when compared with ISO-ineffective patients.

Conclusion

Oral ISO treatment could temporarily alter gene expression in acne patients. ISO therapeutic mechanisms were involved, not only in regulating the inflammatory reaction but also in the process of DNA repair. GATA2, C4BPA, CCR5, DEFA3, ELANE, MMP9, and RPS4Y1 might be susceptible to genes that could participate in the ISO-induced aggravation of acne.

Putative Genes and Pathways Involved in the Acne Treatment of Isotretinoin via Microarray Data Analyses

Acne is the eighth most common disease worldwide. Disease burden of acne such as anxiety, reduced self-esteem, and facial scarring lowers the life quality of acne patients. Isotretinoin is the most potent treatment for moderate-severe acne. However, the adverse events of isotretinoin especially teratogenicity limit its use. This study aims at investigating the therapeutical mechanisms of isotretinoin using bioinformatics analysis. Differentially expressed genes (DEGs) were filtered from microarray datasets GSE10432, GSE10433, and GSE11792. Functional and pathway enrichment analyses of DEGs were performed. Protein–protein interaction (PPI) network and module analyses were also conducted based on DEGs. Using isotretinoin for 1 week, LCN2, PTGES, and GDF15 were upregulated and might mediate sebocytes apoptosis and thus decreased sebum production; CCL2 originated from activated TNF signaling pathway and S100A7 could be related with “acne-flare”. While treating with isotretinoin for 8 weeks, key genes were downregulated, including HMGCS1, HMGCR, FDFT1, MVD, IDI1, and FDPS, which may be associated with decreased sebum synthesis; HMGCS1, HMGCR, and FDFT1 also probably associated with apoptosis of sebocytes. There were only two common genes including ACSBG1 and BCAT2 which worked in both 1 week and 8 weeks and could associate with decreased sebum synthesis and apoptosis of sebocytes, respectively. These results indicate potential therapeutics and side effect mechanisms of isotretinoin in the acne treatment and provide a research direction to further investigate the therapeutic mechanism of isotretinoin and thus develop retinoid-like compounds with similar curative effect and without teratogenicity.

Review-Current Concepts in Inflammatory Skin Diseases Evolved by Transcriptome Analysis: In-Depth Analysis of Atopic Dermatitis and Psoriasis

During the last decades, high-throughput assessment of gene expression in patient tissues using microarray technology or RNA-Seq took center stage in clinical research. Insights into the diversity and frequency of transcripts in healthy and diseased conditions provide valuable information on the cellular status in the respective tissues. Growing with the technique, the bioinformatic analysis toolkit reveals biologically relevant pathways which assist in understanding basic pathophysiological mechanisms. Conventional classification systems of inflammatory skin diseases rely on descriptive assessments by pathologists. In contrast to this, molecular profiling may uncover previously unknown disease classifying features. Thereby, treatments and prognostics of patients may be improved. Furthermore, disease models in basic research in comparison to the human disease can be directly validated. The aim of this article is not only to provide the reader with information on the opportunities of these techniques, but to outline potential pitfalls and technical limitations as well. Major published findings are briefly discussed to provide a broad overview on the current findings in transcriptomics in inflammatory skin diseases.

Integrated Analysis of Key Genes and Pathways Involved in Nonalcoholic Steatohepatitis Improvement After Roux-en-Y Gastric Bypass Surgery

BACKGROUND

As the incidence of nonalcoholic fatty liver disease (NAFLD) increases globally, nonalcoholic steatohepatitis (NASH) has become the second common cause of liver transplantation for liver diseases. Recent evidence shows that Roux-en-Y gastric bypass (RYGB) surgery obviously alleviates NASH. However, the mechanism underlying RYGB induced NASH improvement is still elusive.

METHODS

We obtained datasets, including hepatic gene expression data and histologic NASH status, at baseline and 1 year after RYGB surgery. Differentially expressed genes (DEGs) were identified comparing gene expression before and after RYGB surgery in each dataset. Common DEGs were obtained between both datasets and further subjected to functional and pathway enrichment analysis. Protein-protein interaction (PPI) network was constructed, and key modules and hub genes were also identified.

RESULTS

In the present study, GSE106737 and GSE83452 datasets were included. One hundred thirty common DEGs (29 up-regulated and 101 down-regulated) were identified between GSE106737 and GSE83452 datasets. KEGG analysis showed that mineral absorption, IL-17 signaling pathway, osteoclast differentiation, and TNF signaling pathway were significantly enriched. Based on the PPI network, IGF1, JUN, FOS, LDLR, TYROBP, DUSP1, CXCR4, ATF3, CXCL2, EGR1, SAA1, CTSS, and PPARA were identified as hub genes, and three functional modules were also extracted.

CONCLUSION

This study identifies the global gene expression change in the liver of NASH patients before and after RYGB surgery in a bioinformatic method. Our findings will contribute to the understanding of molecular biological changes underlying NASH improvement after RYGB surgery.