Association between nitrogen dioxide and incident breast cancer in Alberta’s tomorrow project

TL;DR

This study examines if nitrogen dioxide pollution is linked to post-menopausal breast cancer in Alberta, but finds no significant association.

Contribution

The study provides new evidence on the relationship between traffic-related air pollution and breast cancer risk in a Canadian population.

Findings

No significant association was found between nitrogen dioxide exposure and post-menopausal breast cancer risk.

Hazard ratios ranged from 1.01 to 1.10 across different models, indicating a risk magnitude similar to other studies.

The lack of significant association may be due to sample size limitations and lower pollution levels in Alberta.

Abstract

Breast cancer is a significant global health concern, with millions of cases diagnosed annually. Risk factors for breast cancer include obesity, physical inactivity, poor diet, and alcohol use. Urbanization introduces exposures such as traffic-related air pollution (TRAP), which may impact breast cancer risk. To assess the impact of TRAP, measured in this study using nitrogen dioxide (NO2), on the risk of post-menopausal breast cancer we conducted a secondary analysis of baseline and follow-up data from the Alberta Tomorrow Project (ATP) cohort, linked to geographic data at the postal code level from the Canadian Urban Environmental Health Research Consortium (CANUE). Cox proportional hazard regression models were built to assess the association between breast cancer risk and NO2. Direct acyclic graphs (DAG) were used to identify and select relevant confounders for adjustment. The study included 15,536 post-menopausal women and 523 incident breast cancer cases over 12.6 years of follow-up. The hazard ratio (HR) for the initial model, which included all potential confounders (age, education, number of births, fruits and vegetables consumption, income, physical activity, BMI, and smoking) was 1.01 per 10-ppb increase in NO2. The adjusted HR for a 10-ppb increase in NO2 measured at baseline (the year participants completed the health and behaviours survey) was 1.10 (95% CI = 0.90, 1.34). Our study found no significant association between NO2 exposure and post-menopausal breast cancer risk. The variation in reported hazard ratios (HRs) — ranging from 1.01 to 1.10 — reflects different models adjusting for different variables. Our findings suggest a risk magnitude similar to other studies, yet the lack of significant association may be due to sample size limitations and lower NO2 levels in Alberta compared to other regions globally.

The online version contains supplementary material available at 10.1038/s41598-025-04373-x.

Figures

PRISMA diagram of study population. Abbreviations: ATP: Alberta’s Tomorrow Project.a.

Full text

Introduction

Breast cancer is a significant disease among women, with approximately 2.3 million cases diagnosed globally in 2020, resulting in 685,000 deaths^1^. In Canada, about 84 women receive a breast cancer diagnosis every day^2^. Established risk factors for post-menopausal breast cancer, which account for around two thirds of all cases, include obesity, physical activity (PA), dietary intake, smoking, family history of breast cancer, and alcohol^3–8^.

A large proportion of breast cancer cases occur in highly developed countries, often associated with urbanization and economic growth^9^. Currently, about half of the global population reside in urban areas, a figure projected to increase to more than 90% by 2100^10^. Living in urban areas offers many benefits, including better healthcare access, but it may also expose individuals to elevated environmental pollution levels^11^, primarily from air pollution and stationary-source fuel combustion^12^. Traffic-related air pollution (TRAP), a mixture of various compounds, such as carbon dioxide (CO_2_), nitrogen dioxide (NO_2_), sulfur dioxide (SO_2_), ozone (O_3_), and particulate matter (PM)^13^, has been identified as a Group 1 carcinogen by the International Agency for Research on Cancer (IARC)^14^. Recent studies suggest a potential causal link between these compounds, specifically NO_2_, and female breast cancer, especially in post-menopausal women^15–18^.

