Skip to main content
  • Research article
  • Open access
  • Published:

Spatio-temporal mapping of breast and prostate cancers in South Iran from 2014 to 2017

Abstract

Background

The most common gender-specific malignancies are cancers of the breast and the prostate. In developing countries, cancer screening of all at risk is impractical because of healthcare resource limitations. Thus, determining high-risk areas might be an important first screening step. This study explores incidence patterns of potential high-risk clusters of breast and prostate cancers in southern Iran.

Methods

This cross-sectional study was conducted in the province of Kerman, South Iran. Patient data were aggregated at the county and district levels calculating the incidence rate per 100,000 people both for cancers of the breast and the prostate. We used the natural-break classification with five classes to produce descriptive maps. A spatial clustering analysis (Anselin Local Moran’s I) was used to identify potential clusters and outliers in the pattern of these cancers from 2014 to 2017.

Results

There were 1350 breast cancer patients (including, 42 male cases) and 478 prostate cancer patients in the province of Kerman, Iran during the study period. After 45 years of age, the number of men with diagnosed prostate cancer increased similarly to that of breast cancer for women after 25 years of age. The age-standardised incidence rate of breast cancer for women showed an increase from 29.93 to 32.27 cases per 100,000 people and that of prostate cancer from 13.93 to 15.47 cases per 100,000 during 2014–2017. Cluster analysis at the county level identified high-high clusters of breast cancer in the north-western part of the province for all years studied, but the analysis at the district level showed high-high clusters for only two of the years. With regard to prostate cancer, cluster analysis at the county and district levels identified high-high clusters in this area of the province for two of the study years.

Conclusions

North-western Kerman had a significantly higher incidence rate of both breast and prostate cancer than the average, which should help in designing tailored screening and surveillance systems. Furthermore, this study generates new hypotheses regarding the potential relationship between increased incidence of cancers in certain geographical areas and environmental risk factors.

Peer Review reports

Background

Cancers are the second leading cause of death worldwide [1], which can partly be explained by the fact that the world’s population is ageing [2]. Furthermore, human exposure to multiple risk factors has increased the cancer burden worldwide [3]. Despite advances in timely diagnosis and medical treatment of neoplasms in recent years, malignancies in middle to low-income countries are expected to almost double by 2030 compared to high-income nations [4]. If cancer is diagnosed promptly, cures can sometimes be found and life prolonged leading to considerably lower disease burdens [5, 6]. However, health systems, particularly in developing countries, are not capable of screening all people to identify patients in the early stages of the disease. Identifying high-risk geographical areas could help decreasing the cost of screening, finding the people at risk and implementing more efficient diagnostic strategies [7]. Investigating high-risk areas should also provide valuable knowledge to scientists about the aetiology of some malignancies [3].

Cancer of the breast and the prostate are the two most common, gender-specific malignancies worldwide [8]. Furthermore, these neoplasms cause high numbers of disability-adjusted life years (DALYs) [9]. Risk factors for these two diseases are diverse and interrelated, as they include genetic [10], social-economic [11] as well as lifestyle and environmental factors [12]. Further, there are interactions between these risk factors, particularly with those involving the environment [13, 14], whose spatial variation may lead to heterogeneity in the pattern of cancers in a given geographic catchment area. Studies by Wang et al. [15, 16] found a significant spatial variation of prostate cancer incidence and prostate cancer-specific mortality in Pennsylvania, USA. They evaluated potential effects of individual and county-level risk factors and found that spatial variations in prostate cancer-specific mortality rates existed in Pennsylvania with a particularly high risk in the Pen State catchment area. County-level health and environmental factors might contribute to spatial heterogeneity in prostate cancer-specific mortality as shown by Olfatifar et al. [17], who examined spatial clustering of breast cancer at the provincial level in Iran between 2004 and 2010. Their results highlight that the breast tumour incidence varied across the provinces [17]. At the same time, Rohani-Rasaf et al. [18] detected some high-risk regions in Tehran, the capital of Iran, both for cancers of the breast and the prostate. Most studies in Iran have applied spatial analyses at a very coarse level (province scale) and the results are therefore not as useful as a finer scale quite naturally.

Geographical information systems (GIS) constitute a set of useful tools for the identification of high-risk areas of cancer occurrence as well as investigation of the environmental effects on cancer incidence [19,20,21]. GIS approaches combine spatial and non-spatial data producing geodatabases that make it possible to perform spatial analyses using this data structure [22, 23]. For example, spatial autocorrelation is a method of exploratory data analysis which allows detecting spatial data dependence [24]. There are two kinds of spatial autocorrelation methods: global and local statistics. Global methods are more sensitive to departures from the null hypothesis, which examine whether data, here patients, are randomly distributed or if there is a spatial pattern. They can identify spatial structures in the pattern of cancer incidence but do not determine where the clusters are. Local cluster statistics, on the other hand, can quantify spatial autocorrelation and clustering, but only in limited areas. These methods may find restricted areas characterized as high-high (HH), high-low (HL), low-low (LL) or low-high (LH) risk of incidence within a region. HH and LL are defined as target areas surrounded by areas with similar incidence rates, while for HL and LH, the target areas are surrounded by areas with dissimilar cancer incidence rates. In other words, HH and LL indicate clusters, while HL and LH point to outliers [25]. This study aimed to identify the spatial patterns of cancer of the breast and the prostate and to investigate the potential clustering in gender-specific patterns of these cancers in southern Iran between 2014 and 2017.

