Abstract
Background: Primary Sjögren's syndrome (pSS) is a chronic systemic autoimmune disease of the exocrine glands characterized by specific pathological features. Previous studies have pointed out that salivary glands from pSS patients express a unique profile of cytokines, adhesion molecules, and chemokines compared to those from healthy controls. However, there is limited evidence supporting the utility of individual markers for different stages of pSS. This study aimed to explore potential biomarkers associated with pSS disease progression and analyze the associations between key genes and immune cells. Methods: We combined our own RNA sequencing data with pSS datasets from the NCBI Gene Expression Omnibus (GEO) database to identify differentially expressed genes (DEGs) via bioinformatics analysis. Salivary gland biopsies were collected from 14 pSS patients, 6 non-pSS patients, and 6 controls. Histochemical staining and transmission electron micrographs (TEM) were performed to macroscopically and microscopically characterize morphological features of labial salivary glands in different disease stages. Then, we performed quantitative PCR to validate hub genes. Finally, we analyzed correlations between selected hub genes and immune cells using the CIBERSORT algorithm. Results: We identified twenty-eight DEGs that were upregulated in pSS patients compared to healthy controls. These were mainly involved in immune-related pathways and infection-related pathways. According to the morphological features of minor salivary glands, severe interlobular and periductal lymphocytic infiltrates, acinar atrophy and collagen in the interstitium, nuclear shrinkage, and microscopic organelle swelling were observed with pSS disease progression. Hub genes based on above twenty-eight DEGs, including MS4A1, CD19, TCL1A, CCL19, CXCL9, CD3G, and CD3D, were selected as potential biomarkers and verified by RT-PCR. Expression of these genes was correlated with T follicular helper cells, memory B cells and M1 macrophages. Conclusion: Using transcriptome sequencing and bioinformatics analysis combined with our clinical data, we identified seven key genes that have potential value for evaluating pSS severity.