Although there has been extensive research on the association between NO_2_ and other chronic conditions, such as lung cancer and cardiovascular diseases, its association with breast cancer is less well studied. While some studies reported positive and significant associations between breast cancer and NO_2_ exposure, others have found non-significant associations. In a large US-based cohort study, no overall increase in breast cancer risk per interquartile range difference (5.8 parts per billion - ppb) in NO_2_ exposure was observed (HR = 1.02, 95% CI: 0.97, 1.07)^18^. In the European ESCAPE cohort study, which combined data from 15 prospective European cohorts, there was a suggestive positive association between exposure to most ambient air pollutants and post-menopausal breast cancer^19^. However, the association between ambient NO_2_ and breast cancer did not reach statistical significance, with a HR of 1.02 per 10 µg/m^3^^19^. Two large cohort studies, in the Catalonia region of Spain and in Canada, reported a modest non-significant increase of breast cancer incidence with greater NO_2_ exposure. In the Canadian National Breast Screening Study, while no association was observed between breast cancer risk and NO_2_ in the post-menopausal cohorts, in pre-menopausal women, using 50- and 52-years to identify pre- and post-menopausal women, a 9.7 ppb increase in NO_2_ exposure (per interquartile range increase) was linked to rate ratios (RR) of 1.13 and 1.17^16^. The Catalonia study, on the other hand, reported a HR of 1.05 (95% CI, 1.02, 1.08) per interquartile range (27.6 µg/m^3^) increase in NO_2_ exposure in post-menopausal women^20^.

In line with these mixed findings, significant uncertainty remains regarding a causal association between air pollution, particularly NO_2_, and post-menopausal breast cancer. Further research is needed to fully elucidate the mechanisms linking NO_2_ to breast cancer.

The objective of the present study was to assess the association between NO_2_ exposure and post-menopausal breast cancer risk in Alberta’s Tomorrow Project (ATP), a longitudinal cohort study that enrolled participants aged 35 to 69 residing in the province of Alberta, Canada. Given the age range of participants in ATP the majority of breast cancer cases were post-menopausal, and the analysis was restricted to the post-menopausal women given the distinct differences between pre- and post-menopausal breast cancer^3^.

Methods

Study population

The ATP is a longitudinal cohort study that enrolled participants through two phases of recruitment. Phase I, conducted between 2000 and 2008, employed a two-stage telephone random digit dialing (RDD) approach to recruit participants. In this phase, the RDD was mapped to Alberta Regional Health Authorities (RHA) to select households with eligible residents. One or two residents from each household were then invited for participation in the study. Phase II, spanning from 2009 to 2015, enrolled participants through two methods: first, by extending invitations to individuals already involved in Phase I to join Phase II, and second, via volunteer sampling. The recruitment process for Phase II did not employ RHA-based RDD. Between 2000 and 2015, ATP recruited 55,000 men and women between the ages of 35 and 69 years, from towns and cities across the province, with the goal of long-term follow-up over 50 years. Baseline information for individuals recruited during Phase I was collected through the Health and Life Questionnaire (HLQ), Canadian Dietary History Questionnaire (CDHQ), and Physical Activity Questionnaire (PAQ). For participants recruited during Phase II, the same baseline data were collected using an Updated Health and Life Questionnaire (UHLQ) and Physical Activity and Nutrition Survey (PANS). Starting in January 2011, the UHLQ and PANS were replaced by the Canadian Partnership of Tomorrow Project (CanPath) CORE Questionnaire. Consequently, participants who joined the study in 2011 or later exclusively completed the CORE Questionnaire. Our study focused on post-menopausal women. Women were considered post-menopausal either through their response to the ATP questionnaire (options included: yes – natural menopause, yes – other reasons, no, don’t know) or their age at enrollment. In the ATP questionnaire, menopausal status was defined as whether an individual’s menstrual period stopped for one year and did not restart. Further, in cases where the participant answered “no” and were over the age of 55, the age cut-off of 60 was used. In cases where the participant answered “don’t know” or did not answer the question, then the age cut-off of 55 years was used to ascertain menopausal status^21^. This study received ethical approval from the University of British Columbia Behavioral Research Ethics Board. We confirm that all methods were performed in accordance with the relevant guidelines and regulations.

