40 years of progress in female cancer death risk: a Bayesian spatio-temporal mapping analysis in Switzerland
BMC Cancer volume 15, Article number: 666 (2015)
In the past decades, mortality of female gender related cancers declined in Switzerland and other developed countries. Differences in the decrease and in spatial patterns within Switzerland have been reported according to urbanisation and language region, and remain controversial. We aimed to investigate geographical and temporal trends of breast, ovarian, cervical and uterine cancer mortality, assess whether differential trends exist and to provide updated results until 2011.
Breast, ovarian, cervical and uterine cancer mortality and population data for Switzerland in the period 1969–2011 was retrieved from the Swiss Federal Statistical office (FSO). Cases were grouped into <55 year olds, 55–74 year olds and 75+ year olds. The geographical unit of analysis was the municipality.
To explore age- specific spatio-temporal patterns we fitted Bayesian hierarchical spatio-temporal models on subgroup-specific death rates indirectly standardized by national references. We used linguistic region and degree of urbanisation as covariates.
Female cancer mortality continuously decreased in terms of rates in all age groups and cancer sites except for ovarian cancer in 75+ year olds, especially since 1990 onwards.
Contrary to other reports, we found no systematic difference between language regions. Urbanisation as a proxy for access to and quality of medical services, education and health consciousness seemed to have no influence on cancer mortality with the exception of uterine and ovarian cancer in specific age groups. We observed no obvious spatial pattern of mortality common for all cancer sites.
Rate reduction in cervical cancer was even stronger than for other cancer sites.
Female gender related cancer mortality is continuously decreasing in Switzerland since 1990. Geographical differences are small, present on a regional or canton-overspanning level, and different for each cancer site and age group. No general significant association with cantonal or language region borders could be observed.
Female gender related cancers, in particular cancer of the breast, corpus uteri, ovary and cervix uteri account for more than 40 % of newly diagnosed cancers and for about 30 % of cancer related deaths in Swiss women . In the past decades, female cancer mortality declined in Switzerland and the more developed countries  mainly due to advances in the understanding of tumour biology and in early detection, as well as the introduction of targeted therapies. However, differences in the decrease within Switzerland have been reported, such as for breast cancer in four selected cantons .
Switzerland is a small, affluent and culturally diverse confederation of 26 relatively autonomous states called cantons. Health care policies are developed at cantonal level resulting in a large geographical variation in health expenditures, control programs and care planning. I.e. population based mammography screening programs were and are implemented at very different time points over a period of more than 20 years in the various cantons. Most studies, including the above, investigated differences on the same regional level –cantons–, but it remained unknown whether these are consistent geographical disparities related to cantonal decisions or artefacts due to the choice of geographical and time units; driven by sub regions or complete region. The only more detailed maps of female cancer mortality rates are those of Schüler and Bopp  depicting geographical variation in mortality during 1970–1990 on the basis of so called MS-regions, 106 ‘unofficial’ regions smaller than cantons defined by mobility considerations. Since they have not applied temporal and geographical smoothing, the results may be distorted especially in areas where the population is small. This makes it difficult to distinguish chance variability from real differences. To our knowledge, covariate-adjusted and smooth, nationwide maps of female cancer mortality depicting the changes over time and space are not available.
Therefore, we studied geographical and temporal trends of breast, ovarian, cervical and uterine cancer mortality in Switzerland, adding 20 years of data to previous work, using state-of-the-art methodology for results with more detail and fewer artefacts, and without prejudice of geographical unit or shape of time trends. Hence, we used the most detailed available data (municipality level) and accounted for non-linear time trends. We hypothesized similar patterns for the different cancer sites and/or age group. Bayesian spatial models are the state-of-the-art modelling approach for assessing spatio-temporal patterns and trends. They “smooth” or improve estimation of an unstable rate by “borrowing” strength from its neighbours . They can also assess the significance of risk factors taking into account the geographical correlation, and are able to show spatial patterns after adjustment for geographical differences in certain risk factors.
Female cancer mortality data was obtained for the period 1969–2011 from death certificates coded centrally by the Swiss Federal Statistical office (FSO). The data include age at death, year of birth and death for each individual, nationality, municipality of residence, the cause of death and co-morbidities. Cause of death and co-morbidities are coded using the 8th revision of the International Classification of Diseases (ICD) until 1994/1995 and afterwards using the 10th revision. The transition to the 10th revision of the ICD-10 was accompanied by changes in death certificate coding practices (priority rules). We used age- and cancer site-specific correction factors as proposed by Lutz et al.  for the death counts. We included all cases coded with main causes of death being cancer of the female breast (ICD-10 C50.0-C50.9), cervix (ICD-10 C53.0- C53.9), corpus uterine (ICD-10 C54.0-C55.9) and ovary (ICD-10 C56.9). According to federal regulations, mortality data excluding any identifiable information can be used in epidemiological studies without additional ethics committee approval.
Detailed population data on municipality level is only available from census that takes place in Switzerland every 10 years with the last one taking place in 2010. We aggregated the mortality data in five 4-year periods around the census years, i.e. 1969–1972, 1979–1982, 1989–1992, 1999–2002 and 2008–2011, in which population was assumed to be constant.
There are around 2,500 municipalities in the country. Over the study period, the number of municipalities has changed due to fusion, separation, deletion or new occurrences. We aligned all data on the 2011 municipality structure using spatial data for 2011 and municipality transition protocols for each year obtained from the FSO. From the same source, we retrieved data on language region (German, French and Italian/Romansh) and urbanisation (Fig. 1). We grouped municipalities classified as central agglomeration city, greater agglomeration and isolated city into “urban” leaving the classification “rural” unchanged.
Age was grouped into three groups (<55, 55–74, 75+ year olds). The geographical unit of analysis was the municipality.
In a preliminary analysis, we investigated SMR ratio values in a non-spatial model. Spatio-temporal Poisson and negative binomial regression models were fitted separately for each age group on the number of deaths aggregated by municipality and year with the mean being equal to the product of the expected death count and age standardised mortality rate. Indirect standardisation used 5 years age intervals. Expected mortality counts for each municipality, year and age group were obtained from the study population using nationwide age-specific mortality rates for all periods.
Space and temporal random effects as well as possible non-linear temporal trends were modelled on the log of the mean standardised mortality rate following model formulations of Jürgens et al.  (cf. Appendix 1). In particular, municipality-specific random effects were modelled via conditional autoregressive (CAR) models to filter out the noise and highlight the observed patterns. The models were formulated as hierarchical Bayesian models with parameter estimation via Markov chain Monte Carlo simulation (MCMC). We used the Deviance Information Criterion (DIC) to select the regression models from Poisson/Negative binomial regression with or without an additional set of unstructured random effects for each municipality.
Data on language and urbanisation were included as covariates in the model. These analyses will indicate whether there are statistically significant differences in the cancer mortality for each one of the above covariates, assessed by 95 % Bayesian Credible Intervals (CI).
From the estimates of the model, we produced smoothed maps displaying geographical patterns of female gender cancer mortality for each age group, cancer site and year since 1969 till recent almost to date.
Table 1 shows the number of female cancer deaths and crude rates per 100,000 person years in Switzerland by age group within the 4-year periods under investigation. Among the cancer sites studied, breast cancer was the most common cause of death, followed by ovarian, uterine and cervical cancer.
Mortality rates continuously decreased for cervical and uterine cancer, and for ovarian cancer in the <55 year olds. For breast cancer and the other age groups of ovarian cancer, mortality rates decreased only as from 1979–1982 and from 1989–1992 for 75+ year olds respectively.
Table 2 shows the results of the spatio-temporal regression analysis by cancer site and age group. With the spatial analysis, we could confirm the time trends observed in the crude rates in Table 1, while only in few cases the covariates had a significant effect on the standardized mortality ratio (SMR). Language region had in none of the models a significant effect on mortality, urbanisation only in 3 models: An urban environment was associated with a significantly lower mortality of 55–74 year olds in uterine cancer and <55 year olds in ovarian cancer, and associated with higher ovarian cancer mortality in 75+ year olds.
The spatial patterns of mortality based on smoothed small area estimates (Figs. 2, 3, 4 and 5, Additional file 1) are different for the female cancers and age groups and not homogenous among the country. No general, significant coincidence with cantonal or language region borders could be observed, with the latter additionally being confirmed by spatial regression for all cancer sites and age groups (Table 2). The spatial patterns form either sub-cantonal areas or canton-overspanning areas.
For all cancer sites and age group combinations the model 1 with Poisson distributed data and only one, spatially structured, random effect was identified as the best model, with lowest DIC (see Table 3). SMR ratios in the non-spatial models were close to the results presented in Table 2, and significance was the same for all but 4 out of 84 coefficients, with their CIs being close to zero in both models.
Using modern Bayesian small area modelling and mapping techniques we have been able to show that all investigated groups of women in Switzerland have benefited from progress in cancer control regardless of place of residence in the past 40 years. We observed only small differences in the geographical variation of mortality.
A factor, which may have contributed to breast and uterine cancer mortality reductions, is the change in the use of hormone replacement therapy (HRT) . After an association of HRT use with breast cancer occurrence was reported , its use declined sharply.
We were also not able to show similar spatial patterns in breast and ovarian cancer mortality although they share several life style related, environmental and genetic risk factors. It should be noted however, that hereditary cancer accounts only for about 5-10 % of the cases in breast cancer  and about 15 % in ovarian cancer . They are shown to occur at younger age and more advanced stage; still, a visible effect on the mortality map may only be seen in areas with ethnic groups or very large families with a highly elevated risk for hereditary cancer. Such a risk has been described for Ashkenazi Jewish women. The BRCA Ashkenazi founder gene mutations are prevalent in approximately 2 % of these women  with communities of Ashkenazi mainly found in urban areas; largest communities are in the cities of Zürich, Geneva and Basel contributing to 1-2 % of the population [12, 13]. However, the breast and ovarian cancer risk in BRCA carriers is affected by genetic modifiers and non-genetic factors, for example, reproductive behaviour, hormonal exposure, lifestyle and risk reduction surgeries . We could not observe an elevated mortality for the three cities in contrast to the surrounding area and it remains unclear to which extent the mortality rates are driven by these hereditary forms of cancer.
Considerable differences in health and health related behaviour have been reported for the Swiss language regions including alcohol intake, smoking and a healthy diet [15, 16] but lacked significance as regression factors in our analysis.
Only for three cancer site-age group combinations was the urbanisation level identified as a significant factor. Urbanisation is serving as a proxy for access to and quality of medical services, education and health consciousness . By our regression with 20 years of new data, we could not formally confirm an urban–rural gradient for breast cancer as described by Schüler & Bopp  as significant.
Overall, no general pattern across age groups or cancer sites was present.
The reduction of mortality was stronger in the younger age groups, which is probably the result of better survival and therefore a shift in the age of death. This would also explain the temporary increase in breast and ovarian cancer death risk around the year 1990 in the 75+ year olds. In addition, in this age group multi-morbid conditions and fewer treatments are common . Sant et al.  noted that poor survival for gynaecological cancers in the elderly could be due to advanced stage at diagnosis, or failure to give adequate treatment, perhaps because of comorbidity. In general, the interpretability of results in this age group is limited due to its small size, more multi-morbid conditions together with possible inconsistencies in death certification over time, because of only allowing one single cause of death.
Strengths and limitations
As cancer deaths are rare events and in order to increase the power, different geographical units have been used when analysing cancer mortality data in the past. Some authors have used selected cantons  and Schüler & Bopp  used for their cancer atlas somewhat smaller mobility regions based on the accessibility to goods and services but which do not take into account population size. As a result, this choice was too aggregated for some urban areas and not aggregated enough for some sparsely populated areas in order to reveal robust, underlying trends. In view that the choice of the geographical unit of analysis may greatly influence results , the combination of small geographical units with a state-of-the art smoothing technique enabled a more detailed analysis. With this analysis, we could additionally show the driving age groups or subareas of elevated or reduced mortality in certain regions, while reducing uncertainties due to small numbers and adding an investigation of non-linear time trends.
In general, smoothing allows an estimation of the underlying risk, in a sort of a long-year average, rather than the actual situation. However, for single municipalities, without fully eliminating it, the use of Bayesian smoothing reduces the probability to detect narrow areas with specifically high or low risk. Municipalities at the country border may not benefit from smoothing to the same extent as municipalities in the interior of the country due to unknown data on the other side of the border. Therefore, in the interpretation of the results emphasis should be given to the broader spatial patterns rather than to single municipalities.
Comparing with the previous work of Schüler & Bopp  our study not only extended their work by 20 more years and corrected for non-linear time effects, more importantly, we were able to correct the foreseen overestimation in mortality numbers until 1994, which could not be adequately addressed earlier. Priority rules in the coding of causes of death led to an overestimation in cancer deaths due to their prioritization over other comorbidities. The applied methodology of age standardisation takes advantage of the actual age structure rather than a standard population.
There are important limitations to our study. Risk factors affect incidence but are not necessarily linked to mortality . The progression stage of the tumours and their histological type could not be taken into account, as the ICD-classification does not include histological type for the sites studied. The regional case mix and its changes over time therefore may have distorted the results.
Further distortions may arise from the uncertainty as to what level the reported main cause of death and comorbidities are comparable in time and between regions, although the central coding speaks in favour of a certain homogeneity in the coding procedure. In the elderly with frequent multi-morbid conditions, the probability of misclassification is higher.
Furthermore, after prior analysis the covariates language region and urbanisation level were fixed in time for the municipalities, so that varying developments therein may have resulted in inaccuracies.
Female gender related cancer mortality continuously decreased in Switzerland. In most age groups, this decline was significant and quite strong in the past decades, resulting in values more than 6 times lower within 40 years. The strongest reduction of mortality was observed for cervical cancer, followed by uterine, ovarian and breast cancer.
Geographical differences are small and do not follow cantonal borders. Spatial patterns were different for each cancer site and age group. The reasons for these differences are manifold, rising awareness, major advances in cancer therapy and ongoing developments in the field had a major impact on the cancer mortality.
Information on the geographical patterns and temporal trends of the disease burden at different regional scales are important for the design, implementation and evaluation of programs for cancer control. Access to specialized medical facilities should be increased especially in high priority areas in order to further reduce disparities. However, existing disparities are small.
GLOBOCAN 2012 v1.0, Cancer Incidence and Mortality Worldwide [http://globocan.iarc.fr]
Bulliard JL, La Vecchia C, Levi F. Diverging trends in breast cancer mortality within Switzerland. Ann Oncol. 2006;17:57–9.
Schüler G, Bopp M. Atlas der Krebsmortalität in der Schweiz 1970–1990. Basel: Birkhäuser Verlag; 1997.
Bernardinelli L, Montomoli C. Empirical Bayes versus fully Bayesian analysis of geographical variation in disease risk. Stat Med. 1992;11:983–1007.
Lutz JM, Pury P, Fioretta G, Raymond L. The impact of coding process on observed cancer mortality trends in Switzerland. Eur J Cancer Prev. 2004;13:77–81.
Jurgens V, Ess S, Phuleria HC, Fruh M, Schwenkglenks M, Frick H, et al. Bayesian spatio-temporal modelling of tobacco-related cancer mortality in Switzerland. Geospat Health. 2013;7:219–36.
Bouchardy C, Usel M, Verkooijen HM, Fioretta G, Benhamou S, Neyroud-Caspar I, et al. Changing pattern of age-specific breast cancer incidence in the Swiss canton of Geneva. Breast Cancer Res Treat. 2010;120:519–23.
Beral V, Million Women Study C. Breast cancer and hormone-replacement therapy in the Million Women Study. Lancet. 2003;362:419–27.
Campeau P, Foulkes W, Tischkowitz M. Hereditary breast cancer: new genetic developments, new therapeutic avenues. Hum Genet. 2008;124:31–42.
Pal T, Permuth-Wey J, Betts JA, Krischer JP, Fiorica J, Arango H, et al. BRCA1 and BRCA2 mutations account for a large proportion of ovarian carcinoma cases. Cancer. 2005;104:2807–16.
Struewing JP, Hartge P, Wacholder S, Baker SM, Berlin M, McAdams M, et al. The Risk of Cancer Associated with Specific Mutations of BRCA1 and BRCA2 among Ashkenazi Jews. N Engl J Med. 1997;336:1401–8.
STAT-TAB: Die interaktive Statistikdatenbank [http://www.pxweb.bfs.admin.ch/]. Access date 20.06.2014
Juden in der Schweiz [http://www.swissjews.ch/de/kultur/juden_in_der_schweiz/index.php]. Access date 20.06.2014
Levy-Lahad E, Friedman E. Cancer risks among BRCA1 and BRCA2 mutation carriers. Br J Cancer. 2007;96:11–5.
Calmonte R, Galati-Petrecca M, Lieberherr R, Neuhaus M, Kahlmeier S: Gesundheit und Gesundheitsverhalten in der Schweiz 1992–2002. Neuchâtel: Schweizerische Gesundheitsbefragung. Bundesamt für Statistik; Neuchâtel 2005.
Lieberherr R, Marquis J-F, Storni M, Wiedenmayer G. Gesundheit und Gesundheitsverhalten in der Schweiz 2007 - Schweizerische Gesundheitsbefragung. Neuchâtel: Bundesamt für Statistik (BFS); 2010.
Joerger M, Thurlimann B, Savidan A, Frick H, Rageth C, Lutolf U, et al. Treatment of breast cancer in the elderly: a prospective, population-based Swiss study. J Geriatr Oncol. 2013;4:39–47.
Sant M, Aareleid T, Berrino F, Bielska Lasota M, Carli PM, Faivre J, et al. EUROCARE-3: survival of cancer patients diagnosed 1990-94--results and commentary. Ann Oncol. 2003;14 Suppl 5:v61–118.
Woods LM, Rachet B, Coleman MP. Choice of geographic unit influences socioeconomic inequalities in breast cancer survival. Br J Cancer. 2005;92:1279–82.
Barnett GC, Shah M, Redman K, Easton DF, Ponder BAJ, Pharoah PDP. Risk Factors for the Incidence of Breast Cancer: Do They Affect Survival From the Disease? J Clin Oncol. 2008;26:3310–6.
This research was co-funded by the Cancer League Eastern Switzerland and an SNF grant, project no. 32003B_135769.
The authors declare that they have no competing interests.
PV, SE conceived of the study. CH carried out the analysis and data acquisition. CH, SE, PV contributed to the analysis of the data and the writing of the manuscript. SE, BT, NP and PV contributed to interpretation of the findings and critically revised the manuscript. All authors read and approved the final manuscript.
Appendix model formulations
Observed age and cancer site-specific counts of deaths Y it in municipality i(i = 1, …, N) in period t to follow a poisson distribution Y it ~ Pois(μ it ). Age and cancer specific random effects as well as possible non-linear trends were modelled on the log of the mean Age Standardized Mortality Ratio (SMR).
Detailed Figures of SMR development by cancer sites and age groups. Development of age standardized breast (Figures S2a-c), cervical (Figures S3a-c), uterine (Figures S4a-c) and ovarian (Figures S5a-c) cancer mortality (SMR) and spatial differences therein among all time periods by age group. (PDF 5957 kb)
Color key for figures 2-5. (PDF 164 kb)
About this article
Cite this article
Herrmann, C., Ess, S., Thürlimann, B. et al. 40 years of progress in female cancer death risk: a Bayesian spatio-temporal mapping analysis in Switzerland. BMC Cancer 15, 666 (2015). https://doi.org/10.1186/s12885-015-1660-8
- Breast cancer
- Ovarian cancer
- Cervical cancer
- Uterine cancer
- Bayesian inference
- Disease mapping
- Time trends