Mechanisms of Berberine in anti-pancreatic ductal adenocarcinoma revealed by integrated multi-omics profiling

Abstract

Abstract This study integrates pharmacology databases with bulk RNA-seq and scRNA-seq to reveal the latent anti-PDAC capacities of BBR. Target genes of BBR were sifted through TargetNet, CTD, SwissTargetPrediction, and Binding Database. Based on the GSE183795 dataset, DEG analysis, GSEA, and WGCNA were sequentially run to build a disease network. Through sub-network filtration acquired PDAC-related hub genes. A PPI network was established using the shared genes. Degree algorithm from cytoHubba screened the key cluster in the network. Analysis of differential mRNA expression and ROC curves gauged the diagnostic performance of clustered genes. CYBERSORT uncovered the potential role of the key cluster on PDAC immunomodulation. ScRNA-seq analysis evaluated the distribution and expression profile of the key cluster at the single-cell level, assessing enrichment within annotated cell subpopulations to delineate the target distribution of BBR in PDAC. We identified 425 drug target genes and 771 disease target genes, using 57 intersecting genes to construct the PPI network. CytoHubba anchored the top 10 highest contributing genes to be the key cluster. mRNA expression levels and ROC curves confirmed that these genes showed good robustness for PDAC. CYBERSORT revealed that the key cluster influenced immune pathways predominantly associated with Macrophages M0, CD8 T cells, and naïve B cells. ScRNA-seq analysis clarified that BBR mainly acted on epithelial cells and macrophages in PDAC tissues. BBR potentially targets CDK1, CCNB1, CTNNB1, CDK2, TOP2A, MCM2, RUNX2, MYC, PLK1, and AURKA to exert therapeutic effects on PDAC. The mechanisms of action appear to significantly involve macrophage polarization-related immunological responses.

Keywords

• Berberine
• Pancreatic ductal adenocarcinoma
• Single-cell RNA sequencing
• Bulk RNA sequencing
• Multi-omics analysis