Outcome classification

In the present study, cases were defined as post-menopausal individuals who experienced incident, in situ, or invasive breast cancer, and were censored at the year of diagnosis. Incident cases were identified through record linkage to the Alberta Cancer Registry, and data on breast cancer diagnosis, such as tumor staging and other pathological characteristics, were not extracted for this study. Post-menopausal participants who had been diagnosed with breast cancer prior to their enrollment in the study were excluded. Non-cases were post-menopausal women who did not develop breast cancer during the follow-up period (follow-up end date: June 30, 2022), with outcomes including mortality, loss to follow-up, or reaching the end of the follow-up period.

Exposure assessment for NO2

The analysis focused on NO_2_ as a marker or proxy of TRAP. We accessed this information through the Canadian Urban Environmental Health Research Consortium (CANUE) dataset, which compiles and generates standardized environmental data at the dissemination area level related to different built environment characteristics. Following data compilation, CANUE facilitates the integration of the data into existing Canadian cohort studies and administrative health databases, including the ATP database. CANUE data was linked to ATP based on participants’ postal code of residence at enrolment. The initial land-use regression (LUR) annual average concentration at postal code in 2006, which represented annual average concentration at baseline, were provided for all participants in parts per billion (ppb).

Analytic population

The ATP recruited both men and women; however, the study’s dataset included only women, as men were excluded based on the applied exclusion criteria before the data was provided. This aligns with the study’s focus on post-menopausal breast cancer in women. Individuals were considered eligible for inclusion if they were either post-menopausal at the time of enrollment or transitioned into the post-menopausal stage during the follow-up period. These menopausal status transitions were identified through the questionnaires that participants completed or based on their age at baseline. In addition, participants with follow-up time of less than one year were excluded from analysis. With this exclusion criteria, in some cases, incident breast cancers diagnosed within one year of enrolment were also eliminated from the study. In instances where participants had two breast cancer diagnoses at different time points, only first diagnosis was retained.

Statistical analysis

In this study, multivariable Cox proportional-hazards regression models were used. We built a model for NO_2_ measured at baseline (the year participant completed the UHLQ or CORE survey) – and estimated HRs and 95% confidence intervals for the association between incident breast cancer and 10 ppb change in NO_2_. Follow-up time was used as the time scale for the Cox model, calculated as the difference between the year of enrollment and the year of diagnosis for cases, or between the year of enrollment and the end of follow-up in 2022. Schoenfeld’s tests assessed Cox proportional-hazards assumptions, while a 5% change-in-estimate threshold guided the selection of potential confounders for final models, alongside a directed acyclic graph (DAG) to identify a minimally sufficient adjustment set.

The covariate adjustment set included body mass index (BMI) (underweight - BMI less than 18.5 kg/m^2, normal weight - BMI between 18.5 and 24.9 kg/m^2, overweight - BMI between 25 and 29.9 kg/m^2, and obese - BMI greater than or equal to 30 kg/m^2), physical activity (met or did not meet World Health Organization PA recommendations), ethnicity (white or other), vegetable and fruit consumption (categories of 0, 1, 2, 3, 4, 5 or more), social and material deprivation (categorized into quintile of 1 – least deprived to 5 – most deprived), smoking (never smokers, daily smoker, and occasional smoker), income (less than 50 K, 50–100 K, 100–150 K, more than 150 K), alcohol consumption (never, less than once a month, once to 3 times a month, once per week, 2 to 5 times a week, and more than 6 times a week), education (category A - no formal education, elementary school, or high school; category B - included trade/technical or vocational school, a diploma from a community college, or a university certificate below bachelor’s level; category C comprised a bachelor’s degree or a graduate degree), contraceptive use (never or ever), birth (2 or less, 3–4, greater than 4), maternal breast cancer history (yes or no), and age (less than 44, 45–54, 55–64, more than 65). Social deprivation represents the deprivation of relationships within family, workplace, and community – the index is calculated based on the proportion of the population who are separated, divorced, or widowed; the proportion living alone; and the proportion who have relocated in the past five years^22^. Material deprivation represents the deprivation of goods and amenities – its index is based on average household income, unemployment rate, and the rate of high school education^22^.

