- State Key Laboratory of Kidney Diseases, Beijing Key Laboratory of Medical Devices and Integrated Traditional Chinese and Western Drug Development for Severe Kidney Diseases, Beijing Key Laboratory of Digital Intelligent TCM for the Prevention and Treatment of Pan-vascular Diseases, Department of Nephrology, First Medical Center of Chinese PLA General Hospital, National Clinical Research Center for Kidney Diseases, Key Disciplines of National Administration of Traditional Chinese Medicine (zyyzdxk-2023310), Beijing, China
Objective: The Global Burden of Disease Study (GBD) 2019 reveals an increasing prevalence of diabetic nephropathy caused by excessive intake of sugar-sweetened beverages in high and low SDI regions from 1990 to 2021.
Methods: This study comprehensively analyzed the burden of DN caused by SSBs in high and low SDI regions from 1990 to 2021 and projected the trends until 2040 using ARIMA models.
Results: The results revealed distinct patterns and trends in age-standardized rates of deaths, DALYs, YLDs, and YLLs between the two regions. In high SDI regions, the age-standardized rates of deaths, DALYs, YLDs, and YLLs exhibited a consistent upward trend from 1990 (0.044, 1.425, 0.537, 0.888 per 100,000) to 2021 (0.096, 2.284, 0.67, 2.154 per 100,000). Projections for the age-standardized rates in high SDI regions from 2021 to 2040 indicate that these upward trends will persist, with deaths, DALYs, YLDs, and YLLs further increasing to 0.127, 3.629, 1.070, and 2.861 per 100,000. Conversely, low SDI regions displayed different trends. Although the age-standardized rates of deaths, DALYs, YLDs, and YLLs also increased from 1990 (0.018, 0.488, 0.071, 0.417 per 100,000) to 2021 (0.020, 0.522, 0.081, 0.441 per 100,000), the rate of increase was relatively slower compared to high SDI regions. Projections for the age-standardized rates in low SDI regions from 2021 to 2040 suggest that these upward trends will continue, with deaths, DALYs, YLDs, and YLLs further increasing to 0.016, 0.498, 0.105, and 0.398 per 100,000.
Conclusion: This study provides a comprehensive analysis of the burden of type 2 diabetes nephropathy caused by excessive intake of sugar-sweetened beverages in high and low SDI regions. The results revealed distinct patterns and trends in age-standardized rates of deaths, DALYs, YLDs, and YLLs between the two regions.
1 Introduction
Globally, the excessive intake of sugar-sweetened beverages (SSBs) has been on the rise, emerging as a significant public health concern. SSBs, which include sodas, energy drinks, and fruit-flavored drinks sweetened with added sugars, have been linked to a host of health problems (1, 2). Among these, their contribution to the development of type 2 diabetes mellitus (T2DM) is well-established. Excessive intake of SSBs leads to weight gain, insulin resistance, and impaired glucose tolerance, all of which are key factors in the pathogenesis of T2DM (3, 4). Diabetic nephropathy (DN), a serious complication of T2DM, is a leading cause of end-stage renal disease (ESRD) worldwide (5). It not only impairs the quality of life of patients but also imposes a heavy economic burden on healthcare systems. The development and progression of DN are influenced by a variety of factors, including genetic predisposition, duration and glycemic control of diabetes, hypertension, and dyslipidemia (6–8). In recent years, there has been growing evidence suggesting that dietary factors, particularly the excessive intake of SSBs, may also play a role in the development and progression of DN (9, 10). The results of the National Health and Nutrition Examination Survey indicate that sugary soft drink consumption is associated with albuminuria (11). Excessive intake of SSBs has been shown to induce rapid spikes in blood glucose and insulin levels from high levels of sugars or high fructose cornsyrup, which in conjunction with the large volumes consumed contribute to a high dietary glycemic load (GL) (12, 13). Fructose has also been shown to promote accumulation of visceral adiposity or ectopic fat deposition (14). Visceral adiposity, rather than BMI per se, appears to be the most significant predictor of the risk factors comprising metabolic syndrome, and a growing body of literature shows that metabolic syndrome is a risk factor for CKD (15, 16). Besides, High-fructose corn syrup-sweetened soft drink consumption increases vascular resistance in the kidneys at rest and during sympathetic activation (17).
The Global Burden of Disease (GBD) study is a comprehensive and systematic effort to quantify the health-related loss due to diseases and risk factors worldwide (18, 19). It provides a unique opportunity to assess the global burden of DN caused by excessive intake of SSBs. By analyzing the GBD database, we can obtain detailed information on the incidence, prevalence, deaths, and disability-adjusted life years (DALYs) attributable to DN due to excessive intake of SSBs in different regions and populations. Diabetic nephropathy is exacerbated by excessive intake of SSBs as a complication of diabetes. The Global Burden of Disease (GBD) Study was initiated in the 1990s for the purpose of assessing health outcomes in a timely, valid, and relevant manner (20–22). This study aims to provide a comprehensive analysis of the global burden of DN caused by excessive intake of SSBs using the GBD database. We will examine the trends in incidence, prevalence, deaths, and DALYs of DN due to excessive intake of SSBs over time and across different regions and age-groups. The results of this study will have important implications for public health policy-making and interventions aimed at reducing the burden of DN caused by excessive intake of SSBs. It will also provide valuable insights for further research on the mechanisms underlying the association between excessive intake of SSBs and DN. The results highlight the significant and increasing public health challenge posed by this condition, particularly in high SDI (Socio-demographic Index) regions.
2 Materials and methods types
2.1 Study data
GBD 2021 represents a multinational collaborative integrated surveillance system, providing comprehensive and systematic estimates of 369 diseases and injuries as well as 87 risk factors from 1990 to 2021 (23–25). We utilized the GBD Results Tool to retrieve data, selecting “Cause of death or injury” under the GBD Estimate directory; “Deaths,” “DALYs,” “YLDs (Years Lived with Disability),” “YLLs (Years of Life Lost)” under the Measure directory; “Number” and “Rate” under the Metric directory; “Excessive intake of sugar-sweetened beverages” under the Cause directory; “All countries and territories” under the Location directory; “1–95 + years” and “ALL-age standardized” under the Age directory; “Both,” “Female,” and “Male” under the Sex directory; and “1990–2021″ under the Year directory. For this study, we acquired estimates of deaths, DALYs, YLDs, and YLLs attributed to DN due to excessive intake of SSBs based on high and low SDI, age, location, and sex from the Global Health Data Exchange query tool (accessed on 14 January 2025)1 (26). Subsequently, the retrieved raw data was exported. It is important to note that the GBD initiative is guided by ethical principles governing data collection and utilization, and the study protocol received approval from the Institutional Review Board at the University of Washington (Seattle, WA, United States).
The SDI serves as a metric gaging per capita income, educational attainment, and total fertility rate in each state, with a higher index signifying improved socio-demographic development (27, 28). In GBD 2021, the SDI was employed as a metric, encompassing aspects such as income per capita distribution, average years of education for individuals aged 15 and above, and fertility rates in females under 25, for each location and year under consideration. Categorization of countries was based on their SDI values in 2019, grouped into five categories: low SDI (<0.455), low-middle SDI (≥0.455 and < 0.608), middle SDI (≥0.608 and < 0.690), high-middle SDI (≥0.690 and < 0.805), and high SDI (≥0.805).
2.2 Statistical analysis
To measure trends related to deaths, Disability-Adjusted Life Years (DALYs), Years Lived with Disability (YLDs), and Years of Life Lost (YLLs) from DN linked to excessive intake of SSBs in regions categorized by Socio-Demographic Index (SDI), we employed the Age-Standardized Rate (ASR) and assessed the annual percentage change (EAPC) (45, 46). Standardizing data is essential for comparing different populations with distinct age structures or for examining the same population across different time periods, given that the age profile evolves (47, 48). The Age-Standardized Incidence Rate (ASIR) and Age-Standardized Death Rate (ASDR) per 100,000 individuals were calculated by aggregating the products of age-specific rates, where dk represents the distribution of the chosen reference population’s age groups, and Bk denotes the age-specific rate. A weighted average of these rates was derived based on the age distribution, with the relevant formulas detailed:
The EAPC and their 95% confidence intervals (Cls) are a comprehensive measure of ASR trends over a specific period, with lower bounds above 0points to an upward trend, whereas upper bounds below 0 points to a downward trend. We use a log-linear regression to calculate (where y = In [ASR], and x = calendar year), as follows:
The analysis of deaths, DALYs, YLDs, and YLLs across high and low SDI regions includes the Age-Standardized Rate (ASR) for the year 2021, the percentage change in cases, and the Estimated Annual Percentage Change (EAPC) from 1990 to 2021. An upward trend in the ASR of deaths, DALYs, YLDs, and YLLs is indicated by positive estimates for both the EAPC and the 95% confidence intervals (CIs), and conversely, a downward trend is indicated by negative values. Furthermore, Pearson correlation analysis was conducted to estimate the correlation coefficients and p-values. All data visualization and statistical analyses were conducted using GraphPad Prism (version 8.02) and R software (version 3.5.2). The ARIMA model is a differential integrated moving average autoregressive model, also known as an integrated moving average autoregressive model, which is one of the time series forecasting analysis methods. In ARIMA (p, d, q), AR is “autoregressive” and p is the number of autoregressive terms. MA is the “moving average,” q is the number of terms in the moving average, and d is the number of differences (order) made to make it a stationary sequence. The model fitting was evaluated by coefficient of, Mean Absolute Error (MAE), and Root Mean Squared Error (RMSE). The lower The MAE and RMSE values, the better the model. The specific values and related images for ARIMA (p, d, q) can be found in the Supplementary Data Sheets 1–9.
3 Results
3.1 Diabetic nephropathy caused by excessive intake of SSBs in high SDI regions
As shown in Figure 1 and Table 1, this image displays the burden analysis results of type 2 DN caused by excessive intake of SSBs in high SDI regions in2021. The results indicate that the number of YLDs exhibit a certain trend of change across different age groups. Starting from the 25- to 29-year-old age group, YLDs values gradually increase with age, reaching to1201 in the 50- to 54-year-old age group. Subsequently, there is some fluctuation in the higher age groups, but overall, the values remain at a relatively high level. There are certain differences in YLDs values between males and females in each age group, with males having slightly higher YLDs values than females in most age groups. Notably, the peaks for the age groups 65 to 69 year old and 75 to 79 year old are 1,286 and 1,060 respectively, indicating that the loss of life span due to this disease is more severe in these age groups in 2021. Moreover, the number of YLLs for males are significantly higher than those for females in most age groups. The number of deaths shows a certain increasing trend across different age groups, especially in the age group of 65 years and above, where the increase in the number of deaths is more evident. Additionally, the number of male deaths is higher than that of females in each age group, demonstrating that this disease poses a greater threat to the lives of people in the older age group and that males face a higher risk of death. The number of DALYs show a similar trend of change to YLDs and YLLs across different age groups. The relative peaks in the 50- to 54-year-old and 65- to 69-year-old age groups are 3,518 and 4,366, respectively. The DALYs values for males are generally higher than those for females in each age group, reflecting that the overall health loss caused by this disease to the male population in high SDI regions is more severe.

