Anatomical Region-Specific Transcriptomic Signatures and the Role of Epithelial Cells in Pterygium Inflammation: A Multi-Omics Analysis

BackgroundThe anatomy of pterygium consists of the head, neck, and body. Prior studies have examined the three parts collectively, yet the treatment methods and their success rates differ based on the anatomical characteristics. In this study, we divide the pterygium tissue into Main (head and neck) and Acc (body) to investigate its pathogenesis across regions. PurposesThis study aims to understand the pathogenesis and identify potential therapeutic targets of pterygium. Through a robust multi-omics analysis of bulk and single-cell RNA sequencing (RNA-seq) data from pterygium patients, we focus on the expression of inflammatory and mitochondrial energy pathways and identify potential genes responsible for the upregulated pathways in pterygium. MethodsWe collected bulk RNA-seq data from six pterygium patients and single-cell RNA-seq data from two pterygium patients. We then investigated the pathway enrichment, pathway correlation, differential gene expressions, protein-protein interactions, and cell-cell communications of pterygium. ResultsFrom the analysis of bulk RNA-seq data, the distribution of sample points from the principal component analysis plot and the pathways enriched from the Gene Ontology analysis showed distinct expression patterns in the Acc group compared to the Main and control groups. This suggested the need to separate the Main and Acc regions within pterygium samples and utilize single-cell RNA-seq data to understand the differences between the Main and control groups that the bulk data could not capture. The annotation of integrated single-cell data revealed a cluster of epithelial cells containing only pterygium samples. Cells from this cluster exhibited significant contributions to the ANGPTL, IL1, and KLK signaling networks in the cell-cell communication analysis. We also observed significant upregulation of the cells in the inflammatory pathways related to integrated stress response and the renin-angiotensin-aldosterone system, both of which showed high correlations with energy metabolism pathways. Significant changes in the expression of multiple pro-inflammatory, antioxidant, and immune-related genes were also identified. ConclusionThe different expression patterns between the Main and Acc groups suggest the need to consider different anatomical regions separately in future studies of pterygium. Additionally, the significant role epithelial cells from the Main group play in the inflammation of pterygium presents a potential clinical approach to the disease.