Skip to main content

ORIGINAL RESEARCH article

Front. Oncol., 24 January 2023
Sec. Breast Cancer

More recent insights into the breast cancer burden across BRICS-Plus: Health consequences in key nations with emerging economies using the global burden of disease study 2019

Sumaira Mubarik&#x;&#x;Sumaira Mubarik1†‡Lisha Luo&#x;Lisha Luo2†Mujahid IqbalMujahid Iqbal3NawsherwanNawsherwan4Jianjun BaiJianjun Bai1Chuanhua Yu*&#x;Chuanhua Yu1*‡
  • 1Department of Epidemiology and Biostatistics, School of Public Health, Wuhan University, Wuhan, China
  • 2Center for Evidence-Based and Translational Medicine, Zhongnan Hospital of Wuhan University, Wuhan, Hubei, China
  • 3Department of Psychology, School of Philosophy, Wuhan University, Wuhan, China
  • 4Xiamen Cardiovascular Hospital of Xiamen University, School of Medicine, Xiamen University, Xiamen, China

Background: Brazil, Russia, India, China, South Africa, and 30 other Asian nations make up the BRICS-Plus, a group of developing countries that account for about half of the world’s population and contribute significantly to the global illness burden. This study aimed to analyzed the epidemiological burden of female breast cancer (BC) across the BRICS-Plus from 1990 to 2019 and studied the associations with age, period, birth cohort and countries’ sociodemographic index (SDI).

Methods: The BC mortality and incidence estimates came from the 2019 Global Burden of Disease Study. We estimated cohort and period effects in BC outcomes between 1990 and 2019 using age-period-cohort (APC) modeling. The maximum likelihood (ML) of the APC-model Poisson with log (Y) based on the natural-spline function was used to estimate the rate ratio (RR). We used annualized rate of change (AROC) to quantify change over the previous 30 years in BC across BRICS-Plus and compare it to the global.

Results: In 2019, there were about 1.98 million female BC cases (age-standardized rate of 45.86 [95% UI: 41.91, 49.76]) and 0.69 million deaths (age-standardized rate of 15.88 [95% UI: 14.66, 17.07]) around the globe. Among them, 45.4% of incident cases and 51.3% of deaths were attributed to the BRICS-Plus. China (41.1% cases and 26.5% deaths) and India (16.1% cases and 23.1% deaths) had the largest proportion of incident cases and deaths among the BRICS-Plus nations in 2019. Pakistan came in third with 5.6% cases and 8.8% deaths. Over the past three decades, from 1990 to 2019, the BRICS-Plus region’s greatest AROC was seen in Lesotho (2.61%; 95%UI: 1.99-2.99). The birth cohort impacts on BC vary significantly among the BRICS-Plus nations. Overall, the risk of case-fatality rate tended to decline in all BRICS-Plus nations, notably in South Asian Association for Regional Cooperation (SAARC) and China-ASEAN Free Trade Area (China-ASEAN FTA) countries, and the drop in risk in the most recent cohort was lowest in China and the Maldives. Additionally, there was a substantial negative link between SDI and case fatality rate (r1990= -0.91, p<0.001; r2019= -0.89, p<0.001) in the BRICS-Plus in both 1990 and 2019, with the Eurasian Economic Union (EEU) nations having the highest case fatality rate.

Conclusions: The BC burden varies remarkably between different BRICS-Plus regions. Although the BRICS’ efforts to regulate BC succeeded, the overall improvements lagged behind those in high-income Asia-Pacific nations. Every BRICS-Plus country should strengthen specific public health approaches and policies directed at different priority groups, according to BRIC-Plus and other high-burden nations.

Introduction

Breast cancer (BC) has now overtaken lung cancer as the most commonly diagnosed malignancy across the globe, with the majority of its health burden occurring in women (1). The GLOBOCAN 2020 estimates indicated that about 2.3 million new incidence cases and 685000 deaths were caused by breast cancer globally in 2020, accounting for 11.7% and 6.9% of total new cancer cases and deaths, respectively (2). The global burden of disease study 2019 (GBD 2019) estimated that the age-standardized incidence and mortality rates for females with breast cancer were 45.86 and 15.88 per 100,000 worldwide in 2019, respectively (3). Without new interventions, the BC burden will increase to more than 3 million new cases and more than 1 million deaths annually by 2040 (1).

Several internationally sustainable efforts have been made to improve the current situation of high burden of breast cancer, such as The Global Breast Cancer Initiative (GBCI), a strategic collaboration established by the World Health Organization in 2021, to control breast cancer burden by global advocacy (4). The GBCI aimed to the shared goal of reducing breast cancer by 2.5% per year, saving 2.5 million lives over 20 years by the three key strategies (health promotion; early detection and timely diagnosis; comprehensive breast cancer management) (4). Therefore, a better understanding of epidemiological characteristics of breast cancer burden across different regions/countries is of great value for providing a scientific foundation and experiences for policymakers to develop targeted strategies.

The previous studies have revealed that marked geographic discrepancies were observed in the BC burden, with the highest burden occurring in countries that underwent economic transition (1, 5). BRICS (Brazil, Russia, India, China, and South Africa) is an informal and financial group comprised of five countries with fast-economic transition, and BRICS-Plus is a new economic group to bring together the regional integration blocks where BRICS economies play a leading role, which was made up by 35 developing countries “in transition” (6). Therefore, it is urgent to investigate the patterns of breast cancer burden in BRICS Plus countries to explore the underlying causes of these disparities using advanced statistical models, thereby helping policymakers tackle this public health challenge and better promote health and economic integration. Additionally, the Socio-Demographic Index (SDI) has been considered an effective measure for predicting and comparing regional and national burdens of disease, evaluating the effectiveness of health policies in countries with socioeconomic transition (7). A previously published review showed there were 33.6% of global new cases and 36.9% of global deaths of breast cancer occurred in BRICS countries, and a marked rising trend in incidence was also observed in these countries, but the corresponding epidemiological variations of its burden and its correlation to age group, birth cohort, and socioeconomic status are not well studied (8).

