Outdoor air pollution and terminal duct lobular involution of the normal breast
Breast Cancer Research volume 22, Article number: 100 (2020)
Exposure to certain outdoor air pollutants may be associated with a higher risk of breast cancer, though potential underlying mechanisms are poorly understood. We examined whether outdoor air pollution was associated with involution of terminal duct lobular units (TDLUs), the histologic site where most cancers arise and an intermediate marker of breast cancer risk.
Pathologist-enumerated TDLUs were assessed in H&E (hematoxylin and eosin)-stained breast tissue sections from 1904 US women ages 18–75 who donated to the Susan G. Komen Tissue Bank (2009–2012). The 2009 annual fine particulate matter < 2.5 μm in diameter (PM2.5) total mass (μg/m3) at each woman’s residential address was estimated from the Environmental Protection Agency’s Downscaler Model combining Community Multiscale Air Quality (CMAQ) System modeling with air quality monitoring data. We secondarily considered CMAQ-modeled components of PM2.5 and gaseous pollutants. We used K-means clustering to identify groups of individuals with similar levels of PM2.5 components, selecting groups via cluster stability analysis. Relative rates (RRs) and 95% confidence intervals (95% CIs) for the association between air pollutants and TDLU counts were estimated from a zero-inflated negative binomial regression model adjusted for potential confounders.
PM2.5 total mass was associated with higher TDLU counts among all women (interquartile range (IQR) increase, RR = 1.06; 95% CI: 1.01–1.11). This association was evident among both premenopausal and postmenopausal women (premenopausal RR = 1.05, 95% CI: 1.00–1.11; postmenopausal RR = 1.11, 95% CI: 1.00–1.23). We identified 3 groups corresponding to clusters that varied geographically and roughly represented high, medium, and low levels of PM2.5 components relative to population mean levels. Compared to the cluster with low levels, the clusters with both high (RR = 1.74; 95% CI: 1.08–2.80) and medium (RR = 1.82; 95% CI: 1.13–2.93) levels were associated with higher TDLU counts; although not significantly different, the magnitude of the associations was stronger among postmenopausal women.
Higher PM2.5 levels were associated with reduced TDLU involution as measured by TDLU counts. Air pollution exposure may influence the histologic characteristics of normal tissue which could in turn affect breast cancer risk.
Terminal duct lobular units (TDLUs), the structures within the breast where most cancers arise , involute (shrink in number and size) with age. Involution is a complex process that is not observable; however, standardized, reproducible measures to quantify involution have been determined including TDLU count, TDLU span, and number of acini (milk producing substructures) per TDLU. Reduced involution (i.e., higher TDLU count, TDLU span, and acini/TDLU) has been associated with an increased risk of breast cancer among women with benign breast disease [2,3,4]. Certain established reproductive breast cancer risk factors have been related to TDLU count; later age (≥ 30) at first term birth and lack of breastfeeding have been associated with higher TDLU count in premenopausal women, and earlier age at menarche has been associated with higher TDLU count in postmenopausal women . Circulating sex hormones and estrogen metabolites have also been associated with reduced TDLU involution as measured by TDLU count [6, 7]. Thus, reduced TDLU involution is an intermediate marker of breast cancer risk.
Studies have suggested that outdoor air pollution exposure is related to an increased risk of breast cancer . The most consistent associations have been found for nitrogen dioxide (NO2) and nitrogen oxides (NOx) [8,9,10,11,12], traffic-derived air pollutants. Although less evidence has been found for fine particulate matter < 2.5 μm in diameter (PM2.5) in relation to breast cancer [9, 11, 13, 14], two recent, large epidemiologic studies observed increased risks associated with higher levels [10, 15]. White et al. reported that PM2.5 was associated with a higher breast cancer risk overall, for ductal carcinoma in situ, and that there was variability by geographic region and PM2.5 component profiles ; Villeneuve et al. reported an increased risk of premenopausal breast cancer with higher PM2.5 exposure . Additionally, studies have reported associations between multiple different air pollutants (PM2.5, airborne metals, polycyclic aromatic hydrocarbons) and increased mammographic density [16, 17], a strong risk factor for breast cancer . Lack of TDLU involution has been associated with higher mammographic density [19,20,21], and it has been suggested that associations between mammographic density and breast cancer may partially reflect the amount of at-risk epithelium [19, 20]. Therefore, there is support for examining the association between air pollution and TDLU involution. Addressing this question may help inform mechanisms underlying associations between outdoor air pollution and breast cancer risk; however, no epidemiologic studies to date have examined this question.
The objective of this study was to determine whether air pollution was associated with measures of TDLU involution (TDLU count, TDLU span, and acini/TDLU) in normal breast tissue samples from healthy volunteers. We considered PM2.5 total mass as the air pollutant of primary interest, but also examined associations for individual components of PM2.5, clusters of participants classified by PM2.5 component profiles, and gaseous pollutants, including NO2 and NO.
The Komen Tissue Bank (KTB) is a biorepository of breast tissue donated by healthy volunteer women. Additional details have been described elsewhere (https://komentissuebank.iu.edu/). A total of 2,197 women who enrolled in the Komen Tissue Bank between 2009 and 2012 were eligible for the original TDLU cohort . Women were excluded if they did not have available breast tissue (n = 10), had a history of cancer (n = 177), were pregnant at the time of donation (n = 18), did not have information on menopausal status (n = 37), or were older than age 75 (n = 17). Breast tissue characteristics were measured in the remaining 1,938 tissue samples. Of these, 10 of the samples were from women who donated tissue twice within a year; we only included information from the first tissue donation (n = 5), which resulted in 1,933 women with breast tissue characteristics. At the time of donation, women filled out a questionnaire that asked about demographics, medical history, lifestyle factors, and reproductive factors.
All women provided written informed consent, and the protocol was approved by the Indiana University Institutional Review Board.
The primary exposure of interest was the 2009 annual PM2.5 level at each woman’s residential address at the time of donation. The Environmental Protection Agency’s Downscaler Model is a Bayesian space-time downscaling fusion model that combines modeled air pollution levels from the Community Multiscale Air Quality (CMAQ) Modeling System with fixed site air pollution measurements from the National Air Monitoring Stations/State and Local Air Monitoring Stations . The 2009 PM2.5 concentrations were determined at a 12 × 12 km grid level and linked to each woman’s geocoded residential address. We also considered components of PM2.5 (sulfate [SO4], nitrate [NO3], ammonium [NH4], elemental carbon [EC], and organic carbon [OC]) and gaseous pollutants (carbon monoxide [CO], nitrogen oxide [NO], nitrogen dioxide [NO2], and sulfur dioxide [SO2]). Concentrations of these pollutants were only available from CMAQ modeling (version 5.0.2) , but not the more complete Downscaler Model. Therefore, the “raw” (i.e., not combined with air monitoring data) concentrations of PM2.5 components and gaseous pollutants were of secondary interest. Although daily measures of these pollutants were also available from the EPA, we considered annual averages because we hypothesized that an annual average would be more biologically relevant for an outcome such as TDLU involution that is likely to reflect cumulative exposure to risk factors over a long latency period. Prior research has shown that former smoking and certain reproductive factors that occurred years prior to breast tissue donation, such as age at menarche and age at first birth, to be associated with TDLU measures , supporting the relevance of long-term exposures.
The procedures for the collection of breast tissue and assessment of TDLU measures in the KTB population have been described in detail . Briefly, up to four breast tissue cores were collected from the upper outer quadrant of the breast with a standard 10-gauge needle. One core was fixed in 10% formalin buffer and stained with hematoxylin and eosin (H&E). A pathologist, blinded to other covariates, reviewed the digitized section to enumerate TDLUs (TDLU count) and the percent of fat on the slide (0–25%, 26–50%, 51–75%, 76–100%). Among women with a TDLU count > 0, up to 10 TDLUs were measured for TDLU span (using an electronic ruler (microns)) and acini counts per TDLU in categories (1 = ≤ 10, 2 = 11–20, 3 = 21–30, 4 = 31–50, 5 = > 50). For both TDLU span and acini/TDLU, the median value was calculated across the TDLUs assessed.
Reproducibility of the TDLU assessment was previously evaluated for a subset of 72 randomly selected images and eight duplicate images by three pathologists . Spearman correlations for intra-observer reproducibility of three reviewers were ρ = 0.99 for TDLU count, ρ = 0.68–0.96 for TDLU span, and ρ = 0.79–0.90 for acini/TDLU. Spearman correlations between the reviewers were ρ = 0.86–0.96 for TDLU count, ρ = 0.66–0.76 for TDLU span, and ρ = 0.71–0.76 for acini/TDLU. We considered TDLU count our primary outcome given the higher reproducibility and larger sample size of that measure. However, we included TDLU span and acini/TDLU as secondary outcomes given that the different measures may pertain to different stages of TDLU involution and together provide a more complete picture of the biological processes [5, 7].
To examine the association between each air pollutant and TDLU count, we used zero-inflated negative binomial regression models to estimate relative rates (RRs) and 95% confidence intervals (CIs). The zero-inflated portion of the model was used to account for the fact that 34% of women had a TDLU count of 0, which was higher than what would be expected in standard Poisson and negative-binomial regression models. Negative binomial regression accounted for the larger variance compared to the mean of TDLU count . Median TDLU span (53–266, > 266–382, > 382 units) and acini count/TDLU (1, > 1–2, > 2) were categorized into tertiles. To examine associations between each air pollutant and TDLU span or acini count/TDLU, we used ordinal logistic regression models to estimate odds ratios (ORs) and 95% CIs. We tested the proportional odds assumption with the chi-square score test; violations were not observed. Each air pollutant was examined both continuously per interquartile range (IQR) increase in pollutant concentration and categorically based on quartiles of the pollutant among the study population.
Potential confounders were selected based on the literature and a directed acyclic graph (DAG) [25, 26]; our adjustment set included age (continuous), smoking status (never, former, current), education (high school graduate or less, vocational/technical/associate’s degree, college degree, graduate/professional degree, other), race/ethnicity (non-Hispanic white, non-Hispanic Black, Asian, Hispanic, other), body mass index (< 25.0, 25.0 to < 30.0, ≥ 30.0 kg/m2), and percent of fat on the slide (0–25%, 26–50%, 51–75%, 76–100%) (Supplemental Fig. 1). We did not adjust for other TDLU risk factors that are not also related to air pollution exposure, such as number of live births or hormone replacement therapy use, that may affect precision but not control confounding. In the binary, logit portion of the zero-inflated negative binomial models to predict the zero counts, we used the same set of covariates as in the negative binomial portion of the model. This set demonstrated a better model fit than including age alone according to the Akaike information criterion (AIC). Given that conclusions were unchanged regardless of whether the exposures were also included in the zero-inflated portion of the model (data not shown), they were only included in the negative binomial, and not the zero-inflated, portion.
We excluded women whose address was not in the USA or did not provide a donation address (and therefore could not be linked to the air pollution data, n = 6) or were missing information on covariates (n = 23), leaving 1,904 (98% of those eligible) women for inclusion in our analysis.
We examined associations for TDLU count, TDLU span, and acini/TDLU for all women combined and stratified by menopausal status, using an interaction term between menopausal status and the air pollutant. Women were classified as postmenopausal (n = 574) if they reported that they had not had a menstrual period in the past 12 months, had a bilateral oophorectomy, had a hysterectomy without a bilateral oophorectomy and were ≥ 55 years old, or had a uterine ablation and were ≥ 55 years old. Women were classified as premenopausal (n = 1,322) if they reported having a period in the past 12 months, no history of bilateral oophorectomy, history of a hysterectomy without a bilateral oophorectomy and were < 55 years old, or a uterine ablation and were < 55 years old. Eight women were missing information on menopausal status; these women were included in the analyses for the overall population but were not included in the analyses stratified on menopausal status.
Given that a majority (89%) of the women in this population lived in Indiana at the time of tissue donation and that outdoor air pollution sources may vary between Indiana and the rest of the USA, we conducted a sensitivity analysis restricted to Indiana residents (TDLU count only). Individual pollutant analyses were conducted in SAS 9.4 (Cary, NC).
We used K-means clustering to identify groups of individuals with similar patterns of PM2.5 component exposure. K-means is an unsupervised clustering method to identify a pre-specified number (k) of representative joint exposure values (centroids) that can be used to classify individuals into k groups according to which centroid is closest to an individual’s exposure . We examined a plot of the proportion of within over total sums of squares by number of clusters and used cluster stability analysis to determine the smallest number of clusters with the highest Rand index . Roughly, the Rand index selects the value of k that maximizes the repeatability of clustering across subsets of the population and is seen as a measure of how robust the clustering algorithm is to variability in the population. The number of clusters identified at the bend in the plot of the within over the total sums of squares was five, but the Rand index identified three clusters as optimal. Given that the difference in the proportion of the within over total sums of squares between three and five clusters was small, we used three clusters. Identification of the clusters was completed in R 3.6.0.
We then described the clusters by comparing the mean concentrations of the PM2.5 components across the individuals in each cluster to the population means. We compared participant characteristics and geographic variation across the clusters. We examined the association between the clusters and TDLU count using the same analytic methods as described above. One of the clusters had a smaller number of women and given that the TDLU span and acini/TDLU analyses were only among those with a TDLU count > 0, we did not have sufficient power to examine the association between cluster and TDLU span and acini/TDLU.
The average age of the 1,904 women included in this study was 41.5 years (Table 1). Overall, 72% of the women were non-Hispanic white, 58% had a college degree or higher, 73% had never smoked, and 37% had a BMI ≥ 30.0 kg/m2. Compared to postmenopausal women, women who were premenopausal at tissue donation were more likely to have a college degree or higher, be never smokers, have a lower percent of fat on the slide, and have a BMI < 25.0 kg/m2. The state of residence of the women at the time of their tissue donation is provided in Supplemental Table 1.
Air pollutant characteristics
The 2009 annual mean residential PM2.5 total mass was 12.5 μg/m3 (standard deviation (SD) = 0.7 μg/m3) (Additional file 1, Supplemental Table 2). SO4 was the PM2.5 component with the highest concentration (mean = 2.4 μg/m3); CO was the gaseous pollutant with the highest concentration (mean = 223.4 ppb). The strongest correlations among the pollutants examined were between EC, NO, and NO2 (r = 0.99 and 0.98), and the weakest correlation was between NO3 and SO4 (r = 0.02) (Additional file 1, Supplemental Fig. 2).
PM2.5 total mass and TDLU involution measures
An IQR increase in PM2.5 total mass (0.5 μg/m3) was associated with higher TDLU count among all women (RR = 1.06; 95% CI: 1.01, 1.11) (Table 2). This association was evident in both premenopausal (RR = 1.05; 95% CI: 1.00, 1.11) and postmenopausal (RR = 1.11; 95% CI: 1.00, 1.23) women. Quartiles of PM2.5 above the referent were associated with an elevated TDLU count among all women, although the association was non-linear with stronger associations in quartile 2 and quartile 3 than in quartile 4 (e.g., quartile 2 vs. 1 RR = 1.17; 95% CI: 1.00, 1.38). The associations between PM2.5 total mass and TDLU count were similar when restricted to Indiana residents among all women (Additional file 1, Supplemental Table 3). In premenopausal women, the quartile-based associations were slightly stronger whereas among postmenopausal women the associations were slightly attenuated among Indiana residents. However, there was considerable overlap in the CIs between the estimates for the overall population and the Indiana residents. PM2.5 total mass was not associated with TDLU span (Additional file 1, Supplemental Table 4). A positive, non-statistically significant association with acini count/TDLU was observed for an IQR increase in PM2.5 total mass (OR = 1.06; 95% CI: 0.96, 1.17), but this association was not observed for PM2.5 categorized in quartiles and did not significantly vary by menopausal status (Additional file 1, Supplemental Table 5).
PM2.5 components and TDLU involution measures
An IQR increase in all components of PM2.5 (SO4, NO3, NH4, EC), except for OC, was associated with higher TDLU count among all women, although these associations were not statistically significant (Table 3). The 4th vs. 1st quartile of NO3 (RR = 1.28; 95% CI: 0.92, 1.77) and NH4 (RR = 1.39; 95% CI: 1.01, 1.91) were associated with higher TDLU counts only among postmenopausal women, although there was no statistically significant heterogeneity in associations by menopausal status. The associations were slightly attenuated in the sensitivity analysis restricted to Indiana residents, with the exception of stronger associations for the 2nd and 3rd quartile of EC and OC overall and in premenopausal women (Additional file 1, Supplemental Table 6).
A 4th vs. 1st quartile increase in NO3 (OR = 0.71; 95% CI: 0.52, 0.97), NH4 (OR = 0.63; 95% CI: 0.47, 0.86), and EC (OR = 0.73; 95% CI: 0.54, 0.99) and an IQR increase in EC (OR = 0.80; 95% CI: 0.67, 0.96) and OC (OR = 0.81; 95% CI: 0.67, 0.97) were inversely associated with TDLU span among all women (Additional file 1, Supplemental Table 7). PM2.5 components were not associated with acini count/TDLU (Additional file 1, Supplemental Table 8).
Gaseous pollutants and TDLU involution measures
All gaseous pollutants associated with higher TDLU count in a non-linear manner, with significant associations in the 2nd vs. 1st quartile for all women and premenopausal women (Table 3). Except for SO2, the associations were slightly stronger when restricted to Indiana residents (Additional file 1, Supplemental Table 6).
IQR increases in gaseous pollutants, except for SO4, were inversely associated with TDLU span, but not acini count/TDLU (Additional file 1, Supplemental Tables 7 and 8).
PM2.5 component clusters and TDLU involution measures
Using K-means, three clusters were identified that roughly represented high, medium, and low levels of all PM2.5 components. These clusters followed distinct geographic patterns (Fig. 1, Supplemental Table 9) and followed patterns of demographic characteristics. Cluster 1 was composed of 31 individuals who only lived in states other than Indiana; cluster 2 included 1,058 individuals who mostly lived in other parts of Indiana besides Indianapolis and also in a few other states, particularly Ohio and Kentucky; and cluster 3 consisted of 815 individuals who were mostly concentrated around Indianapolis, IN, and around urban areas of surrounding states such as Louisville, KY; Cincinnati and Columbus, OH; and Chicago, IL (Fig. 1). Average levels of PM2.5 components among cluster 1 individuals were lower than levels in the study population as a whole (Additional file 1, Supplemental Table 10). These individuals were more likely to have a graduate or professional degree, be past smokers, and have a BMI < 25.0 kg/m2 compared to the other clusters. Levels of PM2.5 components among cluster 2 individuals were similar to the overall population for SO4, NO3, and NH4, but slightly lower for EC and OC. Individuals in this cluster were more likely to be non-Hispanic white and have a high school degree or less. Levels of all PM2.5 components among cluster 3 individuals were higher than the population mean levels, and this group was more likely to be non-Hispanic Black.
Compared to cluster 1, which had lower concentrations of PM2.5 components, clusters 2 (RR = 1.82; 95% CI: 1.13, 2.93) and 3 (RR = 1.74; 95% CI: 1.08, 2.80) were associated with higher TDLU count among all women (Table 4). The magnitude of these associations was stronger among postmenopausal (cluster 3 vs. 1 RR = 3.33; 95% CI: 1.19, 9.35) compared to premenopausal women (cluster 3 vs. 1 RR = 1.56; 95% CI: 0.91, 2.68), although there was no statistically significant heterogeneity by menopausal status.
In this study of women who volunteered to donate healthy breast tissue to the Komen Tissue Bank, we found that living in areas of higher exposure to PM2.5 was associated with reduced involution of the breast as measured by higher TDLU count. Consistent with this, most of the assessed individual components of PM2.5 were non-linearly associated with higher TDLU count and clusters of individuals with levels of all PM2.5 components above or at study population mean levels had higher TDLU counts compared to the cluster with levels below the population mean. Reduced TDLU involution has been associated with an increased risk of breast cancer and may reflect higher amounts of at-risk epithelium. Therefore, our results suggest PM2.5 could impact the histologic characteristics of breast tissue and inform early carcinogenic mechanisms relating air pollution to breast cancer risk.
To our knowledge, this was the first study to examine the association between air pollution and TDLU involution measures. While there is little prior epidemiologic work with which to compare our results, biologic and indirectly related results support our overall finding that PM2.5 is associated with reduced involution of the breast. Reproductive factors, estrogen metabolites, and circulating sex hormones have all been associated with TDLU count, which suggests that involution of the breast may partially occur through a hormone-related pathway [5,6,7]. Estrogenic and antiestrogenic effects of airborne particles were reported in a study of human T47D-KBluc breast cancer cells , particulate matter < 1 μm in diameter is composed of compounds that affect estrogen-regulated pathways in vivo , and genotoxic effects of PM2.5 have been reported in mice . Polycyclic aromatic hydrocarbons and metals, components of PM2.5 that we were unable to evaluate here but may have contributed to our PM2.5 total mass finding, have been shown to be estrogenic and induce mammary tumors in animal models [32, 33]. Additionally, among breast cancer patients, higher levels of multiple pro-inflammatory markers were associated with reduced involution of healthy breast tissue . PM2.5 has been shown to increase systemic inflammation, as measured by C-reactive protein [35, 36].
Although earlier epidemiologic studies of PM2.5 and breast cancer did not report evidence of elevated risk [9, 11, 13, 14], two recent large cohort studies found PM2.5 was associated with an increased risk of breast cancer , and that composition of PM2.5 and geographic variability were important . Further, in a Breast Cancer Surveillance Consortium study of 279,967 women, PM2.5 was positively associated with mammographic breast density (i.e., heterogeneously dense compared to scattered fibroglandular breasts) , although a smaller study reported no significant association . Mammographic density is one of the strongest known breast cancer risk factors , and reduced TDLU involution has been associated with higher breast density [19, 21]. One hypothesis for the relationship between mammographic density and breast cancer risk is that it may be due in part to the amount of at-risk epithelium such as that measured by TDLU involution . Therefore, our results suggest exposure to PM2.5 may impact characteristics of healthy breast tissue, which could influence future breast cancer risk.
In addition to PM2.5 and its components, we found that the gaseous pollutants, CO, NO, NO2, and SO2, were non-linearly associated with elevated TDLU count. CO, NO, and NO2 are common traffic-related air pollutants, and NO and NO2 are the pollutants that have been most consistently associated with an increased risk of breast cancer . While biological mechanisms linking NO and NO2 to TDLU involution are not established, these pollutants may be proxies for other traffic-related air pollutants such as polycyclic aromatic hydrocarbons, which have demonstrated both antiestrogenic and estrogenic activity .
PM2.5 total mass was not consistently associated with TDLU span or acini count/TDLU. However, individual components of PM2.5, except for SO4, and gaseous pollutants were inversely associated with TDLU span but not associated with acini count/TDLU. Differences in the associations across measures may be because TDLU span and acini count/TDLU were only assessed among women who had > 0 TDLUs, so there was reduced statistical power for these analyses. The underlying biologic mechanisms and significance of the inverse associations generally observed in relation to TDLU span, in contrast to the positive associations in relation to TDLU count, are unclear. However, as suggested previously, there is a weak correlation between TDLU count and acini count/TDLU or TDLU span and associations for certain other factors also vary by TDLU measure . For example, parity and higher levels of circulating sex hormones were associated with TDLU counts but not TDLU span or acini/TDLU [5, 6]. Compared to never use of hormone replacement therapy in postmenopausal women, current use was associated with TDLU span but not TDLU count or acini/TDLU while former use was associated with TDLU count and span but not acini/TDLU . Therefore, the different markers of involution may represent distinct processes or stages of involution (i.e., complete disappearance vs reduction in size) [5, 7, 39] and more research is needed to understand how they may contribute to breast cancer etiology.
Our results are based on a large sample size of women who were demographically diverse and donated healthy breast tissue. Given the dearth of research on air pollutants in relation to breast tissue characteristics, a strength of this study was the range of pollutants considered. However, it is important to note that the correlations across certain pollutants are high and independent associations may be difficult to disentangle. Each of the gaseous pollutants, NO2, CO, and SO2, is distinct criteria pollutants with varying relative contribution from different sources and is a regulatory priority of the Environmental Protection Agency, so understanding their individual health effects is important. For the PM2.5 components, we were able to leverage these high correlations in our K-means approach to identify subgroups of women who had similar patterns of exposure. This was an important consideration because PM2.5 is a heterogeneous mixture and studies have noted that health effects may vary by the composition of the PM2.5 [10, 40, 41]. Another strength of this study was the consideration of clusters of individuals based on distinct PM2.5 profiles that had differential associations with TDLU involution. While our clusters generally separated individuals based on high, average, or low levels of all PM2.5 components rather than varying combinations of the components, these findings supported our overall PM2.5 results.
A small breast tissue core has been shown to be generally representative of involution throughout the breast , but we cannot exclude the possibility of non-representative sampling of breast tissue. We adjusted for percent of fat on the slide, which is inversely related to TDLU count , to partially account for this possibility. A strength of this study was the inclusion of multiple measures of TDLU involution in addition to a range of different air pollutants. However, this resulted in many statistical comparisons, and thus, it is possible that some findings may be due to chance. Given this, we focused the interpretation of our results on the magnitude of point estimates, precision of confidence intervals, and trends observed in the data. Women provided their residential address at the time of tissue donation, but we did not have additional information about residential history, such as how long they had lived there or information on past residences. Therefore, we could not evaluate past air pollution exposure. Further, concentrations at a 12 × 12km grid level linked to residences do not fully account for variations in an individual’s daily activities, such as where they work, that could impact their exposure. PM2.5 total mass data came from the EPA’s Downscaler Model that combines monitoring data (direct ambient measurements, but sparse geographic coverage) with CMAQ data (modeled estimates at all 12 × 12km grids across the USA, but subject to calibration/accuracy of the modeling parameters) . While this hinges on both sources’ strengths to increase the accuracy of the exposure assessment, non-differential exposure measurement error is still likely. Further, the components of PM2.5 and gaseous pollutants were available from CMAQ-modeled data alone and are less accurate than the PM2.5 total mass estimates from the Downscaler Model. One study reported that while fused models outperformed raw CMAQ data, the CMAQ predictions of PM2.5 and certain components were within recommended model performance criteria . While we examined five components of PM2.5, information on other components of PM2.5 that may be carcinogenic, such as trace metals and polycyclic aromatic hydrocarbons, was not available from CMAQ. Finally, women in the Komen Tissue Bank volunteered to donate healthy breast tissue and a majority lived in Indiana, which may limit the generalizability of our results.
Our results suggest outdoor air pollution exposure may affect characteristics of normal breast tissue through reduced TDLU involution as measured by TDLU counts. We specifically found that PM2.5 total mass, most of the individual components of PM2.5 that were considered, and gaseous pollutants, including those that are proxies of traffic-related air pollution, were non-linearly associated with higher TDLU counts. Future studies of lobular involution of the breast should continue to examine geographic variations in air pollution and consider other air pollution exposure metrics that we were not able to examine in this study. More research on the relationship between widespread environmental exposures and TDLU involution will enhance our understanding of the mechanisms underlying the role of the environment in breast carcinogenesis.
Availability of data and materials
Requests for Komen Tissue Bank data, including the data used in this manuscript, can be requested through the study website (https://komentissuebank.iu.edu/researchers/).
Russo J, Hu YF, Yang X, Russo IH. Developmental, cellular, and molecular basis of human breast cancer. J Natl Cancer Inst Monogr. 2000;27:17–37.
Milanese TR, Hartmann LC, Sellers TA, Frost MH, Vierkant RA, Maloney SD, Pankratz VS, Degnim AC, Vachon CM, Reynolds CA, et al. Age-related lobular involution and risk of breast cancer. J Natl Cancer Inst. 2006;98(22):1600–7.
Figueroa JD, Pfeiffer RM, Brinton LA, Palakal MM, Degnim AC, Radisky D, Hartmann LC, Frost MH, Stallings Mann ML, Papathomas D, et al. Standardized measures of lobular involution and subsequent breast cancer risk among women with benign breast disease: a nested case-control study. Breast Cancer Res Treat. 2016;159(1):163–72.
Baer HJ, Collins LC, Connolly JL, Colditz GA, Schnitt SJ, Tamimi RM. Lobule type and subsequent breast cancer risk: results from the Nurses’ Health Studies. Cancer. 2009;115(7):1404–11.
Figueroa JD, Pfeiffer RM, Patel DA, Linville L, Brinton LA, Gierach GL, Yang XR, Papathomas D, Visscher D, Mies C et al: Terminal duct lobular unit involution of the normal breast: implications for breast cancer etiology. J Natl Cancer Inst. 2014;106(10).
Khodr ZG, Sherman ME, Pfeiffer RM, Gierach GL, Brinton LA, Falk RT, Patel DA, Linville LM, Papathomas D, Clare SE, et al. Circulating sex hormones and terminal duct lobular unit involution of the normal breast. Cancer Epidemiol Biomarkers Prevention. 2014;23(12):2765–73.
Oh H, Khodr ZG, Sherman ME, Palakal M, Pfeiffer RM, Linville L, Geller BM, Vacek PM, Weaver DL, Chicoine RE, et al. Relation of serum estrogen metabolites with terminal duct lobular unit involution among women undergoing diagnostic image-guided breast biopsy. Hormones Cancer. 2016;7(5–6):305–15.
White AJ, Bradshaw PT, Hamra GB. Air pollution and breast cancer: a review. Current Epidemiol Reports. 2018;5(2):92–100.
Andersen ZJ, Stafoggia M, Weinmayr G, Pedersen M, Galassi C, Jorgensen JT, Oudin A, Forsberg B, Olsson D, Oftedal B, et al. Long-term exposure to ambient air pollution and incidence of postmenopausal breast cancer in 15 European cohorts within the ESCAPE project. Environ Health Perspect. 2017;125(10):107005.
White AJ, Keller JP, Zhao S, Carroll R, Kaufman JD, Sandler DP. Air pollution, clustering of particulate matter components, and breast cancer in the sister study: a U.S.-wide cohort. Environ Health Perspectives. 2019;127(10):107002.
Reding KW, Young MT, Szpiro AA, Han CJ, DeRoo LA, Weinberg C, Kaufman JD, Sandler DP. Breast cancer risk in relation to ambient air pollution exposure at residences in the sister study cohort. Cancer Epidemiol Biomarkers Prevention. 2015;24(12):1907–9.
Goldberg MS, Labreche F, Weichenthal S, Lavigne E, Valois MF, Hatzopoulou M, Van Ryswyk K, Shekarrizfard M, Villeneuve PJ, Crouse D, et al. The association between the incidence of postmenopausal breast cancer and concentrations at street-level of nitrogen dioxide and ultrafine particles. Environ Res. 2017;158:7–15.
Andersen ZJ, Ravnskjaer L, Andersen KK, Loft S, Brandt J, Becker T, Ketzel M, Hertel O, Lynge E, Brauner EV. Long-term exposure to fine particulate matter and breast cancer incidence in the Danish nurse cohort study. Cancer Epidemiol Biomarkers Prevention. 2017;26(3):428–30.
Hart JE, Bertrand KA, DuPre N, James P, Vieira VM, Tamimi RM, Laden F. Long-term particulate matter exposures during adulthood and risk of breast cancer incidence in the nurses’ health study II prospective cohort. Cancer Epidemiol Biomarkers Prevention. 2016;25(8):1274–6.
Villeneuve PJ, Goldberg MS, Crouse DL, To T, Weichenthal SA, Wall C, Miller AB. Residential exposure to fine particulate matter air pollution and incident breast cancer in a cohort of Canadian women. Environ Epidemiol. 2018;2(3):e021.
Yaghjyan L, Arao R, Brokamp C, O'Meara ES, Sprague BL, Ghita G, Ryan P. Association between air pollution and mammographic breast density in the Breast Cancer Surveilance Consortium. Breast Cancer Research. 2017;19(1):36.
White AJ, Weinberg CR, O'Meara ES, Sandler DP, Sprague BL. Airborne metals and polycyclic aromatic hydrocarbons in relation to mammographic breast density. Breast Cancer Research. 2019;21(1):24.
McCormack VA, dos Santos SI. Breast density and parenchymal patterns as markers of breast cancer risk: a meta-analysis. Cancer Epidemiol Biomarkers Prevention. 2006;15(6):1159–69.
Gierach GL, Patel DA, Pfeiffer RM, Figueroa JD, Linville L, Papathomas D, Johnson JM, Chicoine RE, Herschorn SD, Shepherd JA, et al. Relationship of terminal duct lobular unit involution of the breast with area and volume mammographic densities. Cancer prevention research (Philadelphia, Pa). 2016;9(2):149–58.
Sung H, Guo C, Li E, Li J, Pfeiffer RM, Guida JL, Cora R, Hu N, Deng J, Figueroa JD, et al. The relationship between terminal duct lobular unit features and mammographic density among Chinese breast cancer patients. Int J Cancer. 2019;145(1):70–7.
Ghosh K, Hartmann LC, Reynolds C, Visscher DW, Brandt KR, Vierkant RA, Scott CG, Radisky DC, Sellers TA, Pankratz VS, et al. Association between mammographic density and age-related lobular involution of the breast. J Clin Oncol. 2010;28(13):2207–12.
Environmental Protection Agency. Technical Information abous Fused Air Quality Surface Using Downscaling Tool: Metadata Description. Research Triangle Park, NC; 2016. https://www.epa.gov/sites/production/files/2015-09/documents/dsmetadataair_0612_0.pdf.
Environmental Protection Agency. CMAQ: The Community Multiscale Air Quality Modeling System [https://www.epa.gov/cmaq]. Accessed Dec 2019.
Yang S, Harlow L, Puggioni G, Redding C. A comparison of different methods of zero-inflated data analysis and an application in health surveys. J Mod Appl Stat Methods. 2017;16(1):518–43.
Greenland S, Pearl J, Robins JM. Causal diagrams for epidemiologic research. Epidemiology (Cambridge, Mass). 1999;10(1):37–48.
Shrier I, Platt RW. Reducing bias through directed acyclic graphs. BMC Med Res Methodol. 2008;8:70.
Gibson EA, Nunez Y, Abuawad A, Zota AR, Renzetti S, Devick KL, Gennings C, Goldsmith J, Coull BA, Kioumourtzoglou MA. An overview of methods to address distinct research questions on environmental mixtures: an application to persistent organic pollutants and leukocyte telomere length. Environ Health. 2019;18(1):76.
Hubert L, Arabie P. Comparing partitions. J Classif. 1985;2(1):193–218.
Chen ST, Lin CC, Liu YS, Lin C, Hung PT, Jao CW, Lin PH. Airborne particulate collected from central Taiwan induces DNA strand breaks, poly (ADP-ribose) polymerase-1 activation, and estrogen-disrupting activity in human breast carcinoma cell lines. J Environ Sci Health A Tox Hazard Subst Environ Eng. 2013;48(2):173–81.
Wenger D, Gerecke AC, Heeb NV, Schmid P, Hueglin C, Naegeli H, Zenobi R. In vitro estrogenicity of ambient particulate matter: contribution of hydroxylated polycyclic aromatic hydrocarbons. J Applied Toxicol. 2009;29(3):223–32.
de Oliveira AAF, de Oliveira TF, Dias MF, Medeiros MHG, Di Mascio P, Veras M, Lemos M, Marcourakis T, Saldiva PHN, Loureiro APM. Genotoxic and epigenotoxic effects in mice exposed to concentrated ambient fine particulate matter (PM (2.5)) from São Paulo city, Brazil. Part Fibre Toxicol. 2018;15(1):40.
el-Bayoumy K, Chae YH, Upadhyaya P, Rivenson A, Kurtzke C, Reddy B, Hecht SS: Comparative tumorigenicity of benzo [a] pyrene, 1-nitropyrene and 2-amino-1-methyl-6-phenylimidazo [4,5-b] pyridine administered by gavage to female CD rats. Carcinogenesis 1995, 16(2):431–434.
Byrne C, Divekar SD, Storchan GB, Parodi DA, Martin MB. Metals and breast cancer. J Mammary Gland Biol Neoplasia. 2013;18(1):63–73.
Hanna M, Dumas I, Orain M, Jacob S, Têtu B, Sanschagrin F, Bureau A, Poirier B, Diorio C. Association between local inflammation and breast tissue age-related lobular involution among premenopausal and postmenopausal breast cancer patients. PLoS One. 2017;12(8):e0183579.
Liu Q, Gu X, Deng F, Mu L, Baccarelli AA, Guo X, Wu S. Ambient particulate air pollution and circulating C-reactive protein level: a systematic review and meta-analysis. Int J Hyg Environ Health. 2019;222(5):756–64.
Ostro B, Malig B, Broadwin R, Basu R, Gold EB, Bromberger JT, Derby C, Feinstein S, Greendale GA, Jackson EA, et al. Chronic PM2.5 exposure and inflammation: determining sensitive subgroups in mid-life women. Environ Res. 2014;132:168–75.
DuPre NC, Hart JE, Bertrand KA, Kraft P, Laden F, Tamimi RM. Residential particulate matter and distance to roadways in relation to mammographic density: results from the Nurses’ Health Studies. Breast Cancer Research. 2017;19(1):124.
Sievers CK, Shanle EK, Bradfield CA, Xu W. Differential action of monohydroxylated polycyclic aromatic hydrocarbons with estrogen receptors alpha and beta. Toxicol Sci. 2013;132(2):359–67.
Oh H, Bodelon C, Palakal M, Chatterjee N, Sherman ME, Linville L, Geller BM, Vacek PM, Weaver DL, Chicoine RE, et al. Ages at menarche- and menopause-related genetic variants in relation to terminal duct lobular unit involution in normal breast tissue. Breast Cancer Res Treat. 2016;158(2):341–50.
Krall JR, Mulholland JA, Russell AG, Balachandran S, Winquist A, Tolbert PE, Waller LA, Sarnat SE. Associations between source-specific fine particulate matter and emergency department visits for respiratory disease in four U.S. cities. Environ Health Perspect. 2017;125(1):97–103.
Peng RD, Bell ML, Geyh AS, McDermott A, Zeger SL, Samet JM, Dominici F. Emergency admissions for cardiovascular and respiratory diseases and the chemical composition of fine particle air pollution. Environ Health Perspect. 2009;117(6):957–63.
Vierkant RA, Hartmann LC, Pankratz VS, Anderson SS, Radisky D, Frost MH, Vachon CM, Ghosh K, Distad TJ, Degnim AC, et al. Lobular involution: localized phenomenon or field effect? Breast Cancer Res Treat. 2009;117(1):193–6.
Chen G, Li J, Ying Q, Sherman S, Perkins N, Sundaram R, Mendola P. Evaluation of observation-fused regional air quality model results for population air pollution exposure estimation. Sci Total Environ. 2014;485-486:563–74.
Data from the Susan G. Komen Tissue Bank at the IU Simon Cancer Center were used in this study. We thank contributors, including Indiana University who collected data used in this study, as well as donors and their families, whose help and participation made this work possible. We also thank Bob McConnaughey for performing the geocoding of the residential addresses and linkage to the EPA air pollution data.
This work was supported by projects Z01-ES044005 and Z1AES103332-01 of the Intramural Research Program of the NIH, the National Institute of Environmental Health Sciences, and the National Cancer Institute.
Ethics approval and consent to participate
All women provided written informed consent, and the protocol was approved by the Indiana University Institutional Review Board.
Consent for publication
The authors declare they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
. State of residence at time of tissue donation. Supplemental Table 2. Mean concentrations and quartile cut-points for PM2.5 total massa, PM2.5 componentsb, and gaseous pollutantsc. Supplemental Table 3. Associations between PM2.5 total massa and terminal duct lobular unit (TDLU) counts restricted to Indiana residents. Supplemental Table 4. Associations between PM2.5 total massa and TDLU span. Supplemental Table 5. Associations between PM2.5 total massa and acini count/TDLU. Supplemental Table 6. Associations between PM2.5 components and gaseous pollutantsa and terminal duct lobular unit (TDLU) counts restricted to Indiana residents. Supplemental Table 7. Associations between PM2.5 components and gaseous pollutantsa and TDLU span. Supplemental Table 8. Associations between PM2.5 components and gaseous pollutantsa and acini count/TDLU. Supplemental Table 9. State of residence at time of tissue donation by PM2.5 component cluster. Supplemental Table 10. Participant characteristics by PM2.5 component cluster. Supplemental Figure 1. Directed acyclic graph for the relationship between air pollution and TDLU measures. Supplemental Figure 2. Spearman correlations between PM2.5 total mass, PM2.5 components, and gaseous pollutants.
About this article
Cite this article
Niehoff, N.M., Keil, A.P., Jones, R.R. et al. Outdoor air pollution and terminal duct lobular involution of the normal breast. Breast Cancer Res 22, 100 (2020). https://doi.org/10.1186/s13058-020-01339-x