Vegetarian diet, change in dietary patterns, and diabetes risk: a prospective study

Background/objectives Vegetarian diets are inversely associated with diabetes in Westerners but their impact on Asians—whose pathophysiology differ from Westerners—is unknown. We aim to investigate the association between a vegetarian diet, change in dietary patterns and diabetes risk in a Taiwanese Buddhist population. Methods We prospectively followed 2918 non-smoking, non-alcohol drinking Buddhists free of diabetes, cancer, and cardiovascular diseases at baseline, for a median of 5 years, with 183 incident diabetes cases confirmed. Diet was assessed through a validated food frequency questionnaire at baseline and a simple questionnaire during follow-ups. Incident cases of diabetes were ascertained through follow-up questionnaires, fasting glucose and HbA1C. Stratified Cox Proportional Hazards Regression was used to assess the effect of diets on risk of diabetes. Results Consistent vegetarian diet was associated with 35% lower hazards (HR: 0.65, 95% CI: 0.46, 0.92), while converting from a nonvegetarian to a vegetarian pattern was associated with 53% lower hazards (HR: 0.47, 95% CI: 0.30, 0.71) for diabetes, comparing with nonvegetarians while adjusting for age, gender, education, physical activity, family history of diabetes, follow-up methods, use of lipid-lowering medications, and baseline BMI. Conclusion Vegetarian diet and converting to vegetarian diet may protect against diabetes independent of BMI among Taiwanese.


Introduction
The rapid growth of diabetes creates tremendous health and economic burdens worldwide 1 . In Taiwan, diabetes patients incur 2.8 times more medical expense than matched non-diabetes individuals, and used up 29% of total healthcare expenditure 2 , while high blood glucose accounts for the highest adult mortality among major metabolic, lifestyle, and environmental factors 3 . Preventive strategies are desperately needed, particularly in Asia, where the largest number of diabetes cases is expected 1 , and diabetes tend to occur despite lower BMI 4 .
Diabetes is driven by impaired insulin action (insulin resistance) and insulin secretion (islet dysfunction and incretin failure), co-influenced by genetics and environment 5 . While diabetes in Caucasian is highly attributed to obesity and insulin resistance, emerging evidence suggests that β-cell dysfunction may be more predictive of diabetes in East Asians 6 . In fact, Japanese with normal glucose tolerance already showed an insulin secretion ability similar to Caucasian diabetes patients 7 . The potential pathophysiological difference warrants research on population-specific preventive strategies.
Vegetarian diets 8 and plant-based dietary patterns 9 have been shown to dramatically reduce diabetes risks in American populations, but their effect in Asians is less clear. Unlike in Westerns, the effect of meat and fish on diabetes is inconsistent in Asians [10][11][12][13] , and may be modified by BMI 12,14 . Whether complete avoidance of meat and fish is necessary to boost protection for diabetes remains puzzling, as Asians tend to be leaner and consume less meat than Westerners. Our primary aim is to examine the association between a vegetarian diet, converting from a nonvegetarian to a vegetarian dietary pattern, and diabetes risk in a prospective cohort study. Our secondary aim is to explore potential interaction between dietary patterns and baseline metabolic status (BMI, metabolic syndrome, and other metabolic characteristics).

Study design and population
The Tzu Chi Health Study (TCHS) is a prospective cohort study that recruited 4625 Tzu Chi volunteers (age 18 to 87 years) from 2007 to 2009 at the Buddhist Dalin Tzu Chi Hospital. Tzu Chi volunteers are devoted Buddhists who volunteer regularly for the Buddhist Tzu Chi Foundation for a variety of charity works and disaster relief. Many nonvegetarian volunteers converted to vegetarian in the year 2011 due to a religious event.
Because the nonvegetarians consumed only small amount of meat (median intake per day, women: 10 g, men: 22 g) and fish (median intake per day, women: 6 g, men: 15 g) 15 , this cohort is equipped with a unique opportunity to distinguish disease risk between true vegetarians and those with little to modest meat consumption.
Participants were followed from 2010 to 2012 (first follow-up), and from 2013 to mid-2016 (second follow-up). Every three years, a post-card was sent to invite each participant for a follow-up health examination during which participants also fill out a questionnaire on disease diagnosis. Participants who did not return for health examination by the end of 2015 were sent a follow-up questionnaire in May 2016 to assess their dietary practices and disease conditions (for each disease, choices includes: 'no', 'yes', 'not sure', and the time of diagnosis). If the questionnaire was not returned within a month, a research assistant called the participants to interview this questionnaire. The study was approved by the Institutional Review Board for ethics in the Buddhist Dalin Tzu Chi Hospital. All participants gave written informed consents. The study has been registered at ClinicalTrials.gov (ID: NCT: 03204552).

