Incidence trend of breast Cancer in women of eastern Mediterranean region countries from 1998 to 2019: A systematic review and meta-analysis

Background This study was conducted to provide evidence on the current status of breast cancer and its incidence trend in Eastern Mediterranean Region during 1998–2019. Also, this study aimed to investigate the association between the incidence of breast cancer and Human Development Index and some factors related to this index, including total fertility rate, and obesity, using a meta-analysis. Method Data on incidence of breast cancer were collected from various sources, including PubMed, Embase, Web of Science, and WHO, from 1998 to 2019 using systematic review and meta-analysis. Pooled age standardized rate was calculated based on study duration and quality of data using a subgroup analysis and random effect meta-analysis. Results A total of 80 studies (545 data points) were analyzed. Pooled age standardized rate of breast cancer for Eastern Mediterranean Region was 37.1 per 100,000 person-year (95% confidence interval [CI], 34.5, 39.8) during 2011–2019. age standardized rate of breast cancer had an upward trend in Eastern Mediterranean Region from 2005 to 2019. However, the increasing trend was found to be slightly different in various regions based on quality of data. Moreover, pooled age standardized rate had a significant association with Human Development Index [− 89.2 (95% CI, − 119.8, − 58.7)] and obesity [1.2 (95% CI, 0.9, 1.5)]. Conclusion Pooled age standardized rate of breast cancer in Eastern Mediterranean Region was lower than the global average. Also, the age standardized rate value and its incremental trend have been higher in countries with high-quality data than in other countries of this region in recent years. Data quality or physiological factors, such as increase in obesity rates, could be the reasons for this incremental trend.


Background
Breast cancer (BC) is the most prevalent cancer among women worldwide and has also been the fifth leading cause of death among cancers in both sexes globally between 2005 and 2015 [1]. According to previous studies, BC will be one of the most important causes of death in women in reproductive age in developing countries in the future [2]. The average age at presentation of BC was less in Asian women compared to those in Western countries [3][4][5][6][7]. previous studies showed the age standardized rate (ASR) of BC in women increased from 58 to 65.5 per 100,000 person-year worldwide during 2005-2015. A more incremental trend of ASR in BC has been observed in countries with low socioeconomic index than in countries with a high socioeconomic index. In 2015, the lowest ASR of BC belonged to Southeast Asian countries, with 35.8 (95% CI, 27.5-45.4), followed by South Asian countries, with 44.4 per 100,000 personyear (95% CI, 37.1-52.3) [1]. According to GLOBOCAN 2018, the 5-year prevalence of BC was 88.8 per 100,000 women in Eastern Mediterranean Region (EMR) [8].
According to World Health Organization (WHO), EMR comprises of 21 countries, with the population of about 583 million. The countries in this region have diverse economic, social, health indexes, and life expectancy [9].
In some studies, socioeconomic factors have been mentioned as possible factors influencing BC incidence [1,10]. Human Development Index (HDI), a composite index comprised of life expectancy at birth, salary, and education, is one of the indicators used to study the level of welfare in human societies [11]. Considering the lack of proper cancer registration systems and surveillance structure in most of developing countries, particularly in EMR countries [12], and given that EMR includes countries with diverse socioeconomic and health status [9], similar patterns may be found in this region and in other similar countries in the present and future. This study was conducted to provide evidence on the current status of BC and its incidence trend in the EMR during1998-2019. Moreover, this study investigated the association between the incidence of BC and HDI and some factors related to this index, including total fertility rate (TFR) and obesity, using a meta-analysis.

Search strategy
Studies were selected using a systematic and comprehensive search of the literature, review of references, government publications, and recommendations by active researchers in the field. Electronic databases were searched using the following keywords: "breast neoplasms", "breast cancer", "breast tumor", "incidence", "frequency", "distribution", and "epidemiology" by adding the names of all the 21 EMR countries separately (Web table1). Relevant studies were identified by searching the WHO Global Index Medicus Database, Medline, PubMed, Embase, and Web of Science. The grey literature was found on the websites of WHO, IARC, IRCT, Pecos, Google, and Google Scholar. The papers and reports published up to November 2019 were searched.

