- 1Department of Colorectal Surgery, Heping Hospital Affiliated to Changzhi Medical College, Changzhi, Shanxi, China
- 2Department of Pediatric Health Care, Zhangzi County Maternal and Child Health Family Planning Service Center, Changzhi, Shanxi, China
Background: Colorectal cancer (CRC) persists as a major global public health priority, with insufficient whole grain consumption recognized as a clinically significant and modifiable dietary risk factor. However, comprehensive analyses of its spatiotemporal burden and socioeconomic disparities are limited. This study evaluates the global, regional, and national CRC burden attributable to low whole grains intake from 1990 to 2021 and projects trends to 2046.
Methods: This study employed the Global Burden of Disease (GBD) 2021 dataset to quantify the association between low whole grain intake and CRC burden. Analyses were conducted across multiple dimensions: temporal (1990–2021), demographic (stratified by sex and 5-year age intervals), and geographic (categorized by national and regional Socio-demographic Index (SDI) quintiles).
Results: From 1990 to 2021, global CRC deaths and disability-adjusted life years (DALYs) attributable to low whole grains intake increased by 82.94% (101,813 to 186,257) and 70.30% (2,540,867 to 4,327,219), respectively. However, age-standardized mortality rates (ASMR) and age-standardized DALYs rates (ASDR) declined globally (AAPC: −0.73 and −0.74). High-SDI regions showed the steepest reductions (ASMR AAPC: −1.17; ASDR AAPC: −1.15), while low-middle-SDI regions experienced rising trends (ASMR AAPC: 0.42; ASDR AAPC: 0.35). Gender disparities persisted: males had higher absolute burdens (2021 deaths: 104,344 vs. 81,912 for females; 2021DALYs: 2527992 vs. 1,799,227), but female ASMR and ASDR declined faster (AAPC: −1.09, −1.12, respectively, vs. −0.47 and −0.47, respectively, for males). Geographically, East Asia, Western Europe, and High-income North America had the highest absolute burdens, whereas Uruguay (ASMR: 5.18/100,000) and Hungary (ASDR: 113.76/100,000) led in age-standardized rates. Cabo Verde exhibited the sharpest increases (ASMR AAPC: 3.16; ASDR AAPC: 2.80). Frontier analysis identified Uruguay and Hungary as high-SDI countries with the largest gaps from efficiency targets. Projections to 2046 suggest continued ASMR/ASDR declines but persistent SDI-driven disparities.
Conclusion: Despite declining age-standardized rates, absolute CRC burden attributable to low whole grains intake increased due to population growth and aging. Socioeconomic disparities highlight the need for targeted interventions in transitioning regions adopting Western diets. Promoting whole grains consumption and prioritizing high-burden populations could mitigate future CRC burden.
1 Introduction
Colorectal cancer (CRC) is one of the most common malignancies worldwide, with incidence and mortality rates steadily increasing in recent years (1–3) According to GLOBOCAN 2022 data, CRC ranks third among newly diagnosed cancers globally and second in terms of cancer-related deaths, posing significant challenges to global public health (4). The epidemiological characteristics of CRC vary markedly across regions, age groups, genders, and racial populations, partly attributable to differences in exposure to risk factors, demographic shifts, and genetic susceptibility (5–7). With economic growth and Westernization of lifestyles, the incidence of CRC is rapidly increasing in many low- and middle-income countries, particularly in Eastern Europe, Asia, and South America (8–10).
To date, dietary factors have been recognized as important determinants of CRC risk, with growing evidence that low whole-grain intake is associated with higher CRC incidence (11). Whole grains are rich in dietary fiber, vitamins, and minerals, which can improve gut health and reduce inflammation, potentially decreasing the risk of CRC (3, 12). Despite these mechanistic insights, global data show that low whole-grain consumption remains prevalent and is contributing to CRC burden across regions (13).
Although existing studies indicate that low whole grains intake is significantly correlated with elevated CRC incidence, epidemiological research on this dietary pattern on a global scale remains limited. Particularly in low- and middle-income countries, the widespread adoption of Western dietary patterns may exacerbate the burden of CRC as low whole grains intake becomes more common (14–16). Given regional and population differences in dietary patterns, a systematic evaluation of the CRC burden attributable to low whole-grain intake is of substantial public health importance, informing dietary guidelines and targeted prevention efforts.
The Global Burden of Disease (GBD) study covers 371 diseases and 88 related risk factors, providing data support for exploring and elucidating the epidemiological characteristics of CRC (17–19). We obtained detailed figures, percentages, and ratios of CRC deaths and disability-adjusted life years (DALYs) attributable to low whole grains diets from the GBD 2021 database. Our analysis focuses on examining temporal trends in CRC burden and elucidating how gender, age, geographic location, and the Sociodemographic Index (SDI) influence deaths and DALY attributable to low whole grains intake. Additionally, we conducted projections to anticipate future trends, providing valuable insights to strengthen CRC prevention strategies.
2 Methods
2.1 Data sources and GBD framework overview
The GBD 2021 study estimated data for 371 diseases and 88 related risk factors at the global, regional, and national levels (19). In this study, we utilized the GBD 2021 database to identify the most significant dietary risk factor for CRC: “Low whole grains” Subsequently, we selected “Colorectal cancer” from the “Tumors” list. Data on CRC deaths, DALY, age-standardized death rates (ASMR), and age-standardized DALY rates (ASDR) attributable to low whole grains intake were extracted from the GBD Results Tool1. These indicators were analyzed across 21 GBD regions and 204 countries from 1990 to 2021, including information on gender, age, and the SDI (19).
All estimates derived from the GBD study, including death counts, DALYs, and age-standardized rates, are presented with their corresponding 95% uncertainty intervals (UIs). Based on the GBD framework, 95% UI for all estimates were calculated by averaging the data from 1,000 draws, with the lower and upper bounds of the 95% UIs established by the 25th and 975th ranked values among all 1,000 draws (20, 21).
The GBD hierarchy organizes locations into three levels: super regions, regions, and countries. They divided countries into 21 regions on the basis of two criteria: epidemiological homogeneity, and geographical contiguity. For some statistical analyses, they grouped regions into 7 super-regions: High-income, Central Europe Eastern Europe and Central Asia, Latin America and Caribbean, Southeast Asia East Asia and Oceania, South Asia, North Africa and Middle East, and Sub-Saharan Africa (22).
The SDI is a composite indicator reflecting a region or country’s per capita income, fertility rate, and educational attainment. Based on SDI values (ranging from 0 to 1), regions and countries were classified into five groups: high SDI (>0.81), high-middle SDI (0.70–0.81), middle SDI (0.61–0.69), low-middle SDI (0.46–0.60), and low SDI (<0.46) (23–25).
2.2 Definition of disease and risk factors
The CRC was identified using the following international classification of diseases, tenth edition (ICD10) codes: C18C19.0, C20, C21-C21.8, Z12.1-Z12.13, Z85.03-Z85.048, Z86.010 (19).
Whole grains intake refers to the average daily consumption, in grams, of whole grains and their products, including breakfast cereals, bread, rice, cookies, muffins, corn cakes, pancakes, pasta, and other related foods. Diets with daily whole grains intake below 140–160 grams are considered to be low in whole grains (26, 27).
2.3 Statistical analyses
We conducted a descriptive analysis to illustrate the overall burden of CRC attributable to low whole grains intake. The analysis included the number of deaths, DALY, ASMR, and ASDR. This study characterized the burden of CRC associated with low whole grains intake across different SDI levels, geographic regions, countries, genders, and age groups by analyzing death counts and DALY. Trends over time were assessed using ASMR and ASDR.
The Joinpoint software (version 5.3.0.0; National Cancer Institute) was used to calculate the Average Annual Percentage Change (AAPC) to evaluate trends in disease burden. A regression model was applied to the natural logarithm of age-standardized indicators, using the formula ln(y) = α + β x + ε, where y represents the age-standardized metric and x is the calendar year. The AAPC was calculated as 100 × (exp(β) − 1). The 95% confidence intervals (CIs) for the AAPC were derived from the standard error of the estimated regression coefficient β (in multi-segment models, the SE is obtained by the appropriate combining rule) and were computed by the software under an asymptotic normal approximation. (28). An AAPC with a 95% CI above zero indicates an increasing trend; below zero indicates a decreasing trend; and if the interval includes zero, the trend is considered stable. It is important to note that the Joinpoint regression analysis was performed on the age-standardized rates (ASMR and ASDR) point estimates obtained from GBD. The 95% CIs for the AAPC reflect the uncertainty in the trend estimation itself, and are distinct from the 95% UIs provided by GBD for each annual rate estimate, which reflect the uncertainty in the data and modeling of the rates.
Frontier analysis evaluates the performance and efficiency of decision-making units (such as countries or regions) in converting inputs into outputs (29). We used Data Envelopment Analysis (DEA) and Stochastic Frontier Analysis (SFA) to construct efficiency frontiers representing best practices, comparing each decision unit’s performance to these frontiers to determine efficiency levels. The performance of countries and regions with varying SDI levels was assessed, with the lower bounds of ASMR and ASDR representing the minimum achievable values given each SDI level.
The age–period–cohort (APC) model was used to evaluate the independent effects of age, period, and cohort on mortality rates and DALYs rates. The analysis was conducted in five-year intervals, with the indicators including the relative risk (RR) for period and cohort effects. Since the CRC mortality rate and DALY rate attributed to low whole grains intake were zero among populations under 25 years old, this group was excluded from the model analysis in this study.
The Bayesian Age-Period-Cohort (BAPC) model is a Bayesian statistical approach that analyzes and forecasts demographic data by simulating the individual effects of age, period, and birth cohort (30). In the BAPC model, the age effect reflects changes in the risk of events (such as disease onset) as individuals age. Period effects account for external factors such as environmental or policy changes that influence all individuals at specific time points. Cohort effects refer to characteristics inherent to individuals born in a particular period, due to their unique experiences and exposures. We used the BAPC R package to fit the model and project global age-standardized rates (ASRs) by gender up to 2046 (31).
Spearman’s rank correlation test was used to examine the relationship between CRC burden attributable to low whole grains intake (ASMR and ASDR) in 2021 and SDI. Statistical analyses and data visualization were performed using R (version 4.4.1) and Joinpoint (version 5.3.0.0). A p-value < 0.05 was considered statistically significant. Our data were obtained from the GBD 2021 Results website (see Footnote 1), which provides publicly accessible, free data. Since ethical approval was obtained for the original study, our research did not require additional ethical approval.
3 Results
3.1 Global CRC burden attributable to low whole grains intake
Between 1990 and 2021, global CRC deaths linked to low whole grains intake exhibited a substantial increase, rising from 101,813 (95% UI: 42,588 to 151,170) to 186,257 (95% UI: 76,127 to 284,803), representing an increase of 82.94%. Meanwhile, DALYs associated with this dietary risk factor climbed from 2,540,867 (95% UI: 1,050,794 to 3,754,416) to 4,327,219 (95% UI: 1,754,865 to 6,578,232), a rise of 70.30%. During the same period, both ASMR and ASDR showed declines, with AAPCs of −0.73 (95% CI: −0.83 to −0.63) and −0.74 (95% CI: −0.85 to −0.64), respectively. Over the same period, the CRC ASMR and ASDR linked to this dietary risk factor decreased in both males and females globally, with more pronounced declines observed in females; the AAPCs were −1.09 (95% CI: −1.19 to −0.99) and −1.12 (95% CI: −1.21 to −1.03), respectively.
At the SDI regional level, ASMR and ASDR were consistently highest in high and high-middle SDI regions in both 1990 and 2021. From 1990 to 2021, high, high-middle, middle, and low SDI regions showed decreasing trends in CRC ASMR and ASDR attributable to low whole grains intake. The most significant decline was in high SDI region, with AAPCs of −1.17 (95% CI: −1.31 to −1.04) for ASMR and −1.15 (95% CI: −1.28 to −1.01) for ASDR. The middle SDI region experienced the slowest declines, with an ASMR AAPC of −0.08 (95% CI: −0.24 to 0.09) and an ASDR AAPC of −0.19 (95% CI: −0.31 to −0.08). Conversely, the low-middle SDI region showed increasing trends, with an ASMR AAPC of 0.42 (95% CI: 0.33 to 0.51) and an ASDR AAPC of 0.35 (95% CI: 0.29 to 0.40).
From 1990 to 2021, CRC deaths and DALYs attributable to low whole grains intake showed significant increases across 21 GBD regions. The top three regions with the highest CRC deaths and DALYs due to low whole grains intake were East Asia, Western Europe, and High-income North America, while Oceania had the lowest values. Central Sub-Saharan Africa, Southern Sub-Saharan Africa, and Andean Latin America also had notable lower values. The most pronounced increases in ASMR were observed in Southern Sub-Saharan Africa, Central Latin America, and Southeast Asia, with AAPCs ranging from 0.69 to 0.92. Similarly, the top three regions for ASDR increases were Central Latin America, Southern Sub-Saharan Africa, and Tropical Latin America, with AAPCs between 0.73 and 1.11. In contrast, Australasia experienced the largest declines in ASMR and ASDR, with AAPCs of −1.64 (95% CI: −1.83 to −1.44) and −1.76 (95% CI: −1.93 to −1.58), respectively (Table 1).
Table 1. Global burden of CRC attributable to low whole grains intake in 1990 and 2021, with Average Annual Percentage Change (AAPC).
In 2021, China, the United States, and Japan ranked the top three for CRC deaths attributable to low whole grains intake, while China, the US, and India ranked the top three for CRC DALYs. In 2021, Uruguay had the highest ASMR at 5.18 (95% UI: 2.09 to 7.73) per 100,000, and Hungary had the highest ASDR at 113.76 (95% UI: 47.61 to 173.37) per 100,000. However, Cabo Verde experienced the greatest increases from 1990 to 2021, with AAPCs of 3.16 (95% CI: 2.37 to 3.96) for ASMR and 2.80 (95% CI: 2.09 to 3.51) for ASDR (Figure 1).
Figure 1. ASMR (A) and ASDR (B) for CRC attributable to low whole grains intake in 204 countries in 2021. The AAPC of the ASMR (C) and ASDR (D) for CRC attributable to low whole grains intake from 1990 to 2021; ASMR, age-standardized mortality rate; ASDR, age-standardized DALYs rate; AAPC, average annual percentage change; CRC, colorectal cancer.
3.2 Joinpoint regression analysis
Figure 2 illustrates the Joinpoint regression trends in ASMR and ASDR for CRC due to low whole grains intake among global males, females, and five SDI regions from 1990 to 2021. Globally, the ASMR for all patients, male patients, and female patients shows a steady decline, with a smaller decrease in male patients (AAPC: −0.47, 95% CI: −0.59 to −0.34, p < 0.001) and a more significant decrease in female patients (AAPC: −1.09, 95% CI: −1.19 to −0.99, p < 0.001). Regarding SDI regions, different SDI areas exhibit varying trends. In high SDI regions, the ASMR shows a significant declining trend (AAPC: −1.17, 95% CI: −1.31 to −1.04, p < 0.001). In middle SDI regions, the overall trend of ASMR remains stable (AAPC: −0.08, 95% CI: −0.24 to 0.09, p = 0.34), while in low-middle SDI regions, an increasing trend is observed (AAPC: 0.42, 95% CI: 0.33 to 0.51, p < 0.001) (Table 1 and Figures 2A,C).
Figure 2. Temporal trends of CRC ASMR and ASDR attributed to low whole grains intake by joinpoint regression, 1990–2021. * means p-values < 0.05 and significant results. Data was provided for global, male, and female ASMRs (A), global, male, and female ASDRs (B), global and different SDI region ASMRs (C) and global and different SDI region ASDRs (D). CRC, colorectal cancer; ASMR, age-standardized mortality rates; ASDR, age-standardized DALYs rates.
Similar to the trend in ASMR, the ASDR for global patients, as well as male and female patients, shows a significant declining trend. Likewise, regarding SDI regions, the trend in ASDR changes remains similar to that of ASMR, with the most significant decline observed in high SDI regions (AAPC: −1.15, 95% CI: −1.28 to −1.01, p < 0.001), while low-middle SDI regions exhibit a notable increasing trend (AAPC: 0.35, 95% CI: 0.29 to 0.40, p < 0.001) (Table 1 and Figures 2B,D).
3.3 Age and sex CRC burden attributable to low whole grains intake
From 1990 to 2021, the ASMR for CRC due to low whole grains intake showed a significant decline in both males and females, with males decreasing from 3.19 (95% UI: 1.35 to 4.78) per 100,000 to 2.75 (95% UI: 1.11 to 4.2) per 100,000, and females decreasing from 2.47 (95% UI: 1.04 to 3.69) per 100,000 to 1.76 (95% UI: 0.74 to 2.66) per 100,000. However, the number of deaths exhibited an increasing trend year by year. A similar pattern was observed in the ASDR. During this period, the ASDR for male CRC decreased from 72.49 (95% UI: 30.36 to 107.77) per 100,000 to 62.39 (95% UI: 25.20 to 95.64) per 100,000, while the ASDR for females decreased from 55.85 (95% UI: 22.91 to 82.83) per 100,000 to 39.36 (95% UI: 16.10 to 59.45) per 100,000. However, the DALYs exhibited an increasing trend year by year (Table 1 and Figure 3).
Figure 3. Comparison of full-age cases ASMR and ASDR of CRC attributable to low whole grains intake among male and female from 1990 to 2021. CRC, colorectal cancer.
In 2021, the number of deaths from CRC due to low whole grains intake among global males and females showed a trend of initially increasing and then decreasing with age. In all age groups under 84, the number of deaths from CRC due to low whole grains intake in males consistently exceeded that in females, highlighting the higher exposure risk for males under 84. In the age group over 84, females deaths surpassed those of male. Males had the highest number of deaths in the 70–74 age group, while females had the highest number in the 80–84 age group, with males mortality rates consistently higher than those of females. Similarly, DALYs also exhibited a trend of initially increasing and then decreasing with age, highlighting the significant impact of CRC due to low grains diets on the lifespan of middle-aged and elderly populations. In all age groups under 84, males consistently had higher DALYs, further emphasizing the heavy burden of CRC on males. However, there were no significant changes in those under 45. The DALYs for both males and females were most severe in the 65–69 age group.
From a gender perspective, female mortality rates gradually increased starting at age 25, with an acceleration after 40, becoming more pronounced after 60, and peaking beyond age 95. In contrast, male mortality rates began to rise around age 30 and surged significantly after 45, indicating an earlier onset of increased risk compared to females.
Similarly, for females, the rate of DALYs was close to zero at age 25, then began to rise slowly, accelerated after the age of 40, and peaked at the age of 95. The trend for males is similar, but the rate of DALYs began to increase earlier and more rapidly, remaining at a higher level in most age groups, indicating a faster accumulation of health loss due to CRC in males. Compared to males, females had lower mortality and DALYs rates across all age groups, indicating a higher burden of CRC in males. However, with increasing age, the gender gap gradually narrows, especially after the age of 90, where the mortality rate and DALYs for females increase at a higher rate than for males. This may be due to the longer life expectancy of females and specific health risks in later life. The results are illustrated in Figure 4.
Figure 4. Comparison of the burden of CRC attributable to low whole grains intake in different age groups and genders in 2021. CRC, colorectal cancer.
3.4 CRC burden attributable to low whole grains intake association with the SDI
Overall, there is a nonlinear “S” shaped association between the ASMR for CRC due to low whole grains intake and the SDI. The ASMR gradually decreases when SDI is less than 0.38, rapidly increases between 0.38 and 0.74, peaks at an SDI of 0.74, and then rapidly decreases when SDI is greater than 0.74. Among different regions, the highest ASMR associated with low whole grains intake was observed in Central Europe, followed by Southern Latin America and Eastern Europe, while the lowest ASMR associated with low whole grains intake was found in South Asia and Oceania (Figure 5A). A similar association exists between the ASDR and SDI (Figure 5B). In 2021, among 204 countries and regions worldwide, the relationship between the ASMR and ASDR for CRC attributable to low whole grains intake and SDI initially decreased, then slowly increased with rising SDI. It began to increase rapidly when SDI exceeded 0.67, peaked at an SDI of 0.79, and then started to decline (Figures 6A,B).
Figure 5. The correlation between ASMR for CRC attributable to low whole grains intake in 21 GBD regions from 1990 to 2021 and SDI (A). The correlation between ASDR for CRC attributable to low whole grains intake in 21 GBD regions from 1990 to 2021 and SDI (B). ASMR, age-standardized mortality rate; ASDR, age-standardized DALYs rate; SDI, sociodemographic index; CRC, colorectal cancer.
Figure 6. The correlation between ASMR for CRC attributable to low whole grains intake in 204 countries in 2021 and SDI (A). The correlation between ASDR for CRC attributable to low whole grains intake in 204 countries in 2021 and SDI (B). ASMR, age-standardized mortality rate; ASDR, age-standardized DALYs rate; SDI, sociodemographic index; CRC, colorectal cancer.
From 1990 to 2021, a significant negative correlation was observed between the AAPC of ASMR and SDI (R = −0.49, p < 0.001) (Figure 7A). Similarly, there was an inverse relationship between the AAPC of ASDR and the SDI of various countries (R = −0.47, p < 0.001) (Figure 7B).
Figure 7. The correlation between the AAPC of ASMR for CRC attributable to low whole grains intake and SDI in 2021 (A). The correlation between the AAPC of ASDR for CRC attributable to low whole grains intake and SDI in 2021 (B). ASMR, age-standardized mortality rate; ASDR, age-standardized DALYs rate; SDI, sociodemographic index; AAPC, average annual percentage change; CRC, colorectal cancer.
3.5 Frontier analysis of CRC burden attributable to low whole grains intake
To better understand the potential improvements in the burden of CRC due to low whole grains intake, we conducted a frontier analysis based on ASMR, ASDR, and the SDI levels of regions and countries from 1990 to 2021. The frontier line shows the lowest achievable ASRs for each country or region at corresponding SDI levels, becoming stable when SDI surpasses 0.27. The gap between individual points and this black line reflects the potential for reducing future disease burden, with larger gaps indicating greater opportunities. In the four diagrams, these points form an inverted triangle, with lower values on the left and higher values on the right. As SDI increases, the variation in point distribution widens, indicating that countries with higher SDI have greater potential to lower CRC ASMR and ASDR linked to low whole grains intake. The trends in changes of ASMR and ASDR from 1990 to 2021 are shown in Figures 8A,C. For CRC ASMR due to low whole grains intake, Bangladesh, Gambia, and Papua New Guinea are the top three countries with the highest SDI and the closest effective distances, demonstrating the best ASMR control in these countries. Interestingly, we found that the 10 countries with the largest effective distances for ASMR have SDI levels above 0.70. These countries are Uruguay, Hungary, Bulgaria, Monaco, Slovakia, Poland, Croatia, Greenland, Barbados, and Serbia (Figure 8B). The frontier for ASDR also shows a similar trend (Figure 8D).
Figure 8. Frontier analysis results of CRC attributable to low whole grains intake ASMR and ASDR. Frontier analysis of ASMR from 1990 to 2021 (A,B). Frontier analysis of ASDR from 1990 to 2021 (C,D). Black lines represent the lower limits of ASRs achievable at different SDI levels, with points representing different countries and regions. The 15 countries and regions with the largest effective differences globally are labeled in black font, the 5 countries and regions with the smallest effective differences among low SDI countries are labeled in blue font, and the 5 countries and regions with the largest effective differences among high SDI countries are labeled in red font. In (A,C), the blue dots represent the ASRs of CRC from 1990 to 2021, with darker shades indicating later years. In (B,D), the dots represent changes in CRC ASR from 1990 to 2021. Blue dots indicate countries and territories where the ASR increased from 1990 to 2021, while red dots indicate countries and territories where the ASR decreased. SDI, Socio-demographic index; ASMR, age-standardized mortality rate; ASDR, age-standardized DALYs rate; CRC, colorectal cancer; ASR, age-standardized rate.
3.6 APC analysis of CRC burden attributable to low whole grains intake
Under the assumption of controlling for period and cohort effects, from 1990 to 2021, the globally CRC deaths rates and DALYs rates due to low whole grains intake increased with age, reaching a peak in the age group over 95 years (Figures 9A,D). After controlling for age and cohort effects, from 1990 to 2021, the period RR values for CRC mortality and DALYs rates attributable to low whole grains intake globally showed a decreasing trend over time, using the 2007–2011 period as the reference group (RR = 1). The highest mortality RR was observed in 1992–1996 at 1.11 (95% CI 1.10–1.12) (Figure 9B); similarly, the DALYs rate RR peaked at 1.10 (95% CI 1.09–1.12) during 1992–1996 (Figure 9E). After controlling for age and period effects, the cohort RR values for CRC mortality and DALYs rates linked to low whole grains intake globally from 1990 to 2021 initially increased, then fluctuated with a declining trend over time. Using the 1955–1959 birth cohort as the reference group (RR = 1), the highest mortality RR was observed in the 1910–1914 cohort at 1.33 (95% CI 1.28–1.39) (Figure 9C); similarly, the DALYs rate RR peaked at 1.32 (95% CI 1.23–1.41) during this period (Figure 9F).
Figure 9. Age-period-cohort model analysis of CRC attributable to low whole grains intake. (A–C) Age period and cohort effect of deaths rate. (D–F) Age period and cohort effect of DALYs rate; CRC, colorectal cancer.
3.7 Projections of ASRs of CRC attributable to low whole grains intake to 2046
Using the BAPC model, we predicted the trends of ASMR and ASDR for CRC due to low whole grains intake from 2022 to 2046 by gender. The predictions indicate a reduction in ASMR to varying degrees in both sexes, with ASDR expected to decrease in both male and female groups, showing similar trends. At the same time, the predictive analysis shows that both ASMR and ASDR decrease at a significantly faster rate in females compared to males (Figure 10).
Figure 10. Global trends in ASRs of CRC attributable to low whole grains intake by gender from 1990 to 2046. (A–C) ASMRs for all sexes, males, and females. (D–F) ASDRs for all sexes, males, and females; ASRs, age-standardized rates; CRC, colorectal cancer; ASMRs, age-standardized mortality rates; ASDR, age-standardized DALYs rate.
4 Discussion
The CRC ranks third among newly diagnosed cancers worldwide, imposing a significant disease burden on the global population (4). Dietary factors play a crucial role in its etiology and progression; for instance, low whole grains intake, high-fat diets, and high consumption of processed meats all promote the development of CRC (32–34). Insufficient whole grains intake has been the leading cause of CRC-related deaths since 1990. However, studies examining the long-term trends of the burden of CRC attributable to low whole grains diets remain limited. This study systematically analyzed the global changes in the burden of CRC due to low whole grains intake from 1990 to 2021, examined its associations with the SDI, gender, age, and geographic regions, and conducted frontier and predictive analyses. The aim of this study is to elucidate the evolving trends in the impact of low whole grains intake on CRC burden, providing evidence to inform targeted nutritional interventions and public health strategies.
The study found that, despite significant increases in the absolute numbers of CRC deaths and DALYs attributable to low whole grains intake globally—rising by 82.94 and 70.30%, respectively—ASMR and ASDR showed decreasing trends. This apparent contradiction may be due to population growth and aging offsetting improvements in age-standardized rates. Notably, significant differences exist among SDI regions and genders: high-SDI regions experienced the most pronounced declines in ASMR and ASDR, while middle-to-low SDI regions showed increasing trends. Additionally, worldwide, the disease burden is higher in males than in females; however, females exhibit faster decreases in ASMR and ASDR, and among the elderly (>84 years), female deaths and DALYs surpass those of males. Geographically, the highest absolute burdens are observed in East Asia, Western Europe, and high-income North America, while countries such as Uruguay and Hungary have the highest age-standardized rates globally.
The findings of this study are consistent with previous epidemiological research linking diet and CRC, indicating a negative association between whole grains intake and CRC risk. Some studies have demonstrated a dose–response relationship, whereby each additional serving of whole grains is associated with approximately a 17.00% reduction in CRC risk (12, 35–37). Moreover, this study further reveals heterogeneity in temporal trends: high-SDI countries may have reduced disease burden through improved dietary patterns, increased screening, and advances in medical technology, while middle- and low-SDI countries may experience increasing burdens due to Westernized diets, increased processed food consumption, and insufficient public health resources. Unlike studies focusing solely on absolute case numbers, this study employed age-standardized rates to eliminate the effects of population structure, thereby providing a more accurate reflection of disease prevention and control effectiveness. Additionally, the faster decline in ASMR and ASDR among females compared to males may be related to increased health awareness, the protective effects of estrogen, and differences in gut microbiota; however, further validation is needed (38–41).
The study suggests that whole grains diets may inhibit CRC development through multiple mechanisms. First, whole grains are rich in insoluble fiber, which increases stool bulk, shortens intestinal transit time, and reduces contact between carcinogens (such as secondary bile acids) and the intestinal mucosa (42, 43). Second, components such as polyphenols, vitamin E, and selenium in whole grains can suppress chronic inflammation (such as NF-κB pathway) and oxidative stress, both of which are key drivers of CRC (44–47). Finally, higher whole grains intake is associated with lower obesity rates and insulin resistance, both of which are independent risk factors for CRC (47–50). In summary, low whole grains intake may promote CRC development by reducing dietary fiber intake, enhancing chronic inflammation and oxidative stress, altering the gut microbiota, and exacerbating metabolic conditions—these factors collectively contribute to carcinogenesis.
This study provides critical evidence to guide public health strategies aimed at mitigating the CRC burden attributable to low whole grains intake. The revealed nonlinear, “S”-shaped relationship between SDI and disease burden offers a predictive framework for policymakers. It indicates that low-middle and middle SDI countries, currently experiencing the steepest rise in burden, stand at a critical juncture where proactive intervention during this nutritional transition phase is essential to avert a future surge in CRC cases. Conversely, high-SDI countries must sustain and intensify their efforts to maintain the observed declining trends. To translate these findings into action, a multi-tiered strategy is proposed. Firstly, at the national policy level, governments should implement fiscal policies—such as subsidizing whole grain production or taxing highly refined foods—to enhance affordability and accessibility. Complementing this, mandatory front-of-pack labeling and clear national standards for “whole-grain” products are crucial to empower informed consumer choices and prevent misleading marketing. National dietary guidelines and public awareness campaigns must also be updated to explicitly promote daily consumption targets (e.g., 140–160 g/day) and highlight the specific role of whole grains in cancer prevention. Secondly, healthcare systems must deliver targeted interventions. Regions with the highest burden, namely Eastern Europe and South America, should integrate dietary risk assessment and counseling into routine primary care and CRC screening programs. Within these regions, men and the elderly, identified as the highest-risk demographics, should be prioritized for screening and educational initiatives, with physicians trained to provide specific, actionable dietary advice. Finally, global coordination is key. International bodies like the WHO should incorporate whole grains intake targets into global non-communicable disease (NCD) frameworks and assist lower-SDI countries in building capacity for robust dietary surveillance. Furthermore, facilitating knowledge transfer from high-SDI countries with successful CRC prevention experiences to those earlier in their development trajectory could accelerate progress. Implementing these evidence-based, targeted policies can reshape food environments, guide individual choices, and ultimately contribute to a sustained reduction in the global CRC burden.
Previous studies have explored the relationship between low whole grain diet and the burden of CRC (27). Building on this prior work, our research employed advanced analytical methods, including the AAPC, Joinpoint regression analysis, Frontier analysis, APC models, and BAPC prediction models to provide a more comprehensive and nuanced understanding of the trends in disease burden. Despite these methodological improvements, our study has certain limitations. For instance, the data relied primarily on publicly available statistics, which may be subject to biases and incompleteness. Additionally, the quality of dietary intake data and cancer registry data at the national level varies, which may introduce estimation bias. Third, the analysis did not fully account for interactions with other dietary factors (such as red meat consumption) and lifestyle behaviors (such as smoking and physical activity). Lastly, ecological analysis cannot establish a direct causal relationship between whole grains intake and CRC. Collectively, these factors may contribute to biases in the results.
Future research should focus on: (1) conducting multicenter cohort studies integrating genomics and metabolomics to elucidate molecular mechanisms of whole grains’ protective effects; (2) developing regional dietary intervention models to evaluate economic costs and health benefits; (3) exploring synergistic effects between whole grains and other dietary components (such as prebiotics) to optimize nutritional recommendations. In summary, this study offers valuable insights into global CRC prevention and control, emphasizing that improving whole grains intake should be a central component of multidimensional public health strategies, while also drawing attention to nutritional inequalities during socioeconomic transitions. It also provides health implications for certain regions and countries, offering a theoretical foundation for implementing targeted measures.
5 Conclusion
This study systematically analyzed the impact of low whole grains intake on the burden of CRC at global, regional, and national levels, revealing complex spatiotemporal dynamics in the association between CRC and low whole grains intake as populations change and dietary habits evolve. The results indicate that, despite an overall downward trend in ASMR and ASDR globally, some regions, particularly low- and middle-income countries, have experienced increases. These findings suggest that low whole grains intake remains a significant modifiable risk factor for CRC and holds important implications for future public health interventions.
These findings highlight the potential preventive value of improving whole grains intake levels. Nutritional intervention strategies should be tailored to local contexts. Policymakers should enhance dietary health promotion and encourage the consumption of whole grains to mitigate the future burden of CRC. It is recommended to promote whole grains food production through tax incentives, require mandatory labeling of the proportion of whole grains on food packaging, and include whole grains intake as part of the national chronic disease surveillance indicators. Additionally, regional disparities suggest that differentiated public health measures are necessary for countries at various stages of development.
Future research should focus on evaluating the effectiveness of intervention strategies, exploring the mechanisms linking whole grains intake to CRC pathogenesis, and utilizing micro-level data to uncover individual differences. Through multidisciplinary collaborations, efforts can be made to optimize prevention strategies, thereby providing a solid scientific foundation for effective control of CRC and the promotion of overall population health.
Data availability statement
Publicly available datasets were analyzed in this study. This data can be found here: https://ghdx.healthdata.org/gbd-2021.
Ethics statement
Ethical approval was not required for the study involving humans in accordance with the local legislation and institutional requirements. Written informed consent to participate in this study was not required from the participants or the participants’ legal guardians/next of kin in accordance with the national legislation and the institutional requirements.
Author contributions
KW: Methodology, Software, Supervision, Visualization, Writing – original draft, Writing – review & editing. XL: Methodology, Resources, Software, Writing – original draft. ZG: Methodology, Software, Writing – original draft. JC: Methodology, Software, Writing – original draft. YL: Methodology, Resources, Writing – original draft. HL: Methodology, Supervision, Writing – original draft. SG: Supervision, Validation, Writing – review & editing.
Funding
The author(s) declare that no financial support was received for the research and/or publication of this article.
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.
Generative AI statement
The authors declare that no Gen AI was used in the creation of this manuscript.
Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.
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.
Footnotes
References
1. Morgan, E, Arnold, M, Gini, A, Lorenzoni, V, Cabasag, CJ, Laversanne, M, et al. Global burden of colorectal cancer in 2020 and 2040: incidence and mortality estimates from GLOBOCAN. Gut. (2023) 72:338–44. doi: 10.1136/gutjnl-2022-327736,
2. Kazemi, E, Zayeri, F, Baghestani, A, Bakhshandeh, M, and Hafizi, M. Trends of colorectal cancer incidence, prevalence and mortality in worldwide from 1990 to 2017. Iran J Public Health. (2023) 52:436–45. doi: 10.18502/ijph.v52i2.11897,
3. Wu, X, Xu, J, Yang, X, Wang, D, and Xu, X. Integrating transcriptomics and metabolomics to explore the novel pathway of fusobacterium nucleatum invading colon cancer cells. Pathogens. (2023) 12:201. doi: 10.3390/pathogens12020201,
4. Bray, F, Laversanne, M, Sung, H, Ferlay, J, Siegel, RL, Soerjomataram, I, et al. Global cancer statistics 2022: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clinicians. (2024) 74:229–63. doi: 10.3322/caac.21834,
5. Wong, MCS, Huang, J, Lok, V, Wang, J, Fung, F, Ding, H, et al. Differences in incidence and mortality trends of colorectal cancer worldwide based on sex, age, and anatomic location. Clin Gastroenterol Hepatol. (2021) 19:955–966.e61. doi: 10.1016/j.cgh.2020.02.026,
6. Douaiher, J, Ravipati, A, Grams, B, Chowdhury, S, Alatise, O, and Are, C. Colorectal cancer—global burden, trends, and geographical variations. J Surg Oncol. (2017) 115:619–30. doi: 10.1002/jso.24578,
7. Dong, W, Kim, U, Rose, J, Hoehn, RS, Kucmanic, M, Eom, K, et al. Geographic variation and risk factor association of early versus late onset colorectal cancer. Cancer. (2023) 15:1006. doi: 10.3390/cancers15041006,
8. Khan, SZ, and Lengyel, CG. Challenges in the management of colorectal cancer in low- and middle-income countries. Cancer Treat Res Commun. (2023) 35:100705. doi: 10.1016/j.ctarc.2023.100705,
9. Abreu Lopez, BA, Pinto-Colmenarez, R, Caliwag, FMC, Ponce-Lujan, L, Fermin, MD, Granillo Cortés, AV, et al. Colorectal cancer screening and management in low- and middle-income countries and high-income countries: a narrative review. Cureus. (2024) 16:e70933. doi: 10.7759/cureus.70933,
10. Nawi, AM. Public health: prevention In: A. Vannelli (ed.) Colorectal cancer. London, UK: IntechOpen (2020)
11. Gaesser, GA. Whole grains, refined grains, and cancer risk: A systematic review of meta-analyses of observational studies. Nutrients. (2020) 12:3756. doi: 10.3390/nu12123756,
12. Aune, D, Chan, DSM, Lau, R, Vieira, R, Greenwood, DC, Kampman, E, et al. Dietary fibre, whole grains, and risk of colorectal cancer: systematic review and dose-response meta-analysis of prospective studies. BMJ. (2011) 343:d6617. doi: 10.1136/bmj.d6617,
13. Abdullah, MMH, Hughes, J, and Grafenauer, S. Whole grain intakes are associated with healthcare cost savings following reductions in risk of colorectal cancer and total cancer mortality in Australia: A cost-of-illness model. Nutrients. (2021) 13:2982. doi: 10.3390/nu13092982,
14. Pasdar, Y, Shadmani, FK, Fateh, HL, Soleimani, D, Hamzeh, B, Ghalandari, M, et al. The burden of colorectal cancer attributable to dietary risk in middle east and north african from 1990 to 2019. Sci Rep. (2023) 13:20244. doi: 10.1038/s41598-023-47647-y,
15. Clemente-Suárez, VJ, Beltrán-Velasco, AI, Redondo-Flórez, L, Martín-Rodríguez, A, and Tornero-Aguilera, JF. Global impacts of western diet and its effects on metabolism and health: A narrative review. Nutrients. (2023) 15:2749. doi: 10.3390/nu15122749,
16. Zhou, E, and Rifkin, S. Colorectal cancer and diet: risk versus prevention, is diet an intervention? Gastroenterol Clin N Am. (2021) 50:101–11. doi: 10.1016/j.gtc.2020.10.012,
17. GBD 2021 Causes of Death Collaborators. Global burden of 288 causes of death and life expectancy decomposition in 204 countries and territories and 811 subnational locations, 1990-2021: a systematic analysis for the global burden of disease study 2021. Lancet. (2024) 403:2100–203. doi: 10.1016/S0140-6736(24)00367-2
18. GBD 2021 Risk Factors Collaborators. Global burden and strength of evidence for 88 risk factors in 204 countries and 811 subnational locations, 1990-2021: a systematic analysis for the global burden of disease study 2021. Lancet. (2024) 403:2162–2132. doi: 10.1016/S0140-6736(24)00933-4
19. Zhang, X, Xu, Z, Shang, L, Yang, Q, Ye, H, Liu, H, et al. Global burden of colorectal cancer attributable to metabolic risks from 1990 to 2021, with predictions to 2046. BMC Cancer. (2025) 25:228. doi: 10.1186/s12885-025-13643-w,
20. Wu, Z, Xia, F, and Lin, R. Global burden of cancer and associated risk factors in 204 countries and territories, 1980–2021: A systematic analysis for the GBD 2021. J Hematol Oncol. (2024) 17:119. doi: 10.1186/s13045-024-01640-8,
21. GBD 2019 Mental Disorders Collaborators. Global, regional, and national burden of 12 mental disorders in 204 countries and territories, 1990–2019: a systematic analysis for the global burden of disease study 2019. Lancet Psychiatry. (2022) 9:137–50. doi: 10.1016/S2215-0366(21)00395-3
22. Murray, CJ, Ezzati, M, Flaxman, AD, Lim, S, Lozano, R, Michaud, C, et al. GBD 2010: design, definitions, and metrics. Lancet. (2012) 380:2063–6. doi: 10.1016/S0140-6736(12)61899-6,
23. Global Burden of Disease 2019 Cancer CollaborationKocarnik, JM, Compton, K, Dean, FE, Fu, W, Gaw, BL, et al. Cancer incidence, mortality, years of life lost, years lived with disability, and disability-adjusted life years for 29 cancer groups from 2010 to 2019: a systematic analysis for the global burden of disease study 2019. JAMA Oncol. (2022) 8:420–44. doi: 10.1001/jamaoncol.2021.6987,
24. Zhang, Q, Feng, J, Xu, Z, Guo, Y, Zhu, B, and Qian, P. Burden of colorectal cancer attributable to high body-mass index in 204 countries and territories, 1990–2021: results from the global burden of disease study 2021. Public Health. (2025) 242:388–98. doi: 10.1016/j.puhe.2025.02.040,
25. Wei, J, Fan, L, He, Z, Zhang, S, Zhang, Y, Zhu, X, et al. The global, regional, and national burden of type 2 diabetes mellitus attributable to low physical activity from 1990 to 2021: A systematic analysis of the global burden of disease study 2021. Int J Behav Nutr Phys Act. (2025) 22:8. doi: 10.1186/s12966-025-01709-8,
26. Deng, Y, Wei, B, Zhai, Z, Zheng, Y, Yao, J, Wang, S, et al. Dietary risk-related colorectal cancer burden: estimates from 1990 to 2019. Front Nutr. (2021) 8:690663. doi: 10.3389/fnut.2021.690663,
27. Ma, Y, Ni, J, Mei, P, Chen, Y, and Guo, X. The burden of colorectal cancer attributable to diet low in whole grains from 1990 to 2021: A global, regional and national analysis. Front Nutr. (2025) 12:1527522. doi: 10.3389/fnut.2025.1527522,
28. Kim, H-J, Fay, MP, Feuer, EJ, and Midthune, DN. Permutation tests for joinpoint regression with applications to cancer rates. Statist Med. (2000) 19:335–51. doi: 10.1002/(SICI)1097-0258(20000215)19:3<335::AID-SIM336>3.0.CO;2-Z,
29. Sultana, S, Hossain, MM, and Haque, MN. Estimating the potato farming efficiency: A comparative study between stochastic frontier analysis and data envelopment analysis. PLoS One. (2023) 18:e0284391. doi: 10.1371/journal.pone.0284391,
30. Jürgens, V, Ess, S, Cerny, T, and Vounatsou, P. A bayesian generalized age–period–cohort power model for cancer projections. Stat Med. (2014) 33:4627–36. doi: 10.1002/sim.6248,
31. Knoll, M, Furkel, J, Debus, J, Abdollahi, A, Karch, A, and Stock, C. An R package for an integrated evaluation of statistical approaches to cancer incidence projection. BMC Med Res Methodol. (2020) 20:257. doi: 10.1186/s12874-020-01133-5,
32. Ahmed, FE. Effect of diet, life style, and other environmental/chemopreventive factors on colorectal cancer development, and assessment of the risks. J Environ Sci Health C. (2004) 22:91–148. doi: 10.1081/LESC-200038263,
33. Baena, R, and Salinas, P. Diet and colorectal cancer. Maturitas. (2015) 80:258–64. doi: 10.1016/j.maturitas.2014.12.017,
34. Vargas, AJ, and Thompson, PA. Diet and nutrient factors in colorectal cancer risk. Nutr Clin Pract. (2012) 27:613–23. doi: 10.1177/0884533612454885,
35. Zhang, B, Zhao, Q, Guo, W, Bao, W, and Wang, X. Association of whole grain intake with all-cause, cardiovascular, and cancer mortality: A systematic review and dose–response meta-analysis from prospective cohort studies. Eur J Clin Nutr. (2018) 72:57–65. doi: 10.1038/ejcn.2017.149,
36. Wei, H, Gao, Z, Liang, R, Li, Z, Hao, H, and Liu, X. Whole-grain consumption and the risk of all-cause, CVD and cancer mortality: A meta-analysis of prospective cohort studies. Br J Nutr. (2016) 116:514–25. doi: 10.1017/S0007114516001975,
37. Zong, G, Gao, A, Hu, FB, and Sun, Q. Whole grain intake and mortality from all causes, cardiovascular disease, and cancer. Circulation. (2016) 133:2370–80. doi: 10.1161/CIRCULATIONAHA.115.021101,
38. Rebersek, M. Gut microbiome and its role in colorectal cancer. BMC Cancer. (2021) 21:1325. doi: 10.1186/s12885-021-09054-2,
39. Das, B, and Nair, GB. Homeostasis and dysbiosis of the gut microbiome in health and disease. J Biosci. (2019) 44:117. doi: 10.1007/s12038-019-9926-y,
40. Zackular, JP, Rogers, MAM, Ruffin, MT, and Schloss, PD. The human gut microbiome as a screening tool for colorectal cancer. Cancer Prev Res Phila. (2014) 7:1112–21. doi: 10.1158/1940-6207.CAPR-14-0129,
41. Yadav, D, Sainatham, C, Filippov, E, Kanagala, SG, Ishaq, SM, and Jayakrishnan, T. Gut microbiome–colorectal cancer relationship. Microorganisms. (2024) 12:484. doi: 10.3390/microorganisms12030484,
42. Tullio, V, Gasperi, V, Catani, MV, and Savini, I. The impact of whole grain intake on gastrointestinal tumors: A focus on colorectal, gastric, and esophageal cancers. Nutrients. (2021) 13:81. doi: 10.3390/nu13010081,
43. Hullings, AG, Sinha, R, Liao, LM, Freedman, ND, Graubard, BI, and Loftfield, E. Whole grain and dietary fiber intake and risk of colorectal cancer in the NIH-AARP diet and health study cohort. Am J Clin Nutr. (2020) 112:603–12. doi: 10.1093/ajcn/nqaa161,
44. Luo, H, Fang, Y-J, Zhang, X, Feng, X-L, Zhang, N-Q, Alinuer, A, et al. Association between dietary zinc and selenium intake, oxidative stress-related gene polymorphism, and colorectal cancer risk in Chinese population – a case-control study. Nutr Cancer. (2021) 73:1621–30. doi: 10.1080/01635581.2020.1804950
45. Khan, J, Gul, P, Rashid, MT, Li, Q, and Liu, K. Composition of whole grain dietary fiber and phenolics and their impact on markers of inflammation. Nutrients. (2024) 16:1047. doi: 10.3390/nu16071047,
46. Munteanu, C, and Schwartz, B. Interactions between dietary antioxidants, dietary fiber and the gut microbiome: their putative role in inflammation and cancer. Int J Mol Sci. (2024) 25:8250. doi: 10.3390/ijms25158250,
47. Slavin, JL. Mechanisms for the impact of whole grain foods on cancer risk. J Am Coll Nutr. (2000) 19:300S–7S. doi: 10.1080/07315724.2000.10718964,
48. Malin, SK, Kullman, EL, Scelsi, AR, Haus, JM, Filion, J, Pagadala, MR, et al. A whole-grain diet reduces peripheral insulin resistance and improves glucose kinetics in obese adults: A randomized-controlled trial. Metabolism. (2018) 82:111–7. doi: 10.1016/j.metabol.2017.12.011,
49. Lutsey, PL, Kori, S, Mayer-Davis, E, Shea, S, Steffen, LM, Szklo, M, et al. Whole grain intake and its cross-sectional association with obesity, insulin resistance, inflammation, diabetes and subclinical CVD: the MESA study. Br J Nutr. (2007) 98:397–405. doi: 10.1017/S0007114507700715
Keywords: colorectal cancer (CRC), global burden of disease, low whole grains intake, death, disability-adjusted life years (DALYs), age-period-cohort (APC) model
Citation: Wang K, Li X, Guo Z, Chen J, Li Y, Liu H and Guo S (2025) Trend and age-period-cohort analysis of colorectal cancer burden attributable to low whole grains intake: a global, regional, and national analysis with predictive modeling (1990–2046). Front. Nutr. 12:1632103. doi: 10.3389/fnut.2025.1632103
Edited by:
Emmanouella Magriplis, Agricultural University of Athens, GreeceCopyright © 2025 Wang, Li, Guo, Chen, Li, Liu and Guo. 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: Kai Wang, d2FuZ2thaTE5OTEwMjE4QDE2My5jb20=; Shuwei Guo, Z3Vvc2h1d2VpMzEyODQwOEAxNjMuY29t
Xiaodan Li2