First, we constructed the initial model (model I), which included all covariates. Next, based on a 5% cut-off for change-in-estimate procedure, variables including age, education, number of births, fruit and vegetable consumption, income, physical activity, BMI, and smoking met the criteria for inclusion (model II). However, guided by the DAG and existing literature (see supplementary Figure S1), only social deprivation, material deprivation, income, education, and age were classified as preliminary confounders. However, in model III, age, education, and social deprivation were retained. Social deprivation was included instead of income and material deprivation due to redundancy in the effects captured by these variables. We further constructed three models to assess the association between NO_2_ quartiles and postmenopausal breast cancer. Model I included all covariates. In the change-in-estimate model (model II), only age, number of births, and material deprivation model were retained. In the DAG-informed model (model III), the same confounders as previously identified (age, education, and social deprivation) were included, with social deprivation retained over material deprivation and income due to redundancy in their effects.

Physical activity had the highest missing value rate at 69%, while education and contraceptive use had the lowest at less than 1%. Most other variables had missing value rates below 15%, except for ethnicity, maternal history of breast cancer, and smoking, which all had rates around 30–40%. Notably, age was the only variable devoid of any missing values. To handle high missingness, missing data was categorized into a “not reported” category.

There were 24,889 women who met all the specified criteria for this study. Among them, 12 participants had missing menopausal status data. From the 24,899 participants, 1,942 individuals had prevalent cancers at study baseline, excluding nonmelanoma skin cancer, which resulted in a study population of 22,947 after exclusion. An additional 7,411 participants who were premenopausal throughout the follow-up period were excluded from the analysis. This retained a final study population of 15,536 participants, see Fig. 1. During the follow-up period, there were 523 incident breast cancer diagnoses, accounting for 3.4% of the participants. The remaining 15,013 individuals did not receive a breast cancer diagnosis and were classified as non-cases. For the analysis of baseline NO_2_ concentration a further 485 observations (including 9 breast cancer cases) were excluded due to missing NO_2_ data. The proportional hazard assumption for all models was met. In this study, NO_2_ concentration measurements were utilized as a proxy for TRAP due to data availability.

Fig. 1PRISMA diagram of study population. Abbreviations: ATP: Alberta’s Tomorrow Project.a.

Results

Descriptive statistics

The study population included 21% of individuals under the age of 44 who were post-menopausal, with 12% of cases in this age group identified as post-menopausal. In addition, 29% were aged 45–54, 34% aged 55–64, and 17% were over 65. Among the 17,001 participants who completed the CORE questionnaire, 94% identified as White. Participants were followed for an average of 12.6 years (SD = 4.2), with cases followed for 7.9 years on average and non-cases for 12.8 years. The mean age at baseline for the study population was 53.5 years, with a standard deviation of 9.6 years. Detailed descriptive statistics among cases and non-cases are presented in Table 1. Variables that were statistically significantly different in the two groups included age, ethnicity, maternal history of breast cancer, physical activity, and BMI.