Assessment of metabolic characteristics and covariates
At baseline, trained research assistants interviewed all participants on basic demographics, medical history and relevant medications, family history of diseases, diet, use of cigarette and alcohol, and leisure time physical activities (LTPA). Fasting glucose and blood lipids, including triglyceride (TG), high density lipoprotein cholesterol (HDL-C) and other cholesterol were assessed by the Integra 800 system (Roche Diagnostics, Indianapolis, IN) at baseline and Dimension RxL Max (Siemens Healthcare Diagnostics Inc, Newark, NJ) at follow-ups. HbA1C was assessed by Variant Turbo (BIO-RAD, Hercules, CA).
Heights and weights were measured using an electronic scale while participants stood in an upright position without shoes. Waist circumference was measured at navel.
Blood pressures were assessed using VP1000 system (Colin, Komaki, Japan). Metabolic syndrome was defined by the National Cholesterol Education Program, Adult Treatment Panel (ATP III) but using the Taiwanese cut points for waist circumference (≥90 cm for men, ≥80 cm for women). Fatty liver was assessed using ultrasonography by hospital gastroenterologists.

Assessment of diet
At baseline, participants were interviewed on a 64-item quantitative food frequency questionnaire (FFQ) that had been validated within our cohort participants, and showed good validity and reliability when compared with repeated dietary records and biomarkers 16 . At follow-ups, all participants answered a simple questionnaire that asked whether they are vegetarians (choices including: not vegetarian, breakfast vegetarian, vegetarian on 1st and 15th day of each lunar month (a practice for many Buddhists), irregular dates of vegetarian diets, full time vegetarian), and the type of vegetarian diet (vegan, lacto-ovo vegetarian, lacto-vegetarian, ovo-vegetarian). Only full-time vegetarians were considered vegetarians in our analysis.
Dietary patterns are divided into 4 types: (1) vegetarians are defined as those who have been following vegetarian diets at baseline and all the follow-ups; (2) the reverted are those who were initially vegetarians but became nonvegetarians at one of the follow-ups; (3) the converted are those who were initially nonvegetarians but became vegetarians later; and (4) the nonvegetarians are those who had consistently reported eating nonvegetarian diet in baseline and follow-ups. Diets assessed after diabetes diagnosis were discarded from the analysis.

Case ascertainment
Diabetes cases were identified if participants reported diabetes diagnosis at follow-up questionnaires, or if their HbA1C ≥6.5%. Participants with only one fasting blood glucose ≥7.0 mmol/L but otherwise normal HbA1C were identified as possible diabetes cases. For these individuals, a physician coauthor (MNL) further reviewed their medical records to check if they have additional blood tests or prescription of diabetes medication to determine their diabetes status. Participants without further tests or available medical records were considered unconfirmed diabetes events (n = 25) and were excluded, but included as cases in a sensitivity analysis.

Statistical analysis
Baseline characteristics between different dietary patterns were compared using ANOVA test (for continuous variables) and Chi-square test (for categorical variables).
Cox regression was used to analyze the association between dietary patterns and risk of diabetes, with stratification for LTPA and follow-up methods (questionnaire only vs health examination) as the interaction term of time with these variables violated proportional hazard assumption. Model 1 adjusted for age, sex, education, family history of diabetes, LTPA, lipid-lowering medication, methods of follow-up. Model 2 additionally adjusted for BMI to estimate the protective association independent of BMI (a mediator). Time of disease occurrence was set to be the time that the first abnormal glucose was identified (HbA1c ≥6.5% or fasting blood glucose ≥7.0 mmol/L). For participants who reported diagnosis of diabetes at questionnaire but could not remember the time of diabetes diagnosis, censor time was set to be halfway between the previous known disease-free time point and the follow-up time in which diabetes was reported. For those who did not report having diabetes in the questionnaire, but were found to have diabetes during health examination, the date of health examination was used as the date of disease onset. Subgroup analyses by gender, baseline metabolic syndrome, fatty liver, BMI (<24 or ≥24 kg/m 2 ), fasting glucose (<5.56 or ≥5.56 mmol/L), TG (<1.69 or ≥1.69 mmol/L), and HDL-C (men: <1.03 or ≥1.03 mmol/L, women: <1.29 or ≥1.29 mmol/L) levels were performed. Test of interaction was performed by including a cross-product term in the regression model. Several sensitivity analyses were performed: (1) 25 unconfirmed diabetes event were treated as diabetes cases. (2) To ensure our result was not affected by detection bias from different follow-up methods (health examination vs questionnaire-only), we performed another sensitivity analysis in which only self-reported diabetes were counted as cases. (3) We adjusted for metabolic syndrome in addition to Model 2. (4) Age was used as the underlying time scale for Cox regression, with left truncation at study entry. All analysis were conducted using SAS Statistical Software (version 9.4, SAS Institute, Cary, NC).