Figure 1. Distribution of deaths, DALYs, YLDs, and YLLs attributable to diabetic nephropathy caused by excessive sugary drink intake in high SDI regions in 2021 by age and sex.

Table 1. All-age cases and age-standardized deaths, DALYs, YLDs, and YLLs rates in 2021 for DN caused by excessive sugary drink intake in high SDI regions.
3.2 Diabetic nephropathy caused by excessive intake of SSBs in low SDI regions
As depicted in Figure 2 and Table 2, this image presents the burden analysis results of type 2 DN caused by excessive intake of SSBs in low SDI regions in 2021. The findings reveal the following that from the 25- to 29-year-old age group to the 95 year old and above age group, the number of YLDs show an overall upward trend. There is a more pronounced increase in the 55 to 59 year old and 70- to 74-year-old age groups. The differences in YLDs values between males and females across various age groups are relatively small, but males have slightly higher values than females in some age groups. Relative peaks occur in the 60 to 64 year old and 75- to 79-year-old age groups, which are 47 and 8. The YLLs values for males are higher than those for females in most age groups, indicating that the impact of this disease on life span loss in low SDI regions is more significant for males. The number of deaths shows an increasing trend across different age groups, with a more evident increase in the age group of 65 years and above. The number of male deaths is higher than that of females in each age group, suggesting that the threat of this disease to male life and health in low SDI regions is more prominent. The number of DALYs exhibit a similar trend to YLDs and YLLs across different age groups. Relative peaks appear in the 55 to 59 year old and 65- to 69-year-old age groups which are 178 and 158. The number of DALYs for males are generally higher than those for females in each age group, reflecting that the overall health loss caused by this disease to the male population in low SDI regions is more severe.

