- Research article
- Open Access
- Open Peer Review
Outdoor air pollution, green space, and cancer incidence in Saxony: a semi-individual cohort study
BMC Public Healthvolume 18, Article number: 715 (2018)
There are a few epidemiological studies that (1) link increased ambient air pollution (AP) with an increase in lung cancer incidence rates and (2) investigate whether residing in green spaces could be protective against cancer. However, it is completely unclear whether other forms of cancer are also affected by AP and if residential green spaces could lower cancer incidence rates in general. Therefore, the objective was to estimate whether AP and green space are associated with several cancer types.
The analysis was based on routine health care data from around 1.9 million people from Saxony who were free of cancer in 2008 and 2009. Incident cancer cases (2010–2014) of mouth and throat, skin (non-melanoma skin cancer - NMSC), prostate, breast, and colorectum were defined as: (1) one inpatient diagnosis, or (2) two outpatient diagnoses in two different quarters within one year and a specific treatment or death within two quarters after the diagnosis. Exposures, derived from freely available 3rd party data, included particulate matter with aerodynamic diameter of less than 10 μm (PM10) and nitrogen dioxide (N02) as well as green space (Normalized Difference Vegetation Index - NDVI). Associations between air pollutants, green space, and cancer incidence were assessed by multilevel Poisson models. Age, sex, physician contacts, short- and long-term unemployment, population density, and having an alcohol-related disorder were considered as potential confounders.
Three thousand one hundred seven people developed mouth and throat cancer, 33,178 NMSC, 9611 prostate cancer, 9577 breast cancer, and 11,975 colorectal cancer during the follow-up period (2010–2014). An increase in PM10 of 10 μg/m3 was associated with a 53% increase in relative risk (RR) of mouth and throat cancer and a 52% increase in RR of NMSC. Prostate and breast cancer were modestly associated with PM10 with an increase in RR of 23 and 19%, respectively. The associations with N02 were in the same direction as PM10 but the effect estimates were much lower (7–24%). A 10% increase in NDVI was most protective of mouth and throat cancer (− 11% RR) and of NMSC (− 16% RR). Colorectal cancer was not affected by any of the exposures.
In addition to the studies carried out so far, this study was able to provide evidence that higher ambient AP levels increase the risk of mouth and throat cancer as well as of NMSC and that a higher residential green space level might have a protective effect for NMSC in areas with low to moderate UV intensity. Nevertheless, we cannot rule out residual confounding by socioeconomic or smoking status.
Outdoor air pollution (AP) is ubiquitous with its exposure having effects on a large proportion of the world population . There is strong evidence from experimental and epidemiological studies that AP such as particulate matter (PM), nitrogen dioxide (NO2), and ozone (O3) are major risk factors for cardiovascular and cardiopulmonary diseases, and potentially cancer [2,3,4,5,6,7,8,9]. While the trigger function of AP on cardiac and pulmonary events is considered as causal, the role in cancer onset is only suggestive. In the case of PM and cardiovascular mortality, there is a solid association that fulfils both a temporal and a close exposure-response relationship. There is coherence of results between several scientific disciplines, including experimental studies offering plausible biological mechanisms reviewed in . Both long-term and short-term effects are involved. An increase in PM2.5 long-term exposure per 10 μg/m3 increased the risk for cardiovascular mortality by 11% in a recent meta-analysis . Even a daily rise in mean PM2.5 level per 10 μg/m3 increased the risk for cardiovascular mortality by approximately 0.4 to 1.0% . In the case of cancer, it is much more difficult to prove causality. One reason is the long latency between exposure to carcinogens and the development of cancer. Another reason is that doses are typically low in the environment, and therefore direct causal inferences are hardly realized. Nevertheless, there is suggestive epidemiological evidence that outdoor AP increases incidence rates of some cancer types with strongest evidence for lung cancer [11,12,13,14,15,16,17,18,19]. Molecular epidemiological studies [20,21,22] have shown that the biological mechanisms causing cancer from outdoor AP involve genotoxic effects of the chemical compounds that accumulate over time, including PAH-DNA adducts, chromosome aberrations, sister chromatid exchanges, ras oncogene overexpression, and radically induced (oxidative) DNA damage. PM is listed as Group I carcinogen by the International Agency for Research on Cancer (IARC). Sources of PM, other than household, are numerous including traffic, agricultural and industrial emissions as main anthropogenic origins [23, 24]. We use PM10 (particles with aerodynamic diameter of less than 10 μm) instead of PM2.5, because we have more variation in PM10 data. PM2.5 is far more evenly distributed throughout Saxony. We use NO2as marker for traffic-related AP to investigate the effect of urban agglomeration in Saxony on cancer incidence [25, 26]. We also want to narrow the research gap on the impact of green spaces on health, especially cancer. While numerous studies have analyzed associations between AP and lung cancer, few studies have addressed associations between green space and cancer incidence [27,28,29,30]. Green space is known to be inversely related to AP due to lack of emission sources in green places . Additionally, people residing in greener places might be more motivated to conduct physical activity , which is suggested to decrease cancer risk  through reductions of adipose tissue volume and endocrine activity (e.g. sex hormones) [33, 34]. On the other hand, some studies have also reported adverse associations with green space that may be context-specific, like a negative effect on the human skin through increased sunlight exposure, while spending time in green space . Since a handful of studies have reported a link between lower all-cause mortality among people residing in green spaces [28, 35,36,37], it is a valid assumption that residential green space might decrease risk of some cancer types.
The aim of the study was to analyze the association between residential exposure to PM10 and NO2 as well as green space on different cancer types. We restricted the analysis to those cancer types, which are less strongly associated with smoking as lung cancer is, but which might be affected by AP exposure either by inhalation or dermal route of exposure to PM10 particles and NO2. Further, we needed a sufficient number of cases to ensure that we have enough power for our analyses. Therefore, we had to use cancers that are frequent in the population. We identified mouth and throat cancer, non-melanoma skin cancer (NMSC) and colorectal cancer in both sexes, as well as prostate cancer in men and breast cancer in women, which fulfilled these requirements. At the beginning of the study, we also wanted to examine the effects of pesticides on various forms of cancer in addition to AP and green spaces. According to the literature [15, 16], glandular tissues like the prostate or the breast may be susceptible to some agricultural used pesticides (hormone disruptive agents). Unfortunately, the data situation was not sufficient and therefore this investigation was not continued. Nevertheless, prostate and breast cancer were still used as outcomes.
AOK PLUS study population and case definitions
In Germany, approximately 90% of the population is covered by statutory health insurances. We used routine health care data from AOK PLUS, a large statutory health insurance in Saxony (area ~ 18,000 km2, population ~ 4 Mio), which covers almost half of the local general population. The data include information from inpatient and outpatient care with respect to diagnosis, procedures and prescriptions as well as socio-demographic information of the insured population such as age, sex and residential district (first four digits of the 5-digit postal code of the residential address). Age distribution and sex-ratio of the AOK PLUS beneficiaries in Saxony are comparable to the Germany-wide population . We used data for the years 2007–2014 from the outpatient as well as from the inpatient sector. All beneficiaries were allocated to 186 four digit postal code districts based on their residential address in 2007. We excluded all cancer cases from the year 2007 for the determination of prevalent cases, because we had only outpatient data for this year. Further, beneficiaries that were diagnosed (outpatient or inpatient) in the years 2008 and 2009 with one of the analyzed cancer types were excluded, to estimate incidence rates for the years 2010 until 2014 (see Additional file 1: Figure S3). Age groups in increments of 10 years were built based on the age of every beneficiary in the year 2012. Therefore, the newborns of the years 2013 and 2014 had to be excluded as well. Further, the ages from 0 to 49 years were collapsed to one group. See Table 1 for an exact partitioning of ages.
In accordance with the Good epidemiological practice for secondary data analyses , cancer diseases were defined as one diagnosis of the corresponding ICD-10 code (International Classification of Diseases) for the inpatient data, while two diagnoses in different quarters within one year were necessary for all outpatient data plus a prescription of a specific treatment (e.g. radiotherapy, cytostatic medication) or death within two quarters after the second diagnosis (see Additional file 2: Table S1, Additional file 3: Table S2, Additional file 4: Table S3). We used the codes C00-C14 for mouth and throat cancer, C44 or L57.0 for NMSC, C61 for prostate cancer, C50 for breast cancer, and C18-C21 for colorectal cancer.
Ethics approval and consent to participate
The present analysis is based on secondary data from the health insurance company AOK PLUS, which were collected for the purpose of billing medical services. The study is supported by the AOK PLUS with which a data use and transfer agreement exists. The data was available for us only in anonymous form, so that no conclusions could be drawn about the individuals. Personal data of participants were anonymized through AOK PLUS before data sharing. Personal identifiers were masked or deleted (clear name to pseudonym; no social security number provided). Quasi-identifiers were generalized (only year of birth used; dropping of last digit of the zip code). Unfortunately, the data are not publicly accessible. There was no influence whatsoever on the policyholders and no intervention was carried out. It is therefore a purely observational study. According to paragraph 75 SGB X (Zehntes Buch Sozialgesetzbuch - German federal law) it is not reasonable to get consent for data sharing and analysis from around 2 million people, like in our investigation. Therefore, we submitted and get granted an application to the Saxon State Ministry for Social Affairs and Consumer Protection for obtaining consent to data transmission and analysis on behalf of the insured.
This study was conducted in accordance with the Helsinki Declaration  and follows the principles of Good Epidemiological Practice and Good Practice in Secondary Data Analysis . The study was also registered in the database “Versorgungsforschung Deutschland” under the number VfD_ECo_epi_16_003770.
Air pollution and green space exposure assessment
Annual NO2 and PM10 concentrations for the year 2007 were derived from freely available maps in resolution of 100 m developed for Western Europe . These maps were created by land use regression models based on more than 1500 EuroAirnet monitoring sites. Predictor variables for land use regression models included land use characteristics, population density, road length, altitude, distance to sea, and satellite-derived NO2 and PM10 data .
Green spaces were defined by the Normalized Difference Vegetation Index (NDVI), which was derived from freely available MODerate-resolution Imaging Spectroradiometer (MODIS) satellite images at the resolution of 250 m . Briefly, NDVI is a commonly used indicator of vegetation level, ranging from − 1 (water) to + 1 (absolutely vegetated area). The algorithm for NDVI is based on two vegetation-informative bands: near-infrared (841 nm to 876 nm) and visible red (from 620 nm to 670 nm). For current analysis, we averaged 115 16-day composite NDVI images for the years 2005 to 2009 [45, 46].
These PM10, NO2, and NDVI estimates for 386 five digit postal codes were weighted by the population number from the German census 2011 and then averaged to 186 four digit postal code districts, because the address information were only available for the four digit postal codes due to privacy regulations. Freely available postal code vector data were obtained from postleitzahl.org.
Geographic data management and calculations were conducted using the ArcGIS 10.1 Geographical Information System (GIS) (ESRI, Redlands, CA, USA) software program.
As control variables for the individual data analysis, age, sex, alcohol-related disorder, absolute number of physician contacts in the four digit postal code districts (AOK data) and the proportion of short- (up to one year) and long-term unemployment (more than 1 year; statistical office Saxony), ranging from 0 to 1, were considered. Additionally, for each four digit postal code district, population density per km2, all-cause mortality and proportion of persons with an alcohol-related disorder were considered in regression analysis on aggregated data.
Alcohol-related disorder was defined as one inpatient F10 diagnosis without F10.0 (i. e. acute intoxication) or three F10 prescriptions within four quarters of a year in the outpatient sector and was used to adjust the models to correct for its influence on cancer incidence [47, 48].
All the analyses were conducted with the program R, version 3.3.2 (Vienna, Austria) R Core Team , but data preprocessing was done with the software Stata (StataCorp. 2013. Stata Statistical Software: Release 13. College Station, TX: StataCorp LP). The individual associations between residential air pollutants, green space and cancer incidence rates were assessed by multilevel Poisson models with 95% confidence intervals with the R software package lme4. Due to the fact that not all individuals were fully insured or alive in the five year-long study period (2010–2014), we used the exact observation time in days as offset in the models. Because of high correlations (Pearson correlations above |0.75|) between PM10, NO2, and NDVI, associations with each exposure were analyzed individually (see Additional file 5: Figure S1).
For the aggregated data, cancer incidence rates of the 186 four digit postal codes were age-standardized by the European standard population  and adjusted for measurement errors due to the limited observation period using conditional autoregressive models (CAR), which use the first order spatial dependencies (shared borders) with the software BayesX for R [51, 52]. These associations were analyzed with linear models.
We used the best subset approach according to the lowest Akaike Information Criterion (AIC) for variable selection, implemented in the R package glmulti , to filter out unassociated variables for each model (adjusted models). We computed crude and adjusted models for each outcome and exposure pair (Additional file 6: Table S5).
We conducted a sensitivity analysis for the influence of the change of residence on the effect estimates by exclusion of movers and by comparing these results with the primary analysis where both, movers and non-movers, were included (Additional file 7: Table S4).
Baseline characteristics of the study population
Of more than 1.9 million initially cancer-free persons included into the analysis (Table 1), 3107 people developed mouth and throat cancer, 33,178 NMSC, 9611 prostate cancer, 9577 breast cancer, and 11,975 colorectal cancer during the follow-up period (2010–2014). Although there were slightly more females in the total population (53%), more males were affected by mouth and throat (74 vs. 26%) and colorectal cancer (53 vs. 47%), respectively. The number of incident cancer cases increased with age as expected, but decreased for the oldest age-groups (80–89, 90+). Mean age for developing cancer of interest ranged from 65 to 76 years. In the mouth and throat cancer subgroup, the proportion of people with an alcohol-related disorder reached almost 36%. Altogether, 69,722 people had a prevalent alcohol-related disorder (~ 3.6% from total) and 302,818 people (~ 15.8% from total) changed their place of residence between different four digit postal code districts at least once within the study period.
Air pollution and green space
PM10 and NO2 levels were differentially distributed over the 186 four digit postal code districts in Saxony (Fig. 1). While NO2 showed highest concentrations in agglomeration areas with high population density (cities Dresden, Leipzig, Chemnitz, and Zwickau), PM10 was much smoother distributed, even in rural areas, with highest concentrations in the north-west part of Saxony around Leipzig and in Middle Saxony to Dresden. Lowest concentrations of PM10 were observed in the Erzgebirge region in the south-west of Saxony. On the contrary, green space, as defined by NDVI, was higher in rural areas and the lowest in the metropolitan areas of Leipzig, Dresden, Chemnitz, and Görlitz and their vicinities. PM10 concentrations ranged from 15.47 to 26.30 μg/m3 (mean: 20.89), NO2 concentrations ranged from 9.32 to 31.55 μg/m3 (mean: 20.44), and NDVI varied between 0.38 and 0.64 (mean: 0.51). The admissible annual averages for PM10 and NO2 in accordance with the Federal Immission Control Act (39. BImSchV) and the European Directive on Air Quality (2008/50/EG) are 40 μg/m3 each.
Additional file 5: Figure S1 shows Pearson correlation coefficients for associations between AP and NDVI, as well as physician contacts (PC). Most correlations were at least 75%, except for NDVI and PC (~ − 58%) or NO2 and PC (~ 61%).
Association between air pollution, green space, and cancer incidence
Semi-individual data modeling
An increase in PM10 of 10 μg/m3 was associated with a 53% increase in relative risk (RR) of mouth and throat cancer and a 52% increase in RR of NMSC (Table 2). Prostate and breast cancer were modestly associated with PM10 with an increase in RR of 23 and 19%, respectively. The associations with N02 were in the same direction as PM10 but the effect estimates were much weaker (7–24%) and did not reach statistical significance for prostate cancer. An increase in NDVI by 10% revealed associations with mouth and throat cancer with an 11% decrease in RR and with NMSC with a 16% decrease in RR. No associations were found for prostate or breast cancer and NDVI. Colorectal cancer was not affected by any of the exposures (Table 2).
Alcohol-related disorder was associated with a more than 9 times increased risk of mouth and throat cancer a 50% increase of colorectal cancer and a 22% increase of breast cancer. Men were also more affected by any of the cancers which occur in both sexes (RR 1.6 to 2.7×; Table 2).
When only the people who did not change the place of residence were considered, similar results were observed except for the effect of PM10 on mouth and throat cancer which increased by 8% (Additional file 7: Table S4).
Aggregated data analysis
Cancer incidence rates generally showed a scattered distribution over Saxony (Fig. 1). The range of differences between the four digit postal code districts was small for colorectal cancer (86 to 94 cases per 100.000 insured persons) and breast cancer (144 to 158 cases per 100.000 insured persons), but high for mouth and throat cancer (22 to 36 cases), prostate cancer (140 to 240 cases), and especially for NMSC (150 to 400 cases). High correlation of physician contacts with PM10 (~ 75%) prevented using this variable of personal demand as a confounder in the adjusted models.
Crude linear regression analyses showed associations between PM10 and mouth and throat cancer (R2 0.203) or NMSC (R2 0.144). Associations with breast- or prostate cancer were weak (Fig. 2). NO2 was associated with NMSC (R2 0.164), but asssociations with breast-, and prostate cancer, or colorectal cancer were weak (Additional file 8: Figure S2 and Additional file 6: Table S5 in supplement). Moderate negative correlations were observed for NDVI and mouth and throat cancer (R2 0.114) or NMSC (R2 0.137), but again, the association with prostate cancer was weak (Additional file 8: Figure S2 and Additional file 6: Table S5).
Adjusted linear regression models were in the same direction as the crude analysis. Also PM10 and mouth and throat cancer (R2 0.331) or NMSC (R2 0.258) showed associations, but effect estimates were lower compared to the crude analyses (Additional file 6: Table S5). NO2showed also some association with NMSC (R2 0.259) in the adjusted analysis, while associations to prostate-, breast-, and colorectal cancer remained weak. Negative correlations between NDVI and mouth and throat cancer (R2 0.275) or NMSC (R2 0.250) were still present in the adjusted analysis, but again, effect estimates were lower than in the crude models.
To the best of our knowledge, this study is novel and substantially adds to the current knowledge, because first of all the modeling is based on a very large cohort (1.9 million beneficiaries), which is needed to study associations between environmental exposures and cancer incidence. While earlier work has exclusively focused on the association between AP and lung cancer [1, 12,13,14,15,16,17,18,19, 54,55,56], we observed that PM10 is also related to mouth and throat cancer, NMSC, and maybe also to prostate and breast cancer. It is not so surprising that other types of cancer such as mouth and throat cancer or NMSC are also associated with AP, since chemical compounds are inhaled from the tidal air or even passed through the skin into the human body and accumulate there in the respective tissue over time. Some of these compounds could induce different forms of DNA damage [20,21,22], which could lead to tumor development. On the other hand, no association was found with colorectal cancer. Perhaps the carcinogens of the ambient air cannot accumulate in the intestine as they are already stopped by the lung or skin tissue beforehand. Associations between cancer incidence and green space (NDVI) or NO2 were present, but much smaller than for PM10.
We speculate that the high regional variation in NMSC (Fig. 1) could be due to differences in health care utilization in the population, whereas this is unlikely for more severe cancer types. One Australian study reported a higher skin cancer risk for regions with more vegetation , but exposure to UV light is certainly significantly greater in Australia compared to our study setting. Our results indicate that green space might be protective against NMSC in Germany. This observed protective effect might be secondary as the distribution of green spaces and PM10 were highly inversely correlated, and, as aforementioned, there are fewer emission sites around green spaces . Thus, in areas where the UV intensity is much weaker than in Australia like in Germany, residential green space could lower the skin cancer risk. However, more studies are needed to clarify whether residing in greener places is beneficial or detrimental for NMSC and if this effect is causal or only indirectly transmitted through a lack of AP at these sites.
In the age group 18 until 64 years 5.0% of the sample population had an alcohol-related disorder. Survey data from 2012 for Germany and this age group found 3.1% with alcohol abuse and additionally 3.4% with alcohol dependency according to DSM IV. Our 5.0% estimate lies between 3.1 and 6.5% . We expect therefore, that not all cases of alcohol abuse or dependency have already been recognized by physicians, but yet we consider it a valuable confounding variable. Further, a recent meta-analysis  and other studies [48, 59] found positive associations with high alcohol consumption on colorectal, mouth and throat, as well as breast cancer, and limited evidence for prostate cancer. Our effect estimate of alcohol on mouth and throat cancer is possibly a bit overestimated, as part of the people with an alcohol-related disorder are also smokers and these risk factors interact with each other .
We speculate that the high correlation of physician contacts and PM10 (~ 75%) suggests that PM10 is a wide-ranging risk factor for many diseases and therefore reducing it should be a central focus for community health.
Strength and limitations
One of the strengths of the presented analysis is that the dataset contains a huge population. The population covers almost half of the residential population of Saxony and is regarding age and sex distribution very similar to the general population of Saxony, which allows us to investigate cancer incidence for the federal state of Saxony. Further, the study offers complete information about out- and inpatient treatments and is not affected by sampling or nonresponse bias. Selection bias because of changes of the health insurance due to the cancer disease is also unlikely, since no differences in cancer treatment between statutory health insurances are present in Germany. We have invested a lot of work in exact case definitions and can therefore build our analyses on valid cancer diagnoses.
Nevertheless, several important limitations should be acknowledged. Address information up to four digit postal code districts was available, which did not allow a high spatial resolution of individual outcome and exposure assessments. Therefore, we could only use the exposure of the four digit postal code district as a surrogate for the individual exposure (semi-individual). Additional information on some potentially important individual confounders, such as smoking, socioeconomic status, diet, or physical activity was unfortunately not available, which could have led to residual confounding. Further confounding could occur due to genetic predisposition of patients for certain cancer types, or due to other environmental pollutants, like pesticides. But farm workers, with a potential high exposition with pesticides, are not included in our data set, since they have their own statutory health insurance in Germany. Unfortunately, we had no information about occupation, indoor radon exposure of patients and virus infections possible related to cancer which could also confound our analyses.
The observed associations could therefore be overestimated or be biased due to exposure misclassification, for example, but the magnitude should not be too high. In one study on lung cancer, the hazard ratio of PM10 decreased by only 11% if it was additionally controlled for smoking status, smoking intensity, square of smoking intensity, smoking duration, time since quitting smoking, environmental tobacco smoke, occupation, fruit intake, marital status, education level, and employment status . Our models assumed that the residential population did not move and AP was constant over time. Both assumptions, especially the first, are unlikely. Nevertheless, it is reassuring that our findings remained robust, when analyzing only the population that did not move outside their four digit postal code district during the observation period. Still, it is unknown what their movement history was until 2007. From the literature it is known that migrants have better health because people with low socioeconomic status or severe diseases are less mobile . This could bias our results and therefore, we trust our primary analysis more than the results of the sensitivity analysis with excluded migrants.
Beyond the current study situation we found some evidence that higher ambient AP levels increase the risk of mouth and throat cancer and of NMSC, while a higher residential green space level might have a protective effect. In summary, we assume that our effect estimates are not strongly biased by residual confounding, but we cannot exclude that for sure. Further research should try to measure the environmental exposures through the life course and focus on the effects of relocation together with environmental factors and extent the analyses to other cancer types.
Akaike information criterion
Diagnostic and Statistical Manual of Mental Disorders
Normalized Difference Vegetation Index
Non-melanoma skin cancer
- NO2 :
Uniform Value Scale
German modification of the International Classification of Procedures in Medicine
Pharmaceutical registration numbers
Cohen AJ, Brauer M, Burnett R, Anderson HR, Frostad J, Estep K, Balakrishnan K, Brunekreef B, Dandona L, Dandona R, et al. Estimates and 25-year trends of the global burden of disease attributable to ambient air pollution: an analysis of data from the global burden of diseases study 2015. Lancet (London, England). 2017;389(10082):1907–18.
Brook RD, Newby DE, Rajagopalan S. The global threat of outdoor ambient air pollution to cardiovascular health: time for intervention. JAMA cardiology. 2017;2(4):353–4.
Brook RD, Rajagopalan S, Pope CA 3rd, Brook JR, Bhatnagar A, Diez-Roux AV, Holguin F, Hong Y, Luepker RV, Mittleman MA, et al. Particulate matter air pollution and cardiovascular disease: an update to the scientific statement from the American Heart Association. Circulation. 2010;121(21):2331–78.
Cosselman KE, Navas-Acien A, Kaufman JD. Environmental factors in cardiovascular disease. Nat Rev Cardiol. 2015;12(11):627–42.
Hoek G, Krishnan RM, Beelen R, Peters A, Ostro B, Brunekreef B, Kaufman JD. Long-term air pollution exposure and cardio- respiratory mortality: a review. Environmental health : a global access science source. 2013;12(1):43.
Lee BJ, Kim B, Lee K. Air pollution exposure and cardiovascular disease. Toxicological research. 2014;30(2):71–5.
Li J, Sun S, Tang R, Qiu H, Huang Q, Mason TG, Tian L. Major air pollutants and risk of COPD exacerbations: a systematic review and meta-analysis. International journal of chronic obstructive pulmonary disease. 2016;11:3079–91.
Faustini A, Rapp R, Forastiere F. Nitrogen dioxide and mortality: review and meta-analysis of long-term studies. Eur Respir J. 2014;44(3):744–53.
Atkinson RW, Butland BK, Dimitroulopoulou C, Heal MR, Stedman JR, Carslaw N, Jarvis D, Heaviside C, Vardoulakis S, Walton H, et al. Long-term exposure to ambient ozone and mortality: a quantitative systematic review and meta-analysis of evidence from cohort studies. BMJ Open. 2016;6(2):e009493.
Pope CA 3rd, Dockery DW. Health effects of fine particulate air pollution: lines that connect. Journal of the Air & Waste Management Association (1995). 2006;56(6):709–42.
WHO. An overview of the evidence on environmental and occupational determinants of Cancer. In. 2011:1–7.
Perera FP. Environment and cancer: who are susceptible? Science. 1997;278(5340):1068–73.
Hamra GB, Laden F, Cohen AJ, Raaschou-Nielsen O, Brauer M, Loomis D. Lung Cancer and exposure to nitrogen dioxide and traffic: a systematic review and meta-analysis. Environ Health Perspect. 2015;123(11):1107–12.
Loomis D, Grosse Y, Lauby-Secretan B, El Ghissassi F, Bouvard V, Benbrahim-Tallaa L, Guha N, Baan R, Mattock H, Straif K, et al. The carcinogenicity of outdoor air pollution. The Lancet Oncology. 2013;14(13):1262–3.
Clapp RW, Jacobs MM, Loechler EL. Environmental and occupational causes of cancer: new evidence 2005-2007. Rev Environ Health. 2008;23(1):1–37.
Blackadar CB. Historical review of the causes of cancer. World journal of clinical oncology. 2016;7(1):54–86.
Pope CA 3rd, Burnett RT, Thun MJ, Calle EE, Krewski D, Ito K, Thurston GD. Lung cancer, cardiopulmonary mortality, and long-term exposure to fine particulate air pollution. Jama. 2002;287(9):1132–41.
Perez L, Rapp R, Kunzli N. The year of the lung: outdoor air pollution and lung health. Swiss Med Wkly. 2010;140:w13129.
Fajersztajn L, Veras M, Barrozo LV, Saldiva P. Air pollution: a potentially modifiable risk factor for lung cancer. Nat Rev Cancer. 2013;13(9):674–8.
Farmer PB, Singh R, Kaur B, Sram RJ, Binkova B, Kalina I, Popov TA, Garte S, Taioli E, Gabelova A, et al. Molecular epidemiology studies of carcinogenic environmental pollutants. Effects of polycyclic aromatic hydrocarbons (PAHs) in environmental pollution on exogenous and oxidative DNA damage Mutation research. 2003;544(2–3):397–402.
Perera FP, Hemminki K, Gryzbowska E, Motykiewicz G, Michalska J, Santella RM, Young TL, Dickey C, Brandt-Rauf P, De Vivo I, et al. Molecular and genetic damage in humans from environmental pollution in Poland. Nature. 1992;360(6401):256–8.
Perera FP, Jedrychowski W, Rauh V, Whyatt RM. Molecular epidemiologic research on the effects of environmental pollutants on the fetus. Environ Health Perspect. 1999;107(Suppl 3):451–60.
Council NR. Global sources of local pollution: an assessment of long-range transport of key air pollutants to and from the United States. Washington, DC: the National Academies Press; 2010.
Weichenthal S, Farrell W, Goldberg M, Joseph L, Hatzopoulou M. Characterizing the impact of traffic and the built environment on near-road ultrafine particle and black carbon concentrations. Environ Res. 2014;132:305–10.
Parent ME, Goldberg MS, Crouse DL, Ross NA, Chen H, Valois MF, Liautaud A. Traffic-related air pollution and prostate cancer risk: a case-control study in Montreal, Canada. Occup Environ Med. 2013;70(7):511–8.
Al-Ahmadi K, Al-Zahrani A. NO(2) and cancer incidence in Saudi Arabia. Int J Environ Res Public Health. 2013;10(11):5844–62.
Astell-Burt T, Feng X, Kolt GS. Neighbourhood green space and the odds of having skin cancer: multilevel evidence of survey data from 267072 Australians. J Epidemiol Community Health. 2014;68(4):370–4.
Gascon M, Triguero-Mas M, Martinez D, Dadvand P, Rojas-Rueda D, Plasencia A, Nieuwenhuijsen MJ. Residential green spaces and mortality: a systematic review. Environ Int. 2016;86:60–7.
Richardson EA, Mitchell R. Gender differences in relationships between urban green space and health in the United Kingdom. Soc Sci Med. 2010;71(3):568–75.
Mitchell R, Popham F. Effect of exposure to natural environment on health inequalities: an observational population study. Lancet. 372(9650):1655–60.
Markevych I, Schoierer J, Hartig T, Chudnovsky A, Hystad P, Dzhambov AM, de Vries S, Triguero-Mas M, Brauer M, Nieuwenhuijsen MJ, et al. Exploring pathways linking greenspace to health: theoretical and methodological guidance. Environ Res. 2017;158:301–17.
Kerr J, Anderson C, Lippman SM. Physical activity, sedentary behaviour, diet, and cancer: an update and emerging new evidence. The Lancet Oncology. 2017;18(8):e457–71.
Friedenreich CM, Shaw E, Neilson HK, Brenner DR. Epidemiology and biology of physical activity and cancer recurrence. Journal of molecular medicine (Berlin, Germany). 2017;
McTiernan A. Mechanisms linking physical activity with cancer. Nat Rev Cancer. 2008;8(3):205–11.
James P, Hart JE, Banay RF, Laden F. Exposure to greenness and mortality in a Nationwide prospective cohort study of women. Environ Health Perspect. 2016;
Toftager M, Ekholm O, Schipperijn J, Stigsdotter U, Bentsen P, Gronbaek M, Randrup TB, Kamper-Jorgensen F. Distance to green space and physical activity: a Danish national representative survey. J Phys Act Health. 2011;8(6):741–9.
van den Berg M, Wendel-Vos W, van Poppel M, Kemper H, van Mechelen W, Maas J. Health benefits of green spaces in the living environment: a systematic review of epidemiological studies. Urban For Urban Green. 2015;14(4):806–16.
Trautmann F, Schuler M, Schmitt J. Burden of soft-tissue and bone sarcoma in routine care: estimation of incidence, prevalence and survival for health services research. Cancer Epidemiol. 2015;39(3):440–6.
Swart E, Gothe H, Geyer S, Jaunzeme J, Maier B, Grobe TG, Ihle P. German Society for Social M, prevention, German Society for E: [good practice of secondary data analysis (GPS): guidelines and recommendations]. Gesundheitswesen. 2015;77(2):120–6.
World Medical Association Declaration of Helsinki. ethical principles for medical research involving human subjects. Jama. 2013;310(20):2191–4.
Swart E, Gothe H, Geyer S, Jaunzeme J, Maier B, Grobe TG, Ihle P. Good practice of secondary data analysis (GPS): guidelines and recommendations. Gesundheitswesen. 2015;77(2):120–6.
Small Area Health Statistics Unit [http://www.sahsu.org/content/data-download].
Vienneau D, de Hoogh K, Bechle MJ, Beelen R, van Donkelaar A, Martin RV, Millet DB, Hoek G, Marshall JD. Western European land use regression incorporating satellite- and ground-based measurements of NO2 and PM10. Environmental science & technology. 2013;47(23):13555–64.
Moderate Resolution Imaging Spectroradiometer Land Product Subsets [http://daacmodis.ornl.gov/cgi-bin/MODIS/GLBVIZ_1_Glb/modis_subset_order_global_col5.pl].
Tucker CJ. Red and photographic infrared linear combinations for monitoring vegetation. Remote Sens Environ. 1979;8(2):127–50.
Weier J, Herring D: Measuring vegetation (NDVI & EVI). 2000.
IARC_Working_Group. Personal habits and indoor combustions. Volume 100. E A review of human carcinogens IARC Monographs on the Evaluation of Carcinogenic Risks in Humans. 2012;100(Pt E):373–472.
IARC_Working_Group. Alcohol consumption and ethyl carbamate. IARC Monographs on the Evaluation of Carcinogenic Risks in Humans. 2010;96:3–1383.
R Development Core team. R: a language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing; 2010.
Eurostat: Revision of the European Standard Population - Report of Eurostat's task force. 2013.
Umlauf N, Adler D, Kneib T, Lang S, Zeileis A. Structured Additive Regression Models. An R Interface to BayesX: Faculty of Economics and Statistics, University of Innsbruck; 2012.
Brezger A, Kneib T, Lang S. BayesX: analyzing Bayesian structural additive regression models. J Stat Softw. 2005;014(i11)
Calcagno V, de Mazancourt C. Glmulti: an R package for easy automated model selection with (generalized) linear models. J Stat Softw. 2010;34(12):1–29.
Raaschou-Nielsen O, Andersen ZJ, Beelen R, Samoli E, Stafoggia M, Weinmayr G, Hoffmann B, Fischer P, Nieuwenhuijsen MJ, Brunekreef B, et al. Air pollution and lung cancer incidence in 17 European cohorts: prospective analyses from the European study of cohorts for air pollution effects (ESCAPE). The Lancet Oncology. 2013;14(9):813–22.
Raaschou-Nielsen O, Beelen R, Wang M, Hoek G, Andersen ZJ, Hoffmann B, Stafoggia M, Samoli E, Weinmayr G, Dimakopoulou K, et al. Particulate matter air pollution components and risk for lung cancer. Environ Int. 2016;87:66–73.
Vineis P, Husgafvel-Pursiainen K. Air pollution and cancer: biomarker studies in human populations. Carcinogenesis. 2005;26(11):1846–55.
Kraus L, Pabst A, Piontek D, EGd M. Substanzkonsum und substanzbezogene Störungen: Trends in Deutschland 1980 – 2012. Sucht. 2013;59(6):333–45.
Bagnardi V, Rota M, Botteri E, Tramacere I, Islami F, Fedirko V, Scotti L, Jenab M, Turati F, Pasquali E, et al. Alcohol consumption and site-specific cancer risk: a comprehensive dose-response meta-analysis. Br J Cancer. 2015;112(3):580–93.
Purdue MP, Hashibe M, Berthiller J, La Vecchia C, Dal Maso L, Herrero R, Franceschi S, Castellsague X, Wei Q, Sturgis EM, et al. Type of alcoholic beverage and risk of head and neck cancer--a pooled analysis within the INHANCE consortium. Am J Epidemiol. 2009;169(2):132–42.
Eronen J, von Bonsdorff M, Rantakokko M, Portegijs E, Viljanen A, Rantanen T. Socioeconomic status and life-space mobility in old age. J Aging Phys Act. 2016;24(4):617–23.
We thank the health insurance AOK PLUS for cooperation in data utilization and Stefan Baßler for proofreading the manuscript. We thank Beate Brecht for conducting data queries. We acknowledge support by the German Research Foundation and the Open Access Publication Funds of the TU Dresden.
The project was carried out without funding.
Availability of data and materials
The data that support the findings of this study are available from AOK PLUS, but restrictions apply to the availability of these data which were used under license for the current study and so are not publicly available. Data are however available from the authors upon reasonable request and with permission of AOK PLUS.
Ethics approval and consent to participate
The consultation of an ethics committee and consent to participate is not required for analyses based solely on secondary data.
Jochen Schmitt reports institutional funding for IITs from Sanofi, Novartis, ALK, and Pfizer. The other authors do not indicate any financial or non-financial competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Figure S3. STROBE diagram; Flowchart showing selection of incident cancer cases following the guidelines of the STROBE initiative (STrengthening the Reporting of OBservational studies in Epidemiology). (PNG 26412 kb)
Table S1. ICD-10 Codes; Analyzed cancer entities and corresponding ICD-10-GM codes are shown. (DOCX 12 kb)
Table S2. Validation of outpatient cancer cases other than NMSC; Used OPS-, EBM codes and prescribed medications in outpatient cancer care. (DOCX 14 kb)
Table S3. Validation of outpatient cancer cases (NMSC); Used OPS- and EBM codes in outpatient cancer care of NMSC. (DOCX 14 kb)
Figure S1. Correlation plot; Correlation between AP (PM10; NO2), NDVI, and physician contacts. Numbers are Pearson correlation coefficients. (PNG 1050 kb)
Table S5. Crude and adjusted linear regression of aggregated data; Results of the linear regression analysis are shown. Estimators, 95% confidence intervals (CI) and coefficients of determination (R2) are given for crude and adjusted models. (DOCX 19 kb)
Table S4. Non-mover multilevel Poisson regression (CIs are Wald confidence intervals); Relative risk estimates for an increase of environmental exposures (10 μg/m3 for PM10 / N02; 10% for NDVI) on cancer incidence in Saxony for patients with no relocation within study period between 2010 until 2014. (DOCX 15 kb)
Figure S2. Scatter plots of the crude linear regression analysis; Associations between PM10, NDVI, and NO2 and different cancer types. For each diagram, coefficients of determination (R2) are given. We see positive but weak associations between NO2 and mouth and throat cancer, prostate cancer and breast cancer, but an elevated positive association with NMSC exists (R2 = 0.164). With increasing vegetation level (NDVI) of the neighborhood, cancer incidence rates decrease for breast cancer, prostate cancer, mouth and throat cancer, and NMSC (in increasing order). (PNG 1295 kb)
Supplemental introduction text
All case definitions were based on the respective coding systems for diagnosis (ICD-10-GM), procedures (Uniform Value Scale (EBM) and German modification of the International Classification of Procedures in Medicine (OPS)), as well as prescriptions (Anatomical Therapeutic Chemical code - ATC and pharmaceutical registration numbers - PZN).