Method

Study area and time

This study was conducted in the province of Kerman, located in southern Iran (Fig. 1). The first administrative level of Iran subdivisions is the province, each of which is further divided into counties that are in turn divided into districts. Our study area contained the 22 counties and 58 districts of Kerman Province, which covers an area of 183,285 km2 and has, according to the National Census of 2015, a population of 3,164,718 people [26]. The study covered the time span of March 2014–March 2017.

Fig. 1
figure 1

Map of Kerman counties and districts including the distribution of cancer (prostate and breast) patients during 2014–2017. The figure was created by authors

Data sources

Data were obtained from two different sources with three different spatial scales (individual, county and district). The individual patient data (supplementary files 12) were obtained through the population-based cancer registry of Kerman. They were geocoded and aggregated to both county and district. The digital maps (county and district) were obtained through the mapping organisation of the country.

We used a crude incidence per 100,000 people and the age-standardised rate (ASR) per 100,000 people for the descriptive statistics.

Inclusion and exclusion criteria

We included all residents of Kerman diagnosed with either cancer of the breast or the prostate. Individuals who had come to the province for cancer treatment but lived outside of the province were excluded.

Spatial analyses

For the thematic maps, we used the natural-break classification with five classes. This approach is a data-clustering method designed to determine the best arrangement of values into different classes. This is done by seeking to minimise each class’s average deviation from the class mean, while maximising each class’s deviation from the means of the other groups. In other words, the method seeks to reduce the variance within classes and maximise the variance between them [27]. For spatial visualisation, the crude incidence per 100,000 people and the ASR per 100,000 people were used.

Spatial Cluster Analysis

Incidence rates of the two target cancers were calculated using total population and number of cases in each county and district of the province. The Local Moran’s I statistic [28] was performed to quantify spatial autocorrelation of cancers frequency at county and district level. This test calculates a z-score and p-value to determine whether the apparent similarity (spatial clustering of either high or low values) or dissimilarity (presence of spatial outliers) is more pronounced than one would expect in a random distribution. The null hypothesis states that the cancers are randomly distributed across the study area. A high positive z-score for a feature indicates that the surrounding features have similar values (either high values or low values). However, a low negative z-score for a feature indicates a statistically significant spatial data outlier [28]. We used a 95% Confidence Level (CL), and all clusters and outliers found in this study were significant at this CL.