Inclusion and exclusion criteria
All national and international studies and reports on the incidence of BC in EMR countries (Iran, Afghanistan, Pakistan, Qatar, Kuwait, Egypt, Lebanon, Oman, Jordan, Yemen, Iraq, Bahrain, Libya, Morocco, Saudi Arabia, Sudan, Syria, Tunisia, Djibouti, United Arab Emirates (UAE), and Somalia; based on WHO classification) during 1998-2019 were included in this study.
Review studies and abstract papers of conferences and congresses that did not have full-texts and studies conducted on specific age groups or on a specific histology of BC were excluded. When studies with similar data were found, the one with higher quality was selected. Quality assessment of studies was done using JBI (The Joanna Brigge Institute) checklist [13] and risk of bias assessment checklist [14]. Each checklist has 10 items, with the scores ranging from 0 to 10. The items included representativeness of the target population, recruitment methods, adequacy of sample size, reporting details of study methods, appropriateness of measurements, and analysis and reporting methods. These 2 checklists were filled up by RZ&FT independently to decrease risk of bias. Moreover, the validity and reliability of quality assessment checklists between 2 investigators were checked through coefficient agreement. The average scores were calculated using the checklists, and studies with a score below 5 were excluded from the study, based on the average score of both checklists.

Extracted data
Extracted data included title of the study, DOI (digital object identifier) of the paper, writer's name, year of the study, location of the study, name of the journal, year of publication, and writer's address. Specific information included study period, sampling method, sample size, number of BC cases, data collection source, and study results. Also, crude incidence rate and ASR of BC were standardized based on WHO population and standard error (SE), standard deviation (SD), or confidence interval (CI) of these indexes.

Data manipulation
Crude incidence rate and sample size were used to calculate the number of BC cases when their number was not reported. However, when the number was reported but no information was provided on dispersion index (SD, SE, or CI), Keyfitz formula was used to estimate SE [15]. When neither the dispersion index (such as SE) nor the number of BC cases was reported, SE was obtained using MICE (Multiple Imputation via Chained Equations) [16]. SE was calculated using this method in approximately 13% of the cases.
In cases where a study was conducted during several years and the annual ASR was not reported, but the average of ASR was reported for that period, the mid-point year of the study was considered as a data point. Quality assessment of the studies and data extraction were done by 2 investigators (ZR & TF) independently. Disagreement between the investigators was resolved by discussion and review or by referring to the third author (ZF).

Final data used in the analysis
In addition to the data extracted from the systematic review of the literature, other information, such as TFR, was extracted from HDI of United Nations Development Program during 1998-2018 (UNDP) [11] and from Index Mundi website, from which data on the prevalence of obesity in every EMR country were extracted [17].

Data analysis
All ASR index analyses were run based on SQRT (square root) and according to Poisson distribution due to the positive skewness of the studied variables [18]. Considering the high values of I square index (68.4%, p value < 0.001) and Tau 2 = 2.2(Tau 2 is the estimated standard deviation of underlying effects across studies), and according to the results of the initial meta-regression analyses, the variables of time, location, and quality of cancer registry system were significant. Therefore, subgroup analysis was done to reduce heterogeneity based on 4 separate time periods (1998-2000; 2000-2005; 2006-2010; and 2011-2019) and also based on 3 categories of the quality of cancer registry system. As the major part of data were reported based on the data of the cancer registry system, the extracted data were categorized according to the quality of cancer registry system in each country. of data, Group 1 had high-quality data and included countries in which the coverage of populationbased registry of cancer was over 50% (Qatar, Kuwait, Bahrain, and Oman); group 2 had medium-quality data and included countries in which the coverage of population-based registry of cancer was lower than 50% (Iran, Morocco, Saudi Arabia, Libya, Tunisia, Egypt, and Jordan); and group 3 had low-quality data based on pathology, treatment centers, or survey (Afghanistan, Pakistan, Iraq, Yemen, Djibouti, Somalia, Syria, Sudan, UAE, Lebanon).
In this study, a meta-analysis was done in 3 parts:

Determining the current status of breast Cancer
The data on the ASR of BC during 2011-2019 were used to determine the current status of BC. The pooled ASR of BC was calculated separately for each country in EMR using random-effect meta-analysis method.