Exclusion and loss to follow-up
To calculate the incidence of diabetes and avoid bias related to comorbidities, we excluded participants with baseline diabetes or fasting blood glucose ≥7.0 mmol/L (n = 322), reported history of cancer (n = 172), coronary heart disease (n = 194), stroke (n = 26). We also excluded those who have ever used cigarette (n = 691) or alcohol (n = 606), as cigarette smoking may modify the effect of meat on diabetes 17 and alcohol is frequently associated with smoking. This left a total of 3185 participants. We further excluded those lost to follow-up (n = 210), those who did not know about their diabetes status in follow-up questionnaire (n = 8), those missed the item on diabetes diagnosis in the questionnaire (n = 42, likely without diabetes as they responded to a few diseases they had a diagnosis and left all other diseases blank), and unconfirmed diabetes events (n = 25, some may be included as non-cases up until the time that normal blood glucose and HbA1C were observed), leaving a total of 2918 participants in the final analysis. A flow chart summarizing number of participants at each stage is available at Supplemental Figure. Of those who were followed through health examination, 1902 returned for first follow-up (2010 to 2012), 1739 returned for second follow-up (2013 to mid-2016), and 1247 returned for both. The baseline characteristics by follow-up status and methods were compared at Supplemental Table 1. Men were less likely to be followed (either through health examination or mailed questionnaire), while women with lower education were more likely to return for health examination. No significant difference in diet, lifestyle and metabolic risk factors were observed among those returned for health examination, those responded to the follow-up questionnaire, and those lost to follow-up. Table 1 shows the baseline characteristics of different dietary patterns. Nonvegetarians tend to have higher BMI, waist circumference (among female), and fasting blood glucose. Female were more likely to consume vegetarian diets at baseline or switch to a vegetarian diet later. The converted had the lowest proportion with metabolic syndrome and elevated TG.

Dietary patterns and diabetes risks
In the median 5.2 years of follow-up, 183 individuals developed incident diabetes. The association between dietary patterns and diabetes in all participants are shown in Table 2. Consistent vegetarians and the converted showed about 40-60% reduction in risk of diabetes, compared with nonvegetarians. Those reverting from vegetarian diet to nonvegetarian diet (Reverted) also experienced a lower risk of diabetes compared with nonvegetarians (as indicated by the magnitude of hazard ratio) but the association did not reach statistical significance due to small sample size (only 6 cases in 583 person-years of follow-up).
Additional subgroup analyses and the test for interactions are presented in Table 3. All the analysis in this table were adjusted for BMI (same as Model 2 in Table 2). The association between dietary patterns and diabetes did not differ significantly between subgroups by gender, metabolic syndrome, fasting glucose, and HDL at baseline, as indicated by the non-significant test of interactions and similar trends and magnitudes in hazard ratios. The test of interaction between dietary patterns and TG was significant (P-interaction = 0.039) and that the hazard ratio of diabetes in vegetarians and the converted (vs nonvegetarians) are significant mainly in those with normal TG but not those with elevated TG. Overall, converting to vegetarian is strongly associated with protection among those with healthier baseline metabolic characteristics (normal TG: HR = 0.32, without metabolic syndrome: HR = 0.40, and BMI < 24: HR = 0.24). Due to limited case numbers in the reverted, all the results in the subgroup analysis for the reverted (versus nonvegetarians) were statistically insignificant and some subgroup analyses were not possible (data not shown).

Sensitivity analyses
Our sensitivity analyses (Supplemental Table 2) showed similar findings as our main analysis: (1) When unconfirmed diabetes were counted as diabetes cases,

Discussion
In this prospective cohort study, both vegetarian diet and nonvegetarians switching to vegetarian diets are associated with~50% reduction in risk of diabetes. This trend did not differ across gender, metabolic syndrome, HDL-C, and fasting glucose status. The potential protective effect of converting to vegetarian appeared to be particularly pronounced in those with initially healthier metabolic status. To the best of our knowledge, this is the first study that (1) examined the prospective association between a vegetarian diet and diabetes in high risk Asians, and (2) addressed how switching from a nonvegetarian diet to a vegetarian diet may influences diabetes risk.