Table 1. Baseline participant characteristics.VariablesOverallCasesNon-casesPN or mean% or SDN or mean% or SDN or mean% or SDNumber of participants15,536100%5233.40%15,01396.60%Follow-up time (years)12.64.27.94.912.84.1Age at baseline53.59.656.98.753.49.6< 0.05Alcohol consumptionNever11277.7397.810887.60.6Less than once per month296420.110721.4285720.11–3 times per month365224.811422.8353824.9Once per week162311479.41576112–5 times per week422928.715130.1407828.7More than 6 times per week11327.7438.510897.7Unknown/Missing8095.2224.27875.2BirthTwo or less433432.314632.7418832.30.193–4789358.826158.4763258.8Greater than 411968.9408.911568.9Unknown/Missing211313.67614.5203713.6RaceWhite10,14094.931259.7982894.90.002Other5405.1122.35285.1Unknown/Missing485631.319938465731EducationHigh school or below369723.813125.1356623.80.2Diploma or trade556535.819837.9536735.8University-level education626640.419337607340.4Unknown/Missing------------Contraceptive useEver13,8258946789.313,358890.87Never1704115610.7164811Unknown/Missing------------Deprivation score - socialQ1 (least deprived)269818.29218.9260618.20.25Q22809199719.9271219Q32956209519.5286120Q4333722.611523.6322222.5Q5 (most deprived)298520.28818.1289720.3Unknown/Missing7514.8366.97154.8Deprivation score – materialQ1(least deprived)432129.213327.3418829.30.31Q2357224.212425.5344824.1Q3283219.29519.5273719.1Q4251016.98317242717Q5 (most deprived)155010.55210.7149810.5Unknown/Missing7514.8366.97154.8Fruit consumption, servings/dayZero5783.72145573.70.6One313720.39518.2304220.4Two51073319036.4491732.9Three360623.312123.2348523.3Four191512.45811.1185712.4Five or more11197.3377.110827.3Unknown/Missing------------Vegetable consumption, servings/dayZero------------0.33One186212.17414.2178812Two421027.212824.6408227.3Three391825.414127.1377725.3Four284218.49518.2274718.4Five or more245015.88115.5236915.9Unknown/Missing------------IncomeLess than 50,000304621.212325.5292321.10.09$50,000-$100,000488934.116935472034$100,000-$150,000312521.89920.5302621.8Greater than $150,000329022.99219319823.1Unknown/Missing11867.6407.611467.6Maternal history of breast cancerYes402341.117146.7385241< 0.05No575958.919553.3556459Unknown/Missing57543715730559737.3Physical activityMet51410.81512.649910.7< 0.05Not Met426389.210487.4415989.3Unknown/Missing10,75969.340477.210,35569SmokingNever767487.227989.1739587.20.14Occasionally------------Daily8709.9309.68409.9Unknown/Missing673843.421040.2652843.5BMIUnderweight------------< 0.05Healthy596941.616836.2580141.8Overweight469732.715232.7454532.7Obese352324.614230.5338124.4Unknown/Missing11887.65811.111307.5Abbreviations: BMI: body mass index; p: p-value; PA: physical activity; MET: metabolic equivalent time; SD: standard deviation; For all variables, proportions were calculated excluding missing values from the denominator, and proportions of missing values were calculated using the total sample size as the denominator. --: participant count < 10 suppressed.

Distribution of NO2 exposure estimated by LUR

The mean concentration of baseline NO_2_ was 10.1 ppb, with a standard deviation (SD) of 4.4. The highest and lowest NO_2_ concentrations were 34.8 ppb and 2.1 ppb, respectively (see supplementary Figure S2). The average baseline NO_2_ concentration for cases and non-cases were 10.3 ppb (SD = 4.4) and 10.1 ppb (SD = 4.4), respectively.

NO2 exposure (continuous) at baseline and breast cancer

The HR for a 10-ppb increase in NO_2_ was 1.14 (95% CI = 0.90, 1.45) in model I. Further adjustments for breast cancer risk factors (model II) yielded an HR of 1.01 (95% CI = 0.82, 1.24) for a 10-ppb increase in NO_2_. In the final model the HR for a 10 ppb increases in NO_2_ was 1.10 (95% CI = 0.90, 1.34) (Table 2). This model had a concordance value of 0.62. Based on the Schoenfeld’s test, the proportional hazard assumption was met. The global p-value was 0.08, and the p-value for NO_2_ was 0.20. There was a total of 15,051 observations across all models, of which included 514 breast cancer cases.