Studying the ASR trend of breast Cancer
To find the trend of BC, data were analyzed based on the data of the whole region from 1998 to 2019 and then based on the quality of data in the corresponding period. Data were analyzed using metan command in randomeffects meta-analysis and fractional polynomial regression method for the subgroups based on quality of data.

Studying the association between related factors influencing the ASR of breast Cancer
The association between the ASR of BC and HDI, TFR, and obesity were investigated using unadjusted (univariate) and adjusted (multivariate) meta-regression method.
In addition, meta-regression method was used to investigate the effect of time, location, and quality of cancer registry system on the heterogeneity of the extracted ASR of BC. As about 60% of data points were related to Iran, sensitivity analysis was used to estimate the pooled ASR of BC in EMR. In the sensitivity analysis, the results were reported by including the studies related to Iran and then by excluding them. All the analyses were done using STATA 12.0 software (StataCorp, College Station, TX, USA). All p values were 2-sided, and significance level was set at less than 0.05.

Results
In the initial search, 4235 papers, international, and national reports were identified, among which 80 studies met the criteria to be included in the analysis [1,6,8,12,. The search processes are illustrated in Fig. 1, and the summarized characteristics of the included studies are presented in supplementary file (Table S2). The average agreement between the 2 researchers (ZR&TF) to select the studies, conduct quality assessment, and data extraction was good (Cohen's unweighted κ = 0.87). A total of 545 data points of BC incidence were extracted from 80 studies, among which 55% (298 data points) belonged to Iran. In this study, most of the data (77%) belonged to the national cancer registry system and 40% to the third period of the study (2006)(2007)(2008)(2009)(2010). Also, 72% of the data belonged to those EMR countries that had medium-quality of data (group 2) ( Table 1). Countries with medium quality of data (group 2) Unlike the group with high data quality, a downward trend, followed by an increasing trend was observed in the group with medium data quality; the corresponding figures were 21.2 (95% CI, 16.8, 26.01) to 18.5 (95% CI, 16.1, 21.2) and to 33.6 (95% CI, 30.3, 37.21) per 100,000

The association between the ASR of breast Cancer and developmental risk factors in EMR
The results of the initial analysis showed that location of the study, year of study, and quality of data had significant effects on the heterogeneity of data. The results also revealed that tau 2 [from 2.2 to 1.01] was decreased significantly after insertion of these variables in the metaregression model (Table 4).

Discussion
The current status of the ASR of breast cancer The pooled ASR of BC occurrence was estimated to be about 37 per 100,000 person-year in EMR during 2011-2019. However, after excluding Iran, the results of sensitivity analysis showed that ASR of BC had a significant increase. The global average of BC in women was 65.5 in 2015 according to Global Burden of Disease [41] and it was 41.9 in EMR in 2012 according to GLOBOCAN [5]. The results of the present study were more consistent with GLOBOCAN 2012.
In Lebanon, the ASR of BC was almost 3 times more than that of Syria during 2011-2019. This was confirmed by another study conducted in 2004, in which the ASR of BC was higher than other Arab countries and Iran [6].

Trend of breast Cancer incidence
There was an upward trend for the ASR of BC in EMR from 1998 to 2019. However, based on the quality of data during 1998-2019, the pattern of BC incidence has been somewhat different in EMR countries. The studies performed on Asian and EMR countries in recent years have shown an increasing trend in the occurrence of BC in these countries [48,59,85,[87][88][89]95]. Also, WHO has predicted that the increasing trend of BC will continue in EMR countries in the next decade [84].
The highest value for the ASR of BC in the last period and the highest increase in the incidence of BC from first period to last periods of the study belonged to countries with high-quality data. In fact, only 4 countries (Kuwait, Oman, Jordan, and Tunisia) have had a population-based cancer registry prior to 1998 (11). Cancer registration systems were changed from pathological data collection to a population-based registration system during 1998-2005 in Bahrain, Qatar, Morocco, Egypt, Iran, and Libya [47]. Moreover, in addition to increase in westernized lifestyle of people living in the EMR region, cancer registry has also improved overtime [85,96,97]. These factors can be explained somewhat by the fluctuation in the trend of ASR of BC, with a noticeable increase in the last period. Moreover, the diversity in the increasing patterns of BC in various EMR countries may be due to the difference in the time of upgrading the cancer registry systems in these countries.

Effects of developmental factors on breast Cancer incidence
HDI showed a negative association with the incidence of BC in the adjusted and unadjusted models, which was significant in the adjusted model; the negative association of HDI with BC incidence contradicted the findings of other studies [1,48,85,98]. HDI cannot cause an increase in the incidence of BC directly; however, increased longevity, decreased fertility, increased obesity, or changes in lifestyle can increase the incidence of BC [85,96].The highest HDI in EMR belonged to Arab countries, such as Saudi Arabia, Oman, Kuwait, Qatar, UAE, and Bahrain. The ASR of BC in these countries was lower than that of some other countries in EMR and European countries. The lower ASR of BC reported by Arab countries, compared to the Western countries, may be due to women's higher fertility and breastfeeding. Therefore, the reverse association between HDI and ASR of BC in this study, compared to other studies, may be due to the reproductive pattern and lifestyle variety in highest HDI countries in EMR than the Western countries [99]. Lebanon, Pakistan, and Iraq, which had the highest ASR of BC during the recent years belonged to the middle, low, and middle HDI countries, respectively [11]. Several studies suggested that the higher incidence of BC in Pakistan and Iraq may be less affected by hormonal-and parity-related factors and more affected by genetic and environmental factors [100,101]. The ASR of BC showed a positive association with obesity, as this association was confirmed by other studies [1,48,85,96]. The association between ASR of BC and TFR was positive in the unadjusted and negative in the adjusted models, but it was not significant, indicating that the positive association in the unadjusted model may be due to the confounder effect of other variables, such as quality of data, location, time of the study, HDI, and obesity. Other studies have reported a significant negative association between TFR and ASR of BC [102,103]. The difference between the findings of this study and those of others may be due to the decline in TFR in all EMR countries [17].

Limitations of the study
Considering that no comprehensive study has been conducted on the incidence rate of BC and its trend in EMR, the results of this study can help policymakers of these countries to develop and implement programs to reduce BC incidence and improve the data registry system. However, this study had some limitations that should be taken into consideration when interpreting the results: the small number of studies; lack of availability of adequate data in some counties, especially in recent years; and the mere use of English papers. Nevertheless, in this study, it was not possible to investigate some other risk factors of BC, such as breastfeeding, diet, hormone therapy, and physical activity, which could justify the incremental trend of BC incidence. Thus, when interpreting the results of studies that have used ecological methods to investigate the association between BC incidence and risk factors, their limitations should also be taken into account.

Conclusion
Based on the results of the study, pooled ASR of BC in EMR was lower than its global average. Also, the highest and lowest value of ASR belonged to countries with lowquality data during 2011-2018. Furthermore, there was an increasing trend in the ASR of BC in EMR in recent years, especially in the low-quality data group. The possible reasons for the incremental trend of BC incidence may be data availability and quality, or physiological factors, such as increase in the rate of obesity. Although the quality of data in cancer registry systems has improved, the published data on the incidence of BC in EMR countries have been limited in recent years. Some EMR countries still lack a national registry system or a population-based system. Thus, after development and improvements in the registry systems of EMR countries, an upward trend of BC is highly expected in this region in the future.
According to these results, it is highly recommended to investigate the factors associated with the lower incidence of breast cancer in EMR, especially in higher HDI countries in this region.
Additional file 1: Table S1. Strategy for systematic searches of the published literature. Table S2. Summarized characteristics of included studies. Graph S1. The ASR trend of breast cancer in Eastern Mediterian Region by qulity of data registration throughout 1998-2019. Graph S2. The ASR trend of breast cancer in Eastern Mediterian Region by qulity of data registration throughout 1998-2019. Figure S1. Geographic Distribution of ASR of breast cancer in Eastern Mediterranean Countries from 1998 to 2018.

Availability of data and materials
The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.
Ethics approval and consent to participate The research was approved by ethical committee of Kerman University of Medical Sciences. The Ethic approval Cod is IR.KMU.REC.1396.1617.

Consent for publication
Not applicable.