Figure 2. Distribution of deaths, DALYs, YLDs, and YLLs attributable to diabetic nephropathy caused by excessive sugary drink intake in low SDI regions in 2021 by age and sex.

Table 2. All-age cases and age-standardized deaths, DALYs, YLDs, and YLLs rates in 2021 for DN caused by excessive sugary drink intake in low SDI regions.
3.3 Age-specific rates of deaths, DALYs, YLDs, and YLLs in high SDI regions
The age-specific rates illustrated in Figure 3 reveal pronounced sex and age-related differences across deaths, DALYs, YLDs, and YLLs for type 2 DN caused by excessive intake of SSBs in high SDI regions in 2021. The age-specific rate of YLDs in high SDI regions shows a certain trend of change across different age groups. Starting from the 25- to 29-year-old age group, YLDs values gradually increase with age, reaching to 69 in the 50- to 54-year-old age group. Subsequently, there is some fluctuation in the higher age groups, but overall, the values remain at a relatively high level. There are certain differences in YLDs values between males and females in each age group, with males having slightly higher YLDs values than females in most age groups. The age-specific rate of YLLs is widely distributed across various age groups, with relatively high values ranging from the 25- to 29-year-old age group to the 95 year old and above age group. Notably, there are relatively pronounced peaks in the 65 to 69 year old and 75- to 79-year-old age groups, which are 7.203 per 100,000 and 12.337 per 100,000, indicating that the loss of life span due to this disease is more severe in these age groups. Moreover, the YLLs values for males are significantly higher than those for females in most age groups. The age-specific death rate shows a certain increasing trend across different age groups, especially in the age group of 65 years and above, where the increase is more evident. The age-specific rate of DALYs shows a similar trend of change to YLDs and YLLs across different age groups. The relative peaks in the 50 to 54 year old and 65- to 69-year-old age groups which are 6 per 100,000 and 9 per 100,000. The DALYs values for males are generally higher than those for females in each age group, reflecting that the overall health loss caused by this disease to the male population in high SDI regions is more severe.

Figure 3. Age-specific rates of deaths, DALYs, YLDs, and YLLs attributable to diabetic nephropathy caused by excessive sugary drink intake in high SDI regions in 2021 by age and sex.
3.4 Age-specific rates of deaths, DALYs, YLDs, and YLLs in low SDI regions
The age-specific rates illustrated in Figure 4 reveal pronounced sex and age-related differences across deaths, DALYs, YLDs, and YLLs for type 2 DN caused by excessive intake of SSBs in low SDI regions in 2021. The age-specific rate of YLDs in low SDI regions shows an upward trend from the 25- to 29-year-old age group to the 95 year old and above age group. There is a more noticeable increase in the 55 to 59 year old and 70- to 74-year-old age groups, which are, respectively, 0.234 per 100,000 and 0.263 per 100,000. Relative peaks occur in the 60 to 64 year old and 75- to 79-year-old age groups, which are 0.24 per 100,000 and 0.29 per 100,000. The YLLs values for males are higher than those for females in most age groups, indicating that the impact of this disease on life span loss in low SDI regions is more significant for males. The age-specific rate of death shows an increasing trend across different age groups, with a more evident increase in the age group of 65 years and above. The number of male deaths is higher than that of females in each age group, suggesting that the threat of this disease to male life and health in low SDI regions is more prominent. The age-specific rate of DALYs exhibits a similar trend to YLDs and YLLs across different age groups. Relative peaks appear in the 55 to 59 year old and 65- to 69- year-old age groups, which are 1.651 per 100,000 and 2.143 per 100,000. The DALYs values for males are generally higher than those for females in each age group, reflecting that the overall health loss caused by this disease to the male population in low SDI regions is more severe.

Figure 4. Age-specific rates of deaths, DALYs, YLDs, and YLLs attributable to diabetic nephropathy caused by excessive sugary drink intake in low SDI regions in 2021 by age and sex.
3.5 Trends in age-standardized rates of deaths, DALYs, YLDs, and YLLs in high SDI regions, 1990–2021
The age-standardized rates illustrated in Figure 5 reveal pronounced sex- and age-related differences in deaths, DALYs, YLDs, and YLLs due to DN resulting from excessive intake of SSBs in high SDI regions from 1990 to 2021. In these regions, the ASR of the burden of chronic kidney disease attributed to type 2 diabetes mellitus, caused by excessive intake of SSBs, exhibited an increasing trend from 2019 to 2021. Specifically, the ASR for deaths rose from 0.044 (95% CI: 0.023–0.069) in 2019 to 0.096 (95% CI: 0.051–0.144) in 2021. The ASR for DALYs increased from 1.425 (95% CI: 0.715–2.284) to 2.824 (95% CI: 1.462–4.174). For YLDs, the ASR increased from 0.537 (95% CI: 0.239–0.908) to 0.670 (95% CI: 0.312–1.103). Additionally, the ASR for YLLs climbed from 0.888 (95% CI: 0.458–1.390) to 2.154 (95% CI: 1.123–3.189). Among females, the ASR for deaths was 0.042 (95% CI: 0.021–0.068) in 2019 and increased to 0.088 (95% CI: 0.044–0.136) in 2021. For males, the ASR was 0.047 (95% CI: 0.024–0.074) in 2019 and rose to 0.106 (95% CI: 0.055–0.158) in 2021. This indicates a significant increase in the disease burden in high SDI regions over these 2 years. From 1990 to 2021, the age-standardized rate of YLDs in high SDI regions exhibited an upward trend. The YLDs for both males and females have increased, with a more pronounced rise observed in females compared to males. Similarly, the age-standardized rate of YLLs exhibited an upward trend from 1990 to 2021. The YLLs for both sexes have risen, with females experiencing a higher increase than males. The age-standardized rate of deaths has gradually ascended from 1990 to 2021. In 1990, the deaths was notably low; however, after more than 30 years of development, it had increased substantially by 2021. The deaths for both males and females have risen, with the increase in female mortality rates occurring at a relatively faster pace. Moreover, the age-standardized rate of DALYs has continued to rise from 1990 to 2021. The DALYs was low in 1990, but by 2021, it had significantly increased, indicating a continuous rise in overall health loss due to this disease in SDI regions. The DALYs for both males and females have increased, with a more significant rise in females, particularly among middle-aged and older adult populations.