Table 2. Association for breast cancer risk with baseline NO_2_ (continuous).Regression modelHR95% CIP-valueModel I1.140.90–1.450.30Model II1.010.82–1.240.90Model III1.100.90–1.340.36Abbreviations: NO_2_: nitrogen dioxide; HR: hazard ratio; CI: confidence interval; *Adjusted for all covariates; **Adjusted for age, education, number of births, fruit and vegetable consumption, income, physical activity, BMI, and smoking; ***Adjusted for education, social deprivation, and age.

NO2 exposure (quartiles) at baseline and breast cancer

The HR for the second quartile (Q2) were slightly above 1 in all models (ranging from 1.12 to 1.16) with p-values and 95% CIs indicating statistical non-significance (Table 3). For the third quartile (Q3) there was a borderline significant association in model I (HR = 1.27, 95% CI: 0.97–1.66), which further strengthened in model II (HR = 1.29, 95% CI: 0.99–1.67), and became statistically significant in model III (HR = 1.31, 95% CI: 1.01–1.68). Compared to Q1, no association was observed for the highest quartile (Q4) of NO_2_, with HRs close to 1 across all models. The global p-value for NO_2_ was 0.07.

Table 3. Association for breast cancer risk with baseline NO_2_ (quartile).NO_2_ QuartileModel IHR (95% CI), p-valueModel IIHR (95% CI), p-valueModel IIIHR (95% CI), p-valueQ1 (lowest)1.001.001.00Q21.12 (0.85–1.45), 0.411.14 (0.88–1.48), 0.311.16 (0.90–1.50), 0.25Q31.27 (0.97–1.66), 0.081.29 (0.99–1.67), 0.061.31 (1.01–1.68), 0.04Q40.96 (0.70–1.30), 0.790.94 (0.71–1.24), 0.650.99 (0.75–1.29), 0.91Abbreviations: NO_2_: nitrogen dioxide; HR: hazard ratio; CI: confidence interval; *Adjusted for all covariates; **Adjusted for age, number of births, material deprivation; ***Adjusted for education, social deprivation, and age.

Discussion

Findings of the present study and summary of literature

In this study, we assessed the association between NO_2_ and incident post-menopausal breast cancer. The adjusted HR for baseline NO_2_ was 1.10 (95% CI = 0.90, 1.34) for every 10-ppb increase in NO_2_. In our assessment of the functional form of NO_2_, we found that the highest quartile did not show an increased risk which may suggest a non-linear relationship, or there may be other factors that influence risk at higher exposure. However, when a restricted cubic spline model was fitted, we found no significant improvement in the model fit compared to the linear model. This suggests that the data does not support a nonlinear relationship between breast cancer risk and NO_2_ exposure. A significant association between increased NO_2_ and breast cancer risk has been reported in several studies, mostly nested case-control studies, or cross-sectional analyses. While some studies have reported positive and significant associations^15,20,23–25^, others have found positive but non-significant associations, consistent with our results^17,19,26,27^. Our findings are in line with two other studies. A Canadian case-control study found a positive but non-significant association between NO_2_ exposure and incident breast cancer among post-menopausal women, with an odds ratio (OR) of 1.07 (95% CI: 0.86, 1.32) for every 10-ppb increase^26^. Similarly, a study that assessed long-term exposure to air pollution and breast cancer risk among European cohorts, with 12 years of follow-up and 3,600 women developing breast cancer from the total of 75,000 post-menopausal women included in the study population, found a positive but non-significant association^19^. The HR was 1.02 per 10 µg/m3 increase (95% CI: 0.98, 1.07)^19^. The magnitude of effect in these two studies was similar to the present findings.