We used ArcGIS, v. 10.5 (ESRI, Redlands, CA, USA) and GeoDa (https://spatial.uchicago.edu/geoda) for spatial analyses and Microsoft Excel 2016 for the descriptive analyses.

Results

There were 1350 breast cancer patients (including, 42 male cases) and 478 prostate cancer patients in the province of Kerman, Iran during the period March 2014–March 2017. Table 1 shows the crude incidence per 100,000 people and ASR per 100,000 people of these cancers. The ASR of both breast and prostate cancers increased by 29.93 to 32.27 and by 13.93 to 15.47 from 2014 to 2017, respectively.

Table 1 Breast and prostate cancer incidence per 100,000 people during different time periods in the province of Kerman, Iran

As shown in Fig. 2, the number of women (after 25 years of age) who developed breast cancer increased rapidly. Furthermore, the highest incidence occurred in the 50–54 age group in the 2014–2015 period, in the 65–69 age group in the 2015–2016 period and in the 75–79 age group in the 2016–2017 period.

Fig. 2
figure 2

Age distribution of breast cancer patients in the province of Kerman, Iran

After 45 years of age, the number of men with diagnosed prostate cancer increased similarly to that of breast cancer for women after 25 years of age. In contrast, however, the ascent of this cancer in relation to seniority was considerable in all years from 2014 to 2017 (Fig. 3).

Fig. 3
figure 3

Age distribution of prostate cancer patients in the province of Kerman, Iran

The descriptive maps of Fig. 4 reveal that the breast cancer incidence was highest in the north-eastern part of the province from 2014 to 2017. However, as the cluster maps show, there were HH clusters of breast cancer in the north-western part of the province from 2014 to 2017; furthermore, there was a LL cluster in the South-East in 2014–2015 and in 2016–2017, but not in the time between these periods. All clusters and outliers were statistically significant (p < 0.05).

Fig. 4
figure 4

Breast cancer incidence map at the county level in the province of Kerman, Iran. The figure was created by authors

The descriptive maps of Fig. 5 reveal that the breast cancer incidence was highest in the North stretching towards the centre of the province. However, as the cluster maps show, there were HH clusters of breast cancer in the north-western part of the province from 2014 to 2015 and 2016 to 2017; furthermore, there were LL clusters in the north-eastern part of the province from 2014 to 2015 and from 2016 to 2017. All clusters and outliers in the figure were statistically significant (p < 0.05).

Fig. 5
figure 5

Breast cancer incidence map at the district level in the province of Kerman, Iran. The figure was created by authors

The descriptive maps of Fig. 6 reveal that the prostate cancer incidence was highest in the North-East. However, as the cluster maps show, there were HH clusters of prostate cancer in the north-western part of the province from 2014 to 2016; furthermore, there were LL clusters of prostate cancer in the South-East and East of the province from 2014 to 2017. All clusters and outliers in the figure were statistically significant (p < 0.05).

Fig. 6
figure 6

Prostate cancer incidence map at the county level in the province of Kerman, Iran. The figure was created by authors

The descriptive maps of Fig. 7 reveal that the prostate cancer incidence was highest in the North-East of the province from 2015 to 2016. However, there was a HH cluster of prostate cancer in the North-West from 2014 to 2015 and some HH clusters in the North-West stretching towards the centre from 2016 to 2017. All clusters and outliers in the figure were statistically significant (p < 0.05).

Fig. 7
figure 7

Prostate cancer incidence map at the district level in the province of Kerman, Iran. The figure was created by authors

The four categories shown in Figs. 4, 5, 6, 7 correspond to the four quadrants in the Moran scatter plots as shown in Fig. 8. If nearby or neighbouring areas are more alike, this is understood as positive spatial autocorrelation. Negative autocorrelation describes patterns in which neighbouring areas are unlike. For example, the first scatter plot (which shows the county level breast cancer for 2014–2015) has a Moran’s value of 0.271 that should be interpreted as a spatial pattern with a cluster tendency (HH and LL). Note that the incidence values have been standardized and are given in standard deviational units (the mean is zero and the standard deviation is 1). Similarly, the spatial lag was computed for those standardized values.

Fig. 8
figure 8

Moran’s scatter plots for breast and prostate cancers cluster maps

Discussion

The main aim of this study was to explore the spatio-temporal patterns of the incidence of breast and prostate cancer at a high-resolution geography scale. To the best of our knowledge, this is the first study to assess spatial variations in the incidence pattern of breast and prostate cancers in Iran. The study area in the South of Iran revealed a high incidence rate of both cancers in north-western Kerman, while it was low in the south-eastern part of the province. The number of people with breast or prostate cancer increased considerably after patients reached 25 and 45 years of age, respectively. Further investigations are needed to assess the drivers in the high-risk areas identified in north-western Kerman. They might be associated with environmental factors and lifestyles [12], poor access to cancer-specific services [29], hereditary reasons [10] and/or socio-economic inequalities [30, 31].

Environmental risk factors such as air pollution [32,33,34] and presence of heavy metals [35, 36] could be linked to the geographic outcome disparities for both cancer forms. We found high-risk areas for the whole study period of breast and prostate cancers in the north-western part of Kerman, an area extending southeast of the Iranian volcano-plutonic copper belt [37] where arsenic contamination is one of the most significant environmental concerns [38]. This is the location of the Sarcheshmeh copper industrial plant, the biggest copper mine of Iran, and it is likely that the neighbourhood is contaminated also by other heavy metals. Field studies report widely distributed travertine rocks north of the Sarcheshmeh copper mine indicating the presence of a highly concentrated range of arsenic compounds [38], which could seep into the water system and contaminate the drinking water in nearby urban and rural communities [38, 39]. Indeed, the arsenic concentration in the water in these areas is higher than the limit recommended by the World Health Organisation (WHO) for drinking water [38, 40, 41]. Arsenic has been categorised as a Group 1 carcinogen factor by the International Agency for Research on Cancer (IARC) [42] and various studies associate arsenic and breast cancer [43,44,45,46]; its presence in the study area is thus a potential explanation for the increased incidence of breast cancer found. However, it should also be mentioned that other studies do not show any significant association between breast cancer and arsenic [47, 48]. However, it is conceivable that he power of this association can change due to local and individual diversities [45].

Previous studies indicate a significant association between arsenic-enriched water and prostate cancer incidence [46, 49, 50], while increased levels of copper has been linked with the initiation of prostate cancer [35]. Copper smelting and toxic discharges have led to soil pollution, especially in the region of the smelting plant in Sarcheshmeh Copper Complex. Importantly, the most contaminated areas are located in the most common wind directions [51], and it is particularly disturbing that the polluted areas are also used as grazing land enabling toxic elements from this soil to enter the food chain. These elements include various heavy metals in addition to copper and arsenic, e.g., lead, molybdenum and cadmium [51]. Therefore, soil, water and nutrition of Rafsanjan and the adjacent townships, located in the north-western part of the province are subject to these potential negative effects. Indeed, previous studies have found associations between heavy metals and both breast and prostate cancers [52,53,54]. The current study strongly recommends examining the hypothesis that exposure to heavy metals, especially arsenic and copper, may be associated with high incidence of gender-specific cancers. In fact, the high incidence of both breast and prostate cancers in north-western Kerman may be associated with these contaminants and this need to be investigated in future studies.

Air pollution, such as that due to particulate matter, has been shown to be associated with breast cancer [33, 34]. Further studies are suggested to confirm the effects of airborne pollution on the risk of breast cancer, especially particulate matter. Fazzo et al. (2016) used a spatial approach to estimate the industrial air pollution impact on 17 selected neoplasms in the territory around the industrial Sicilian area of Priolo, Italy. This area has been defined as a contaminated site of national priority for remediation because of diffuse environmental contamination caused by large industrial settlements, and their study found a higher incidence of breast cancer in the contaminated area compared to the rest of the province [55].

Previous studies highlight that poor access to health care services, such as increased incidence of cancer [29, 56], lead to poor health outcomes [57, 58]. The high incidence of gender-specific cancers in some regions of the study area may be due to their considerable distance from the provincial capital with limited cancer screening programmes. On the other hand, parts of the study area in the South had the lowest incidence of both cancers investigated here although those affected were located even further away from the provincial capital. However, proximity to health care services does not directly translate into access because of potential factors such as poor socio-economic status and low level of education that also are associated with poor access [57]. GIS enable researchers to assess the revealed access to cancer services through combing spatial and non-spatial factors [58,59,60] and the results suggest measuring access to cancer prevention programmes should be the first step when examining this hypothesis. Previous studies have highlighted the impact of the socio-economic status on the differences in the incidence of cancers of the breast and prostate [61,62,63]. Assessing the impact of socio-economic status on the geographic disparities of the gender-specific cancers incidence in the study area can be done by analysing the overall spatial structure or identifying high-risk areas. This also warrants further studies.

Hereditary cancer syndromes, a type of inherited disorder in which there is a higher-than-normal risk of certain types of cancer, are caused by mutations in certain genes passed from parents to children [64,65,66,67,68]. Certain such family-related cancers are well-known, e.g., hereditary breast cancer [67] and Lynch syndrome, which is a hereditary non-polyposis colorectal cancer [65]. Hereditary cancer screening programmes [69,70,71] have made it possible to detect many of the approximately 5–10% of breast cancers caused by a genetic predisposition [72, 73] thus making it possible to prevent them before they occur. There are also studies assessing the risk of prostate cancer associated with hereditary cancer syndromes. This highlights the risk of prostate cancer in members of families associated with early-onset diseases of various kinds [68]. We strongly recommend researchers and policymakers to perform hereditary cancer screening and genetic testing in areas of the province with a high incidence rate of cancer of either breast or prostate.

Spatio-temporal cluster analysis plays a significant role in visualising and quantifying geographical variation in patterns of disease incidence. Global Moran’s I and Getis-Ord General G statistic are both global cluster methods which can be used to investigate the level of spatial autocorrelation of disease patterns, while Local Moran’s I and Getis-Ord Gi* are indicate the location of the clusters. Although Getis-Ord Gi* statistic is used for identifying hotspots and the opposite, Local Moran’s I is also effective for assessing statistically significant spatial outliers [74] and has therefore been predominantly applied [28, 75] and successfully assessed the hotspots [76, 77]. Those using these methods for analysis of the spatial pattern of incident data should consider aggregating the incident data into polygons. The main question here is the geographical scale that should be used for aggregation because it could affect the results. In this study, we conducted the analyses both at the county and the district level, which are geographical scales providing practically useful information. We recommend other researchers performing cohort studies examining the hypotheses proposed here.

Limitations and future course of the study

We included all patients with breast and prostate cancer in the province of Kerman during 2014–2017. However, we did not have the migration data and used the current address of patients at the time of cancer diagnosis as the patients’ residence. However, some patients might have lived in other regions in the years that might have affected their health.

Conclusions

We identified a great deal of spatial variations with significant clusters in the patterns of cancer of the breast and the prostate. This suggests that policymakers need to develop prevention strategies tailored to areas where the risk of these conditions are greater than elsewhere. Further, there is a need to conduct further research to test the causal relationship between environmental risk factors and cancer incidence.

Availability of data and materials

The cancers data have been uploaded as supplementary files. However, due to protect the patients’ data, the latitude/longitude of patients’ location has been removed.

Abbreviations

GIS:

Geographical information systems

HH:

High-high

LL:

Low-low

HL:

High-low

LH:

Low-high

WHO:

World health Organization

IARC:

International agency for research on cancer

References

  1. Wang H, Naghavi M, Allen C, Barber RM, Bhutta ZA, Carter A, et al. Global, regional, and national life expectancy, all-cause mortality, and cause-specific mortality for 249 causes of death, 1980–2015: a systematic analysis for the global burden of disease study 2015. Lancet. 2016;388(10053):1459–544.

    Article  Google Scholar 

  2. Kennedy BK, Berger SL, Brunet A, Campisi J, Cuervo AM, Epel ES, et al. Geroscience: linking aging to chronic disease. Cell. 2014;159(4):709–13.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  3. Danaei G, Vander Hoorn S, Lopez AD, Murray CJ, Ezzati M. Group CRAc. Causes of cancer in the world: comparative risk assessment of nine behavioural and environmental risk factors. Lancet. 2005;366(9499):1784–93.

    Article  PubMed  Google Scholar 

  4. World Health Organization. Cancer key statistics. Cited on 10 November 2020. Available from https://www.who.int/cancer/resources/keyfacts/en/.

  5. Sankaranarayanan R, Swaminathan R, Brenner H, Chen K, Chia KS, Chen JG, et al. Cancer survival in Africa, Asia, and Central America: a population-based study. Lancet Oncol. 2010;11(2):165–73.

    Article  PubMed  Google Scholar 

  6. Rajaraman P, Anderson BO, Basu P, Belinson JL, D'Cruz A, Dhillon PK, et al. Recommendations for screening and early detection of common cancers in India. Lancet Oncol. 2015;16(7):e352–e61.

    Article  PubMed  Google Scholar 

  7. Lofters AK, Gozdyra P, Lobb R. Using geographic methods to inform cancer screening interventions for south Asians in Ontario, Canada. BMC Public Health. 2013;13(1):395.

    Article  PubMed  PubMed Central  Google Scholar 

  8. Ghoncheh M, Pournamdar Z, Salehiniya H. Incidence and mortality and epidemiology of breast cancer in the world. Asian Pac J Cancer Prev. 2016;17(S3):43–6.

    Article  PubMed  Google Scholar 

  9. Murillo-Zamora E, Mendoza-Cano O, Ríos-Silva M, Sánchez-Piña RA, Higareda-Almaraz MA, Higareda-Almaraz E, et al. Disability-adjusted life years for Cancer in 2010–2014: a regional approach in Mexico. Int J Environ Res Public Health. 2018;15(5):864.

    Article  PubMed Central  Google Scholar 

  10. Tedaldi G, Tebaldi M, Zampiga V, Danesi R, Arcangeli V, Ravegnani M, et al. Multiple-gene panel analysis in a case series of 255 women with hereditary breast and ovarian cancer. Oncotarget. 2017;8(29):47064.

    Article  PubMed  PubMed Central  Google Scholar 

  11. Hemminki K, Zhang H, Czene K. Socioeconomic factors in cancer in Sweden. Int J Cancer. 2003;105(5):692–700.

    Article  CAS  PubMed  Google Scholar 

  12. Kulhánová I, Znaor A, Shield KD, Arnold M, Vignat J, Charafeddine M, et al. Proportion of cancers attributable to major lifestyle and environmental risk factors in the eastern Mediterranean region. Int J Cancer. 2020;146(3):646–56.

    Article  PubMed  CAS  Google Scholar 

  13. Coyle YM. The effect of environment on breast cancer risk. Breast Cancer Res Treat. 2004;84(3):273–88.

    Article  CAS  PubMed  Google Scholar 

  14. Gann PH. Risk factors for prostate cancer. Rev Urol. 2002;4(Suppl 5):S3.

    PubMed  PubMed Central  Google Scholar 

  15. Wang M, Wasserman E, Geyer N, Carroll RM, Zhao S, Zhang L, et al. Spatial patterns in prostate Cancer-specific mortality in Pennsylvania using Pennsylvania Cancer registry data, 2004–2014. BMC Cancer. 2020;20(1):394.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Wang M, Matthews SA, Iskandarani K, Li Y, Li Z, Chinchilli VM, et al. Spatial-temporal analysis of prostate cancer incidence from the Pennsylvania Cancer registry, 2000-2011. Geospat Health. 2017;12(2):611.

    PubMed  PubMed Central  Google Scholar 

  17. Olfatifar M, Karami M, Moghimbeigi A, Motlagh A, Rooshanaee G, et al. Spatial clustering of breast cancer: an epidemiological analysis of Iranian Women. Int J Cancer Manag. 2017;10(1):e5402. https://doi.org/10.17795/ijcp-5402.

  18. Rohani-Rasaf M, Rohani-Rasaf M R, Hashemi Nazari S S, Mohammadian-Hafshejani A, Asadi-Lari M. Cancer Incidence Rate in the Elderly Inhabitants of Tehran: Is there Really any Cluster? Int J Cancer Manag. 2017;10(7):e5753. https://doi.org/10.5812/ijcm.5753.

    Article  Google Scholar 

  19. Sahar L, Foster SL, Sherman RL, Henry KA, Goldberg DW, Stinchcomb DG, et al. GIScience and cancer: state of the art and trends for cancer surveillance and epidemiology. Cancer. 2019;125(15):2544–60.

    PubMed  Google Scholar 

  20. Graves BA. Integrative literature review: a review of literature related to geographical information systems, healthcare access, and health outcomes. Perspect Health Inf Manag. 2008;5:11. PMID: 18698429.

  21. Goshayeshi L, Pourahmadi A, Ghayour-Mobarhan M, Hashtarkhani S, Karimian S, Shahhosein Dastjerdi R, et al. Colorectal cancer risk factors in north-eastern Iran: A retrospective cross-sectional study based on geographical information systems, spatial autocorrelation and regression analysis. Geospat Health. 2019;14(2):219–28.

    Article  Google Scholar 

  22. Hashtarkhani S, Kiani B, Bergquist R, Bagheri N, VafaeiNejad R, Tara M. An age-integrated approach to improve measurement of potential spatial accessibility to emergency medical services for urban areas. Int J Health Plann Manag. 2020;35(3):788–98.

    Article  Google Scholar 

  23. Kiani B, Bagheri N, Tara A, Hoseini B, Tara M. Haemodialysis services in the northeastern region of Iran. Geospat Health. 2017;12(1):130–6.

    Article  Google Scholar 

  24. Halimi L, Bagheri N, Hoseini B, Hashtarkhani S, Goshayeshi L, Kiani B. Spatial analysis of colorectal Cancer incidence in Hamadan Province, Iran: a retrospective cross-sectional study. Appl Spat Anal Policy. 2020;13(2):293–303.

    Article  Google Scholar 

  25. Bagheri N, Furuya-Kanamori L, Doi SA, Clements AC, Sedrakyan A. Geographical outcome disparities in infection occurrence after colorectal surgery: an analysis of 58,096 colorectal surgical procedures. Int J Surg. 2017;44:117–21.

    Article  PubMed  Google Scholar 

  26. Iran SCo. Statistical Center of Iran 2018 [cited 2018 16/6/2018]. Available from: https://www.amar.org.ir/english.

  27. McMaster R. In memoriam: George F. Jenks (1916-1996). Cartography Geogr Inf Syst. 1997;24(1):56–9.

    Article  Google Scholar 

  28. Anselin L. Local indicators of spatial association—LISA. Geogr Anal. 1995;27(2):93–115.

    Article  Google Scholar 

  29. Bhat D, Heiman AJ, Talwar AA, Dunne M, Amanjee K, Ricci JA. Access to breast Cancer treatment and reconstruction in rural populations: do women have a choice? J Surg Res. 2020;254:223–31.

    Article  PubMed  Google Scholar 

  30. Akinyemiju TF, Demb J, Izano MA, Rehkopf DH, Fang ML, Hiatt RA, et al. The association of early life socioeconomic position on breast cancer incidence and mortality: a systematic review. Int J Public Health. 2018;63(7):787–97.

    Article  PubMed  Google Scholar 

  31. Meijer M, Bloomfield K, Engholm G. Neighbourhoods matter too: the association between neighbourhood socioeconomic position, population density and breast, prostate and lung cancer incidence in Denmark between 2004 and 2008. J Epidemiol Community Health. 2013;67(1):6–13.

    Article  PubMed  Google Scholar 

  32. Hwang J, Bae H, Choi S, Yi H, Ko B, Kim N. Impact of air pollution on breast cancer incidence and mortality: a nationwide analysis in South Korea. Sci Rep. 2020;10(1):5392.

    Article  PubMed  PubMed Central  Google Scholar 

  33. White AJ, Keller JP, Zhao S, Carroll R, Kaufman JD, Sandler DP. Air pollution, clustering of particulate matter components, and breast cancer in the sister study: A U.S.-wide cohort. Environ Health Perspect. 2019;127(10):107002.

    Article  PubMed  PubMed Central  Google Scholar 

  34. Cheng I, Tseng C, Wu J, Yang J, Conroy SM, Shariff-Marco S, et al. Association between ambient air pollution and breast cancer risk: the multiethnic cohort study. Int J Cancer. 2020;146(3):699–711.

    Article  CAS  PubMed  Google Scholar 

  35. Saleh SAK, Adly HM, Abdelkhaliq AA, Nassir AM. Serum levels of selenium, zinc, copper, manganese, and Iron in prostate Cancer patients. Curr Urol. 2020;14(1):44–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Lim JT, Tan YQ, Valeri L, Lee J, Geok PP, Chia SE, et al. Association between serum heavy metals and prostate cancer risk - a multiple metal analysis. Environ Int. 2019;132:105109.

    Article  CAS  PubMed  Google Scholar 

  37. Nedimovic R. Exploration for ore deposits in Kerman region. Geological Survey of Iran Report 53. 1973:247.

  38. Khorasanipour M, Esmaeilzadeh E. Geo-genic arsenic contamination in the Kerman Cenozoic magmatic arc, Kerman, Iran: implications for the source identification and regional analysis. Appl Geochem. 2015;63:610–22.

    Article  CAS  Google Scholar 

  39. Abbasnejad A, Mirzaie A, Derakhshani R, Esmaeilzadeh E. Arsenic in groundwaters of the alluvial aquifer of Bardsir plain, SE Iran. Environ Earth Sci. 2013;69(8):2549–57.

    Article  CAS  Google Scholar 

  40. Ebrahimi M. Evaluation of Arsenic dispersivity and the survey of the source of Arsenic in Groundwaters of Rafsanjan plain and represent the suitable ways of removal arsenic: M. Sc. Thesis (in Persian), Shahid Bahonar University of Kerman, p 203; 2009.

  41. Williams M. Arsenic in mine waters: an international study. Environ Geol. 2001;40(3):267–78.

    Article  CAS  Google Scholar 

  42. Straif K, Benbrahim-Tallaa L, Baan R, Grosse Y, Ba S, El Ghissassi F, et al. A review of human carcinogens--part C: metals, arsenic, dusts, and fibres. Lancet Oncol. 2009;10(5):453.

    Article  PubMed  Google Scholar 

  43. López-Carrillo L, Hernández-Ramírez RU, Gandolfi AJ, Ornelas-Aguirre JM, Torres-Sánchez L, Cebrian ME. Arsenic methylation capacity is associated with breast cancer in northern Mexico. Toxicol Appl Pharmacol. 2014;280(1):53–9.

    Article  PubMed  CAS  Google Scholar 

  44. Pineda-Belmontes CP, Hernández-Ramírez RU, Hernández-Alcaraz C, Cebrián ME, López-Carrillo L. Genetic polymorphisms of PPAR gamma, arsenic methylation capacity and breast cancer risk in Mexican women. Salud Publica Mex. 2016;58:220–7.

    Article  PubMed  Google Scholar 

  45. Khanjani N, Jafarnejad A-B, Tavakkoli L. Arsenic and breast cancer: a systematic review of epidemiologic studies. Rev Environ Health. 2017;32(3):267–77.

    Article  CAS  PubMed  Google Scholar 

  46. Bardach AE, Ciapponi A, Soto N, Chaparro MR, Calderon M, Briatore A, et al. Epidemiology of chronic disease related to arsenic in Argentina: a systematic review. Sci Total Environ. 2015;538:802–16.

    Article  CAS  PubMed  Google Scholar 

  47. Aballay LR, Díaz MP, Francisca FM, Muñoz SE. Cancer incidence and pattern of arsenic concentration in drinking water wells in Córdoba, Argentina. Int J Environ Health Res. 2012;22(3):220–31.

    Article  CAS  PubMed  Google Scholar 

  48. Garland M, Morris JS, Colditz GA, Stampfer MJ, Spate VL, Baskett CK, et al. Toenail trace element levels and breast cancer: a prospective study. Am J Epidemiol. 1996;144(7):653–60.

    Article  CAS  PubMed  Google Scholar 

  49. Bulka CM, Jones RM, Turyk ME, Stayner LT, Argos M. Arsenic in drinking water and prostate cancer in Illinois counties: an ecologic study. Environ Res. 2016;148:450–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  50. Roh T, Lynch CF, Weyer P, Wang K, Kelly KM, Ludewig G. Low-level arsenic exposure from drinking water is associated with prostate cancer in Iowa. Environ Res. 2017;159:338–43.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  51. Rastmanesh F, Moore F, Kopaei MK, Keshavarzi B, Behrouz M. Heavy metal enrichment of soil in Sarcheshmeh copper complex, Kerman. Iran Environ Earth Sci. 2011;62(2):329–36.

    Article  CAS  Google Scholar 

  52. Julin B, Wolk A, Johansson J-E, Andersson S-O, Andrén O, Åkesson A. Dietary cadmium exposure and prostate cancer incidence: a population-based prospective cohort study. Br J Cancer. 2012;107(5):895–900.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  53. Gallagher CM, Chen JJ, Kovach JS. Environmental cadmium and breast cancer risk. Aging (Albany NY). 2010;2(11):804.

    Article  CAS  Google Scholar 

  54. Romaniuk А, Lyndin M, Sikora V, Lyndina Y, Romaniuk S, Sikora K. Heavy metals effect on breast cancer progression. Journal of Occupational Medicine and Toxicology. 2017;12(1):32.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  55. Fazzo L, Carere M, Tisano F, Bruno C, Cernigliaro A, Cicero MR, et al. Cancer incidence in Priolo, Sicily: a spatial approach for estimation of industrial air pollution impact. Geospat Health. 2016;11(1):43–55.

    Article  Google Scholar 

  56. Morris LG, Sikora AG, Tosteson TD, Davies L. The increasing incidence of thyroid cancer: the influence of access to care. Thyroid. 2013;23(7):885–91.

    Article  PubMed  PubMed Central  Google Scholar 

  57. Hoseini B, Bagheri N, Kiani B, Azizi A, Tabesh H, Tara M. Access to dialysis services: a systematic mapping review based on geographical information systems. Geospat Health. 2018;13(1):577.

    Article  PubMed  Google Scholar 

  58. McGrail MR, Humphreys JS. Spatial access disparities to primary health care in rural and remote Australia. Geospat Health. 2015;10(2):358.

    Article  PubMed  Google Scholar 

  59. Kiani B, Bagheri N, Tara A, Hoseini B, Hashtarkhani S, Tara M. Comparing potential spatial access with self-reported travel times and cost analysis to haemodialysis facilities in North-eastern Iran. Geospat Health. 2018;13(2):240–6.

    Article  Google Scholar 

  60. Kiani B, Bagheri N, Tara A, Hoseini B, Tabesh H, Tara M. Revealed access to haemodialysis facilities in northeastern Iran: factors that matter in rural and urban areas. Geospat Health. 2017;12(2):584.

    PubMed  Google Scholar 

  61. Mohaghegh P, Yavari P, Akbari ME, Abadi A, Ahmadi F, Shormeij Z. Relationships between family levels of socioeconomic status and distribution of breast Cancer risk factors. Iran J Cancer Prev. 2015;8(1):53–9.

    PubMed  PubMed Central  Google Scholar 

  62. Vicens GR, Zafra MS, Moreno-Crespi J, Ferrer BC, Marcos-Gragera R. Incidence variation of prostate and cervical cancer according to socioeconomic level in the Girona health region. BMC Public Health. 2014;14:1079.

    Article  PubMed  PubMed Central  Google Scholar 

  63. Shafique K, Oliphant R, Morrison DS. The impact of socio-economic circumstances on overall and grade-specific prostate cancer incidence: a population-based study. Br J Cancer. 2012;107(3):575–82.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  64. Goshayeshi L, Ghaffarzadegan K, Khooei A, Esmaeilzadeh A, Rahmani Khorram M, Mosannen Mozaffari H, et al. Prevalence and clinicopathological characteristics of mismatch repair-deficient colorectal carcinoma in early onset cases as compared with late-onset cases: a retrospective cross-sectional study in northeastern Iran. BMJ Open. 2018;8(8):e023102.

    Article  PubMed  PubMed Central  Google Scholar 

  65. Goshayeshi L, Khooiee A, Ghaffarzadegan K, Rahmani Khorram M, Bishehsari F, Hoseini B, et al. Screening for Lynch syndrome in cases with colorectal carcinoma from Mashhad. Arch Iran Med. 2017;20(6):332–7.

    PubMed  Google Scholar 

  66. Khorram MR, Goshayeshi L, Maghool F, Bergquist R, Ghaffarzadegan K, Eslami S, et al. Prevalence of mismatch repair-deficient colorectal adenoma/polyp in early-onset, advanced cases: a cross-sectional study based on Iranian hereditary colorectal Cancer registry. J Gastrointest Cancer. 2020.

  67. Ozanne EM, Loberg A, Hughes S, Lawrence C, Drohan B, Semine A, et al. Identification and management of women at high risk for hereditary breast/ovarian cancer syndrome. Breast J. 2009;15(2):155–62.

    Article  PubMed  Google Scholar 

  68. Beebe-Dimmer JL, Kapron AL, Fraser AM, Smith KR, Cooney KA. Risk of prostate Cancer associated with familial and hereditary Cancer syndromes. J Clin Oncol. 2020;38(16):1807–13.

    Article  PubMed  PubMed Central  Google Scholar 

  69. Thompson DJ, Genovese G, Halvardson J, Ulirsch JC, Wright DJ, Terao C, et al. Genetic predisposition to mosaic Y chromosome loss in blood. Nature. 2019;575(7784):652–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  70. Miki Y, Swensen J, Shattuck-Eidens D, Futreal PA, Harshman K, Tavtigian S, et al. A strong candidate for the breast and ovarian cancer susceptibility gene BRCA1. Science. 1994;266(5182):66–71.

    Article  CAS  PubMed  Google Scholar 

  71. Wooster R, Bignell G, Lancaster J, Swift S, Seal S, Mangion J, et al. Identification of the breast cancer susceptibility gene BRCA2. Nature. 1995;378(6559):789–92.

    Article  CAS  PubMed  Google Scholar 

  72. 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(12):2807–16.

    Article  CAS  PubMed  Google Scholar 

  73. Claus EB, Schildkraut JM, Thompson WD, Risch NJ. The genetic attributable risk of breast and ovarian cancer. Cancer. 1996;77(11):2318–24.

    Article  CAS  PubMed  Google Scholar 

  74. Goovaerts P, Jacquez GM. Accounting for regional background and population size in the detection of spatial clusters and outliers using geostatistical filtering and spatial neutral models: the case of lung cancer in Long Island, New York. Int J Health Geogr. 2004;3(1):14.

    Article  PubMed  PubMed Central  Google Scholar 

  75. Zhang C, Luo L, Xu W, Ledwith V. Use of local Moran's I and GIS to identify pollution hotspots of Pb in urban soils of Galway, Ireland. Sci Total Environ. 2008;398(1–3):212–21.

    Article  CAS  PubMed  Google Scholar 

  76. Li K, Lin G-Z, Li Y, Dong H, Xu H, Song S-F, et al. Spatio-temporal analysis of the incidence of colorectal cancer in Guangzhou, 2010–2014. Chin J Cancer. 2017;36(1):60.

    Article  PubMed  PubMed Central  Google Scholar 

  77. Goovaerts P, Jacquez GM. Detection of temporal changes in the spatial distribution of cancer rates using local Moran’s I and geostatistically simulated spatial neutral models. J Geogr Syst. 2005;7(1):137–59.

    Article  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

We would like to thank Kerman University of Medical Sciences for funding this study.

Funding

This study received funding from Kerman University of Medical Sciences (number = 97000230). The funder provided the cost of geocoding and cleaning the data.

Author information

Authors and Affiliations

Authors

Contributions

M.M and L.A provided the data and funding. B.K and N.F conducted spatial analyses. H.R.M and A.B and F.K geocoded the data. A.D drafted the introduction part of the manuscript. V.Z calculated the age-specific incidence rates and cleaned the data. R.B and N.B and S.E revised the text. B.H drafted the discussion part of the manuscript. B.K drafted the method and result sections of the manuscript, designed the study and reviewed the literature and was the research leader. The author(s) read and approved the final manuscript.

Authors’ information

B.K is an assistant professor at Mashhad University of Medical Sciences. He has been working on health geographical information systems and different aspects of geo-spatial data sciences.

N.B is a Senior Research Fellow and spatial epidemiologist at the Australian National University (ANU) with a particular interest in geo-spatial analysis and modelling.

R.B is the editor-in-chief of the Geospatial Health Journal and has much experience in spatial epidemiology.

B.H is a faculty member at Neyshabur University of Medical Sciences and junior researcher at Mashhad University of Medical Sciences. He has been working on the health informatics projects and also has much experience in the health outcome registries.

Corresponding author

Correspondence to Behzad Kiani.

Ethics declarations

Ethics approval and consent to participate

This study has been approved by Ethical Committee of Kerman University of Medical Sciences (number = IR.KMU.REC.1397.290). For retrospective study, the need for consent to participate was deemed unnecessary according to Ethical Committee of Kerman University of Medical Scicnes.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Montazeri, M., Hoseini, B., Firouraghi, N. et al. Spatio-temporal mapping of breast and prostate cancers in South Iran from 2014 to 2017. BMC Cancer 20, 1170 (2020). https://doi.org/10.1186/s12885-020-07674-8

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12885-020-07674-8

Keywords