Identification of Novel Pleiotropic SNPs Associated with Osteoporosis and Rheumatoid Arthritis.

Genome-wide association studies (GWASs) have identified hundreds of genetic loci for osteoporosis (OP) and rheumatoid arthritis (RA), individually, however, a large proportion of the total trait heritability remains unexplained. Previous studies demonstrated that these two diseases may share some common genetic determination and risk factors, but they were generally focused on individual trait and failed to identify the common variants that play key functional roles in the etiology of these two diseases. Here, we performed a conditional false discovery rate (cFDR) analysis to identify novel pleiotropic variants shared between them by integrating two independent GWASs with summary statistics for total body bone mineral density (TB-BMD, a major risk factor for osteoporosis) (n = 66,628) and RA (n = 58,284). A fine-mapping approach was also applied to identify the most probable causal variants with biological effects on both TB-BMD and RA. As a result, we found 47 independent pleiotropic SNPs shared between TB-BMD and RA, and 40 of them were validated in heel ultrasound estimated BMD (eBMD), femoral neck BMD (FN-BMD) or lumbar spine (LS-BMD). We detected one SNP (rs13299616) was novel and not identified by previous BMD or RA-related studies. Combined with fine-mapping and GWAS-eQTL colocalization analyses, our results suggested that locus 1p13.2 (including PTPN22, MAGI3, PHTF1, and RSBN1) was an important region to regulate TB-BMD and RA simultaneously. These findings provide new insights into the shared biological mechanisms and functional genetic determinants between OP and RA, and novel potential targets for treatment development.

66,628 European ancestry individuals with bone mineral density measurement, 58,284 European ancestry individuals with rheumatoid arthritis