Figure 5. Trends in age-standardized rates of deaths, DALYs, YLDs, and YLLs attributable to diabetic nephropathy caused by excessive sugary drink intake in high SDI regions, 1990–2021.
3.6 Trends in age-standardized rates of deaths, DALYs, YLDs, and YLLs in low SDI regions, 1990–2021
The age-standardized rates illustrated in Figure 6 reveal pronounced sex and age-related differences in deaths, DALYs, YLDs, and YLLs due to DN resulting from excessive intake of SSBs in low Socio-Demographic Index (SDI) regions from 1990 to 2021. Overall, from 2019 to 2021, the age-standardized rates (ASR) for deaths, DALYs, YLDs, and YLLs attributed to this cause in low SDI regions demonstrated an upward trend. Specifically, the ASR for deaths increased from 0.018 (95% CI: 0.009–0.03) in 2019 to 0.02 (95% CI: 0.01–0.032) in 2021; the ASR for DALYs rose from 0.488 (95% CI: 0.242–0.818) to 0.522 (95% CI: 0.254–0.857); the ASR for YLDs increased from 0.071 (95% CI: 0.032–0.123) to 0.081 (95% CI: 0.036–0.144); and the ASR for YLLs grew from 0.417 (95% CI: 0.206–0.706) to 0.441 (95% CI: 0.211–0.725). The increase was relatively more pronounced among females, with the ASR for deaths rising from 0.017 (95% CI: 0.008–0.027) to 0.02 (95% CI: 0.009–0.033), whereas the change among males was comparatively smaller. From 1990 to 2021, the YLDs for both males and females increased, with a more substantial rise observed in females, particularly in younger age groups. Similarly, YLLs for both sexes also increased, with females exhibiting a higher increase than males, particularly in younger cohorts. In 1990, the deaths was low; however, by 2021, it had risen significantly, indicating a growing threat posed by this disease to individuals in low SDI regions. The age-standardized DALY rate continued to rise from 1990 to 2021, reflecting a significant increase in overall health loss attributed to this disease in low SDI regions. Both males and females experienced an increase in DALY rates, with the rise being more pronounced in females than in males.

Figure 6. Trends in age-standardized rates of deaths, DALYs, YLDs, and YLLs attributable to diabetic nephropathy caused by excessive sugary drink intake in low SDI regions, 1990–2021.
3.7 The EAPC in age-standardized rates of DALYs in countries and territories of high and low SDI regions from 1990 to 2021
As illustrated in Figure 7, the number of deaths has increased in most countries and regions during this period. For instance, in the United States, the number of deaths rose from 223 in 1990 to 1,300 in 2021, representing nearly a five-fold increase. Similar upward trends are observed in other countries such as Canada, Australia, and the United Kingdom. The DALYs in high SDI regions have generally increased. Taking the United States as an example, its DALYs increased from 7,508 in 1990 to 33,065 in 2021, an increase of approximately 340%. Significant growth trends are also evident in other countries like Canada, Australia, and the United Kingdom. The YLDs have similarly increased across most countries and regions. For example, in the United States, YLDs grew from 2,573 in 1990 to 5,785 in 2021, an increase of about 125%. Comparable growth trends are observed in other countries such as Canada, Australia, and the United Kingdom. Additionally, YLLs have significantly increased in most countries and regions. For instance, YLLs in the United States surged from 4,935 in 1990 to 27,281 in 2021, representing an increase of approximately 453%. Other countries, including Canada, Australia, and the United Kingdom, also exhibit significant upward trends. As depicted in Figure 8, the number of deaths, DALYs, YLDs, and YLLs have all increased in most countries. The deaths has risen in many nations during this period, particularly in several African countries, including Burkina Faso, Chad, and the Central African Republic. In terms of DALYs, values have increased in nearly all countries, reflecting the growing overall impact of the disease on population health. Notably, the rise in DALYs in Niger and Burundi indicates significant challenges these countries face in addressing the disease. The trends observed in YLDs and YLLs are similar; the increase in YLDs signifies a growing number of individuals living with disabilities due to illness, while the rise in YLLs suggests a worsening situation regarding premature deaths caused by illness. Particularly alarming is the substantial increase in YLLs in countries such as Guinea-Bissau and Madagascar, highlighting the severe impact of the disease on life expectancy.

Figure 7. The EAPC in age-standardized rates of DALYs in countries and territories of high SDI regions from 1990 to 2021.

Figure 8. The EAPC in age-standardized rates of DALYs in countries and territories of low SDI regions from 1990 to 2021.
3.8 Predicted trends of age-standardized rates of deaths, DALYs, YLDs, and YLLs in high and low SDI regions over 20 years (2021–2041)
Figure 9 presents the actual data and projected trends of age-standardized rates in high and low SDI regions from 1990 to 2020, along with projections for 2021 to 2040. In high SDI regions, the mortality rate exhibits a consistent upward trend, increasing from 0.0978 per 100,000 in 2022 to 0.1270 per 100,000 in 2040. This trend indicates that the number of deaths attributable to DN resulting excessive intake of SSBs will continue to rise in high SDI regions in the future. DALYs also demonstrate a year-on-year increase, rising from 2.8610 per 100,000 in 2022 to 3.6290 per 100,000 in 2040, suggesting a growing health burden associated with this disease in these regions. In high SDI regions with advanced healthcare, longer life expectancy increases the proportion of older adult individuals. This aging population experiences physical decline, raising the risk of chronic diseases like cardiovascular illnesses, diabetes, and cancers, which are key contributors to and DALY increases. YLDs increased from 0.6865 per 100,000 in 2022 to 1.0704 per 100,000 in 2040, reflecting an annual increase in YLLs. The YLLs rose from 2.1831 per 100,000 in 2022 to 2.8615 per 100,000 in 2040, indicating a rise in YLLs due to premature death. Conversely, in low SDI regions, the deaths shows a consistent downward trend, decreasing from 0.0194 per 100,000 in 2022 to 0.0161 per 100,000 in 2040. This suggests a decline in the number of deaths due to type 2 diabetic nephropathy linked to excessive intake of SSBs in low SDI regions in the future. DALYs in these regions also exhibit a decreasing trend, declining from 0.5205 per 100,000 in 2022 to 0.4984 per 100,000 in 2040, indicating a gradual reduction in the health burden from the disease. YLDs increased from 0.0819 per 100,000 in 2022 to 0.1052 per 100,000 in 2040, demonstrating an annual increase in YLLs due to disability. In contrast, YLLs decreased from 0.4389 per 100,000 in 2022 to 0.3984 per 100,000 in 2040, reflecting a reduction in YLLs due to premature death.