The current study aimed to compare the long-term trends in incidence, mortality and case fatality rates of breast cancer in BRICS -Plus countries and at the global level, investigating the independent effects of age, period and cohort using data from the Global Burden of Disease Study 2019 (GBD 2019) within the age-period-cohort (APC) framework, and exploring the potential causes of the associations between SDI and these health outcomes. The findings of this study give a combined comparison of breast cancer burden across BRICS Plus countries, which may provide references for the rational allocation of healthcare resources and help policymakers to specify prevention and control strategies according to local conditions.

Material and methods

All anonymized data are accessible online at http://ghdx.healthdata.org/gbd-results-tool, the Institute for Health Metrics and Evaluation. The deidentified, compiled data was used in the GBD investigation. As a result, the University of Washington Institutional Review Board examined and approved a waiver of informed consent.

Data source, estimation of study variables, and study population

The most comprehensive and systematic epidemiological study to date is called the Global Burden of Diseases (GBD). It offers a comprehensive analysis of 359 illnesses and injuries, 84 risk factors, and 282 recorded causes of mortality across 195 nations, 21 regions, and seven super-regions. The GBD uses novel estimation techniques and supporting sources specified in earlier research to give the 2019 estimates and revisions for 1990-2019 data (912). The Global Health Data Exchange (GHDx) online data collecting tool (http://ghdx.healthdata.org/gbd-results-tool) has been used to gather data on breast cancer (BC). The “Availability of data and materials” section at the article’s end defines all the data sources’ specifics. From 1990 to 2019, gathered information on BC from 35 BRICS -Plus nations’ yearly death and incidence rates and age-standardized rates (ASRs) per 100k person-years. The sociodemographic index (SDI), which is based on national per capita income, the average number of years of education for those over 15 years old, and the overall fertility rate, was developed to divide countries into five quintiles (high SDI, high middle SDI, middle SDI, low middle SDI, low SDI). From 0 to 1, or from least to most developed, is the range of this index.

The GBD study states that estimates of BC incidence are based on integrated cancer registry databases or individual cancer registries. The percentage of data gathered from the systematic literature study was subjected to four separate DisMod-MR2.1 inputs (13). As the primary method for estimating the prevalence and incidence of 354 injuries and diseases in 195 nations and territories, DisMod-MR2.1 is based on the Bayesian meta-regression tool.

The BRICS-Plus countries were divided into five regions, including Mercosur (core members as well as acceding members), South African Customs Union (SACU), SAARC (SAFTA members), China-ASEAN FTA and Eurasian Economic Union (EEU), and the specific subdivision list of 35 countries was as follows: Mercosur (Brazil, Argentina, Paraguay, Uruguay, Bolivia, Venezuela); SACU (Botswana, Lesotho, Namibia, South Africa, Swaziland); SAARC (Afghanistan, Bangladesh, Bhutan, India, the Maldives, Nepal, Pakistan, Sri Lanka); China-ASEAN FTA (China, Indonesia, Malaysia, Philippines, Singapore, Thailand, Brunei, Vietnam, Laos, Myanmar, Cambodia); EEU (Russia, Kazakhstan, Belarus, Armenia, Kyrgyzstan).

This study’s female population was diagnosed with BC using the ICD-10 code C50 for the disease (14). The population estimates from 1990–2019 from the online Global Health Data Exchange (GHDx) were used to calculate the matching population estimates for the BRICS-Plus Countries by age group (20-24 to 80-84) and year (1990 to 2019) (15). For instance, the original data for Brazil came primarily from the Ministry of Health’s Mortality Information System (16); New Economic School’s Russian Centre for Demographic Research provides mortality data by age, sex, location, and cause of death (17); The mortality database for India is managed by the Registrar General of India and the Indian Sample Registration System (18); The Chinese Center for Disease Control and Prevention, Disease Surveillance Points, Maternal and Child Surveillance System, and Cause of Death Reporting System are in charge of maintaining the health statistics database for the Chinese population (19). The Department of Home Affairs South Africa Vital Registration system is responsible for keeping track of the country’s overall mortality statistics by age, sex and region and cause of death. Other BRICS-Plus countries used their current health database registration systems to estimate the BC figures. Due to the lack of information on various diseases, injuries, and risk factors from various countries, estimates from GBD rely on complex statistical modeling to handle sparse and usually inconsistent data (20). CODEm (Death Ensemble Model) was run using mortality projections and vital registration system data for BC (21). Based on information and factors, such as education, SDI, lagging distribution income, smoking habit and alcohol consumption, the CODEm calculates mortality rates. The single-cause estimates have been modified to suit all-cause mortality separately using the CodCorrect technique (22, 23). The series of publications written by GBD include specific computation methodologies. In the collection of articles released by GBD, specific calculating techniques are mentioned (11, 2426).

The GBD approach has the advantage of using consistent methods to critically evaluate the currently available data on each condition, making this data comparable and systematic; estimating results from countries with incomplete data; and reporting on the burden of disease with standardized metrics.

Statistical analysis

We investigated the geographic assessment of BC changes by adopting a graphical display and stratifying by time and age group. We presented annualized rates of change (AROC) for changes over time as the proportional difference between the rate’s natural logarithms in 1990 and 2019 divided by 29 (i.e., 100*[ln(2019 Rate/1990 Rate)/(29)]). AROC (%) is a crude trend measure over a 29-year period. An increasing trend/slope or acceleration of mortality/incidence over 29 years in BC is indicated by a positive AROC. In contrast, a decreasing trend/slope or deceleration is shown by a negative AROC. Additionally, based on one thousand bootstrap draws from the posterior distribution, uncertainty intervals (UIs) were used to quantify uncertainty for each outcome (21, 27). We calculated Point estimates from the mean, and UIs were calculated using the 25th and 975th ordered values of the posterior distribution of the 1000 drawings. When the 95% UI of the percentage change did not contain zero, changes over time were deemed statistically significant. For 35 BRICS-Plus nations and globally, we estimated the percent change in age-standardized BC mortality and incidence rates between 1990 and 2019. We gave the 2019 estimates with 95% uncertainty intervals (UI) for those 35 countries. By dividing the age-standardized death rate (ASDR) by the age-standardized incidence rate (ASIR), then multiplying the result by 100, the case-fatality percent (CFP) was determined (28). Spearman correlation coefficient(r) was used to examine the relationship between SDI and BC outcomes.

The relationship between several BC-related indicators, including incidence, mortality, and case fatality rate, and each of the three primary sources of spatial and temporal variability—age, period (year), and cohort (year of birth)—was examined for each BRICS-Plus region/country as well as globally within the framework of the age-period-cohort model. Period effects show population-wide exposure at a specific moment, and the period was defined as the survey year. The term “cohort effects” refers to variations in hazards among birth cohorts. Age-specific rates of BC from 20 to 84 years with subsequent 5-year age intervals, calendar time including a consecutive period from 1990 to 2019, and subsequent cohort (period-age) from 1905 to 1999 were taken into consideration for statistical analysis in the age-period-cohort study. Using the R package (Epi, version 2.44) created by Carstensen et al.,2021, we carried out the APC analysis (29). Maximum Likelihood (ML) of the APC-model Poisson with log (Y) based on the natural-spline function was used to estimate the rate ratio. The reference cohort and reference period were respectively chosen based on the median date of birth and the median date of diagnosis among cases. The estimated models’ deviation table was used to assess the goodness of fit. To assess the study hypothesis, two-sided statistical tests were taken into consideration. The threshold for statistical significance of the findings was set at 1% (level of significance).

Results

Burden and trends in female BC by BRICS-Plus regions and globally

Figure 1 and Supplementary Figures S1, S2 show the burden and trends in female breast cancer (BC) mortality and incidence across the BRICS-Plus countries and globally from 1990 to 2019.

FIGURE 1
www.frontiersin.org

Figure 1 (A) Relative proportion of BRICS-Plus to Global burden of breast cancer incidence and death in year 2019, (B) Country specific breast cancer burden within BRICS-Plus; percent cases (incidence and death) indicate proportion of country specific cases to total cases in BRICS-Plus.

In 2019, there were about 1.98 million female BC cases (age-standardized rate of 45.86 [95% UI: 41.91, 49.76]) and 0.69 million deaths (age-standardized rate of 15.88 [95% UI: 14.66, 17.07]) around the globe. Among these, 45.4% of incident cases and 51.3% of deaths contributed to overall BRICS -Plus. Across the BRICS -Plus countries, there were the highest proportion of incident cases and deaths from China (41.1%cases and 26.5%deaths) and India (16.1%cases and 23.1%deaths), followed by Pakistan (5.6% cases and 8.8%deaths) in 2019. Overall the trends in case fatality percent (CFP) declined across the BRICS -Plus from 1990 to 2019, but these trends were higher than global CFP. The highest CFP was observed in Lesotho during the entire duration from 1990 to 2019 than in other BRICS-Plus countries (Figures 1, S1, S2).

The annualized rate of change in BC over 30 years across BRICS-Plus and globally

Compared to the other BRICS-Plus regions, the SACU’s BC age-standardized mortality rate (ASMR) grew by 7.38% (%AROC Range: 0.29-2.61) from 1990 to 2019. Lesotho had the highest AROC (%) among the BRICS-Plus countries from 1990 to 2019 (2.61% (95%UI: 1.99,2.99). Thirteen of the 35 BRICS-Plus nations, divided into the SACU, SAARC, China-ASEAN FTA, EEU, and Mercosur regions, had a decline in ASMR over the same period. Only 4 of these nations’ ASMR changes were less than the worldwide AROC. In addition, all of the BRICS-Plus countries saw a rise in the AROC for age-standardized incidence rate (ASIR), except for Myanmar and Kyrgyzstan, where ASIR decreased by 1.98% (95%UI: -1.53, -2.29) and 1.19% (95%UI: -1.48, -0.94) over the past three decades from 1990 to 2019. From 1990 to 2019, all of the BRICS-Plus countries had great changes relative to the global AROC in ASIR (Table 1 and Figure 2).

TABLE 1
www.frontiersin.org

Table 1 Age-standardised breast cancer mortality rate and incidence rate (per 100 000 person-years) and Annualized rate of change (AROC%), 1990 to 2019 across BRICS-Plus regions.

FIGURE 2
www.frontiersin.org

Figure 2 Comparison of annualized rate of change (AROC) BRICS-Plus to Global burden of breast cancer. IR, Incidence rate; MR, Mortality rate.

BC age-specific mortality rates across BRICS-Plus and globally

We organized the population and BC (mortality and incidence) data into 19 consecutive cohorts, including those born from 1905 to 1909 (median, 1907) through 1995 to 1999 (median, 1997) and consecutive 5-year periods from 1990 to 1994 (median, 1992) to 2015 to 2019 (median, 2017). Figure 3 displays the BC mortality rate trends by age group over time in the BRICS-Plus and globally. The findings show that the BC mortality rate has risen over time in most of the BRICS-Plus regions, with age group. At the same time, recent years saw the most significant reductions in BC mortality in Myanmar, the Philippines, Singapore, Thailand, Armenia, Belarus, Kazakhstan, and Kyrgyzstan (P<0.01 for all) (Figure 3).

FIGURE 3
www.frontiersin.org

Figure 3 Age-specific mortality rates of breast cancer by period across 35 BRICS-Plus countries and Globally between 1990 to 2019.

From Figure 4, we can observe that globally there was a decreasing trend in BC mortality across the birth cohorts. In contrast, most of the BRICS-Plus countries showed throughout increasing trend of BC mortality rate across birth cohorts. Eswatini, South Africa, Bangladesh, India, Maldives, Sri Lanka, Malaysia, Philippines, Singapore, Thailand, Armenia, Belarus, Kazakhstan, Kyrgyzstan, Russia, Argentina, and Uruguay all displayed an increase in BC mortality rates at first, followed by a decline across all age groups (p<0.001 for all). For instance, those under 50 had consistently lower BC mortality over subsequent birth cohorts than those over 50 (Figure 4).

FIGURE 4
www.frontiersin.org

Figure 4 Cohort-specific mortality rates of breast cancer by age group across 35 BRICS-Plus countries and Globally between 1990 to 2019.

Age period cohort effects on BC across BRICS-Plus and globally

To account for the age, period, and cohort influences of BC outcome in each BRICS-Plus region and country, we computed the APC trends in BC mortality, incidence, and case fatality rates. Figure 5 illustrates these effects while keeping all other variables fixed at their sample mean values.

FIGURE 5
www.frontiersin.org

Figure 5 Age period cohort related female breast cancer trends in (A) death rate (B) incidence rate and (C) case fatality rate (CFR) from 1990-2019 with ages 20 to 84. Rate ratio was estimated using ML of APC-model Poisson with log(Y) based on natural-spline function, for each BRICS-Plus region separately. ML, maximum likelihood; APC, age-period-cohort; Reference cohort for age-effects was chosen as the median date of birth among cases; and Median date of diagnosis among cases was selected as reference period; SACU, South African Customs Union; SAARC, South Asian Association for Regional Cooperation; China-ASEAN FTA, China-ASEAN Free Trade Area; EEU, Eurasian Economic Union; Mercosur, core members as well as acceding members.

For each BRICS-Plus region, age-related effects can be seen to have a significant impact on BC mortality. The SACU, SAARC, and Mercosur regions show the most striking increases in BC mortality risk, with Botswana, Eswatini, South Africa, Pakistan, Uruguay, and Argentina having the highest age-related risk countries.

We found birth cohort impacts vary significantly among the BRICS-Plus nations. In the majority of BRICS-Plus countries, the study revealed an increase in mortality risk trend when birth cohorts increased from the cases’ median birthdate. Except for the EEU and a few other countries from other regions, the risk of death was rising in all BRICS-Plus regions. Lesotho, Pakistan, Burnie, and Venezuela stood out as high cohort-risk nations (Figure 5A). The risk of BC mortality rate tends to rise with age and cohort, and the risk value was more significant than the global effect.

The age period cohort effects on the incidence rate of BC are shown in Figure 5B. Except for the EEU, where a downward trend was seen after the birth cohort of 1945, all BRICS-Plus regions saw an increase in overall age and cohort-associated incidence risks. Most of the nations in the SACU, SAARC, China-ASEAN FTA, and Mercosur regions had higher risk values than the globe. In particular, Lesotho, Sri Lanka, Vietnam, and Venezuela exhibited high incidence cohort-related risks. For cohorts older than 50 and those born between 1940 and after, the effect of separation across different BRICS-Plus regions was more pronounced. In SAARC and Mercosur, the period impact persisted over the entire period. In contrast, SACU, China-ASEAN FTA, and EEU showed moderate changes during the whole study period from 1990 to 2019 (Figure 5).

To fully understand the BC burden, we also did an APC analysis on case-fatality rates (CFR) by BRICS-Plus regions and nations. The Maldives and China experienced the smallest decrease in the risk of case-fatality rate in the most recent cohort, as illustrated in Figure 5C. In the BRICS-Plus region as a whole, case fatality rates tend to drop, particularly in the SAARC and China-ASEAN FTA countries. Between different BRICS-Plus regions, the effect of separation was more evident for cohorts beginning in 1939 and later for those who were greater than 40 years older. In the majority of BRICS-Plus nations and throughout the world, the period impact persisted for the whole time (Figure 5C).

Relationship between sociodemographic index and BC incidence, mortality and case-fatality across BRICS-Plus

Figure 6 illustrates the bivariate relationship between country-wise age-standardized mortality rates, incidence and case fatality percent (CFP) and countries’ SDI for the years 1990 and 2019. We can see that while ASMR and SDI showed mixed patterns throughout BRICS-Plus, for example, the nations that improved their SDI in 2019 also decreased ASMR, ASIR established a positive and statistically significant link with SDI in both 1990 and 2019. Additionally, there was a substantial negative link between SDI and CFP (r1990= -0.91, P<0.001; r2019= -0.89, P<0.001) in the BRICS-Plus region in both 1990 and 2019, with the EEU nations having the highest case fatality rates.

FIGURE 6
www.frontiersin.org

Figure 6 Country wise correlation between age-standardised BC rate and country’s sociodemographic index (SDI) across BRICS-Plus; (A) mortality rate in 1990 (B) incidence rate in 1990 (C) Case fatality percent (CFP) in 1990; (D) mortality rate in 2019 (E) incidence rate in 2019 (F) Case fatality percent (CFP) in 2019; Case-fatality percent was calculated by dividing age-standardised death rate by age-standardised incidence rate and multiplied by 100; LSE, Least Square Error fit; LAE, Least Absolute Error fit; SDI ranges from 0 (less developed) to 1 (most developed).

Additionally, different patterns of epidemiological variations were seen in Figure S3 when age-standardized mortality rate (ASMR), age-standardized incidence rate (ASIR), and case fatality percent (CFP) in each BRICS-Plus country from 1990 to 2019 were plotted against an index of that country’s sociodemographic status in the same year. Along with SACU, SAARC, China-ASEAN FTA, and Mercosur, the ASIR of BC increased with rising SDI in different years, but the ASMR has shown notable variations. As SDI rises, ASIR gradually grows, with more rapid growth for the nations with the greatest levels of SDI. Most of the BRICS-Plus countries showed a downward trend in ASMR, especially those with high SDIs, which had a progressive decline in ASMR as SDI rose. Lesotho, Pakistan, and Bolivia, in contrast to certain BRICS-Plus nations with poor SDI, have had strong or consistent trends in ASMR over time. Additionally, from 1990 to 2019, across the entire BRICS-Plus region, a significant negative association between case fatality percentage and SDI was seen (Figure S4).

Discussion

Breast carcinogenesis has been widely considered a significant global public health concern that continues to increase due to genetic, environmental, lifestyle and socioeconomic risk factors (30). The APC framework allowed estimations of when and how three time-dependent parameters would affect the variations of disease burden, and its result was helpful for theorizing the etiology of the observed trends (31). The current study has examined the age, period and cohort effects of breast cancer incidence, mortality and case fatality rates among females and evaluated the association of these health outcomes with sociodemographic transition across 35 BRICS-Plus countries based on GBD 2019 data from 1990 to 2019. The results revealed that nearly half of global incidence and deaths of breast cancer contributed from BRICS-Plus countries, and its case fatality percent was also higher than the worldwide level, which was consistent with previous studies indicating a shift in the BC burden from high-income to middle- and low-income countries (30, 32, 33). The high proportion suggested that the breast cancer burden was still a serious public health concern in BRICS-Plus countries, so more attention should be paid to exploring the underlying causes of this situation.

Within the BRICS-Plus, striking disparities in breast cancer burden were also observed among these 35 developing countries, with a higher proportion of incidence and deaths in China, India and Pakistan, mainly due to the large population. China, India and Pakistan are among the top five countries in the world with the largest population, with China home to about 1.45 billion, India to 1.40 billion and Pakistan to 0.23 billion in 2022 (34). The total population of China and India account for 36% of the world population, and Pakistan has the highest population growth between 2000 and 2022 (34). Population awareness, healthcare conditions, lifestyles, and socioeconomic factors also played an essential role in the higher cases in these three countries (5, 31). The case fatality percent (CFP) of breast cancer showed decreasing trends across the BRICS-Plus from 1990 to 2019, which indicated the improvement in the quality of healthcare, better prognosis (earlier detection and advanced therapy) and increasing survival rates of cancer cases in these countries (35). Across the BRICS-Plus countries, the highest CFP of breast cancer occurred in Lesotho, a small landlocked country in Southern Africa. It was highly associated with its backward healthcare systems, shortage of health resources and unequal service delivery (36). As we know, Lesotho is one of the most underdeveloped countries declared by the United Nations, which suffers from higher poverty rates, with 72% of the population living in rural areas far away from healthcare services and having experienced an increasing double disease burden from communicable diseases as well as non-communicable diseases.

Additionally, Lesotho was one of 28 African countries that didn’t provide radiotherapy, so many patients cannot be diagnosed and treated early. Although some policies have been made for its health systems challenges in Lesotho, such as the “5-year National Strategic Development Plan (NSDP)”, more effective and targeted measures are also urgent to curb the severe situations in this country (36). Specifically, the most urgent need for Lesotho is adopting radiotherapy for cancer treatment and providing mammographic screening for females above the age of 50 to improve the early diagnosis and treatment of breast cancer and effectively reduce the high case fatality of breast cancer in this country.

The significantly increasing trends in incidence and mortality of breast cancer were observed in most of the 35 BRICS-Plus countries over 30 years, which indicated the increasing burden of breast cancer from developed countries to developing countries confirmed in the prior studies (1, 30). From our results, the SACU region had the largest increase in breast cancer mortality and incidence over the last three decades, which were also consistent with other studies and mainly attributable to the late-stage prognosis and inadequate access to high-quality healthcare resources (1, 37). Additionally, genetic factors also played a key role in the growing incidence of breast cancer in the SACU region, and it’s confirmed the variants in the two major genes, BRCA1 and BRCA2, and in five other genes in the patients in South Africa (38). The incidence increase in most of the BRICS-Plus countries was also associated with mammographic screening, diet changes, improved health awareness, and reduced physical activity (30). The stage at diagnosis was a key prognostic factor for mortality of breast cancer, and a previous meta-analysis had revealed that about 77% of all cases were stage III and IV at diagnosis across 17 sub-Saharan African countries, while only 15% of cases were diagnosed at stages III and IV in some high-income countries (1, 39). Even though the incidence and mortality showed decreasing in a few countries, the AROC in most BRICS-Plus countries were still higher than the global level, so the overall condition of breast cancer in BRICS-Plus still lagged behind those in developed Asia Pacific countries, especially in SACU region. Recognizing this, the BRICS-Plus countries, especially high-risk regions, should learn and translate the experiences from developed countries into strategies suited to their characteristics, including early screening and treatment, improvement of population awareness, and healthier lifestyle, to strengthen the prevention and intervention of breast cancer in high-risk regions, and reduce the public health burden caused by breast cancer.

Age has been considered one of the most important and independent risk factors for developing breast cancer (40), and its incidence and mortality varied across the different age groups in BRICS-Plus countries. After adjusting the influence of period and cohort effect using the APC analysis, the BC mortality and incidence risk showed increasing trends with advancing age, and the older women suffered higher chances, which was consistent with the age trends of breast cancer in the previous studies (31, 41, 42). DNA methylation might be one of the important reasons behind the increased burden with age; additionally, biological factors and hormone levels were also well-proven factors for breast cancer (42, 43). The cohort effects revealed various exposures to behavioral, environmental and socioeconomic factors among different birth cohorts, and our results indicated that the incidence and mortality risk increased with the birth cohorts in most of the BRICS-Plus countries, which indicated that the later birth cohorts had higher risks of breast cancer than earlier cohorts in these countries. On the contrary, the risk of case-fatality rate showed decreasing trends revealing earlier birth cohorts suffered from a higher risk of the case-fatality rate. The previous study also revealed that the risk for breast cancer increased slightly in low-middle and low-socioeconomic regions. At the same time, decreasing trends were observed in the high-socioeconomic areas, which was consistent with our results (30). The possible reasons behind the increasing cohort effects in most countries were westernized lifestyles, inequality of diagnosis, prognosis and treatment of breast cancer across BRICS-Plus countries, and it also highlighted that the improvement of diagnostic tools and easy access to treatment helped decrease breast cancer burden (30, 44). The period effects remained consistent for BC case-fatality rates across most of the BRICS-Plus countries during the entire duration. In contrast, incidence and mortality remained constant in SAARC and Mercosur regions. In contrast, they showed moderate fluctuations in SACU, China-ASEAN FTA and EEU regions, which reflected the immediate effects of the risk factors on breast cancer burden adjusted for age and cohort effects.

In our study, an increase of SDI tended to be positively associated with an increase of breast cancer incidence, which was also confirmed in the previous studies (30, 45). The positive association between socioeconomic status and breast cancer incidence might be caused by easier access to early detection and screening measures in countries with higher socioeconomic status; besides, drastic changes in lifestyle, environmental factors and the use of hormone replacement therapy also played a crucial role (45). The lower incidence of breast cancer in countries with lower socioeconomic status might be underestimated because of unreliable disease surveillance and reporting systems and lower screening rates (45). Therefore, less developed regions should strengthen screening efforts and improve disease surveillance systems to reflect local disease burden levels more accurately. The declining trends of the relationship between case fatality burden and SDI in 1990 and 2019 across the BRICS-Plus were also consistent with the results conducted in 21 world regions, which indicated better treatment facilities in high-income regions and limited treatment resources in less developed regions (30). The mixed patterns of the relationship between mortality and SDI revealed the shift in the death burden of breast cancer from developed to less developed countries (30).

The study also had several limitations. Firstly, all data about breast cancer in our study were obtained from GBD 2019. Some data were not directly measured when the reliably original data was lacking or reporting lag. Hence, it didn’t remain very objective or ambiguities (46). Secondly, interpreting results from the age-period-cohort model at the population level doesn’t necessarily hold for individuals, so there might be an ecological fallacy. Therefore, more individual-based epidemiological studies should be conducted to confirm our results in the future. Thirdly, our study did not analyze the temporal trends of breast cancer burden attributable to risk factors, so further analyses are needed to investigate the variations of attributable burden across the BRICS-Plus countries over the past three decades.

In conclusion, breast cancer remains a serious health challenge among women across BRICS-Plus, and nearly half of the breast cancer incidence and death cases in the world contributed from overall BRICS-Plus, of which China, India and Pakistan were the highest. The age-standardized incidence and mortality rates increased in most of the BRICS-Plus countries over the past three decades, and the highest increase occurred in the SACU region. By APC analysis, breast cancer incidence, mortality and case fatality rates increased with advancing age in most of the BRICS-Plus countries, indicating older women might be the vulnerable population for breast cancer in these regions. The increasing cohort effect on incidence mortality risk revealed that the later birth cohorts were high-risk groups for the development of breast cancer. The rising burden of breast cancer in countries with lower socioeconomic status highlighted that more targeted and cost-effective measures should be strengthened for vulnerable populations in underdeveloped countries to screen and treat breast cancer.

Data availability statement

Publicly available datasets were analyzed in this study. This data can be found here: The dataset analyzed during the current study are available in the Institute for Health Metrics and Evaluation (IHME): http://ghdx.healthdata.org/gbd-results-tool.

Author contributions

CY supervised the study. SM and CY conceptualized the analysis. SM did the data analysis. SM and LL wrote the first draft of the paper. MI, N and JB reviewed and provided comments on the first draft. All authors reviewed the final draft, contributed to the article and approved the submitted version.

Funding

This research was funded by the National Natural Science Foundation of China (Grant No. 82173626) and Health Commission of Hubei Province Scientific Research Project (Grant No. WJ2019H304). The funders had no role in the study design, data collection, analysis, decision to publish or preparation of the manuscript.

Acknowledgments

We appreciate the contributions made by GBD study in Data.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fonc.2023.1100300/full#supplementary-material

Abbreviations

BC, breast cancer; GBD, global burden of diseases; CFR, case-fatality rates; CFP, case-fatality-percent; SDI, sociodemographic index; DR, death rates; IR: incidence rates; ASIR, age-standardized incidence rates; ASMR, age-standardized mortality rates; APC, age-period-cohort, BRICS, Brazil, Russia, India, China, South Africa.

References

1. Arnold M, Morgan E, Rumgay H, Mafra A, Singh D, Laversanne M, et al. Current and future burden of breast cancer: Global statistics for 2020 and 2040. Breast. (2022) 66:15–23. doi: 10.1016/j.breast.2022.08.010

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Sung H, Ferlay J, Siegel RL, Laversanne M, Soerjomataram I, Jemal A, et al. Global cancer statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin (2021) 71(3):209–49. doi: 10.3322/caac.21660

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Diseases GBD, Injuries C. Global burden of 369 diseases and injuries in 204 countries and territories, 1990-2019: A systematic analysis for the global burden of disease study 2019. Lancet. (2020) 396(10258):1204–22. doi: 10.1016/S0140-6736(20)30925-9

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Anderson BO, Ilbawi AM, Fidarova E, Weiderpass E, Stevens L, Abdel-Wahab M, et al. The global breast cancer initiative: a strategic collaboration to strengthen health care for non-communicable diseases. Lancet Oncol (2021) 22(5):578–81. doi: 10.1016/S1470-2045(21)00071-1

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Li N, Deng Y, Zhou L, Tian T, Yang S, Wu Y, et al. Global burden of breast cancer and attributable risk factors in 195 countries and territories, from 1990 to 2017: results from the global burden of disease study 2017. J Hematol Oncol (2019) 12(1):140. doi: 10.1186/s13045-019-0700-2

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Lissovolik Y. BRICS-plus: Alternative globalization in the making. Valdai papers (2017) 69:1–2.

Google Scholar

7. Collaborators GS. Measuring the health-related sustainable development goals in 188 countries: a baseline analysis from the global burden of disease study 2015. Lancet. (2016) 388(10053):1813–50. doi: 10.1016/S0140-6736(16)31467-2

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Basu P, Zhang L, Hariprasad R, Carvalho AL, Barchuk A. A pragmatic approach to tackle the rising burden of breast cancer through prevention & early detection in countries in transition. Indian J Med Res (2020) 152(4):343–55. doi: 10.4103/ijmr.IJMR_1868_19

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Roth GA, Abate D, Abate KH, Abay SM, Abbafati C, Abbasi N, et al. Global, regional, and national age-sex-specific mortality for 282 causes of death in 195 countries and territories, 1980–2017: A systematic analysis for the global burden of disease study 2017. Lancet (2018) 392(10159):1736–88. doi: 10.1016/S0140-6736(18)32203-7

PubMed Abstract | CrossRef Full Text | Google Scholar

10. James SL, Abate D, Abate KH, Abay SM, Abbafati C, Abbasi N, et al. Global, regional, and national incidence, prevalence, and years lived with disability for 354 diseases and injuries for 195 countries and territories, 1990–2017: a systematic analysis for the global burden of disease study 2017. Lancet (2018) 392(10159):1789–858. doi: 10.1016/S0140-6736(18)32279-7

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Kyu HH, Abate D, Abate KH, Abay SM, Abbafati C, Abbasi N, et al. Global, regional, and national disability-adjusted life-years (DALYs) for 359 diseases and injuries and healthy life expectancy (HALE) for 195 countries and territories, 1990–2017: a systematic analysis for the global burden of disease study 2017. Lancet (2018) 392(10159):1859–922. doi: 10.1016/S0140-6736(18)32335-3

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Gakidou E, Afshin A, Abajobir AA, Abate KH, Abbafati C, Abbas KM, et al. Global, regional, and national comparative risk assessment of 84 behavioural, environmental and occupational, and metabolic risks or clusters of risks, 1990–2016: a systematic analysis for the global burden of disease study 2016. Lancet (2017) 390(10100):1345–422. doi: 10.1016/S0140-6736(17)32366-8

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Vos T, Abajobir AA, Abate KH, Abbafati C, Abbas KM, Abd-Allah F, et al. Global, regional, and national incidence, prevalence, and years lived with disability for 328 diseases and injuries for 195 countries, 1990–2016: A systematic analysis for the global burden of disease study 2016. Lancet (2017) 390(10100):1211–59. doi: 10.1016/S0140-6736(17)32154-2

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Organization WH. International classification of diseases, 10th revision (ICD-10) . Available at: https://www.who.int/standards/classifications/classification-of-diseases (Accessed 10 June 2021).

Google Scholar

15. Wang H, Abbas KM, Abbasifard M, Abbasi-Kangevari M, Abbastabar H, Abd-Allah F, et al. Global age-sex-specific fertility, mortality, healthy life expectancy (HALE), and population estimates in 204 countries and territories, 1950–2019: A comprehensive demographic analysis for the global burden of disease study 2019. Lancet (2020) 396(10258):1160–203. doi: 10.1016/S0140-6736(20)30977-6

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Ribeiro ALP, Duncan BB, Brant LC, Lotufo PA, Mill JG, Barreto SM. Cardiovascular health in Brazil: Trends and perspectives. Circulation. (2016) 133(4):422–33. doi: 10.1161/CIRCULATIONAHA.114.008727

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Starodubov VI, Marczak LB, Varavikova E, Bikbov B, Ermakov SP, Gall J, et al. The burden of disease in Russia from 1980 to 2016: A systematic analysis for the global burden of disease study 2016. Lancet (2018) 392(10153):1138–46. doi: 10.1016/S0140-6736(18)31485-5

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Jha P, Gajalakshmi V, Gupta PC, Kumar R, Mony P, Dhingra N, et al. Prospective study of one million deaths in India: rationale, design, and validation results. PLos Med (2006) 3(2):e18. doi: 10.1371/journal.pmed.0030018

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Zhao D, Liu J, Wang M, Zhang X, Zhou M. Epidemiology of cardiovascular disease in China: Current features and implications. Nat Rev Cardiol (2019) 16(4):203–12. doi: 10.1038/s41569-018-0119-4

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Murray CJ, Lopez AD. Measuring the global burden of disease. New Engl J Med (2013) 369(5):448–57. doi: 10.1056/NEJMra1201534

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Foreman KJ, Lozano R, Lopez AD, Murray CJ. Modeling causes of death: An integrated approach using CODEm. Population Health metrics (2012) 10(1):1. doi: 10.1186/1478-7954-10-1

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Wang H, Abajobir AA, Abate KH, Abbafati C, Abbas KM, Abd-Allah F, et al. Global, regional, and national under-5 mortality, adult mortality, age-specific mortality, and life expectancy, 1970–2016: A systematic analysis for the global burden of disease study 2016. Lancet (2017) 390(10100):1084–150. doi: 10.1016/S0140-6736(17)31833-0

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Zhou L, Deng Y, Li N, Zheng Y, Tian T, Zhai Z, et al. Global, regional, and national burden of Hodgkin lymphoma from 1990 to 2017: estimates from the 2017 global burden of disease study. J Hematol Oncol (2019) 12(1):107. doi: 10.1186/s13045-019-0799-1

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Murray CJ, Acharya AK. Understanding DALYs. J Health economics (1997) 16(6):703–30. doi: 10.1016/S0167-6296(97)00004-0

CrossRef Full Text | Google Scholar

25. Murray CJ. Quantifying the burden of disease: The technical basis for disability-adjusted life years. Bull World Health Organization (1994) 72(3):429.

Google Scholar

26. Salomon JA, Haagsma JA, Davis A, de Noordhout CM, Polinder S, Havelaar AH, et al. Disability weights for the global burden of disease 2013 study. Lancet Global Health (2015) 3(11):e712–e23. doi: 10.1016/S2214-109X(15)00069-8

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Naghavi M, Abajobir AA, Abbafati C, Abbas KM, Abd-Allah F, Abera SF, et al. Global, regional, and national age-sex specific mortality for 264 causes of death, 1980–2016: A systematic analysis for the global burden of disease study 2016. Lancet (2017) 390(10100):1151–210. doi: 10.1016/S0140-6736(17)32152-9

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Heer E, Harper A, Escandor N, Sung H, McCormack V, Fidler-Benaoudia MM. Global burden and trends in premenopausal and postmenopausal breast cancer: a population-based study. Lancet Global Health (2020) 8(8):e1027–e37. doi: 10.1016/S2214-109X(20)30215-1

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Plummer M, Carstensen B. Lexis: An R class for epidemiological studies with long-term follow-up. Journal of Statistical Software. (2011) 38:1–2. doi: 10.18637/jss.v038.i05

CrossRef Full Text | Google Scholar

30. Mubarik S, Yu Y, Wang F, Malik SS, Liu X, Fawad M, et al. Epidemiological and sociodemographic transitions of female breast cancer incidence, death, case fatality and DALYs in 21 world regions and globally, from 1990 to 2017: An age-Period-Cohort analysis. J Adv Res (2022) 37:185–96. doi: 10.1016/j.jare.2021.07.012

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Wang Z, Bao J, Yu C, Wang J, Li C. Secular trends of breast cancer in China, south Korea, Japan and the united states: Application of the age-Period-Cohort analysis. Int J Environ Res Public Health (2015) 12(12):15409–18. doi: 10.3390/ijerph121214993

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Ginsburg O, Rositch AF, Conteh L, Mutebi M, Paskett ED, Subramanian S. Breast cancer disparities among women in low- and middle-income countries. Curr Breast Cancer Rep (2018) 10(3):179–86. doi: 10.1007/s12609-018-0286-7

CrossRef Full Text | Google Scholar

33. Yip CH, Buccimazza I, Hartman M, Deo SV, Cheung PS. Improving outcomes in breast cancer for low and middle income countries. World J Surg (2015) 39(3):686–92. doi: 10.1007/s00268-014-2859-6

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Stats IW. The world population and the top ten countries with the highest population (2021). Available at: https://internetworldstats.com/stats8.htm.

Google Scholar

35. Beiki O, Hall P, Ekbom A, Moradi T. Breast cancer incidence and case fatality among 4.7 million women in relation to social and ethnic background: a population-based cohort study. Breast Cancer Res (2012) 14(1):R5. doi: 10.1186/bcr3086

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Ramathebane MM, Sooro MA, Kabuya RM, Sayed AR. Knowledge and attitudes relating to cervical and breast cancer among women in Maseru, Lesotho. African Journal of Primary Health Care & Family Medicine (2022) 14(1):8. doi: 10.4102/phcfm.v14i1.3459

CrossRef Full Text | Google Scholar

37. Verhoeven D, Allemani C, Kaufman C, Mansel R, Riesling S, Anderson B. Breast cancer: Global quality care optimizing care delivery with existing financial and personnel resources. ESMO Open (2020) 4(Suppl 2):e000861. doi: 10.1136/esmoopen-2020-000861

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Abbad A, Baba H, Dehbi H, Elmessaoudi-Idrissi M, Elyazghi Z, Abidi O, et al. Genetics of breast cancer in African populations: A literature review. Glob Health Epidemiol Genom (2018) 3:e8. doi: 10.1017/gheg.2018.8

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Jedy-Agba E, McCormack V, Adebamowo C, dos-Santos-Silva I. Stage at diagnosis of breast cancer in sub-Saharan Africa: a systematic review and meta-analysis. Lancet Global Health (2016) 4(12):e923–e35. doi: 10.1016/S2214-109X(16)30259-5

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Luan HH, Luo LS, Lu ZY. Historical trends in incidence of breast cancer in shanghai, Hong Kong and Los Angeles, 1973-2012: A joinpoint and age-Period-Cohort analysis. Int J Public Health (2021) 66:603810. doi: 10.3389/ijph.2021.603810

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Ataeinia B, Saeedi Moghaddam S, Shabani M, Gohari K, Sheidaei A, Rezaei N, et al. National and subnational incidence, mortality, and years of life lost due to breast cancer in Iran: Trends and age-Period-Cohort analysis since 1990. Front Oncol (2021) 11:561376. doi: 10.3389/fonc.2021.561376

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Mubarik S, Malik SS, Wang Z, Li C, Fawad M, Yu C. Recent insights into breast cancer incidence trends among four Asian countries using age-period-cohort model. Cancer Manag Res (2019) 11:8145–55. doi: 10.2147/CMAR.S208323

PubMed Abstract | CrossRef Full Text | Google Scholar

43. Akoko LO, Rutashobya AK, Lutainulwa EW, Mwanga AH, Kivuyo SL. The effect of reproductive, hormonal, nutritional and lifestyle on breast cancer risk among black Tanzanian women: A case-control study. PLos One (2022) 17(2):e0263374. doi: 10.1371/journal.pone.0263374

PubMed Abstract | CrossRef Full Text | Google Scholar

44. Mubarik S, Liu X, Malik SS, Wang L, Yu Y, Yu C. Evaluation of lifestyle risk factor differences in global patterns of breast cancer mortality and DALYs during 1990-2017 using hierarchical age-period-cohort analysis. Environ Sci pollut Res Int (2021) 28(36):49864–76. doi: 10.1007/s11356-021-14165-1

PubMed Abstract | CrossRef Full Text | Google Scholar

45. Allahqoli L, Mazidimoradi A, Momenimovahed Z, Rahmani A, Hakimi S, Tiznobaik A, et al. The global incidence, mortality, and burden of breast cancer in 2019: Correlation with smoking, drinking, and drug use. Front Oncol (2022) 12:921015. doi: 10.3389/fonc.2022.921015

PubMed Abstract | CrossRef Full Text | Google Scholar

46. Zou Z, Liu G, Hay SI, Basu S, Belgaum UI, Dhali A, et al. Time trends in tuberculosis mortality across the BRICS: an age-period-cohort analysis for the GBD 2019. EClinicalMedicine. (2022) 53:101646. doi: 10.1016/j.eclinm.2022.101646

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: breast cancer, mortality, case fatality, APC, SDI, BRICS-Plus

Citation: Mubarik S, Luo L, Iqbal M, Nawsherwan, Bai J and Yu C (2023) More recent insights into the breast cancer burden across BRICS-Plus: Health consequences in key nations with emerging economies using the global burden of disease study 2019. Front. Oncol. 13:1100300. doi: 10.3389/fonc.2023.1100300

Received: 16 November 2022; Accepted: 04 January 2023;
Published: 24 January 2023.

Edited by:

Nosheen Masood, Fatima Jinnah Women University, Pakistan

Reviewed by:

Ijaz Haq, The University of Haripur, Pakistan
Abbas Khan, University of Home Economic Lahore, Pakistan

Copyright © 2023 Mubarik, Luo, Iqbal, Nawsherwan, Bai and Yu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Chuanhua Yu, yuchua@whu.edu.cn

These authors have contributed equally to this work and share first authorship

ORCID: Sumaira Mubarik, orcid.org/0000-0001-6041-1061, Chuanhua Yu, orcid.org/0000-0002-5467-2481

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.