Background Psoriasis is a chronic and hyperproliferative skin disease featured by hyperkeratosis with parakeratosis, Munro micro-abscess, elongation of rete pegs, granulosa thinning, and lymphocyte infiltration. We previously profiled gene expression and chromatin accessibility of psoriatic skins by transcriptome sequencing and ATAC-seq. However, integrating both of these datasets to unravel gene expression regulation is lacking. Here, we integrated transcriptome and ATAC-seq of the same psoriatic and normal skin tissues, trying to leverage the potential role of chromatin accessibility and their function in histopathology features. Results By inducing binding and expression target analysis (BETA) algorithms, we explored the target prediction of transcription factors binding in 15 psoriatic and 19 control skins. BETA identified 408 upregulated genes (rank product < 0.01) and 133 downregulated genes linked with chromatin accessibility. We noticed that cumulative fraction of genes in upregulation group was statistically higher than background, while that of genes in downregulation group was not significant. KEGG pathway analysis showed that the upregulated 408 genes were enriched in TNF, NOD, and IL-17 signaling pathways. In addition, the motif module in BETA suggested the 57 upregulated genes are targeted by transcription factor AP-1, indicating that increased chromatin accessibility facilitated the binding of AP-1 to the target regions and further induced expression of relevant genes. Among these genes, SQLE, STRN, EIF4, and MYO1B expression was increased in patients with hyperkeratosis, parakeratosis, and acanthosis thickening. Conclusions In summary, with the advantage of BETA, we identified a series of genes that contribute to the disease pathogenesis, especially in modulating histopathology features, providing us with new clues in treating psoriasis.

IntroductionMunro's microabscess is a typical pathological feature in the early psoriatic lesion, mainly characterized by the accumulation of neutrophils in the epidermis. DNA methylation microenvironment of Munro's microabscess and the crosstalk with transcription and its effect on neutrophils have not yet been revealed.MethodsPerformed genome-wide DNA methylation analysis and further differential methylation analysis of psoriatic skin lesions with and without Munro's microabscess from two batch samples consisting of 114 former samples in the discovery stage and 21 newly-collected samples in the validation stage. Utilized GO, MEME, and other tools to conduct downstream analysis on differentially methylated sites (DMSs). Correlation analysis of methylation level and transcriptome data was also conducted.ResultsWe observed 647 overlapping DMSs associated with Munro's microabscess. Subsequently, GO pathway analysis revealed that DNA methylation might affect the physical properties associated with skin cells through focal adhesion and cellsubstrate junction and was likely to recruit neutrophils in the epidermis. Via the MEME tool, used to investigate the possible binding transcription factors (TFs) of 20 motifs around the 647 DMSs, it was found that DNA methylation regulated the binding of AP1 family members and the recruitment of neutrophils in the epidermis through the TGF-beta pathway and the TH17 pathway. Meanwhile, combined with our earlier transcriptome data, we found DNA methylation would regulate the expressions of CFDP, SIRT6, SMG6, TRAPPC9, HSD17B7, and KIAA0415, indicating these genes would potentially promote the process of Munro's microabscess.DiscussionIn conclusion, DNA methylation may affect the course of psoriasis by regulating the progression of Munro's microabscess in psoriatic skin lesions.

Introduction: Systemic lupus erythematosus (SLE) is an autoimmune disorder characterized by the production of autoantibodies, immune complex deposition, and tissue/organ damage. In this study, we aimed to identify molecular features and signaling pathways associated with SLE severity using RNA sequencing (RNA-seq), single-cell RNA sequencing (scRNA-seq), and clinical parameters.Methods: We analyzed transcriptome profiles of 45 SLE patients, grouped into mild (mSLE, SLEDAI ≤ 9) and severe (sSLE, SLEDAI &gt; 9) based on SLE Disease Activity Index (SLEDAI) scores. We also collected clinical data on anti-dsDNA, ANA, ESR, CRP, snRNP, AHA, and anti-Smith antibody status for each patient.Results: By comparing gene expression across groups, we identified 12 differentially expressed genes (DEGs), including 7 upregulated (CEACAM6, UCHL1, ARFGEF3, AMPH, SERPINB10, TACSTD2, and OTX1) and 5 downregulated (SORBS2, TRIM64B, SORCS3, DRAXIN, and PCDHGA10) DEGs in sSLE compared to mSLE. Furthermore, using the CIBERSORT algorithm, we found that Treg cells were significantly decreased in sSLE and negatively correlated with AMPH expression, which was mainly expressed in Treg cells from SLE patients according to public scRNA-seq data (GSE135779).Discussion: Overall, our findings shed light on the molecular mechanisms underlying SLE severity and provide insight into potential therapeutic targets.