Systemic autoimmune rheumatic diseases and multiple industrial air pollutant emissions: A large general population Canadian cohort analysis

Past investigations of air pollution and systemic autoimmune rheumatic diseases (SARDs) typically focused on individual (not mixed) and overall environmental emissions. We assessed mixtures of industrial emissions of fine particulate matter (PM2.5), nitrogen dioxide (NO2), and sulfur dioxide (SO2) and SARDs onset in Ontario, Canada. We assembled an open cohort of over 12 million adults (without SARD diagnoses at cohort entry) based on provincial health data for 2007–2020 and followed them until SARD onset, death, emigration, or end of study (December 2020). SARDs were identified using physician billing and hospitalization diagnostic codes for systemic lupus, scleroderma, myositis, undifferentiated connective tissue disease, and Sjogren’s. Rheumatoid arthritis and vasculitis were not included. Average PM2.5, NO2, and SO2 industrial emissions from 2002 to one year before SARDs onset or end of study were assigned using residential postal codes. A quantile g-computation model for time to SARD onset was developed for the industrial emission mixture, adjusting for sex, age, income, rurality index, chronic obstructive pulmonary disease (as a proxy for smoking), background (environmental overall) PM2.5, and calendar year. We conducted stratified analyses across age, sex, and rurality. We identified 43,931 new SARD diagnoses across 143,799,564 person-years. The adjusted hazard ratio for SARD onset for an increase in all emissions by one decile was 1.018 (95% confidence interval 1.013–1.022). Similar positive associations between SARDs and the mixed emissions were observed in most stratified analyses. Industrial PM2.5 contributed most to SARD risk. Industrial air pollution emissions were associated with SARDs risk.