A US-based cohort study using LUR modeling, a method that estimates air pollution levels based on geographic and environmental data such as traffic density and land use, assessed NO_2_ exposure and observed positive associations with breast cancer risk, though these associations were not significant for all participants. However, in stratified analysis by distance to major roads, a positive and significant association was found for participants (n = 17,075) who lived within 500 m of major roads (combination of < 200 m and 200-<500 m categories) with an HR of 1.26 [95% CI: 1.0–1.59]^24^. In addition, a Canadian case-control study in Montreal found a significant association between air pollution and incident breast cancer, with an OR of 1.31 (95% CI: 1.00, 1.71) per 5 ppb increase in NO_2_^15^. Their study population consisted of 383 incident invasive breast cancer cases and 416 controls^15^. It is important to note the difference in sample size and exposure levels of TRAP between our study and other studies where significant association between TRAP and breast cancer risk have been reported. In our study, the mean TRAP level at baseline was 10.1 ppb. The US-based study that identified significant associations between NO_x_ and NO_2_ and breast cancer risk, reported an average NO_x_ level reported was 70 ppb in 1993, 57 ppb in 2003, and 35 ppb in 2010^24^. The absence of statistically significant associations in this study may be attributed to very low TRAP levels in Alberta^28^.

The study findings underscore the importance of considering TRAP, specifically NO_2_, in breast cancer risk management. Given the consistent positive associations observed, policymakers should recognize the potential role of NO_2_ on breast cancer incidence, while acknowledging the limitations of the current findings. Incorporating such evidence into policy considerations may involve revisiting air quality standards and urban planning regulations to mitigate NO_2_ exposure. Mechanistic studies exploring the biological pathways linking NO_2_ exposure to breast cancer development could further inform targeted policies and interventions. While the current study primarily reported no significant association, the implications for public health are substantial when considering evidence from other studies^15,19,20,24,26^,warranting a proactive approach to reducing NO_2_ exposure to potentially decrease breast cancer risk. Policymakers should consider these findings as part of a broader strategy aimed at creating healthier environments and minimizing potential environmental contributors to breast cancer incidence.

Strengths and limitations

This study had several strengths and limitations. Unlike many previous studies where the association between TRAP and risk of breast cancer were investigated via case-control or cross-sectional studies, our study benefited from a longitudinal cohort study design. This design enabled the identification of incident breast cancer over an extended period which is more likely to capture the latency period of cancer. Furthermore, a key asset of our study was the integration of ATP data with environmental exposure information, precisely linked to each participant’s residential postal code at baseline. Moreover, the incident cancer cases were linked to the Alberta Cancer Registry, ensuring the accuracy and validity of case measurements.

In this study, we were unable to comprehensively assess participants’ total environmental exposure because NO_2_ was only measured at participants’ residence at baseline. Exposure at workplaces, schools, or places of recreation where participants may potentially spend a significant portion of time, could not be accounted for. It is possible that a participant lives in an area with low TRAP but spends time in an environment with high TRAP, and vice versa. Unfortunately, these dynamics could not be accounted for and could potentially introduce bias into our findings. However, this exposure misclassification is likely to be non-differential, which can lead to a bias towards the null, attenuating our results. Another potential limitation was our inability to accurately classify menopausal status of participants. The questionnaire that included menopause status was only administered once at the baseline and was not reassessed during the follow-up period. However, any misclassification related to menopausal status is also likely to be non-differential and may bias our results towards the null hypothesis. The potential impact of missing data, especially for important risk factors like maternal history of breast cancer, may introduce bias and limit the interpretation of our results. In addition, the small sample size further restricts the statistical power of the study, which may introduce bias and limit our ability to detect significant associations. Another limitation of this study was the lack of data on tumor characteristics, such as hormone receptor status (ER/PR).

Conclusion

The magnitude and direction of association for NO_2_ and breast cancer risk in our study was similar to findings established from other studies, although our results were not statistically significant. Considering the wide variation in findings from the current literature, it is premature to definitively conclude a causal association between NO_2_ and breast cancer. Our result necessitate replication in larger populations, particularly populations from diverse climate conditions and variations in air pollution levels, with a specific focus on levels similar to those observed in Alberta.

Electronic supplementary material

Below is the link to the electronic supplementary material.

Supplementary Material 1