Chromatin Accessibility and Transcriptome Integrative Analysis Revealed AP-1-mediated Genes Potentially Modulate Histopathology Features in Psoriasis

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.