Plant-based diets
The magnitude of protective association between a vegetarian diet and diabetes in our study (HR before BMI adjustment = 0.54, HR after BMI adjustment = 0.65) is comparable to the Adventist Health Study -2 (AHS-2; lacto-ovo vegetarians: OR before BMI adjustment = 0.46, OR after BMI adjustment = 0.62) 8 , and the three Harvard cohorts of US nurses and health professionals (highest vs lowest deciles of healthy plant-based diet index: HR before BMI adjustment = 0.55, HR after BMI adjustment = 0.66) 9 . In all these studies, BMI appears to account for only a small portion of the protective association, as the results in all these studies are only slightly attenuated after BMI adjustment. The BMI-independent protective effect is also confirmed in a trial using the Mediterranean diet 18 -another healthy plant-based diet. These studies, together, suggest that vegetarian or healthy plant-based diets may exert protective effect beyond just BMI. All Models adjusted for age, gender, education, leisure time physical activities, family history of diabetes, follow-up methods (health examination or questionnaire only), lipid medication, and BMI (same as Model 2 in Table 2). Converted, diet changed from nonvegetarian to vegetarian; IFGimpaired fasting glucose, MSmetabolic syndrome, TGtriglyceride, HDL-Chigh-density lipoprotein cholesterol. Normal glucose: <5.56 mmol/L; IFG: ≥5.56 mmol/L; normal TG: <1.69 mmol/L; elevated TG: ≥1.69 mmol/L; normal HDL-C: ≥1.03 mmol/L for men, ≥1.29 mmol/L for women; low HDL-C: <1.03 mmol/L for men, <1.29 mmol/L for women. Data are hazard ratios (and 95% CI) estimated using Cox regression stratified by follow-up methods and leisure time physical activities. Model 1 adjusted for age, gender, education, leisure time physical activities, family history of diabetes, follow-up methods (health examination or questionnaire only), and lipid medication. Model 2 additionally adjusted for BMI Vegetarians in our cohort and in AHS-2 consumed more whole grains and vegetables than nonvegetarians 15,19 , and these may protect against diabetes through higher fiber and magnesium 20 , as deficiency in magnesium may impair insulin signaling 21 . In addition, soy is a major source of protein for Taiwanese vegetarians, and soy has been shown to improve insulin resistance when replacing meat in randomized controlled trials 22,23 . Soy and legume are inversely associated with risk of diabetes in a Chinese cohort 24 . A vegetables-fruits-soy dietary pattern is inversely associated with diabetes incidence in Singaporean Chinese 17 .

Meat and fish
Although the protective association is likely caused by various plant components, it may have resulted partly from the simultaneous elimination of meat, as components of meat and fish may affect both insulin secretion and insulin resistance. Meat is high in saturated fats, which have been shown to trigger human β-cell apoptosis 25 . Fatty acids from meat have been adversely associated with insulin secretion and Disposition Index (β-cell function accounting for insulin sensitivity) 26 . A trial showed that while plant polyphenol improves glucose metabolism, fish omega-3 fatty acids decreases insulin secretion and postprandial GLP-1 27 . Besides insulin secretion, high intakes of animal protein may induce insulin resistance through 3-Hydroxyinsbutyrate (a valine metabolite) and fibroblast growth factor 21 28 . The positive association between animal protein and type 2 diabetes has also been confirmed in a recent prospective cohort study and an updated meta-analysis 29 .
However, the association between meat and diabetes is equivocal among Asian women: meat was not associated with diabetes risk in Japanese women 13 and was associated with protection among normal weight Chinese women in Shanghai 12 . None of these studies actually included a diet range of complete meat avoidance, and it is possible that even the lowest quantile in these cohorts did not consume low enough meat (and high enough healthy plant foods) to attain protection. The inverse association between meat intake and diabetes in the Shanghai Women's Health Study may potentially be confounded by unmeasured social economic factors, early life food insecurity, and cohort effect, as the lowest meat consuming quantile appeared to be older and from lower socioeconomic classes 12 . Lower socioeconomic sectors may have experienced early life under nutrition, which could potentially trigger epigenetic changes that induce diabetes risk 30 .

