Causal associations of air pollution with rheumatoid arthritis: A transethnic Mendelian randomization study

Abstract

Background: Rheumatoid arthritis is a common rheumatic disease, and its onset is closely related to genetic and environmental factors, however, the relationship between air pollution and RA is still hotly debated. Further investigation of the relationship between air pollution and rheumatoid arthritis is conducive to a comprehensive understanding of the risk factors of the disease, providing certain value for the clinical prevention and treatment of RA.

Methods: We used a Two-Sample Mendelian Randomization approach, integrating the large-scale public genomewide association study, to assess the genetically predicted causal effect of air pollution (including: PM2.5, PM2.5-10, PM10, nitrogen dioxide, nitrogen oxides) on RA in European and European East Asian populations, respectively. Indicators related to air pollution (2,505 individuals to 423,796 individuals), including European and East Asian populations were obtained from the Integrative Epidemiology Unit open GWAS project. Published East Asian RA data were also obtained from the IEU open GWAS project (212,453 individuals), while large-scale publicly available European RA data were obtained from finngen R10 (13,621 cases and 262,844 controls). Inverse variance weighting was used as the primary analytical method, complemented by MR-egger, Weighed median, and Weighted mode results. Cochran Q tested for heterogeneity, and MR-Egger regression analyses were performed to test for multiplicity. leave-one-out analysis allowed for the robustness and reliability were assessed.

Results: No statistically significant effects of PM2.5, PM2.5-10, PM10, nitrogen dioxide, nitrogen oxides and RA were observed in either European or East Asian populations. Results from European data: PM2.5 (IVW OR: 0.71; 95% CI: 0.27-1.91; p = 0.498; number of SNPs: 5), PM2.5-10 (IVW OR: 1.20; 95% CI: 0.61-2.40; p = 0.596; number of SNPs: 15), PM10 (IVW OR: 1.69; 95% CI: 0.84-3.39; p = 0.142; number of SNPs: 9), nitrogen dioxide (IVW OR: 3.88; 95% CI: 0.19-77.77; p = 0.375; number of SNPs: 2), nitrogen oxides (IVW OR: 0.51; 95% CI: 0.16-1.67; p = 0.268; number of SNPs: 4). East Asian data results: PM2.5 (IVW OR: 1.16; 95% CI: 0.98-1.38; p = 0.086; number of SNPs: 4), PM2.5-10 (IVW OR: 1.14; 95% CI: 0.95-1.38; p = 0.166; number of SNPs: 2), PM10 (IVW OR: 0.95; 95% CI: 0.81-1.11; p = 0.503; number of SNPs: 3), nitrogen dioxide (IVW OR: 0.87; 95% CI: 0.76-1.00; p = 0.051; number of SNPs: 6), nitrogen oxides (IVW OR: 0.96; 95% CI: 0.82-1.14; p = 0.671; number of SNPs: 3). No signs of pleiotropy or heterogeneity were observed in the MR-Egger intercept, MR-PRESSO and Cochrane's Q (p>0.05). In addition, no outliers were found in the MR-PRESSO analysis. The results were further validated by leave-one-out tests, confirming the robustness of the findings.

Conclusions: We performed transethnic MR analysis suggesting that there may not be a genetically predicted causal relationship between air pollution and RA.

Fig 1. Flowchart of Mendelian randomization analyses conducted in this study.

LD, linkage disequilibrium; r², linkage disequilibrium coefficient; kb, linkage disequilibrium block size; SNPs, Single nucleotide polymorphisms; Assumption 1, Correlation hypothesis: There is a strong correlation between SNPs and exposure factors; Assumption 2, Independence assumption: SNPs and confounding factors are independent; Assumption 3, Exclusion hypothesis: SNPs can only affect outcome through exposure factors; MR, Mendelian randomization; IVW, Inverse variance weighted; MR-PRESSO, MR Pleiotropy RESidual Sum and Outlier.

Fig 2. The forest plot displays the associations of genetically predicted levels of PM2.5, PM2.5–10, PM10, nitrogen dioxide, and nitrogen oxide with the risk of rheumatoid arthritis.

SNP: Single-nucleotide polymorphism; CI: confidence interval; OR: odds ratio. P-value: P-value of the IVW method.