Figure 9. Predicted trends of age-standardized rates of deaths, DALYs, YLDs, and YLLs in high and low SDI regions over 20 years (2021–2041).
3.9 Predicted trends of age-standardized rates of DALYs in low SDI and high SDI countries over 20 years (2021–2041)
As illustrated in Figures 10, 11, the burden of DN attributable to excessive intake of SSBs in high SDI regions has exhibited divergent trends from 1990 to 2021. The age-standardized rate of DALYs reveals that certain countries, such as Slovakia, Australia, and Canada, have experienced an overall upward trend during this period, indicating an increased disease burden. In contrast, countries like Iceland and Japan have demonstrated relatively stable trends, with minimal changes in disease burden. Specifically, Slovakia reported a DALYs rate of 0.6247 in 2021, projected to rise to 0.9993 by 2040, indicating a significant upward trajectory. Conversely, Australia’s DALYs was 0.7853 in 2021, projected to remain unchanged by 2040, reflecting a relatively stable disease burden. Additionally, Canada recorded a DALYs of 1.0912 in 2021, with expectations of a decrease to 0.9993 by 2040, suggesting a downward trend. As shown in Figure 11, significant disparities exist in DALYs among various low SDI countries. For instance, Somalia’s DALYs was 0.6113 in 2022, anticipated to decrease slightly each year, reaching 0.5528 by 2040. A similar trend is observed in Niger, where the DALYs was 0.2649 in 2022 and is projected to decline to 0.1975 by 2040. This suggests a potential alleviation of the burden of DN due to excessive intake of SSBs in these countries over the coming decades. However, the outlook for some nations remains less optimistic. For example, Chad’s DALYs was 0.2129 in 2022, expected to decrease only slightly to 0.1888 by 2040, indicating a relatively heavy disease burden. In contrast, Ethiopia’s DALYs is projected to rise from 0.7271 in 2022 to 1.0392 in 2040, signaling greater health challenges for the country in the future.

Figure 10. Predicted trends of age-standardized rates of DALYs in high SDI countries over 20 years (2021–2041).

