Introduction In the past decades, mortality due to breast cancer has declined considerably in Switzerland and other developed countries. The reasons for this decline remain controversial as several factors occurred almost simultaneously, including important advances in treatment approaches, breast cancer awareness and the introduction of mammography screening programmes in many European countries. In Switzerland, mammography screening programmes (MSPs) have existed in some regions for over 20 years but do not yet exist in others. This offers the possibility to analyse its effects with modern spatiotemporal methodology. We aimed to assess the spatiotemporal patterns and the effect of MSPs on breast cancer mortality.
Participants The study covers breast cancer deaths of the female population of Switzerland during the period 1969–2012. We retrieved data from the Swiss Federal Statistical Office aggregated on a small-area level.
Design We fitted Bayesian hierarchical spatiotemporal models on death rates indirectly standardised by national references. We used linguistic region, degree of urbanisation, duration of population-based screening programmes and socioeconomic index as covariates.
Results In Switzerland, breast cancer mortality in women slightly increased until 1989–1992 and declined strongly thereafter. Until 2009–2012, the standardised mortality ratio declined to 57% (95% CI 54% to 60%) of the 1969–1972 value. None of the other coefficients of the spatial regressions had a significant effect on breast cancer mortality. In 2009–2012, no region had significantly elevated or reduced breast cancer mortality at 95% credible interval level compared with the national mean.
Conclusion There has been a strong reduction of breast cancer mortality from the 1990s onwards. No important spatial disparities were observed. The factors studied (urbanisation, language, duration of population-based MSP and socioeconomic characteristics) did not seem to have an influence on them. Low participation rates and opportunistic screening use may have contributed to the low impact of MSPs.
- breast cancer
- bayesian disease mapping
This is an Open Access article distributed in accordance with the Creative Commons Attribution Non Commercial (CC BY-NC 4.0) license, which permits others to distribute, remix, adapt, build upon this work non-commercially, and license their derivative works on different terms, provided the original work is properly cited and the use is non-commercial. See: http://creativecommons.org/licenses/by-nc/4.0/
Statistics from Altmetric.com
Strengths and limitations of this study
A modern Bayesian spatial model was used to improve estimation of an unstable rate by ‘borrowing’ strength from its neighbours.
The model is capable of assessing the significance of risk factors while also taking the geographical correlation into account.
Switzerland with its homogeneous health system and different regional screening policies provides an ideal setting for assessing the impact of population-based mammography screening programmes.
Data on the geographical differences in opportunistic screening use and therefore overall screening participation are not available, where opportunistic screening use is estimated to be high and programme participation less than 50%.
The ecological study design does not allow an assessment of the combined impact of participation in and type (programme vs opportunistic) of mammography screening.
In Switzerland, breast cancer is the most frequently diagnosed cancer in women,1 it is the leading cause of cancer-related deaths2 and of premature mortality for Swiss women.3 Mortality due to breast cancer has declined considerably in the past decades in Switzerland and other developed countries.4 The reasons for the decline remain controversial because several factors including important advances in treatment approaches, breast cancer awareness and the introduction of mammography screening programmes in many European countries occurred almost simultaneously.
Some randomised controlled studies5 have demonstrated a breast cancer mortality reduction of 20% for women invited for breast cancer screening. However, they were conducted in the 1970s to 1980s. Since then, many advances in therapies have been made and adopted6 so that some authors doubt that the difference would persist under present conditions. Therefore, often used historical prescreening control groups are not best suited to disentangle these effects. Autier et al 7 compared countries in Europe but a criticism was that different countries may have different health systems. Kalager et al 8 used comparison groups in Norway and showed that only a third of the total mortality reduction could be attributed to mammography screening. However, a short observation period was used. Olsen et al 9 confirmed these results in principle with the same data but with a somewhat longer follow-up duration. In addition, in a setting where voluntary screening is assumed to be high, it is unknown what the effect an organised screening programme would be for the population as a whole.
In Switzerland, with its homogenous health system, these pitfalls can be avoided. Switzerland is a small confederation of 26 relatively autonomous states called cantons with somewhat low inequalities10 and many health and cancer-related resources.11–13 Although the healthcare system is homogeneous in providing universal and rapid access to and use of almost unlimited healthcare resources, some healthcare policies are developed at the cantonal level; in particular, the decision to initiate a population-based mammography screening programme. These programmes were implemented in Switzerland at different times over the past two decades. The first Swiss mammography pilot programme was established in 1993 in the French-speaking canton of Vaud. However, it was only in 2010 that the first organised programme in a German-speaking canton (St Gallen) started.
In breast cancer incidence, cantonal differences are well known and have been attributed to the differential use of opportunistic or organised mammography screening.14 In addition, considerable differences in health and health-related behaviour that affect the risk of breast cancer, including alcohol intake and a healthy diet, have been reported for the Swiss language regions,15 16 as well as differences in the age at first childbirth and number of children born to a mother.17 Differences in access to mammography screening and in lifestyle may be reflected in spatiotemporal differences in both breast cancer incidence and mortality, whereas only the latter will reflect the management of breast cancer.
In contrast, breast cancer mortality studies in Switzerland showed contradictory results. Bulliard et al 18 observed a steeper decrease from 1980 to 2002 in 55–74-year olds in French-speaking regions where population-based mammography screening started earlier. In a recent study,19 we presented the spatiotemporal trends of female gender-related cancer mortality in Switzerland by age group. The geographical differences found were small. We observed a differential decline in breast cancer mortality by age. The decline was highest in women younger than 50 and lower in women 75 or older. A similar pattern was observed in other European countries4 and attributed to early detection by mammography and to improved treatment.20–22 However, it was not clear to what extent improvements in survival could have affected the age at death. It was difficult to evaluate a shift of deaths into the next higher age group and the influence of screening programmes based on using fixed age groups rather than cohorts.
In the present study, we aimed to assess the spatiotemporal patterns in breast cancer mortality and specifically the effect of population-based mammography screening programmes on it. We corrected for urbanisation for which a mortality gradient was described23 and additionally for area-based socioeconomic factors, which may have influenced results in the previous study.
The Swiss Federal Statistical Office provided data on female breast cancer mortality, electronically available for the period of 1969–2012. The anonymised data included sex, age, year of birth and death for each individual, nationality, municipality of residence, the cause of death and comorbidities. The cause of death and comorbidities were coded centrally from death certificates using the eighth revision of the International Classification of Diseases (ICD) for deaths until 1994, and the 10th revision for deaths that occurred afterwards. The transition to the 10th revision of the ICD 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 24 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). According to federal regulations, mortality data excluding a person’s identifying information can be used in epidemiological studies without additional ethics committee approval.
The administrative borders of Swiss municipalities define the smallest geographical unit for which data were available. There are around 2500 municipalities in the country with a median population of 740 inhabitants in 1970 and 1150 in 2010.
Aggregated population data by age and area unit were extracted from the census that takes place in Switzerland every 10 years. The last census was conducted in 2010. Because of missing detailed intercensal population data, we aggregated the mortality data in five 4-year periods around the census years, that is, 1969–1972, 1979–1982, 1989–1992, 1999–2002 and 2009–2012, in which population was assumed to be constant and identical to the census year.
From the same source, we retrieved data on language region (German, French and Italian and Romansh) and urbanisation (rural/urban). We obtained information on population-based screening programmes from the Swiss Federation of Cancer Screening Programmes25 and categorised their duration in the census years into ‘no programme’, ‘0–4 years’ and ‘5+ years’. Data on socioeconomic position by municipality were provided by the Swiss National Cohort26 based on the census data of 2000.
Table 1 shows the observed number of deaths and mortality rates for each of the covariates.
As a small area geographical unit, we used the municipality borders as of 2012. We used municipality transition protocols from the Federal Statistical Office to align all data to this structure.
We investigated mortality for all ages combined in a spatial and a non-spatial model, one time for the five time periods from 1969 to 2012 to assess possible non-linear time trends and another time only for the period of 2009–2012.
For the spatial model, we used the Bayesian hierarchical spatiotemporal Poisson model formulations as described in Herrmann et al,19 fitted on the number of deaths aggregated by small area and year, with the mean being equal to the product of the expected death count and age-standardised mortality rate. The indirect standardisation used 5-year age intervals. Expected mortality counts for each small area and year were obtained from the study population using nationwide age-specific mortality rates, once for all periods and again only for the period of 2009–2012. The small area-specific random effects were modelled via conditional autoregressive models to filter out the noise and highlight the observed patterns. The deviance information criterion (DIC) was used to select the regression model from Poisson, zero-inflated Poisson and negative binomial regression models. The DIC was lowest with the Poisson regression model.
We accounted for differences that were influenced by linguistic region, life in rural or urban areas, screening programme duration and socioeconomic position. These analyses are used to indicate whether there are significant differences in cancer mortality for each of the above covariates, assessed by 95% Bayesian Credible Intervals.
No patients were involved in this study.
In Switzerland, more than 61 000 women died from breast cancer between 1969 and 2012. Table 2 presents the results of the regressions including all time periods and time trends. In Switzerland, breast cancer mortality in women slightly increased until the 1989–1992 period and has declined strongly since. Until the most recent period (2009–2012), the SMR has fallen to 57% of the 1969–1972 period’s value, both in the non-spatial and the spatial models. The trends and geographical differences are visualised in figure 1.
From the covariates studied, only the year of death and the urbanisation level in the non-spatial model had a significant impact when investigating all periods. An urban environment was associated with a 5% elevated SMR (3% in the spatial model) compared with a rural environment.
Limiting the analysis to the period of 2009–2012, none of the regression factors had a significant effect on breast cancer mortality (table 3).
Most SMR ratios of the non-spatial and the spatial model showed nearly identical values. The length of a screening programme and the French language region showed slightly higher values, but the differences were not significant.
In the 2009–2012 period, no region had a significantly higher or lower breast cancer mortality rate at 95% CI level compared with the national mean (figure 2). A map with covariate-adjusted smoothed SMR values is not shown because there was no information gain. The covariates are not significant and the geographical patterns are the same as for the smoothed SMR values.
The socioeconomic index value for the municipalities ranged from 28 to 85, where 25% of municipalities were below 55% and 25% above 66.
In the past decades, breast cancer mortality has nearly halved in Switzerland when considering all ages together. This trend, including the shift from increasing to decreasing rates around the period of 1989–1992, has been observed in several other European countries.4 Although significant spatial differences in breast cancer incidence are well described for Switzerland, we have not found any significant differences in breast cancer mortality in any of the periods studied. We have not observed any general significant differences between regions classified by duration of screening programmes, urbanisation, language and socioeconomic position. In addition, when limiting the analysis to the most recent period (2009–2012), none of the factors are significant. In fact, at 95% CI, none of the regions have a significantly elevated or reduced breast cancer mortality compared with the national mean.
There are several factors that explain why the significant differences in incidence do not translate into corresponding mortality differences. Most importantly, risk factors such as health and health-related behaviour that are reported to be different for the language regions16 affect incidence but are not necessarily linked to mortality.27 That is, while a temporary increase in the use of hormone replacement therapy has led to an increase in breast cancer incidence, many of those tumours have a favourable prognosis and might have influenced breast cancer mortality only marginally.28 Accordingly, the French language region, despite earlier implementation of mammography screening programmes, does not show a relevant impact on breast cancer mortality in our study.
Because screening has been identified as a potential source of mortality reduction,21 we also included data on population-based screening programme duration. However, our study did not show a significant effect on mortality on the population level. The reasons for this are probably manifold and may include factors such as screen-detected cancers being mainly of low stage, many women having not participated in the screening programmes or having chosen to undergo opportunistic screening. In addition, the effect of advances in diagnosis and therapy on mortality is quite strong and may have outweighed benefits from population-based screening programmes, as suggested by Autier et al.29 Moreover, the level of opportunistic screening in Switzerland has been described to be quite high,30 but data on the geographical differences in opportunistic screening use, and therefore overall screening participation, are not available. Data on participation in population-based screening programmes are published in a national monitoring report showing that participation rates of the programmes are close to the combined mean of 47.8%.31The ecological study design does not allow the assessment of the combined impact of participation in and type (programme vs opportunistic) of mammography screening or the impact of stage of tumour at diagnosis and mortality at an individual level. For the above reasons, the interpretability with regard to screening is limited. In addition, we had to group into 0–4 years and 5+ years of screening, which was done to avoid overfitting issues. There are only a few regions that are in close proximity to each other with 10+ years of screening in the 2009–2012 period only (see online supplementary figure A1).
Supplementary file 1
The present study is an in-depth analysis of our previous study,19 focusing on breast cancer mortality using an additional year of more recent data. We were also interested in the effects on the population as a whole. The applied methodology of age standardisation suits this by taking advantage of the actual age structure rather than of a standard population.
The non-significant fixed effect of socioeconomic position is in line with the results of Panczak et al.32 The additional correction served the disentanglement of affluence from the urbanisation parameter—which is connected with access to medical services—and further possible distortions.33
A strength of Bayesian spatial models is their ‘smoothing’ or improvement of estimation of an unstable rate by ‘borrowing’ strength from its neighbours.34 These models can also assess the significance of risk factors, taking into account the geographical correlation and are able to show spatial patterns after adjusting for geographical differences in certain risk factors. By adding a time dimension, Bayesian spatiotemporal models indicate changes of geographical patterns over time and determine how a disease evolves in different regions and different groups of the population (age, language or affluence groups). These models have provided a state-of-the-art modelling approach over the last 15 years for assessing spatio-temporal patterns and trends. We have not observed that coefficients in our analysis have shrunk towards zero when including geographical correlation as hypothesised by Hodges and Reich.35 In fact, in the spatial model for the 2009–2012 period, the impact of the French language region is 1.03 compared with 1.00 in the non-spatial model. However, we have included the results of the non-spatial models as well.
There has been a strong reduction of breast cancer mortality from the 1990s onwards. Geographical differences are present, but at a moderate level with no significant differences in the overall mean. In addition, they are not explained by the duration of population-based screening programmes, socioeconomic position, urbanisation and language region. Low participation rates and opportunistic screening use may have contributed to the low impact of mammography screening programmes. Continuous evaluation of geographical patterns of breast cancer mortality using modern spatiotemporal methodology is necessary for evaluating the efficacy of programmes.
Contributors 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. CH, PV, BT, NPH, CR and SE contributed to interpretation of the findings and critically revised the manuscript. All authors read and approved the final manuscript.
Funding CH was supported by the Cancer League Eastern Switzerland and CH and PV were supported by a grant of the Swiss National Science Foundation with the number 32003B_135769.
Disclaimer The funders had no role in the study design, data collection and analysis, decision to publish or preparation of the manuscript.
Competing interests None declared.
Patient consent Detail has been removed from this case description/these case descriptions to ensure anonymity. The editors and reviewers have seen the detailed information available and are satisfied that the information backs up the case the authors are making.
Provenance and peer review Not commissioned; externally peer reviewed.
Data sharing statement All data are publically available from the sources stated in the Methods section. Statistical code is available from the corresponding author.
If you wish to reuse any or all of this article please use the link below which will take you to the Copyright Clearance Center’s RightsLink service. You will be able to get a quick price and instant permission to reuse the content in many different ways.