RESUMO
The aryl hydrocarbon receptor (AHR) is an evolutionary conserved environmental sensor identified as an indispensable regulator of epithelial homeostasis and barrier organ function. Molecular signaling cascade and target genes upon AHR activation and their contribution to cell and tissue function are however not fully understood. Multiomics analyses using human skin keratinocytes revealed that upon ligand activation, AHR binds open chromatin to induce expression of transcription factors, for example, TFAP2A, as a swift response to environmental stimuli. The terminal differentiation program, including upregulation of barrier genes, FLG and keratins, was mediated by TFAP2A as a secondary response to AHR activation. The role of AHR-TFAP2A axis in controlling keratinocyte terminal differentiation for proper barrier formation was further confirmed using CRISPR/Cas9 in human epidermal equivalents. Overall, the study provides additional insights into the molecular mechanism behind AHR-mediated barrier function and identifies potential targets for the treatment of skin barrier diseases.
RESUMO
The recent development of single-cell techniques is essential to unravel complex biological systems. By measuring the transcriptome and the accessible genome on a single-cell level, cellular heterogeneity in a biological environment can be deciphered. Transcription factors act as key regulators activating and repressing downstream target genes, and together they constitute gene regulatory networks that govern cell morphology and identity. Dissecting these gene regulatory networks is crucial for understanding molecular mechanisms and disease, especially within highly complex biological systems. The gene regulatory network analysis software ANANSE and the motif enrichment software GimmeMotifs were both developed to analyse bulk datasets. We developed scANANSE, a software pipeline for gene regulatory network analysis and motif enrichment using single-cell RNA and ATAC datasets. The scANANSE pipeline can be run from either R or Python. First, it exports data from standard single-cell objects. Next, it automatically runs multiple comparisons of cell cluster data. Finally, it imports the results back to the single-cell object, where the result can be further visualised, integrated, and interpreted. Here, we demonstrate our scANANSE pipeline on a publicly available PBMC multi-omics dataset. It identifies well-known cell type-specific hematopoietic factors. Importantly, we also demonstrated that scANANSE combined with GimmeMotifs is able to predict transcription factors with both activating and repressing roles in gene regulation.