Integrating multi-omics data to analyze the potential pathogenic mechanism of CTSH gene involved in type 1 diabetes in the exocrine pancreas

Abstract Type 1 diabetes (T1D) is an autoimmune disease caused by the destruction of insulin-producing pancreatic islet beta cells. Despite significant advancements, the precise pathogenesis of the disease remains unknown. This work integrated data from expression quantitative trait locus (eQTL) studies with Genome wide association study (GWAS) summary data of T1D and single-cell transcriptome data to investigate the potential pathogenic mechanisms of the CTSH gene involved in T1D in exocrine pancreas. Using the summary data–based Mendelian randomization (SMR) approach, we obtained four potential causative genes associated with T1D: BTN3A2, PGAP3, SMARCE1 and CTSH. To further investigate these genes’roles in T1D development, we validated them using a scRNA-seq dataset from pancreatic tissues of both T1D patients and healthy controls. The analysis showed a significantly high expression of the CTSH gene in T1D acinar cells, whereas the other three genes showed no significant changes in the scRNA-seq data. Moreover, single-cell WGCNA analysis revealed the strongest positive correlation between the module containing CTSH and T1D. In addition, we found cellular ligand–receptor interactions between the acinar cells and different cell types, especially ductal cells. Finally, based on functional enrichment analysis, we hypothesized that the CTSH gene in the exocrine pancreas enhances the antiviral response, leading to the overexpression of pro-inflammatory cytokines and the development of an inflammatory microenvironment. This process promotes β cells injury and ultimately the development of T1D. Our findings offer insights into the underlying pathogenic mechanisms of T1D.