The influence of overall lifestyle behaviors on diabetic microvascular complications remains unknown. In addition, the potential mediating biomarkers underlying the association is unclear. This study aimed to examine the associations of the combined lifestyle factors with risks of total and individual microvascular complications among patients with type 2 diabetes (T2D) and to explore the potential mediation effects of metabolic biomarkers.
Methods and findings
This retrospective cohort study included 15,104 patients with T2D free of macro- and microvascular complications at baseline (2006 to 2010) from the UK Biobank. Healthy lifestyle behaviors included noncurrent smoking, recommended waist circumference, regular physical activity, healthy diet, and moderate alcohol drinking. Outcomes were ascertained using electronic health records. Over a median of 8.1 years of follow-up, 1,296 cases of the composite microvascular complications occurred, including 558 diabetic retinopathy, 625 diabetic kidney disease, and 315 diabetic neuropathy, with some patients having 2 or 3 microvascular complications simultaneously. After multivariable adjustment for sociodemographic characteristics, history of hypertension, glycemic control, and medication histories, the hazard ratios (95% confidence intervals (CIs)) for the participants adhering 4 to 5 low-risk lifestyle behaviors versus 0 to 1 were 0.65 (0.46, 0.91) for diabetic retinopathy, 0.43 (0.30, 0.61) for diabetic kidney disease, 0.46 (0.29, 0.74) for diabetic neuropathy, and 0.54 (0.43, 0.68) for the composite outcome (all Ps-trend ≤0.01). Further, the population-attributable fraction (95% CIs) of diabetic microvascular complications for poor adherence to the overall healthy lifestyle (<4 low-risk factors) ranged from 25.3% (10.0%, 39.4%) to 39.0% (17.7%, 56.8%). In addition, albumin, HDL-C, triglycerides, apolipoprotein A, C-reactive protein, and HbA1c collectively explained 23.20% (12.70%, 38.50%) of the associations between overall lifestyle behaviors and total diabetic microvascular complications. The key limitation of the current analysis was the potential underreporting of microvascular complications because the cases were identified via electronic health records.
Adherence to overall healthy lifestyle behaviors was associated with a significantly lower risk of microvascular complications in patients with T2D, and the favorable associations were partially mediated through improving biomarkers of glycemic control, systemic inflammation, liver function, and lipid profile.
Why was this study done?
- Diabetic microvascular complications including diabetic retinopathy, diabetic neuropathy, and diabetic kidney disease have placed a significant disease burden.
- This study asked whether multiple lifestyle behaviors including waist circumference, smoking status, habitual diet, physical activity, and alcohol intake could jointly affect the microvascular complications, and whether the metabolic biomarkers could mediate the association.
What did the researchers do and find?
- In this retrospective study of 15,104 patients with type 2 diabetes, participants adhering 4 to 5 low-risk lifestyle behaviors versus 0 to 1 had a 35% lower risk of diabetic retinopathy, 57% lower risk of diabetic kidney disease, 54% lower risk of diabetic neuropathy, and 46% lower risk of the composite outcome.
- Our study found that albumin, HDL-C, triglycerides, apolipoprotein A, C-reactive protein, and HbA1c collectively explained 23.20% of the associations between overall lifestyle behaviors and diabetic microvascular complications.
What do these findings mean?
- Our findings support the importance of public health programs and interventions targeting at improving health behaviors in combination to ameliorate the risk of diabetic microvascular complications.
Citation: Geng T, Zhu K, Lu Q, Wan Z, Chen X, Liu L, et al. (2023) Healthy lifestyle behaviors, mediating biomarkers, and risk of microvascular complications among individuals with type 2 diabetes: A cohort study. PLoS Med 20(1):
Academic Editor: Weiping Jia, Shanghai Jiao Tong University Affiliated Sixth People’s Hospital, CHINA
Received: June 8, 2022; Accepted: October 26, 2022; Published: January 10, 2023
Copyright: © 2023 Geng et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Data Availability: The UK Biobank data are available on application to the UK Biobank (www.ukbiobank.ac.uk/). This research has been conducted using the UK Biobank Resource under Application Number 68307.
Funding: GL is supported by grants from National Nature Science Foundation of China (82073554), the Hubei Province Science Fund for Distinguished Young Scholars (2021CFA048), and the Fundamental Research Funds for the Central Universities (2021GCRC076). AP is supported by grants from National Nature Science Foundation of China (81930124, 82021005, and 82192902), and the Fundamental Research Funds for the Central Universities (2021GCRC075). The Funder had no role in the design and conduct of the study; collection, management, analysis, and interpretation of the data; preparation, review, or approval of the manuscript; and decision to submit the manuscript for publication.
Competing interests: The authors have declared that no competing interests exist.
alkaline phosphatase; ALT,
alanine aminotransferase; AST,
aspartate aminotransferase; BMI,
body mass index; CI,
confidence interval; CRIC,
Chronic Renal Insufficiency Cohort; CRP,
C-reactive protein; CVD,
cardiovascular disease; DBP,
diastolic blood pressure; DIY,
estimated glomerular filtration rate; GGT,
gamma glutamyltransferase; HDL-C,
high-density lipoprotein cholesterol; HR,
hazard ratio; LDL-C,
low-density lipoprotein cholesterol; MET,
metabolic equivalent of task; ONTARGET,
Ongoing Telmisartan Alone and in Combination with Ramipril Global Endpoint Trial; PAF,
population-attributable fraction; SBP,
systolic blood pressure; TC,
total cholesterol; T2D,
type 2 diabetes; WC,
Diabetes is a global public health crisis affecting greater than 0.5 billion adults worldwide . Diabetic microvascular complications including diabetic retinopathy, diabetic neuropathy, and diabetic kidney disease have placed a significant health and economic burden borne by individuals, families, and health systems [2,3]. For example, diabetic retinopathy, the leading cause of vision loss, is present in nearly 30% of patients with diabetes . Furthermore, both diabetic kidney disease and diabetic neuropathy may develop in approximately 50% of patients with diabetes [5,6]. Therefore, it is paramount to identify cost-effective strategies to prevent and delay the development of microvascular complications in patients with diabetes.
Beyond the glucose control by medications, the American Diabetes Association guideline has highlighted that both caregivers and patients should focus on how to optimize lifestyle behaviors to improve diabetes care . Although lifestyle behaviors that are generally recommended, e.g., normal weight, no smoking, moderate alcohol drinking, healthy diet, and physically active, have been associated with lower risks of microvascular complications [8–14], to our best knowledge, the magnitudes of the joint association of multiple lifestyle factors with the development of microvascular complications in diabetes have not yet been quantified, which may have substantial public health implications on translating epidemiological findings to meaningful public health actions. In addition, several studies have linked lifestyle behaviors with a range of intermediate variables including lipid profile [15,16], liver function biomarkers [15,17–19], renal function biomarkers [20,21], blood pressure indices , glucose metabolism measures , and systemic inflammatory factors [15,16]; however, whether and the extent to which these metabolic biomarkers could mediate the association between lifestyle behaviors and diabetic microvascular complications remains unclear.
To shed light on the potential favorable association of overall lifestyle behaviors on microvascular complications in patients with diabetes, we examined the joint association of multiple lifestyle behaviors, including waist circumference (WC), smoking status, habitual diet, physical activity, and alcohol intake with risks of total microvascular complications, diabetic retinopathy, diabetic neuropathy, and diabetic kidney disease among patients with type 2 diabetes (T2D) who participated in the UK Biobank study. In addition, we also comprehensively evaluated the effect of a series of blood biomarkers on mediating the relationship between lifestyle behaviors and diabetic microvascular complications.
The UK Biobank is a large community-based prospective cohort study for common diseases of middle and older adults including over 500,000 participants aged 37 to 73 years from 22 sites across England, Scotland, and Wales between March 2006 and October 2010. Extensive data were obtained through touchscreen questionnaires, physical measurements, and biological samples at recruitment. Specific methods of data collection have been described previously [24,25].
Our sample of 15,104 was generated by including patients with T2D identified by using the algorithms method developed by the UK Biobank study  and excluding participants with prevalent macro- or microvascular complication cases, had incomplete information on lifestyle behaviors, or withdrawal from the study. The flowchart of patients included in the current study is present in S1 Fig.
The study was approved by the North West Multi-Centre Research Ethics Committee, the National Information Governance Board for Health and Social Care in England and Wales, and the Community Health Index Advisory Group in Scotland. All participants provided written informed consent. In the current analysis, we employed the UK Biobank study to test a priori hypothesis; we did not publish an analysis plan before conducting analyses between January 2022 and March 2022. The associations between lifestyle factors and the risk of microvascular complications in participants without excluding those with macrovascular complications and stratified analysis by preexisting cardiovascular disease (CVD) status were performed in response to peer review in July 2022. This study is reported as per the Strengthening the Reporting of Observational Studies in Epidemiology (STROBE) guideline (S1 Checklist).
Measurements of lifestyle behaviors
Five lifestyle behaviors, namely, WC, smoking status, physical activity, habitual diet, and alcohol intake, were evaluated in the current analysis. We used WC instead of body mass index (BMI) to avoid the potential obesity paradox [27,28] as evidence found an obesity paradox when obesity was measured by BMI but not when measured by WC in patients with diabetes . WC was measured using the Wessex nonstretchable sprung tape measurement, and low-risk WC was defined as <80 cm for women and <94 cm for men [30,31]. Data on smoking status were self-reported, and noncurrent smoking was defined as low-risk behavior. The frequency of all types of alcohol intake was reported using 6 predefined categories, between never to daily or almost daily. For participants who reported to drink alcohol, data on the average monthly or weekly alcohol intake from 6 types of alcohol beverages were collected. We calculated the average units of alcohol intake using the abovementioned information and defined low-risk drinking as moderate drinking (1 to 14 g/day for women or 1 to 28 g/day for men). Data on the type and duration of physical activity were derived from the questionnaire. Leisure-time physical activity score based on the 5 activities undertaken in the last 4 weeks was computed by multiplying the metabolic equivalent of task [MET] score of each activity by the minutes performed [32,33]. Light DIY (do-it-yourself), walking for pleasure, other exercises (e.g., swimming, cycling, keep fit, bowling), heavy DIY, and strenuous sports were given 1.5, 3.5, 4.0, 5.5, and 8.0 METs, respectively . The midpoints of the frequency and duration of physical activities were used to calculate the time spent on each activity. We then classified the top third of the physical activity score as the low-risk group. In addition, we generated a dietary score to reflect the overall diet quality including 10 components, namely, fruits, vegetables, whole grains, fish, dairy, vegetable oils, refined grains, processed meat, unprocessed meat, and sugar-sweetened beverages. Low-risk diet was defined as meeting 5 or more ideal diet components . Participants with each low-risk behavior were assigned 1 point; otherwise, 0 points. The overall healthy lifestyle score was the sum of individual score of the 5 lifestyle behaviors, ranging from 0 to 5, with higher score indicating healthier lifestyle.
Assessment of the circulating biomarkers
Blood samples were collected from consenting participants at recruitment, separated by components and stored at UK Biobank (−80°C and LN2) until analysis. Blood biomarkers were externally validated with stringent quality control in the UK Biobank; full details on assay performance have been given elsewhere . We selected the potential biological biomarkers mediating the association between lifestyle factors and microvascular complications based on knowledge of potential pathways, including glycemic control determined by glycated hemoglobin (HbA1c), lipid profile (total cholesterol [TC], high-density lipoprotein cholesterol [HDL-C], low-density lipoprotein cholesterol [LDL-C], triglycerides, apolipoprotein A, apolipoprotein B, and lipoprotein A), liver function (alanine aminotransferase [ALT], alkaline phosphatase [ALP], aspartate aminotransferase [AST], gamma glutamyltransferase [GGT], total bilirubin, total protein, and albumin), renal function (cystatin C, creatinine, urate, and urea), inflammation (C-reactive protein [CRP], and white blood cell count), and blood pressure indices (systolic blood pressure [SBP] and diastolic blood pressure [DBP]).
Ascertainment of outcomes
Diabetic retinopathy (ICD-10: E113, E143, H280, H360), diabetic neuropathy (ICD-10: E114, E144, G590, G629, G632, G990), and diabetic kidney disease (ICD-10: E112, E142, N180, N181, N182, N183, N184, N185, N188, N189) cases were identified through linking the cohort database with the hospital inpatient admissions and death registries.
Comparisons of baseline characteristics across the categories of the overall healthy lifestyle score were made using ANOVA or chi-squared test. We also compared the differences between patients included in the current analysis and those who were excluded due to missing values. Person-years were calculated from the date of recruitment to the date of death, first endpoint, lost to follow-up, or the end of follow-up, whichever came first. The lost to follow-up variable in the UK Biobank has been created by amalgamating data from 5 possible sources: (1) Death reported to UK Biobank by a relative; (2) NHS records indicate they are lost to follow-up; (3) NHS records indicate they have left the UK; (4) UK Biobank sources report they have left the UK; (5) Participant has withdrawn consent for future linkage. The end of follow-up dates were 1 April 2017, 17 September 2016, and 1 November 2016, for centers in England, Wales, and Scotland, respectively. Cox proportional hazards regression models were used to calculated hazard ratios (HRs) and 95% confidence intervals (CIs) for the associations of individual lifestyle behaviors and overall healthy lifestyle score with risks of total and individual microvascular complications in patients with T2D. We imputed the missing values of covariates (≤7%) using multiple imputations by chained equations with 5 imputations (SAS PROC MI with a fully conditional specification method and PROC MIANALYZE). Linear regression model and logistic regression model with all the covariates in the fully adjusted model were used to impute continuous variables and categorical variables, respectively. The percentage of missing values are present in S1 Table.
Three models were built. In Model 1, we adjusted for age (continuous, years), sex (male, female), Townsend Deprivation Index (continuous), and race/ethnicity (White, others). In Model 2, we further adjusted for education attainment (college or university degree, A/AS levels or equivalent or O levels/GCSEs, NVQ or HND or HNC or equivalent or other professional qualifications, none of the above), sleep duration (<6, 6 to 8, or ≥9 hours/day), family history of CVD (yes, no), family history of hypertension (yes, no), and prevalence of hypertension (yes, no). Finally, in Model 3, diabetes duration (continuous, years), HbA1c (continuous, mmol/mol), use of diabetes medication (none, only oral medicine, insulin, and others), use of antihypertensive medication (yes, no), use of lipid-lowing medication (yes, no), and use of aspirin (yes, no) were additionally adjusted. Further, restricted cubic spline analysis was applied to test dose–response relationships between the healthy lifestyle score and risks of outcomes. We also calculated the population-attributable fractions (PAFs) using the %par SAS Macro (https://www.hsph.harvard.edu/donna-spiegelman/software/par/) to estimate the proportion of microvascular complications that could theoretically be avoided if all participants adhered to 4 or more low-risk lifestyle behaviors.
Mediation effects of biomarkers on the associations of overall lifestyle score with risks of total and individual microvascular complications were evaluated using mediation package in R. Indirect, direct, and total effects for each mediator were computed via combining the mediator and outcome models with the adjustment of all the covariates in Model 3. Nonparametric bootstrap resampling was used to compute the CIs of the proportions of mediations. We selected the available biomarkers from the UK Biobank for the mediation analyses based on knowledge of potential causal pathways to predisposing to microvascular complications or mortality [19,37–40]. The selected biomarkers were considered as potential mediators following two-step analysis. First, we assessed the associations of all biomarkers with the overall lifestyle score using the multivariable-adjusted linear regression models. Second, we evaluated the associations of biomarkers that were significantly associated with the overall lifestyle score, with risks of all the outcomes using the multivariable-adjusted Cox regression model. We then chose the biomarkers significantly associated with each outcome for the mediation analysis accordingly.
In addition, stratified analyses were conducted by age (≤60, >60 years), sex (female, male), education (less than college, college, or above), diabetes duration (≤3, >3 years), use of diabetes medication (yes, no), and HbA1c (≤53, >53 mmol/mol). Interactions between the overall healthy lifestyle score and stratified factors on the risk of outcomes were examined using the likelihood ratio test by adding product terms in the multivariable-adjusted Cox models. Further, we examined the associations of different combinations of low-risk lifestyle behaviors with outcomes.
Several sensitivity analyses were conducted to test the robustness of our results. First, to minimize the potential reverse causation, we performed the analysis among patients with T2D after excluding the cases that occurred within 2 years of follow-up. Second, we generated the overall lifestyle score using low-risk drinking defined as moderate alcohol drinking and never drinking and repeated the main analysis using the new lifestyle score. Third, we constructed the healthy lifestyle score using BMI or waist-to-hip ratio instead of WC. Fourth, we generated a weighted healthy lifestyle score and examined the associations of the weighted healthy lifestyle score with risks of outcomes. Fifth, we investigated the association between the overall lifestyle score and risk of diabetic kidney disease, and mediation analysis for diabetic kidney disease with additional adjustment for kidney function biomarkers. Sixth, we performed the analysis via including the patients with CVD (n = 3,397) at baseline and stratified the associations by preexisting CVD status. Finally, given the potential competing risk of death highlighted during the peer review process, we assessed the associations of healthy lifestyle score with risks of microvascular complications using both the cause-specific hazard model and Fine and Gray subdistribution methods.
We used SAS V.9.4 and R software version 4.0.2 (R Foundation for Statistical Computing) for all statistical analyses. A two-tailed P < 0.05 was considered to be statistically significant.
Among 15,104 participants with T2D (60.3% male; mean age, 59.3 years), there were 3,406 (22.6%), 6,080 (40.3%), 4,062 (26.9%), 1,556 (10.3%) having 0 or 1, 2, 3, and 4 or 5 low-risk lifestyle behaviors, respectively. The baseline characteristics are shown in Table 1. Participants with more low-risk lifestyle behaviors were more likely to be men, White, less deprived, highly educated, sleep recommended hours, have a lower level of HbA1c, and have a lower prevalence of hypertension. They were less likely to use aspirins and medications for diabetes, dyslipidemia, and hypertension. In addition, compared the participants who were excluded due to missing values, those included in the current analysis were more likely to be men, White, less deprived, highly educated, noncurrent smokers, physically active, moderate alcohol drinkers, and eat healthier (S2 Table).
Lifestyle behaviors and outcomes
During 117,445 person-years of follow-up (median 8.1 years; interquartile range 7.3 to 8.8 years; maximum 11.9 years), there occurred 1,639 (10.9%) deaths and 1,296 (8.6%) composite microvascular complications cases, including 558 (3.7%) diabetic retinopathy, 625 (4.1%) diabetic kidney disease, and 315 (2.1%) diabetic neuropathy. Among all the cases, one case of diabetic kidney disease was uniquely identified from death records. S3 Table shows the associations between individual lifestyle behaviors and all the outcomes. Being physically active, with lower WC, and moderate alcohol intake were associated with a lower risk of microvascular complications, while noncurrent smoking and healthy diet were not. The overall healthy lifestyle score was associated with lower risks of all the outcomes in a dose–response manner (all Ps for linear trend ≤0.01; Table 2 and Figs 1 and S2). Compared with participants with 0 to 1 low-risk lifestyle behavior, participants with 4 to 5 low-risk lifestyle behaviors had HRs (95% CIs) of 0.65 (0.46, 0.91) for diabetic retinopathy, 0.43 (0.30, 0.61) for diabetic kidney disease, 0.46 (0.29, 0.74) for diabetic neuropathy, and 0.54 (0.43, 0.68) for the composite microvascular complications, respectively. For each number increment in low-risk lifestyle behavior, there was a 13% lower risk of diabetic retinopathy (HR, 0.87; 95% CI: 0.80, 0.95), 22% lower risk of diabetic kidney disease (HR, 0.78; 95% CI: 0.72, 0.85), 27% lower risk of diabetic neuropathy (HR, 0.73; 95% CI: 0.65, 0.83), and a 18% lower risk of the composite microvascular complications (HR, 0.82; 95% CI: 0.77, 0.87).
Fig 1. Dose–response relationship of the healthy lifestyle score with risk of microvascular complications among individuals with T2D.
X-axis showed the numbers of low-risk lifestyle behaviors, and y-axis showed the HRs of the composite microvascular complications (A), diabetic retinopathy (B), diabetic kidney disease (C), and diabetic neuropathy (D). Black curves were HRs, and grey zones were 95% CIs. Multivariable-adjusted models were adjusted for age (continuous, years), sex (male, female), ethnicity (White, others), education attainment (college or university degree, A/AS levels or equivalent or O levels/GCSEs or equivalent or other professional qualifications, or none of the above), Townsend Deprivation Index (continuous), sleep duration (<6, 6–8, or ≥9 hours/day), family history of CVD (yes, no), family history of hypertension (yes, no), prevalence of hypertension (yes, no), diabetes duration (continuous, years), HbA1c (continuous, mmol/mol), use of diabetes medication (none, only oral medication pills, or insulin or others), use of antihypertensive medication (yes, no), use of lipid-lowing medication (yes, no), and use of aspirin (yes, no). All P-nonlinearity were ≥0.09 and all P for overall association were <0.001 (except for diabetic retinopathy: P for overall association = 0.008). CI, confidence interval; CVD, cardiovascular disease; HR, hazard ratio; T2D, type 2 diabetes.
In addition, the estimated PAFs of nonadherence to 4 or more low-risk lifestyle factors were 39.0% (17.7%, 56.8%) for diabetic kidney disease and 25.3% (10.0%, 39.4%) for the composite microvascular complications (Table 2).
All the biomarkers were significantly associated with the overall lifestyle score except for total protein, lipoprotein A, and SBP (S4 Table). The associations between the selected biomarkers and all outcomes are shown in S5 Table. Six significant mediators were detected on the associations of lifestyle score with risk of the composite microvascular complications and diabetic kidney disease, namely, albumin, HDL-C, triglycerides, apolipoprotein A, CRP, and HbA1c. The relationship between the lifestyle behaviors and risk of diabetic neuropathy was mediated by cystatin C, GGT, total bilirubin, albumin, HDL-C, triglycerides, apolipoprotein A, CRP, and HbA1c with the proportion of mediation effect ranging from 3.22% to 11.35%. Collectively, the mediators explained 23.20%, 24.40%, and 31.90% of the associations of overall lifestyle behaviors with composite microvascular complications, diabetic kidney disease, and diabetic neuropathy, respectively. In addition, our data showed that among all the potential biomarkers, only HbA1c was a significant mediator that explained 15.26% of the relationship between the overall lifestyle score and risk of diabetic retinopathy (Table 3).
Secondary analysis and sensitivity analysis
Consistent results were observed when analyses were stratified by age, sex, education, diabetes duration, use of hypoglycemic medication, and HbA1c level. No significant interaction was observed between the healthy lifestyle score and the stratified factors on the outcomes considering multiple comparisons (S3 Fig). Further, the results of different combinations of low-risk lifestyle factors showed that the increased numbers of low-risk lifestyle factors were associated with graded lower risks of diabetic retinopathy, diabetic kidney disease, diabetic neuropathy, and the composite microvascular complications (S6 Table).
In the sensitivity analyses, the results were generally robust when excluding patients with events that occurred within the first 2 years of follow-up, defining low-risk alcohol intake as moderate drinking and nondrinking, generating the lifestyle score using BMI or waist-to-hip ratio instead of WC, or generating the overall lifestyle score as a weighted score (S7–S10 Tables). The association between overall lifestyle behaviors and risk of diabetic kidney disease was slightly attenuated when estimated glomerular filtration rate (eGFR) was additionally adjusted, and the results of mediation analysis for diabetic kidney disease were largely unchanged with the additional adjustment of eGFR (S11 and S12 Tables). Further, we observed similar results when patients with preexisting CVD were included and in patients with preexisting CVD, although diabetic retinopathy did not reach statistical significance in patients with preexisting CVD probably due to the insufficient power (S13 and S14 Tables). Finally, consistent results were demonstrated when we used 2 competing risk models accounting for the death (S15 Table).
In this retrospective cohort study of patients with T2D, adherence to a greater number of healthy lifestyle behaviors, including recommended WC, noncurrent smoking, physically active, healthy diet, and moderate alcohol drinking, was inversely associated with lower risks of diabetic retinopathy, diabetic kidney disease, diabetic neuropathy, and the composite microvascular complications. For each number increment in low-risk lifestyle behavior, there was an 18% lower risk of developing diabetic microvascular complications. Moreover, the results of PAFs suggested that 25.3% of the diabetic microvascular complications could have been avoided if the patients with T2D had 4 or more healthy lifestyle behaviors. In addition, the mediators collectively explained 23.20% of the associations between the overall healthy lifestyle score and diabetic microvascular complications. Specifically, CRP, albumin, HbA1c, and lipids profile (HDL-C, triglycerides, and apolipoprotein A) could explain 4.44% to 10.69% of the association between overall lifestyle behaviors and the total diabetic microvascular complications.
Our study contributes to the literature regarding the influence of combined healthy lifestyle behaviors on the risk of diabetic microvascular complications. To date, many studies have been performed to evaluate the relationship between individual lifestyle behaviors and risk of diabetic microvascular complications; however, the joint association of multiple lifestyle behaviors with microvascular complications remains unknown. For example, the Irish Longitudinal Study showed that a history of smoking was associated with a higher risk of developing microvascular complications . The Ongoing Telmisartan Alone and in Combination with Ramipril Global Endpoint Trial (ONTARGET) studies demonstrated that adherence to a healthy dietary pattern (the Alternate Healthy Eating Index) , being physically active, and moderate alcohol consumption  were associated with a lower risk of incident chronic kidney disease among patients with T2D. Furthermore, general obesity and abdominal obesity were associated with higher risks of diabetic kidney disease , diabetic retinopathy , and diabetic neuropathy .
However, the results of lifestyle interventions on microvascular complications among patients with diabetes or impaired glucose tolerance in clinical trials were inconsistent. The Steno-2 randomized trial including 160 patients with T2D and persistent microalbuminuria showed pharmacological therapies in combination with lifestyle behavior modifications, including adopting a healthy diet, engaging regular physical activity, and participating in smoking cessation courses, significantly reduced the risk of diabetic nephropathy, retinopathy, and neuropathy . Further, the China Da Qing Diabetes Prevention Study including 577 participants with impaired glucose tolerance reported that healthy diet and exercise interventions in combination resulted in a 47% reduction in the diabetic retinopathy incidence, but no beneficial effects were observed for diabetic nephropathy or neuropathy . In addition, the Look AHEAD trial consisting of 5,145 overweight or obese patients with T2D, which focused on weight management through increased energy deficit and physical activity, resulted in a significant decrease in chronic kidney disease , but not diabetic neuropathy measured by physical examinations . Notably, microvascular complications were not predefined primary outcomes in these trials and small numbers of cases might partially explained the heterogeneities in these findings (e.g., 296 cases of very-high-risk chronic kidney disease in the Look AHEAD trial). Further trials with proper designs are needed to corroborate our findings in the future.
Our mediation analyses contribute to better understanding the lower risk of microvascular complications associated with lifestyle behaviors. Our data showed that the associations of overall lifestyle behaviors with diabetic kidney disease, diabetic neuropathy, and total microvascular complications may be explained by the improvement in glycemic control, liver function, lipid profile, and systemic inflammation, with lifestyle behaviors related lower risk of diabetic neuropathy might be additionally explained by kidney function amelioration. However, our data showed that the association between lifestyle and diabetic retinopathy was mainly through the glycemic control rather than other pathways. Our results corroborate prior findings from the observational studies. For example, intensive lifestyle intervention including physical activity and healthy diet recommendations could benefit glycemic control . Adherence to a combined healthy lifestyle score including healthy diet, physically active, nonsmoking, healthy sleep, and social support were associated with lower concentrations of inflammatory markers . Chronic Renal Insufficiency Cohort (CRIC) Study showed that combined healthy lifestyle characterized as physically active, nonsmoking, and BMI ≥25 kg/m2 were associated with lower risks of atherosclerotic events and kidney function decline among patients with chronic kidney disease . Furthermore, lifestyle modifications including promoting healthy diet, physical activity, and weight loss could significantly improve liver function, renal function, lipid profile, endothelial dysfunction, and reduce systemic inflammation in interventional studies [49–54].
The current study is among the first to investigate the relationship between the overall lifestyle behaviors and diabetic microvascular complications. The strengths of this study included the large sample size, long period of follow-up, and extensive collection of data on clinical biomarkers, which allowed us to comprehensively evaluate the potential mechanisms underlying the observed associations. Despite the strengths, this study should be interpreted in the light of its potential limitations. First, as the microvascular complications were identified via hospital inpatient records and death registries, there might be underreporting of the cases, for example, primary care data were not completely available currently. Second, the self-reported and one-time assessment of lifestyle behaviors data are susceptible to measurement errors. In addition, information on lifestyle behaviors was collected at recruitment and the behaviors may change over time; hence, the observed associations might be attenuated due to nondifferential misclassification bias. Third, mediation analysis assumes causality between lifestyles behaviors and biological biomarkers, although both the lifestyle behaviors and biological mediators were assessed at the same time in the UK Biobank. Future studies with repeatedly measured data are required to replicate our findings. Fourth, our study is limited in terms of ethnic diversity (>85% Whites); our results may not be directly generalized to other ethnic groups. Fifth, our study was based on a retrospective sampling from the UK Biobank study; hence, the causality should be interpreted with caution. Sixth, the UK Biobank is not representative of the general population of the UK, particularly relating to socioeconomic deprivation, lifestyles, and noncommunicable disease, with evidence of the healthy volunteer selection bias. Finally, residual or unknown confounding could not be excluded due to the observational study design, although we have in our effort to adjust for the potential confounding factors.
Our findings suggested that adherence to overall healthy lifestyle behaviors, including abstaining from smoking, abdominal obesity management, adopting a healthy diet, engaging in regular physical activity, and moderate alcohol drinking, was associated with a lower risk of diabetic microvascular complications. Our findings support the importance of public health programs and interventions targeting at improving health behaviors in combination to ameliorate the risk of diabetic microvascular complications. In addition, our data showed significant mediation effect of biomarkers involving of glycemic control, liver function, lipid profile, and systemic inflammation on the associations between lifestyle behaviors and microvascular complications. Further studies are warranted to confirm our findings.
S2 Fig. Dose–response relationship of the healthy lifestyle score with risk of microvascular complications among individuals with T2D (unadjusted model).
X-axis showed the numbers of low-risk lifestyle behaviors, and y-axis showed the HRs of the composite microvascular complications (A), diabetic retinopathy (B), diabetic kidney disease (C), and diabetic neuropathy (D). Black curves were HRs, and grey zones were 95% CIs. All P-nonlinearity were ≥0.10, and all P for overall association were <0.001 (except for diabetic retinopathy: P for overall association = 0.005). CI, confidence interval; HR, hazard ratio; T2D, type 2 diabetes.
S3 Fig. Stratified analyses of the associations of overall healthy lifestyle score with microvascular complications.
Model 1: unadjusted model. Model 2: age (continuous, years), sex (male, female), race (White, others), education attainment (college or university degree, A/AS levels or equivalent or O levels/GCSEs or equivalent or other professional qualifications, or none of the above), Townsend Deprivation Index (continuous), sleep duration (<6, 6–8, or ≥9 hours/day), family history of CVD (yes, no), family history of hypertension (yes, no), prevalence of hypertension (yes, no), diabetes duration (continuous, years), use of diabetes medication (none, only oral medication pills, or insulin or others), HbA1c (continuous, mmol/mol), use of antihypertensive medication (yes, no), use of lipid-lowing medication (yes, no), and use of aspirin (yes, no), with exception of stratifying factors. Interactions were tested based on Model 2.
S7 Table. HRs (95% CIs) of microvascular complications according to the overall lifestyle score among individuals with T2D after excluding the cases that occurred within first 2 years of follow-up.
CI, confidence interval; HR, hazard ratio; T2D, type 2 diabetes.
S8 Table. HRs (95% CIs) of microvascular complications according to the healthy lifestyle score using moderate drinking and nondrinking as the low-risk behavior in individuals with T2D.
CI, confidence interval; HR, hazard ratio; T2D, type 2 diabetes.
S9 Table. HRs (95% CIs) of microvascular complications according to the healthy lifestyle score using optimal BMI or optimal WHR as the low-risk behavior in individuals with T2D.
BMI, body mass index; CI, confidence interval; HR, hazard ratio; T2D, type 2 diabetes; WHR, waist-to-hip ratio.
S11 Table. HRs (95% CIs) of diabetic kidney disease according to the healthy lifestyle score among individuals with T2D with additional adjustment for eGFR.
CI, confidence interval; eGFR, estimated glomerular filtration rate; HR, hazard ratio; T2D, type 2 diabetes.
S12 Table. Possible mediators of the associations between the healthy lifestyle score and diabetic kidney disease among individuals with T2D with additional adjustment for eGFR.
eGFR, estimated glomerular filtration rate; T2D, type 2 diabetes.
S13 Table. HRs (95% CIs) of microvascular complications according to the healthy lifestyle score in 18,501 individuals with T2D including patients with preexisting macrovascular disease.
CI, confidence interval; HR, hazard ratio; T2D, type 2 diabetes.
We are grateful to all the participants of UK Biobank and all the people involved in building the UK Biobank study.
IDF. International Diabetes Federation Diabetes Atlas. 10th edition. Brussels; 2021.
Chatterjee S, Khunti K, Davies MJ. Type 2 diabetes. Lancet. 2017;389(10085):2239–51. Epub 2017/02/14. pmid:28190580.
DeFronzo RA, Ferrannini E, Groop L, Henry RR, Herman WH, Holst JJ, et al. Type 2 diabetes mellitus. Nat Rev Dis Primers. 2015;1:15019. Epub 2015/01/01. pmid:27189025.
Zhang X, Saaddine JB, Chou CF, Cotch MF, Cheng YJ, Geiss LS, et al. Prevalence of diabetic retinopathy in the United States, 2005–2008. JAMA. 2010;304(6):649–56. Epub 2010/08/12. pmid:20699456; PubMed Central PMCID: PMC2945293.
Parving HH, Lewis JB, Ravid M, Remuzzi G, Hunsicker LG. Prevalence and risk factors for microalbuminuria in a referred cohort of type II diabetic patients: a global perspective. Kidney Int. 2006;69(11):2057–63. Epub 2006/04/14. pmid:16612330.
Schmader KE. Epidemiology and impact on quality of life of postherpetic neuralgia and painful diabetic neuropathy. Clin J Pain. 2002;18(6):350–4. Epub 2002/11/21. pmid:12441828.
5. Lifestyle Management: Standards of Medical Care in Diabetes-2019. Diabetes Care. 2019;42(Suppl 1):S46–s60. Epub 2018/12/19. pmid:30559231.
Tracey ML, McHugh SM, Fitzgerald AP, Buckley CM, Canavan RJ, Kearney PM. Risk Factors for Macro- and Microvascular Complications among Older Adults with Diagnosed Type 2 Diabetes: Findings from The Irish Longitudinal Study on Ageing. J Diabetes Res. 2016;2016:5975903. Epub 2016/06/14. pmid:27294152; PubMed Central PMCID: PMC4884580.
Dunkler D, Kohl M, Teo KK, Heinze G, Dehghan M, Clase CM, et al. Dietary risk factors for incidence or progression of chronic kidney disease in individuals with type 2 diabetes in the European Union. Nephrol Dial Transplant. 2015;30 Suppl 4:iv76–85. Epub 2015/07/26. pmid:26209742.
Blomster JI, Chow CK, Zoungas S, Woodward M, Patel A, Poulter NR, et al. The influence of physical activity on vascular complications and mortality in patients with type 2 diabetes mellitus. Diabetes Obes Metab. 2013;15(11):1008–12. Epub 2013/05/17. pmid:23675676.
Blomster JI, Zoungas S, Chalmers J, Li Q, Chow CK, Woodward M, et al. The relationship between alcohol consumption and vascular complications and mortality in individuals with type 2 diabetes. Diabetes Care. 2014;37(5):1353–9. Epub 2014/03/01. pmid:24578358.
Dunkler D, Kohl M, Heinze G, Teo KK, Rosengren A, Pogue J, et al. Modifiable lifestyle and social factors affect chronic kidney disease in high-risk individuals with type 2 diabetes mellitus. Kidney Int. 2015;87(4):784–91. Epub 2014/12/11. pmid:25493953.
Man RE, Sabanayagam C, Chiang PP, Li LJ, Noonan JE, Wang JJ, et al. Differential Association of Generalized and Abdominal Obesity With Diabetic Retinopathy in Asian Patients With Type 2 Diabetes. JAMA Ophthalmol. 2016;134(3):251–7. Epub 2016/01/01. pmid:26720805.
Geng TT, Jafar TH, Yuan JM, Koh WP. The impact of diabetes on the association between alcohol intake and the risk of end-stage kidney disease in the Singapore Chinese Health Study. J Diabetes. 2020;12(8):583–93. Epub 2020/03/07. pmid:32142209.
Nivukoski U, Niemelä M, Bloigu A, Bloigu R, Aalto M, Laatikainen T, et al. Impacts of unfavourable lifestyle factors on biomarkers of liver function, inflammation and lipid status. PLoS ONE. 2019;14(6):e0218463. Epub 2019/06/21. pmid:31220128; PubMed Central PMCID: PMC6586311.
Shi N, Aroke D, Jin Q, Lee DH, Hussan H, Zhang X, et al. Proinflammatory and Hyperinsulinemic Dietary Patterns Are Associated With Specific Profiles of Biomarkers Predictive of Chronic Inflammation, Glucose-Insulin Dysregulation, and Dyslipidemia in Postmenopausal Women. Front Nutr. 2021;8:690428. Epub 2021/10/08. pmid:34616762; PubMed Central PMCID: PMC8488136.
Nivukoski U, Niemelä M, Bloigu A, Bloigu R, Aalto M, Laatikainen T, et al. Combined effects of lifestyle risk factors on fatty liver index. BMC Gastroenterol. 2020;20(1):109. Epub 2020/04/16. pmid:32293287; PubMed Central PMCID: PMC7157978.
Paramastri R, Hsu CY, Chuang YK, Lee HA, Wiratama BS, Chao JC. Synergistic Interaction of Dietary Pattern and Concordance Lifestyle with Abnormal Liver Function among Young Adults in Taiwan: A Population-Based Longitudinal Study. Nutrients. 2021;13(10). Epub 2021/10/24. pmid:34684598; PubMed Central PMCID: PMC8539530.
Geng TT, Talaei M, Jafar TH, Yuan JM, Koh WP. Pulse Pressure and the Risk of End-Stage Renal Disease Among Chinese Adults in Singapore: The Singapore Chinese Health Study. J Am Heart Assoc. 2019;8(23):e013282. Epub 2019/11/27. pmid:31766974; PubMed Central PMCID: PMC6912960.
Ricardo AC, Anderson CA, Yang W, Zhang X, Fischer MJ, Dember LM, et al. Healthy lifestyle and risk of kidney disease progression, atherosclerotic events, and death in CKD: findings from the Chronic Renal Insufficiency Cohort (CRIC) Study. Am J Kidney Dis. 2015;65(3):412–24. Epub 2014/12/03. pmid:25458663; PubMed Central PMCID: PMC4339665.
Geng TT, Jafar TH, Neelakantan N, Yuan JM, van Dam RM, Koh WP. Healthful dietary patterns and risk of end-stage kidney disease: the Singapore Chinese Health Study. Am J Clin Nutr. 2021;113(3):675–83. Epub 2021/01/01. pmid:33381807; PubMed Central PMCID: PMC7948892.
Niu M, Zhang L, Wang Y, Tu R, Liu X, Wang C, et al. Lifestyle Score and Genetic Factors With Hypertension and Blood Pressure Among Adults in Rural China. Front Public Health. 2021;9:687174. Epub 2021/09/07. pmid:34485217; PubMed Central PMCID: PMC8416040.
Tu Z-Z, Lu Q, Zhang Y-B, Shu Z, Lai Y-W, Ma M-N, et al. Associations of Combined Healthy Lifestyle Factors with Risks of Diabetes, Cardiovascular Disease, Cancer, and Mortality Among Adults with Prediabetes: Four Prospective Cohort Studies in China, the United Kingdom, and the United States. Engineering. 2022.
Sudlow C, Gallacher J, Allen N, Beral V, Burton P, Danesh J, et al. UK biobank: an open access resource for identifying the causes of a wide range of complex diseases of middle and old age. PLoS Med. 2015;12(3):e1001779. Epub 2015/04/01. pmid:25826379; PubMed Central PMCID: PMC4380465.
Caleyachetty R, Littlejohns T, Lacey B, Bešević J, Conroy M, Collins R, et al. United Kingdom Biobank (UK Biobank): JACC Focus Seminar 6/8. J Am Coll Cardiol. 2021;78(1):56–65. Epub 2021/07/03. pmid:34210415.
Vetter C, Dashti HS, Lane JM, Anderson SG, Schernhammer ES, Rutter MK, et al. Night Shift Work, Genetic Risk, and Type 2 Diabetes in the UK Biobank. Diabetes Care. 2018;41(4):762–9. Epub 2018/02/15. pmid:29440150; PubMed Central PMCID: PMC5860836.
Carnethon MR, De Chavez PJ, Biggs ML, Lewis CE, Pankow JS, Bertoni AG, et al. Association of weight status with mortality in adults with incident diabetes. JAMA. 2012;308(6):581–90. Epub 2012/08/09. pmid:22871870; PubMed Central PMCID: PMC3467944.
Tobias DK, Pan A, Jackson CL, O’Reilly EJ, Ding EL, Willett WC, et al. Body-mass index and mortality among adults with incident type 2 diabetes. N Engl J Med. 2014;370(3):233–44. Epub 2014/01/17. pmid:24428469; PubMed Central PMCID: PMC3966911.
Dallongeville J, Bhatt DL, Steg PH, Ravaud P, Wilson PW, Eagle KA, et al. Relation between body mass index, waist circumference, and cardiovascular outcomes in 19,579 diabetic patients with established vascular disease: the REACH Registry. Eur J Prev Cardiol. 2012;19(2):241–9. Epub 2011/04/01. pmid:21450609.
WHO. Waist circumference and waist-hip ratio: report of a WHO expert consultation. Geneva; 2008.
IDF. Consensus Worldwide Definition of the Metabolic Syndrome. Brussel; 2020.
Chudasama YV, Khunti K, Gillies CL, Dhalwani NN, Davies MJ, Yates T, et al. Healthy lifestyle and life expectancy in people with multimorbidity in the UK Biobank: A longitudinal cohort study. PLoS Med. 2020;17(9):e1003332. Epub 2020/09/23. pmid:32960883; PubMed Central PMCID: PMC7508366.
Chudasama YV, Khunti KK, Zaccardi F, Rowlands AV, Yates T, Gillies CL, et al. Physical activity, multimorbidity, and life expectancy: a UK Biobank longitudinal study. BMC Med. 2019;17(1):108. Epub 2019/06/13. pmid:31186007; PubMed Central PMCID: PMC6560907.
Ainsworth BE, Haskell WL, Whitt MC, Irwin ML, Swartz AM, Strath SJ, et al. Compendium of physical activities: an update of activity codes and MET intensities. Med Sci Sports Exerc. 2000;32(9 Suppl):S498–504. Epub 2000/09/19. pmid:10993420.
Said MA, Verweij N, van der Harst P. Associations of Combined Genetic and Lifestyle Risks With Incident Cardiovascular Disease and Diabetes in the UK Biobank Study. JAMA Cardiol. 2018;3(8):693–702. Epub 2018/06/30. pmid:29955826; PubMed Central PMCID: PMC6143077.
Biobank UK. UK Biobank Biomarker Project—Companion Document to Accompany Serum Biomarker Data. UK Biobank Showcase. 2019;1(2019):2019.
Bonaccio M, Costanzo S, Di Castelnuovo A, Persichillo M, Magnacca S, De Curtis A, et al. Ultra-processed food intake and all-cause and cause-specific mortality in individuals with cardiovascular disease: the Moli-sani Study. Eur Heart J. 2022;43(3):213–24. Epub 2021/12/02. pmid:34849691.
Figueras-Roca M, Molins B, Sala-Puigdollers A, Matas J, Vinagre I, Ríos J, et al. Peripheral blood metabolic and inflammatory factors as biomarkers to ocular findings in diabetic macular edema. PLoS ONE. 2017;12(3):e0173865. Epub 2017/03/23. pmid:28328965; PubMed Central PMCID: PMC5362077.
Saito H, Tanabe H, Kudo A, Machii N, Higa M, Yamaguchi S, et al. High FIB4 index is an independent risk factor of diabetic kidney disease in type 2 diabetes. Sci Rep. 2021;11(1):11753. Epub 2021/06/05. pmid:34083571; PubMed Central PMCID: PMC8175689.
Honigberg MC, Zekavat SM, Pirruccello JP, Natarajan P, Vaduganathan M. Cardiovascular and Kidney Outcomes Across the Glycemic Spectrum: Insights From the UK Biobank. J Am Coll Cardiol. 2021;78(5):453–64. Epub 2021/05/21. pmid:34015477; PubMed Central PMCID: PMC8324525.
Man REK, Gan ATL, Fenwick EK, Gupta P, Wong MYZ, Wong TY, et al. The Relationship between Generalized and Abdominal Obesity with Diabetic Kidney Disease in Type 2 Diabetes: A Multiethnic Asian Study and Meta-Analysis. Nutrients. 2018;10(11). Epub 2018/11/08. pmid:30400648; PubMed Central PMCID: PMC6266073.
Andersen ST, Witte DR, Dalsgaard EM, Andersen H, Nawroth P, Fleming T, et al. Risk Factors for Incident Diabetic Polyneuropathy in a Cohort With Screen-Detected Type 2 Diabetes Followed for 13 Years: ADDITION-Denmark. Diabetes Care. 2018;41(5):1068–75. Epub 2018/03/01. pmid:29487078.
Gaede P, Vedel P, Larsen N, Jensen GV, Parving HH, Pedersen O. Multifactorial intervention and cardiovascular disease in patients with type 2 diabetes. N Engl J Med. 2003;348(5):383–93. Epub 2003/01/31. pmid:12556541.
Gong Q, Gregg EW, Wang J, An Y, Zhang P, Yang W, et al. Long-term effects of a randomised trial of a 6-year lifestyle intervention in impaired glucose tolerance on diabetes-related microvascular complications: the China Da Qing Diabetes Prevention Outcome Study. Diabetologia. 2011;54(2):300–7. Epub 2010/11/04. pmid:21046360.
Effect of a long-term behavioural weight loss intervention on nephropathy in overweight or obese adults with type 2 diabetes: a secondary analysis of the Look AHEAD randomised clinical trial. Lancet Diabetes Endocrinol. 2014;2(10):801–9. Epub 2014/08/16. pmid:25127483; PubMed Central PMCID: PMC4443484.
Effects of a long-term lifestyle modification programme on peripheral neuropathy in overweight or obese adults with type 2 diabetes: the Look AHEAD study. Diabetologia. 2017;60(6):980–8. Epub 2017/03/30. pmid:28349174; PubMed Central PMCID: PMC5423967.
Johansen MY, MacDonald CS, Hansen KB, Karstoft K, Christensen R, Pedersen M, et al. Effect of an Intensive Lifestyle Intervention on Glycemic Control in Patients With Type 2 Diabetes: A Randomized Clinical Trial. JAMA. 2017;318(7):637–46. Epub 2017/08/16. pmid:28810024.
Sotos-Prieto M, Bhupathiraju SN, Falcon LM, Gao X, Tucker KL, Mattei J. Association between a Healthy Lifestyle Score and inflammatory markers among Puerto Rican adults. Nutr Metab Cardiovasc Dis. 2016;26(3):178–84. Epub 2016/02/04. pmid:26838054; PubMed Central PMCID: PMC4788524.
Goodpaster BH, Delany JP, Otto AD, Kuller L, Vockley J, South-Paul JE, et al. Effects of diet and physical activity interventions on weight loss and cardiometabolic risk factors in severely obese adults: a randomized trial. JAMA. 2010;304(16):1795–802. Epub 2010/10/12. pmid:20935337; PubMed Central PMCID: PMC3082279.
Estruch R, Martínez-González MA, Corella D, Salas-Salvadó J, Ruiz-Gutiérrez V, Covas MI, et al. Effects of a Mediterranean-style diet on cardiovascular risk factors: a randomized trial. Ann Intern Med. 2006;145(1):1–11. Epub 2006/07/05. pmid:16818923.
Esposito K, Marfella R, Ciotola M, Di Palo C, Giugliano F, Giugliano G, et al. Effect of a mediterranean-style diet on endothelial dysfunction and markers of vascular inflammation in the metabolic syndrome: a randomized trial. JAMA. 2004;292(12):1440–6. Epub 2004/09/24. pmid:15383514.
Schrauben SJ, Hsu JY, Amaral S, Anderson AH, Feldman HI, Dember LM. Effect of Kidney Function on Relationships between Lifestyle Behaviors and Mortality or Cardiovascular Outcomes: A Pooled Cohort Analysis. J Am Soc Nephrol. 2021;32(3):663–75. Epub 2021/02/07. pmid:33547215; PubMed Central PMCID: PMC7920187.
King CC, Piper ME, Gepner AD, Fiore MC, Baker TB, Stein JH. Longitudinal Impact of Smoking and Smoking Cessation on Inflammatory Markers of Cardiovascular Disease Risk. Arterioscler Thromb Vasc Biol. 2017;37(2):374–9. Epub 2016/12/10. pmid:27932354; PubMed Central PMCID: PMC5269476.
Franco I, Bianco A, Mirizzi A, Campanella A, Bonfiglio C, Sorino P, et al. Physical Activity and Low Glycemic Index Mediterranean Diet: Main and Modification Effects on NAFLD Score. Results from a Randomized Clinical Trial. Nutrients. 2020;13(1). Epub 2021/01/01. pmid:33379253; PubMed Central PMCID: PMC7823843.