Background: Sjogren’s disease (SjD) is an autoimmune disease characterized by reduced function of exocrine glands (i.e., salivary and lacrimal glands). Epithelial cell damage resulting from lymphocytic infiltration has been implicated in SjD etiology [1]. How genetic and epigenetic changes influence epithelial-immune cell interactions in SjD pathogenesis remain understudied.

Objectives: Evaluate the role of SjD risk loci in salivary gland tissue to gain insights into the potential genes involved in salivary gland dysfunction.

Methods: SNPs from 16 regions with SNP-SjD associations (P<5x10-8) in our GWAS study (3232 SjD cases) and meta-analysis of ImmunoChip data (619 SjD cases) [2] were interrogated for eQTLs using Genotype-Tissue Expression (GTEx) minor salivary gland data. Subsequent analysis identified genes that were both eQTLs in the minor salivary gland and significantly expressed in RNA-seq and ATAC-seq data from the submaxillary salivary gland epithelial cell line, A253. Pathway enrichment analysis was performed using gProfiler on the genes where coalescence of eQTL, RNA-seq, and ATAC-seq data was observed. To further validate the results, we performed transcriptome-wide association study (TWAS) analysis using GWAS summary statistics and minor salivary gland eQTL GTEx data.

Results: In total, 5884 genome-wide significant SNPs from 16 SjD risk loci were identified as potential minor salivary gland eQTLs using two discovery thresholds: p(FDR)<0.05 provided by eQTL study (3566 SNPs) and p(FDR)>0.05 and p<0.05 in eQTL study (2318 SNPs). Further analysis revealed 10 SjD risk loci with SNPs that were minor salivary gland eQTLs for a total of 155 unique genes that had a coalescence of RNA- and ATAC-seq data in A253 cells. Many SNPs altered the expression of the nearest gene to the risk allele (i.e., index gene), such as IRF5 and TNPO3 on chromosome 7 at 128Mb; however, this locus had 12 additional genes that were eQTLs in minor salivary gland. In contrast, other loci had no reported eQTLs for the index gene, but several reported eQTLs for other genes, such TYK2 on chromosome 19 at 10Mb that showed no change in TYK2 expression but eQTLs for 8 distant genes, including ICAM1 . Pathway enrichment analysis revealed an enrichment in Butyrophilin (BTN) family interactions (R-HSA-8851) (PAdj=1.564x10-5), including the BTN2A1 , BTN2A2 , BTN3A1 , BTN3A2 and BTN3A3 gene cluster in the MHC region. In further support, TWAS of the minor salivary gland and the SjD GWAS summary statistics (after Bonferroni correction) showed association between SjD and BTN3A2 (p=1.24x10-42), as well as many other loci in the MHC region. In addition, several long non-coding (lnc) RNAs on chromosome 17 were significant, peaking at RP11-259G18.1 (p=4.43x10-10).

Conclusion: This study shows that SjD-associated risk alleles influence disease by altering gene expression in immune cells and minor salivary glands. Further, our analysis suggests that altered gene expression in the minor salivary gland expands beyond effects on the index gene to several genes on each locus. Interestingly, we observed minor salivary gland eQTLs for several BTN family genes, which act as cell-surface binding partners to regulate cell-cell interactions, including interactions between epithelial cells and activated T cells [3]. Future work will assess chromatin-chromatin-interactions within the 10 SjD risk loci in salivary gland cells and tissues to map local chromatin regulatory networks that regulate gene expression. Additional transcriptional studies of SjD minor salivary gland tissues will provide further insights into how altered gene expression in the salivary gland influences SjD pathology.

REFERENCES:

[1]Verstappen. Nat Rev Rheumatol 2021;17(6):333-348.

[2]Khatri, et al. Annals of Rheumatic Diseases 2020;79:30-31.

[3]Arnett HA, Viney JL. Nature Reviews Immunology 2014;14:559-569.

Disclosure of Interests: Farhang Aghakhanian: None declared, Mandi M Wiley: None declared, Bhuwan Khatri: None declared, Kandice L Tessneer: None declared, Astrid Rasmussen: None declared, Simon J. Bowman Consultant of: Abbvie, Galapagos, and Novartis in 2020-2021., Lida Radfar: None declared, Roald Omdal: None declared, Marie Wahren-Herlenius: None declared, Blake M Warner: None declared, Torsten Witte: None declared, Roland Jonsson: None declared, Maureen Rischmueller: None declared, Patrick M Gaffney: None declared, Judith A. James: None declared, Lars Ronnblom: None declared, R Hal Scofield: None declared, Xavier Mariette: None declared, Marta Alarcon-Riquelme: None declared, Wan Fai Ng: None declared, Kathy Sivils Employee of: Current employee of Janssen, Gunnel Nordmark: None declared, Umesh Deshmukh: None declared, A Darise Farris: None declared, Christopher Lessard: None declared