Abstract
Aim: To characterize the spatiotemporal patterns of breast cancer (BC) incidence in females in the area with the highest incidence rate (IR) of the country in 2005-2012. Materials and Methods: The BC-IR was studied using mapping techniques, analysis of spatiotemporal clusters and analysis of spatial variations in temporal trends. Results: The overall BC-IR was 119.13/105 inhabitants. The annual BC-IRs were 17.7, 156.9, 213.3 and 232.9/105 inhabitants for women diagnosed at <40, 40-49, 50-64 and ≥65 years of age. This IR increased overall (by 4.113%/year) and for the four age groups (by 5.935, 3.833, 4.114 and 2.194%/year, respectively). In patients with locoregional and metastatic disease, the IRs were 93.6 and 7.4/105 inhabitants, increasing by 6.976 and 0.303%/year, respectively. Several spatiotemporal clusters and two spatial-variations in temporal trends were detected. The Lisbon region showed high IR clusters for most groups. Conclusion: This study identified critical areas of high IR and increasing trends for female BC-IR, providing evidence of heterogeneities in this area.
Breast cancer (BC) is the second most common cancer worldwide and the most frequent malignant disease in women (1-3). This age-associated malignancy (2, 4-9) represented 25% of female cancer worldwide in 2012 (2, 3) and in Europe corresponded to 27.4% (10) and 28.9% (11) of cancer in women in 2004 and 2006, respectively. In 2012, the age-adjusted (world population) BC incidence rate (IR) in the world, Europe (EU) and Southern EU were 43.3, 71.1 and 74.5/105 inhabitants, respectively (3), being the most common cancer in women in all European countries (12). The age-adjusted European BC-IR (for the EU population) in the same year was 94/105 inhabitants (4, 12) [EU-27 median incidence: 109/105 inhabitants (12)].
In developed countries, the BC-IR has been rising (1, 2, 4, 6, 13-15), a growth mostly detected until the year 2000 (13, 16) when the use of hormone replacement therapy in older women started to become obsolete (2, 16). This change in clinical practice especially affected the IR in older patients (13, 16) [ages 50-64 years (13, 17, 18) and 65-74 years (13, 17, 18)]. Still, the European BC-IR has been shown to increase with age until age 60-65 years, decreasing then in older women (8). Although some authors have described a stable IR in patients age <40 (16) or <50 (13) years, an increase of BC-IR has been reported in younger and middle-age (15) patients [age 35-49 (17)]. A European study on IR detected a growth (1-3%/year) of BC among women age <40 years in 1995-2006 (19), which might be explained by screening programs and exposure to predisposing factors (19). Additionally, a Swiss registry analysis showed a 1.6-2.7% growth of BC-IR in women <40 years (15), and in France a small increment in this IR also was described for women of this age (18). Conversely, in the UK, the BC-IR increased in 1993-2014 for women ≥65 years but stabilized in patients age 25-49 and 50-64 (6). Notably, when looking at women aged 50-69 years, the BC-IR increased for women in the UK in 1981-2009 (20) and for Italian women in 1985-1994 (21). This growth could potentially reflect screening actions (2, 15, 20) as the BC-IR plateaued in subsequent years (6, 20). Finally, an Italian study indicated an increase in local BC cases between 1985 and 1994 (21), but a decrease in regional and metastatic cases was detected (21).
An increase in BC cases also was described in Portugal (22-24), with 47,868 and 4,300 new cases being diagnosed in 1995-2009 (25) and in 2002 (26), respectively. In 2008, 2010 and 2012, about a third of new Portuguese female neoplasm cases were of BC (30.2, 31.1 and 29.4%, respectively) (27-29). In 2007 and 2009, the Portuguese female BC-IR was 101.80/105 inhabitants (standardized IR: 82.40/105 inhabitants) (30) and 110.12/105 inhabitants (standardized IR for the EU population: 87.58/105 inhabitants) (24). In 2010 and 2012, BC was the leading cancer diagnosis in Portuguese women (12, 22, 28) with a standardized IR (EU) of 93.2/105 inhabitants (28) [unadjusted and standardized IR world population: 118.5 and 69.8/105 inhabitants (28)] and 85.6/105 inhabitants (12) [age-adjusted IR 67.6/105 inhabitants (29)], placing Portugal on the highest quintile of BC-IR according to the World Health Organization (1). A European study also revealed that between 1990 and 2008, the incidence of BC in young Portuguese women (<40 years) grew, reaching an IR of 14.54/105 inhabitants (23). Additionally, evidence suggests that significantly more Portuguese BC cases are being diagnosed at a local stage, mainly due to screening actions (30).
In the periods 2000-2004 and 2005-2009, 59.4% and 63.3% of new Portuguese BC cases among women were diagnosed in the region covered by the South Regional Cancer Registry (ROR-Sul) (25), and in 2008 and 2010, this region registered 51.16% (27) and 49.41% (28) of BC incident cases among women in the country. This registry also showed the highest standardized BC-IR of the four active regional registries [2008: 98.3 vs. 86.6, 69.4 and 83.9/105 inhabitants in North, Centre and Azores (27); 2010: 97.6 vs. 91.7, 83.2 and 96.5/105 inhabitants, respectively (28)].
The study of the spatial distribution of diseases and its determinants, the analysis of spatiotemporal clusters (31, 32) and spatial variations in temporal trends (33) have been used to study geographic patterns of cancer (34) and are fundamental as public health (31, 32) statistical tools. Several studies conducted in China (35-37), Japan (38), Canada (39) and the USA (40, 41) detected intra-regional and temporal variations in BC incidence using these methods. This type of cluster analysis allowed the identification of areas with a higher risk of this cancer (36), provided a basis for the study of local risk factors [example: pollution (35)] and analyzed the role of the health services and lifestyles of the different populations in a defined area or city (35). These results can have public health implications, namely when it comes to resource allocation(36), planning of new epidemiological studies (40), risk assessment and health management (37).
BC is still a pathology with a high variability in terms of incidence evolution in different countries and regions (1, 3, 12, 22, 27, 28). Until now, however, only one study on BC incidence clustering has been conducted in Portugal (42). No further evidence was identified regarding the spatial-temporal variations of BC incidence in this country. Therefore, we conducted a spatiotemporal analysis of BC incidence between 2005 and 2012 in the south of Portugal, the region with the highest incidence and IR for this cancer, with the main goal of characterizing the spatiotemporal patterns of BC incidence among women in the period 2005-2012, considering four age groups and according to the stage at diagnosis.
Materials and Methods
Data and sources. A retrospective and observational study was carried out. The number of new BC cases diagnosed among women in 2005-2012 in the south of Portugal was obtained from the ROR-Sul. This database is registered in the National Commission for Data Protection (No. 1973500003, December 1997). For this specific study, Ethics Committee approval and informed consent were not required because the data were based on an official national surveillance system and had been previously anonymized. Six counties (out of 126) were not included in this study, as they only started reporting cancer cases to ROR-Sul in 2009. The median female population per year and per county was obtained from Statistics Portugal (SP).
Analytical methods. Firstly, a descriptive analysis of the number of incident BC cases/area/year was performed. Patients were divided into four age groups based on the literature: <40, 40-49, 50-64 and ≥65 years (5-7, 13-19, 21, 23). The annual IR was defined as the number of new BC cases among women in a certain area in a determined year, divided by the median population (at 30 June) of female inhabitants in that same area and year (20). Two spatial units were used in this study: the area covered by ROR-Sul, and the mainland counties covered by the registry (n=109). The annual and global BC-IR for southern Portuguese women was calculated as an unadjusted rate and was standardized according to the European and World standard population by indirect standardization.
Four different populations were used for the standardizations – two for Europe [EU1976 (43) and EU2010 (44)] and two for the world [World (43) and World2000-2015 (45)] to allow external comparisons (43) and reflect an age distribution closer to the current demographic pattern (44, 45). Furthermore, the global and yearly IR were determined as unadjusted rates according to the defined four age groups and BC stage at diagnosis (locoregional disease or metastatic disease). The identification of temporal and spatiotemporal clusters plus the detection of spatial variations in temporal trends was conducted using the Software for the Spatial and Space–Time Scan Statistics (SaTScan™, version 9.4.4; Martin Kulldorff, Boston, MA, USA), applying circular windows with a maximum of 20% of the studied population. For these analysis, only the 109 counties in the ROR-Sul mainland were considered. The median population density in the period in analysis was calculated for each of the counties and cut-off values of 150 and 300 inhabitants/km2 were applied for the definition of rural, intermediate and urban areas (46). The significance level was set at 0.05, and the IR was expressed per 105 inhabitants. The data were analyzed using Statistical Package for Social Sciences (SPSS™, version 22 for Windows; IBM Corp., Armonk, NY, USA), and the results were mapped using the software QGIS™ (version 2.18s; Free Software Foundation Inc., Boston, MA, USA).
New female breast cancer cases, unadjusted incidence rate (IR) and standardized incidence rate [(SIR) for the European (EU) and World Population], in the Southern Portugal Cancer Registry (ROR-Sul) per 105 inhabitants, in the period 2005-2012.
Results
ROR-Sul. This study included 23,202 women diagnosed with malignant BC (C50) between January 01, 2005 and December 31, 2012 living in the area covered by the ROR-Sul (n=120 counties). Most of the women were Portuguese (96.2%), and the mean age at diagnosis was 60.84 years [standard deviation (SD) 14.37, range=15-102, median=61.00 years]. At diagnosis, most patients presented locoregional disease (n=18,190, 78.4%), and 6.1% of the patients had metastatic disease (n=1,413). The unadjusted IR for 2005-2012 was 119.126/105 inhabitants, and this rate increased with age, peaking in the 65- to 69-year-old group. The unadjusted BC-IR showed some fluctuation in 2005-2012, ranging between 95.42 and 123.76/105 inhabitants. When standardized for the European and world standard populations, the IR presented lower values than in its unadjusted form, this reduction being more apparent when the new standard populations proposed for 2010 (44) or 2000-2025 (45), respectively, were used (Table I).
When considering the age groups <40, 40-49, 50-64 and ≥65 years, most of the BC cases in the period under analysis were diagnosed in the older group (n=9,555, 41.2%). The year that had the highest number of new BC cases was 2012, matching with the peak of new cases in the older group. The BC-IR was lower in the younger age groups (17.38, 155.60 and 213.65/105 inhabitants, respectively) than in the group aged ≥65 years (232.33/105 inhabitants). The IR of all the age groups showed a tendency to increase in the period under analysis. Additionally, the IR for locoregional BC (93.39/105 inhabitants) showed a growing pattern, whereas that for metastatic BC (7.25/105 inhabitants) seemed to be stable (Figure 1).
Analysis by county (n=109). In the spatiotemporal analyses, only the 109 counties in the mainland area covered by ROR-Sul were considered, corresponding to a total of 22,154 incident BC cases in the period 2005-2012. The Portuguese counties presenting the highest BC-IR differed over time. In each of the years 1.83% to 7.34% of the counties showed no new cases of BC among women; however, no counties presented no cases for the entire period (Table II).
Spatiotemporal analysis (n=109). The IR, temporal and spatiotemporal clusters identified for all the patients, the four age groups and the two groups of stage at diagnosis are further described in Table III and IV and Figure 2.
Spatiotemporal clustering analysis (Figure 2, a-f) identified several areas with high IR with diverse characteristics and dimensions, with Lisbon being the only area always highlighted in all these different settings. Additionally, when considering all patients, two statistically significant spatial variations in temporal trends were detected: Cluster I (7 rural counties) – Inland rural Alentejo: 215 cases; inside vs. outside trend: +19.495% vs. +3.972%/year; annual IR: 88.6/105 inhabitants) and Cluster II (3 rural counties) – inland rural Médio Tejo: 387 cases; inside vs. outside trend: +14.434% vs. +3.937%/year; annual IR: 95.6/105 inhabitants. Considering patients with locoregional disease at diagnosis, two spatial variations in temporal trends were identified: Cluster I (3 rural counties) – inland rural Médio Tejo: 213 cases; inside vs. outside trend: +26.739% vs. +6.753%/year; annual IR: 52.6/105 inhabitants) and Cluster II (2 urban counties) – urban seaside Almada-Seixal: 1,374 cases; inside vs. outside trend: +13.215% vs. +6.458%/year; annual IR: 100.4/105 inhabitants) (Figure 2 g-j). No spatial variations in temporal trends were identified for the remaining groups analyzed.
Discussion
To our knowledge, this study was the first spatiotemporal analysis of female BC incidence conducted in Portugal, providing a starting point for future studies conducted in this area. The BC-IR detected in this analysis was 119.13/105 inhabitants, corresponding to a standardized IR of 93.35 or 69.57/105 inhabitants when considering the EU1976 or EU2010 and 108.64 or 75.73/105 inhabitants when adjusting for the World or World2000-2025 standard populations. This rate is similar to the European median IR reported in 2012 of 94/105 inhabitants (4, 12). The difference in the standardized IR detected with the use of the EU1976 (43) and the EU2010 (44) populations points to a potential misalignment of the traditionally used standard populations with regards to the aging women in this population. The increase of the IR in the period of analysis (from 95.42 in 2005 to 123.76/105 inhabitants in 2012, or +4.113%/year for the mainland southern region) is aligned with the growing BC-IR described for developed countries (1, 2, 4, 6, 13-15) and Portuguese women (22-24).
Female breast cancer incidence rate (IR) unadjusted and standardized for the European population(43) per 105 inhabitants. in the Southern Portugal Cancer Registry region between 2005 and 2012: overall and according to the four age groups considered and stage at diagnosis.
Total and annual number of new breast cancer cases and incidence rate (per 105 inhabitants) in the 109 mainland counties covered by the Southern Portugal Cancer Registry in the period of 2005-2012.
Number of new breast cancer cases, annual incidence rate (IR), time trends and statistically significant temporal clusters (p<0.05) in the mainland counties covered by Southern Portugal Cancer Registry (n=109) between 2005 and 2012: overall, in the four age-groups, and according to the stage at diagnosis.
Breast cancer incidence in spatiotemporal clusters on the mainland counties covered by the Southern Portugal Cancer Registry (n=109) between 2005 and 2012 presenting statistical significance (p<0.05): overall, in the four age groups considered and according to the stage at diagnosis.
In this study, the BC-IR in the mainland counties of ROR-Sul increased with age (17.7, 156.9, 213.3 and 232.9/105 inhabitants respectively, in the four age groups studied), which is concordant with the literature (2, 4-9). All the mainland counties showed cases of BC, with most of the diagnosis occurring at the locoregional stage of the disease (77.57%) and in the elderly population (≥65 years: 41.24%). All age groups showed an increasing IR trend (+5.935%, +3.833%, +4.114% and +2.194%/year, respectively), which is concordant with the published international literature (6, 15, 17-21). This increase in IR points to a potential need for healthcare reorganization (25), especially when considering long-term survival in these patients (4, 14, 20, 25, 47) and their need to be followed-up by the national healthcare system (25). A statistically significant increase of this IR was further described in the years 2005 and 2009-2012 for patients diagnosed with metastatic and locoregional disease, respectively, when compared with the remaining period, matching previously published data for local BC(30). Additionally, locoregional cases of BC have been increasing at a higher rate than metastatic disease (by 6.976 vs. 0.303%/year), a fact than could be a result of screening actions (and early detection of cancer cases) (2, 4-6, 12, 13, 17, 21, 26, 30) or growing awareness of the population of this health problem (4, 6, 12, 13, 17, 18).
Spatiotemporal clusters of high breast cancer incidence rate [(IR) per 105] (a-f) and spatial variations in temporal trends (g-j) in the mainland counties covered by the Southern Portugal Cancer Registry (n=109) between 2005 and 2012: overall (a); at age 40-49 (b), 50-64 (c) and ≥65 (d) years; for those with locoregional disease (e), metastatic disease (f); all patients (g and i) and locoregional disease (h and j). Designations A-E correspond to high IR clusters given in Table IV. Designations I-II correspond to spatial variations in temporal trends clusters described in the text.
Some findings from this analysis should be investigated further, namely the increasing BC-IR trend in young women (<40 years: by 5.935%/year), as these cases of BC usually are more aggressive (18, 23) and are associated with a greater loss of potential years of life (18). Secondly, the Lisbon area was identified as a cluster of high incidence overall, in all the age groups analyzed (except women <40 years) and for both metastatic and locoregional disease at diagnosis. High IR clusters identified in this city potentially could be explained by the urbanity of the area (9, 48-50). On the other hand, the Alentejo and Algarve intermediate/rural regions showed clusters with lower IR than the average of the studied area for all the age groups considered. These findings might point to the presence of environmental (4, 13, 16, 49, 51), genetic (4, 51) or lifestyle factors (4, 13, 49, 51) in the capital that predispose patients to the development of this cancer. This excess IR in urban versus rural settings can also be explained by the number of registry hospitals in the urban centers (49) and the fact that patients prefer to be treated in ‘renowned’ hospitals. This could lead to the use of addresses of family and friends to facilitate access to those facilities (49), inflating the number of cases diagnosed in urban hospitals. The low number of inhabitants could account for the high overall BC-IR cluster detected in isolated rural counties (small numbers), whereas the high IR cluster reported for locoregional BC for the intermediate regions matched the previously described trend for BC-IR in Portuguese regions (22).
The spatial variations in temporal trends in the inland rural cluster Médio Tejo showed a dramatic increase in BC-IR (overall and for those with locoregional disease) in 2005-2012, reaching the values of the regional trend in the last year of the period. This ‘exponential’ growth in overall BC cases may be due to a sub-detection in the beginning of the period, as suggested by the spatiotemporal cluster Médio Tejo (relative risk=0.28) or higher adherence to screening actions in these years. In addition, the high incidence spatiotemporal cluster in Lezíria do Tejo for patients with locoregional disease (detected at the end of the analyzed period) includes the three counties highlighted by this spatial variation in temporal trend analysis, showing the complementary nature of these methods for the detection of critical areas of BC incidence. The Alentejo inland rural cluster also presented a variation in the BC-IR dissimilar to that of the region, with two peaks of incidence being described in 2009 and 2011. This phenomenon could be a result of more intense screening actions in these years.
The urban seaside Almada-Seixal cluster showed a trend significantly different from the remaining area in analysis, presenting a higher incidence in all the years studied (except 2007), a fact that might be explained by environmental factors or lifestyle habits in these counties. These counties were also identified as part of the high incidence spatiotemporal clusters for all patients and for those ≥65 years, demonstrating that this area is critical for BC incidence. Further studies should be carried out to further clarify these findings and to assess if the growing IR trend persisted after the period studied, especially considering that these clusters were detected in the later years of the analysis.
This study presents some limitations, namely the fact that the whole country was not studied and that some cases may have been missed or been recorded in another registry. Additionally, the delay in the public presentation of these data might compromise the relevance of this information for evidence-based public health policies. However, this study provides an innovative, detailed and clear insight into the increase of the BC burden in this population and on the heterogeneous distribution of this type of cancer.
Southern Portuguese women have shown a growing BC-IR; however, the trend for this increase varied according to the age, stage at diagnosis, and area of residence. The reason underlying these different patterns could be related to different distribution (6, 12, 17), screening methods (2, 4, 6, 12, 17-20) and increase of early diagnosed cases (11). Altogether, these data support the uneven distribution of BC incidence reported previously in different Portuguese regions (14, 22, 25, 27, 28) and the fundamental need for further studies on this topic for a clearer evaluation of potential inequalities in BC incidence.
Acknowledgements
The Authors thank the ROR-Sul registry for kindly providing the data for this study.
Footnotes
Funding and Conflicts of Interests
This research did not receive any specific grant from funding agencies in the public, commercial, or not-for-profit sectors. There is no conflict of interests to declare.
- Received December 22, 2017.
- Revision received January 24, 2018.
- Accepted January 30, 2018.
- Copyright© 2018, International Institute of Anticancer Research (Dr. George J. Delinasios), All rights reserved