Transcriptome-wide association study identifies new susceptibility genes and pathways for spondyloarthritis

Abstract

Background: Spondyloarthritis (SpA) is a group of multifactorial bone diseases influenced by genetic factors, the environment and lifestyle. However, current studies have found a limited number of SpA-related genes, and the genetic and pathogenic mechanism of SpA is still unclear. Methods: A tissue-specific transcriptome-wide association study (TWAS) of SpA was performed by using genome-wide association study (GWAS, including 3966 SpA patients and 452264 controls) summary data and gene expression weights of whole blood and skeletal muscle. The SpA-associated genes identified by TWAS were further compared with the differentially expressed genes (DEGs) identified in the SpA gene expression profile acquired from the Gene Expression Omnibus database (GEO, GSE58667). Finally, FUMA and Metascape were used to perform functional enrichment and annotation analysis.Results: The TWAS detected 28 significant genes associated with SpA in whole blood and skeletal muscle, such as CTNNAL1 (PSM=0.0304, PWB=0.0096). The gene expression profile of SpA identified 20 candidate genes that overlapped in the TWAS data, such as MCM4 (PTWAS=0.0132, PDEG=0.0275) and KIAA1109 (PTWAS=0.0371, PDEG =0.0467). Enrichment analysis of the genes identified by TWAS identified 93 significant GO terms 33 and KEGG pathways, such as mitochondrion organization (GO:0007005) and axon guidance (hsa04360).Conclusion: We identified multiple candidate genes that were genetically related to SpA. Our study may provide novel clues of the genetic mechanism, diagnosis, and treatment of SpA.