Potential interactions and effect modifications
In our subgroup analysis by gender, the insignificant association between vegetarian diet and diabetes in the male subgroup was likely due to small sample size. The hazard ratio of 0.56 (clearly trends toward protection) and the insignificant test for interaction (P-interaction = 0.35) suggest that the protective association between a vegetarian diet and diabetes is likely similar for both genders. This surmise is consistent with the results reported for healthy plant-based diets in the three Harvard cohorts of US nurses and health professionals 9 .
In the Shanghai Women's Health Study 12 and in Japanese Americans within the Multiethnic Cohort 13 , meat or a meat-fat dietary pattern increased diabetes risk only among overweight individuals. In our subgroup analysis, however, the association between a vegetarian diet and diabetes did not differ significantly across categories of BMI or metabolic syndrome. It is possible that the effect of a vegetarian dietary pattern is not due solely to either the minimization of animal product or the higher functional plant ingredients, but the additive effect of both. While meat may be harmful among the overweight, plant components may be beneficial even among the normal weight, thus extending its effect across different weight status. This result is also found in the PREDIMED-Reus Trial, in which a Mediterranean diet-also a plant-based dietary pattern-showed similar protective trends for different categories of BMI (HR = 0.54 for BMI ≤ 30, HR = 0.56 for BMI > 30), and fasting glucose (HR = 0.53 for fasting glucose ≤ 6.1 mmol/L, HR = 0.32 for fasting glucose > 6.1 mmol/L) 18 .
The effect of lifestyle interventions tend to be more obvious for individuals at high risk than those at low risk 31 , and it should be no surprise that lifestyle trials aiming to prevent diabetes tend to focus on high risk individuals 18,32 . Interestingly, we found that the association between converting to vegetarian diet and diabetes risk to be quite strong for those with supposedly healthier metabolic status: normal TG (HR = 0.32), without metabolic syndrome (HR = 0.32), and BMI < 24 (HR = 0.24). However, our sample size in this subgroup is limited and residual confounding may still exists. Randomized controlled trials are needed to confirm our findings.

Change in diet
The strong inverse association between converting to vegetarian diet and diabetes suggests that diabetes risk or protection may be influenced by recent diets. In US nurses and health professionals, change in red meat consumption over 4 years has been associated with diabetes risk independent of baseline red meat intake 33 . Trials using vegetarian diets improve glycemic control in weeks 34 . Switching to a complete plant-based diet may increase butyrate production through gut microbiota alteration in as short as 1 day 35 . Butyrate may induce incretin secretion and improve glucose metabolism 36 . The magnitude of protective association of converting to vegetarian diet (53% reduction in diabetes risk compared with a persistent nonvegetarian diet) in our population is comparable to the Diabetes Prevention Trial (58% reduction in diabetes risk through lifestyle changes aiming at weight reduction), which is more effective than use of metformin (31% reduction) among high risk individuals 32 . These findings, together, suggest that a powerful diabetes preventive potential may be achieved by either a healthy plant-based diet (calorically non-restrictive) or a BMIlowering intervention aiming to reduce calories and increase physical activities, offering alternatives strategies for both lean (weight loss may lead to undesirable side effects associated with underweight) and obese individuals.

Study strength and limitations
The prospective design with high follow-up rate (93%) reduces recall and selection biases. The majority (75%) of participants have their diabetes status confirmed by blood tests, which minimize potential undiagnosed diabetes and reduce outcome misclassification. The homogenous population of non-smokers and non-alcohol drinkers from the same religious community may reduce unmeasured confounding and strengthen internal validity, though the generalizability to other population requires further confirmation. Although we attempted to adjust for potential confounders, some residual confounding effect may still remain: (1) Measurement of physical activities through questionnaire may not be accurate or adequate 37 .
(2) Conversion to vegetarian diet may be associated with religious devotion and it is unknown whether religious devotion would impact diabetes risk. In addition, the lack of detail in the follow-up questionnaire prevented us from analyzing detail dietary changes-other than meat and fish intake went from small to zero for the converted. Although imperfect, our verification of dietary pattern at follow-ups allowed us to capture dietary changes pertinent to our study aim (vegetarian vs nonvegetarian dietary patterns), providing more insights than most cohorts that rely solely on one baseline diet assessment. Finally, our sample size may be inadequate in some of the subgroups to draw conclusions on potential effect modification by metabolic status. Nevertheless, our finding shows no harm associated with a vegetarian diet in any subgroup and that vegetarian diets could be safely practiced.
Our study suggests that a vegetarian diet may immensely reduce the risk of diabetes in lean Asian populations. Further researches on how animal and plant components affect β-cell function and insulin resistance will be needed to establish the ideal diet for diabetes prevention. Our consistent finding with Western populations has a farreaching public health and environmental implication. The large and consistent protective association between plant-based diets and diabetes, and the over-consumption of meat by the majority today, suggest enormous population-attributable protection potential of a vegetarian diet. At the same time, shifting toward a plant-based diet is estimated to reduce food-related greenhouse gas emissions by 29-70% 38 . A vegetarian diet may be a stunning solution to the diet-environment-health trilemma that our globe urgently need to tackle, for the welfare, if not the survival, of many who are deeply threatened by climate change and diabetes.