Figure 11. Predicted trends of age-standardized rates of DALYs in low SDI countries over 20 years (2021–2041).
4 Discussion
The studies have shown that in 2020, about 2.2 million new cases of T2DM globally were attributable to excessive intake of SSBs, accounting for 9.8% of all new T2DM cases (29). The link between SSBs and T2DM is well-established, as the high sugar content in these beverages can lead to weight gain, insulin resistance, and ultimately the development of T2DM (30–32). DN, a common and severe microvascular complication of T2DM, is a major cause of end-stage renal disease (ESRD) (33, 34). The development of DN is influenced by a combination of factors, including hyperglycemia, hypertension, dyslipidemia, and genetic predisposition (35, 36). In addition to SSB intake, overall dietary patterns may influence the burden of DN. For example, unhealthy dietary habits such as high-salt, high-fat (especially trans and saturated fats), and low-fiber diets may act in concert with SSB consumption to boost DN risk. Also, low physical activity can cut energy expenditure, raising the risk of weight gain and obesity, which in turn affects glycemic control and insulin sensitivity. Moreover, physical activity is linked to inflammatory and oxidative stress processes, which are also key in DN development. However, the role of SSB consumption in the development and progression of DN has gained increasing attention in recent years (37). Some studies have suggested that the high fructose content in SSBs may have a direct impact on renal function, leading to renal injury and the development of DN (38, 39).
This study comprehensively analyzed the burden of DN caused by SSBs in high and low SDI regions from 1990 to 2021 and projected the trends until 2040 using ARIMA models. The results revealed distinct patterns and trends in age-standardized rates of deaths, DALYs, YLDs, and YLLs between the two regions. In high SDI regions, the age-standardized rates of deaths, DALYs, YLDs, and YLLs showed a consistent upward trend from 1990 to 2021. The deaths increased significantly, indicating a growing number of deaths attributable to DN caused by SSBs. Similarly, the DALYs, which represent the overall health burden, also rose considerably, suggesting an increasing impact of the disease on the population’s health. The YLDs and YLLs increased as well, highlighting the rising burden of disability and premature death associated with this condition. In high SDI regions, the significant upward trend in the burden of DN caused by excessive intake of SSBs can be attributed to a combination of socioeconomic and healthcare factors, for example, people in high-SDI areas have a high standard of living, but health education about sugar-sweetened beverages is not widely available. Economically, as these regions develop and urbanize, there is a marked shift in lifestyle, including dietary changes and reduced physical activity. These changes drive up excessive intake of SSBs, boosting the incidence of diabetes and its complications. While healthcare resources are relatively abundant in high SDI regions, public awareness of SSBs’ harms remains low. Current public health policies also fall short in curbing SSB consumption. In low SDI regions, despite a growing DN burden, the pace of increase is slower and some indicators even show a decline. This might stem from enhanced public health policies in recent years. Some countries have intensified health education and publicity to boost awareness of healthy eating. Alongside, efforts to bolster healthcare service accessibility have improved early diabetes detection and diagnosis, partially mitigating DN’s growth. The projections for 2021 to 2040 indicated that these upward trends would continue, with further increases in deaths, DALYs, YLDs, and YLLs, suggesting a worsening public health challenge in the future. In contrast, low SDI regions exhibited different trends. Although the age-standardized rates of deaths, DALYs, YLDs, and YLLs also increased from 1990 to 2021, the pace of increase was relatively slower compared to high SDI regions. However, projections showed that from 2021 to 2040, the mortality rate and DALYs would start to decline, while YLDs would continue to rise slightly and YLLs would decrease. This suggests that while the burden of DN caused by SSBs remains a concern in low SDI regions, there may be potential for improvement in certain aspects such as reducing mortality and overall health burden through targeted interventions.
Our findings are consistent with previous research that has documented the increasing prevalence and burden of type 2 diabetes and its complications in high SDI regions (40). Studies have shown that excessive intake of SSBs is closely associated with the development of type 2 diabetes and its subsequent complications, including DN (41, 42). The rising trends in mortality, DALYs, YLDs, and YLLs observed in this study further emphasize the need for effective public health measures to address this issue. In low SDI regions, limited data has been available on the burden of DN caused by SSBs. Our study provides valuable insights into the trends and projections in these regions, which can inform future research and policy-making. The results of this study have important implications for public health policy. In high SDI regions, the continuous increase in the burden of DN caused by excessive sugary drink intake calls for urgent action. Policies aimed at reducing sugary drink consumption, such as taxation, marketing restrictions, and public education campaigns, should be strengthened (43). Additionally, improving access to healthcare services for early diagnosis and management of type 2 diabetes and its complications can help mitigate the rising burden of DN. In low SDI regions, while the overall burden may be lower compared to high SDI regions, the increasing trends still warrant attention. Public health interventions should focus on raising awareness about the risks of excessive intake of SSBs and promoting healthy dietary habits. Strengthening healthcare infrastructure and capacity building can also contribute to better prevention and management of DN in these regions (44).
In conclusion, this study provides a comprehensive analysis of the burden of type 2 diabetes nephropathy caused by excessive intake of SSBs in high and low SDI regions. The findings highlight the significant and increasing public health challenge posed by this condition, particularly in high SDI regions. In high SDI regions, higher taxes can be imposed on SSBs, with revenues used for public health education and diabetes prevention/treatment projects. For example, Mexico’s tax policy reduced SSB consumption. Limiting SSB advertising, especially to children/adolescents, in media and public places. Norway strictly restricts such advertising. A multidisciplinary DN management team should be established, including endocrinologists, nephrologists, dietitians, and nurses, to provide personalized treatment plans. In low SDI regions, while the trends are relatively more favorable, continued attention and appropriate measures are still required to prevent and control the burden of DN caused by excessive sugary drink intake, governments can encourage businesses to produce low- or no-sugar beverage alternatives via policies and tax incentives, increasing their market availability. Invest in primary healthcare by upgrading facilities, equipping them better, and training healthcare providers to enhance their DN management capabilities. Schools and families should work together to develop students’ healthy eating habits and reduce their consumption of SSBs.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary material, further inquiries can be directed to the corresponding authors.
Author contributions
XW: Software, Writing – original draft, Formal analysis, Data curation, Methodology, Validation, Project administration. GC: Supervision, Writing – original draft, Investigation, Validation. QO: Funding acquisition, Writing – review & editing, Supervision, Visualization, Conceptualization. XC: Supervision, Writing – review & editing, Visualization, Funding acquisition.
Funding
The author(s) declare that financial support was received for the research and/or publication of this article. This work was supported by the National Natural Science Foundation of China (Nos. 82030025 and 32141005), Beijing Natural Science Foundation (No. L232122.). All data can be extracted from the Global Health Data Exchange query tool (http://ghdx.healthdata.org/gbd-results-tool), accessed on 14 January 2025.
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 author(s) declare that no Gen AI was used in the creation of this manuscript.
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/fpubh.2025.1598278/full#supplementary-material
Footnotes
References
1. 2016 Risk Factors Collaborators. 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:1345–422. doi: 10.1016/S0140-6736(17)32366-8
2. GBD 2015 Risk Factors Collaborators. Global, regional, and national comparative risk assessment of 79 behavioural, environmental and occupational, and metabolic risks or clusters of risks, 1990-2015: a systematic analysis for the global burden of disease study 2015. Lancet. (2016) 388:1659–724. doi: 10.1016/S0140-6736(16)31679-8
3. Stepien, M, Duarte-Salles, T, Fedirko, V, Trichopoulou, A, Lagiou, P, Bamia, C, et al. Consumption of soft drinks and juices and risk of liver and biliary tract cancers in a European cohort. Eur J Nutr. (2016) 55:7–20. doi: 10.1007/s00394-014-0818-5
4. Schulze, MB, Manson, JE, Ludwig, DS, Colditz, GA, Stampfer, MJ, Willett, WC, et al. Sugar-sweetened beverages, weight gain, and incidence of type 2 diabetes in young and middle-aged women. JAMA. (2004) 292:927–34. doi: 10.1001/jama.292.8.927
5. Hu, H, Liang, W, and Ding, G. Ion homeostasis in diabetic kidney disease. Trends Endocrinol Metab. (2024) 35:142–50. doi: 10.1016/j.tem.2023.09.009
6. Watanabe, K, Sato, E, Mishima, E, Miyazaki, M, and Tanaka, T. What's new in the molecular mechanisms of diabetic kidney disease: recent advances. Int J Mol Sci. (2022) 24:570. doi: 10.3390/ijms24010570
7. Limonte, CP, Kretzler, M, Pennathur, S, Pop-Busui, R, and de Boer, IH. Present and future directions in diabetic kidney disease. J Diabetes Complicat. (2022) 36:108357. doi: 10.1016/j.jdiacomp.2022.108357
8. Song, S, Cai, X, Shen, D, Hu, J, Zhu, Q, Ma, H, et al. Relationship between cumulative exposure and time course of plasma aldosterone concentrations and chronic kidney disease in hypertensive patients: a real-world cohort study. Am J Nephrol. (2025) 1:1–22. doi: 10.1159/000545451
9. Barbaresko, J, Neuenschwander, M, Schwingshackl, L, and Schlesinger, S. Dietary factors and diabetes-related health outcomes in patients with type 2 diabetes: protocol for a systematic review and meta-analysis of prospective observational studies. BMJ Open. (2019) 9:e027298. doi: 10.1136/bmjopen-2018-027298
10. Nangaku, M, Izuhara, Y, Usuda, N, Inagi, R, Shibata, T, Sugiyama, S, et al. In a type 2 diabetic nephropathy rat model, the improvement of obesity by a low calorie diet reduces oxidative/carbonyl stress and prevents diabetic nephropathy. Nephrol Dial Transplant. (2005) 20:2661–9. doi: 10.1093/ndt/gfi096
11. Shoham, DA, Durazo-Arvizu, R, Kramer, H, Luke, A, Vupputuri, S, Kshirsagar, A, et al. Sugary soda consumption and albuminuria: results from the National Health and nutrition examination survey, 1999-2004. PLoS One. (2008) 3:e3431. doi: 10.1371/journal.pone.0003431
12. Raben, A, Møller, BK, Flint, A, Vasilaris, TH, Christina Møller, A, Juul Holst, J, et al. Increased postprandial glycaemia, insulinemia, and lipidemia after 10 weeks' sucrose-rich diet compared to an artificially sweetened diet: a randomised controlled trial. Food Nutr Res. (2011) 55:5961. doi: 10.3402/fnr.v55i0.5961
13. Janssens, JP, Shapira, N, Debeuf, P, Michiels, L, Putman, R, Bruckers, L, et al. Effects of soft drink and table beer consumption on insulin response in normal teenagers and carbohydrate drink in youngsters. Eur J Cancer Prev. (1999) 8:289–96.
14. Kuzma, JN, Cromer, G, Hagman, DK, Breymeyer, KL, Roth, CL, Foster-Schubert, KE, et al. Consuming glucose-sweetened, not fructose-sweetened, beverages increases fasting insulin in healthy humans. Eur J Clin Nutr. (2019) 73:487–90. doi: 10.1038/s41430-018-0297-5
15. Thomas, G, Sehgal, AR, Kashyap, SR, Srinivas, TR, Kirwan, JP, and Navaneethan, SD. Metabolic syndrome and kidney disease: a systematic review and meta-analysis. Clin J Am Soc Nephrol. (2011) 6:2364–73. doi: 10.2215/CJN.02180311
16. Balkau, B, Deanfield, JE, Després, JP, Bassand, JP, Fox, KA, Smith, SC Jr, et al. International day for the evaluation of abdominal obesity (IDEA): a study of waist circumference, cardiovascular disease, and diabetes mellitus in 168,000 primary care patients in 63 countries. Circulation. (2007) 116:1942–51. doi: 10.1161/CIRCULATIONAHA.106.676379
17. Chapman, CL, Grigoryan, T, Vargas, NT, Reed, EL, Kueck, PJ, Pietrafesa, LD, et al. High-fructose corn syrup-sweetened soft drink consumption increases vascular resistance in the kidneys at rest and during sympathetic activation. Am J Physiol Ren Physiol. (2020) 318:F1053–65. doi: 10.1152/ajprenal.00374.2019
18. Hashemi, SY, Momenabadi, V, Faramarzi, A, and Kiani, A. Trends in burden of chronic obstructive pulmonary disease in Iran, 1995-2015: findings from the global burden of disease study. Arch Public Health. (2020) 78:45. doi: 10.1186/s13690-020-00426-x
19. Rahmani, H, Sarabi Asiabar, A, Niakan, S, Hashemi, SY, Faramarzi, A, Manuchehri, S, et al. Burden of esophageal cancer in Iran during 1995-2015: review of findings from the global burden of disease studies. Med J Islam Repub Iran. (2018) 32:55. doi: 10.14196/mjiri.32.55
20. Zou, J, Sun, T, Song, X, Liu, YM, Lei, F, Chen, MM, et al. Distributions and trends of the global burden of COPD attributable to risk factors by SDI, age, and sex from 1990 to 2019: a systematic analysis of GBD 2019 data. Respir Res. (2022) 23:90. doi: 10.1186/s12931-022-02011-y
21. Zhang, R, Liu, H, Pu, L, Zhao, T, Zhang, S, Han, H, et al. Global burden of ischemic stroke in young adults in 204 countries and territories. Neurology. (2024) 102:e209328. doi: 10.1212/WNL.0000000000209328
22. Dong, C, Bu, X, Liu, J, Wei, L, Ma, A, and Wang, T. Cardiovascular disease burden attributable to dietary risk factors from 1990 to 2019: a systematic analysis of the global burden of disease study. Nutr Metab Cardiovasc Dis. (2022) 32:897–907. doi: 10.1016/j.numecd.2021.11.012
23. GBD 2021 US Burden of Disease and Forecasting Collaborators. Burden of disease scenarios by state in the USA, 2022-50: a forecasting analysis for the global burden of disease study 2021. Lancet. (2024) 404:2341–70. doi: 10.1016/S0140-6736(24)02246-3
24. GBD 2021 Fertility and Forecasting Collaborators. Global fertility in 204 countries and territories, 1950-2021, with forecasts to 2100: a comprehensive demographic analysis for the global burden of disease study 2021. Lancet. (2024) 403:2057–99. doi: 10.1016/S0140-6736(24)00550-6
25. Nahvijou, A, Arab, M, Faramarzi, A, Hashemi, Y, and Javan-Noughabi, J. Burden of esophageal cancer according to World Health Organization regions: review of findings from the global burden of disease study 2015. Health Scope. (2019) 8:55. doi: 10.5812/jhealthscope.64984
26. Murray, CJL. The global burden of disease study at 30 years. Nat Med. (2022) 28:2019–26. doi: 10.1038/s41591-022-01990-1
27. Rostami, R, Moradinazar, M, Moradi, S, Samannejad, B, Cheshmeh, S, Saber, A, et al. Impact of dietary risk on global ischemic heart disease: findings from 1990-2019. Sci Rep. (2024) 14:18012. doi: 10.1038/s41598-024-69089-w
28. Cai, Q, Chen, Y, Xin, S, Zhang, D, Pan, J, Xie, Z, et al. Temporal trends of bladder cancer incidence and mortality from 1990 to 2016 and projections to 2030. Transl Androl Urol. (2020) 9:153–65. doi: 10.21037/tau.2020.02.24
29. Lara-Castor, L, O'Hearn, M, Cudhea, F, Miller, V, Shi, P, Zhang, J, et al. Burdens of type 2 diabetes and cardiovascular disease attributable to sugar-sweetened beverages in 184 countries. Nat Med. (2025) 31:552–564. doi: 10.1038/s41591-024-03345-4
30. Tseng, TS, Lin, WT, Gonzalez, GV, Kao, YH, Chen, LS, and Lin, HY. Sugar intake from sweetened beverages and diabetes: a narrative review. World J Diabetes. (2021) 12:1530–8. doi: 10.4239/wjd.v12.i9.1530
31. Henriksen, HB, and Kolset, SO. Sugar intake and public health. Tidsskr Nor Laegeforen. (2007) 127:2259–62.
32. Xu, F, Park, S, and Siegel, KR. Factors associated with frequency of sugar-sweetened beverage consumption among US adults with diabetes or prediabetes. Am J Health Promot. (2018) 32:1489–97. doi: 10.1177/0890117117746187
33. Jaber, M, Sharabati, A, Nofal, K, Hassan, M, Hamdan, Z, and Nazzal, Z. Decline in eGFR and mortality among type II diabetic patients: a 3-year prospective cohort study from Palestine. BMC Nephrol. (2025) 26:64. doi: 10.1186/s12882-025-03947-1
34. Oh, YJ, An, JN, Kim, CT, Yang, SH, Lee, H, Kim, DK, et al. Circulating tumor necrosis factor α receptors predict the outcomes of human IgA nephropathy: a prospective cohort study. PLoS One. (2015) 10:e0132826. doi: 10.1371/journal.pone.0132826
35. Piko, N, Bevc, S, Ekart, R, Petreski, T, Vodošek Hojs, N, and Hojs, R. Diabetic patients with chronic kidney disease: non-invasive assessment of cardiovascular risk. World J Diabetes. (2021) 12:975–96. doi: 10.4239/wjd.v12.i7.975
36. Wascher, TC, Paulweber, B, Toplak, H, Saely, CH, Drexel, H, Föger, B, et al. Lipids-diagnosis and therapy in diabetes mellitus (update 2019). Wien Klin Wochenschr. (2019) 131:136–8. doi: 10.1007/s00508-019-1444-2
37. Audain, K, Levy, L, and Ellahi, B. Sugar-sweetened beverage consumption in the early years and implications for type-2 diabetes: a sub-Saharan Africa context. Proc Nutr Soc. (2019) 78:547–53. doi: 10.1017/S0029665118002860
38. Bjornstad, P, Lanaspa, MA, Ishimoto, T, Kosugi, T, Kume, S, Jalal, D, et al. Fructose and uric acid in diabetic nephropathy. Diabetologia. (2015) 58:1993–2002. doi: 10.1007/s00125-015-3650-4
39. Meléndez-Salcido, CG, Ramírez-Emiliano, J, and Pérez-Vázquez, V. Hypercaloric diet promotes metabolic disorders and impaired kidney function. Curr Pharm Des. (2022) 28:3127–39. doi: 10.2174/1381612829666221020162955
40. Sun, H, Liu, Y, Xu, Y, Gong, X, Zhai, C, Hu, W, et al. Global disease burden attributed to high sugar-sweetened beverages in 204 countries and territories from 1990 to 2019. Prev Med. (2023) 175:107690. doi: 10.1016/j.ypmed.2023.107690
41. Malik, VS, Popkin, BM, Bray, GA, Després, JP, Willett, WC, and Hu, FB. Sugar-sweetened beverages and risk of metabolic syndrome and type 2 diabetes: a meta-analysis. Diabetes Care. (2010) 33:2477–83. doi: 10.2337/dc10-1079
42. Drouin-Chartier, JP, Zheng, Y, Li, Y, Malik, V, Pan, A, Bhupathiraju, SN, et al. Changes in consumption of sugary beverages and artificially sweetened beverages and subsequent risk of type 2 diabetes: results from three large prospective U.S. cohorts of women and men. Diabetes Care. (2019) 42:2181–9. doi: 10.2337/dc19-0734
43. Mangera, KAS, and Adams, OP. Knowledge, attitudes and practices with regard to sugar sweetened beverages and taxation among people with type 2 diabetes mellitus in the Caribbean island of Barbados - a cross sectional survey in primary care. Prim Care Diabetes. (2021) 15:69–73. doi: 10.1016/j.pcd.2020.04.002
44. Zhang, XX, Kong, J, and Yun, K. Prevalence of diabetic nephropathy among patients with type 2 diabetes mellitus in China: a Meta-analysis of observational studies. J Diabetes Res. (2020) 2020:1–11. doi: 10.1155/2020/2315607
45. GBD 2017 Causes of Death Collaborators. Health effects of dietary risks in 195 countries, 1990-2017: a systematic analysis for the global burden of disease study 2017. Lancet. (2019) 393:1958–72. doi: 10.1016/S0140-6736(19)30041-8
46. GBD 2017 Causes of Death Collaborators. 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:1736–88. doi: 10.1016/S0140-6736(18)32203-7
47. Li, X, Xiao, X, Wu, Z, Li, A, Wang, W, and Lin, R. Global, regional, and national burden of early-onset colorectal cancer and projection to 2050: An analysis based on the global burden of disease study 2021. Public Health. (2025) 238:245–53. doi: 10.1016/j.puhe.2024.12.011
Keywords: diabetic nephropathy, excessive intake of sugar-sweetened beverages, socio-demographic index, deaths, DALYs, YLDs, YLLs, ARIMA
Citation: Wang X, Cai G, Ouyang Q and Chen X (2025) Burden analysis of diabetic nephropathy caused by excessive intake of sugar-sweetened beverages in high and low SDI regions. Front. Public Health. 13:1598278. doi: 10.3389/fpubh.2025.1598278
Edited by:
Jianbo Xiao, University of Vigo, SpainReviewed by:
Seyed Yaser Hashemi, Tehran University of Medical Sciences, IranHuimin Ma, People’s Hospital of Xinjiang Uygur Autonomous Region, China
Copyright © 2025 Wang, Cai, Ouyang and Chen. 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: Qing Ouyang, bm9saW1pdGh5ZXIzMTY5QG91dGxvb2suY29t; Xiangmei Chen, eG1jaGVuMzAxQDEyNi5jb20=