Genome-wide identification of cell type-specific susceptibility genes for Juvenile dermatomyositis through the analysis of N6-methyladenosine-associated SNPs

Genome-wide association studies (GWASs) have pinpointed genetic loci associated with juvenile dermatomyositis (JDM). Functional genes within the GWAS loci may be cell type-specific, but their identity remains largely unknown. N6-methyladenosine (m6A) plays a pivotal role in regulating various cellular processes and is linked to autoimmune diseases. This study aimed to underscore the potential functional genes within the GWAS loci through the analysis of m6A-associated SNPs (m6A-SNPs), specifically within relevant cell types. JDM-associated m6A-SNPs were identified from the GWAS summary dataset. The correlation between m6A-SNPs and gene expression was assessed through bulk tissue and single-cell eQTL analyses. To further investigate the relationship between gene expression and JDM, Mendelian randomization analysis was employed. Additionally, differential expression analyses were conducted on bulk tissues, as well as single-cell transcriptomic data comprising 6 JDM patients and 11 juvenile controls (99,396 cells). Seven m6A-SNPs associated with JDM were identified. Bulk tissue analysis revealed differential expression of HLA-DPA1, HLA-DPB1, MICB, HLA-A, HLA-F, HLA-DQB2, HLA-DRB5, TAP2, PSMB9, MICA, AIF1, and DDX39B influenced by m6A-SNPs, all showing associations with JDM in both differential expression and Mendelian randomization analyses. In single-cell analysis, the six m6A-SNPs within the HLA locus acted as cell-type-specific eQTLs, correlating with the expression of HLA-A, HLA-B, HLA-C, HLA-DPB1, HLA-DQA1, HLA-DQB1 and HLA-DRB1 in myeloid, T or B cells. Notably, these genes displayed abnormal expression in T, B, and myeloid cells of JDM patients. The present study identified m6A-SNPs within JDM susceptibility genes, shedding light on the intricate interplay between m6A-SNPs, gene expression, and JDM.