White rice consumption and risk of type 2 diabetes: meta-analysis and systematic review
BMJ 2012; 344 doi: https://doi.org/10.1136/bmj.e1454 (Published 15 March 2012) Cite this as: BMJ 2012;344:e1454- Emily A Hu, research assistant1,
- An Pan, research fellow1,
- Vasanti Malik, research fellow1,
- Qi Sun, instructor in medicine12
- 1Department of Nutrition, Harvard School of Public Health, 655 Huntington Avenue, Boston, MA 02115, USA
- 2Channing Laboratory, Department of Medicine, Brigham and Women’s Hospital and Harvard Medical School, 181 Longwood Avenue, Boston
- Correspondence to: Qi Sun qisun{at}hsph.harvard.edu
- Accepted 19 January 2012
Abstract
Objectives To summarise evidence on the association between white rice consumption and risk of type 2 diabetes and to quantify the potential dose-response relation.
Design Meta-analysis of prospective cohort studies.
Data sources Searches of Medline and Embase databases for articles published up to January 2012 using keywords that included both rice intake and diabetes; further searches of references of included original studies.
Study selection Included studies were prospective cohort studies that reported risk estimates for type 2 diabetes by rice intake levels.
Data synthesis Relative risks were pooled using a random effects model; dose-response relations were evaluated using data from all rice intake categories in each study.
Results Four articles were identified that included seven distinct prospective cohort analyses in Asian and Western populations for this study. A total of 13 284 incident cases of type 2 diabetes were ascertained among 352 384 participants with follow-up periods ranging from 4 to 22 years. Asian (Chinese and Japanese) populations had much higher white rice consumption levels than did Western populations (average intake levels were three to four servings/day versus one to two servings/week). The pooled relative risk was 1.55 (95% confidence interval 1.20 to 2.01) comparing the highest with the lowest category of white rice intake in Asian populations, whereas the corresponding relative risk was 1.12 (0.94 to 1.33) in Western populations (P for interaction=0.038). In the total population, the dose-response meta-analysis indicated that for each serving per day increment of white rice intake, the relative risk of type 2 diabetes was 1.11 (1.08 to 1.14) (P for linear trend<0.001).
Conclusion Higher consumption of white rice is associated with a significantly increased risk of type 2 diabetes, especially in Asian (Chinese and Japanese) populations.
Introduction
Humans have a long history of cultivating rice crops; rice was first domesticated approximately 8000 to 9000 years ago by people living in the region of the Yangtze River valley in China.1 2 3 Rice is now grown worldwide and provides food for more than half of the world’s population, especially those living in some of the most populous countries, such as China, India, and Japan. Polished rice or white rice, which primarily consists of starch, is produced through a series of mechanised processes including hulling and milling,4 and it is the predominant type of rice consumed worldwide.5 6 7 Although the glycaemic index value of a specific white rice variety depends on the degree of processing, cooking time, and amylose content, the glycaemic index values of white rice are higher on average than those of whole grains.8 For example, the mean glycaemic index values were 64 (SD 7) for white rice, 55 (5) for brown rice, 41 (3) for whole wheat, and 25 (1) for barley in a previous meta-analysis.8 In addition, white rice is the primary contributor to dietary glycaemic load for populations that consume rice as a staple food.5 7
In large scale human observational studies among various populations, diets with a high glycaemic index or glycaemic load were associated with increased risk of developing type 2 diabetes.7 9 10 11 A significant positive association between white rice consumption and risk of diabetes was observed among two cohorts of Chinese and Japanese women,5 7 although the association was not significant for Japanese men.5 Two investigations in Western populations with much lower consumption levels than Asian populations also generated mixed results.12 13 These studies were heterogeneous with respect to sample size, white rice intake levels, and other characteristics that may contribute to inconsistencies in the literature. In addition, whether any dose-response relation exists between white rice consumption and risk of diabetes is unknown. Therefore, we did a meta-analysis on all published prospective cohort studies evaluating white rice intake and incidence of type 2 diabetes and quantified dose-response relations between intake of white rice and risk of type 2 diabetes.
Methods
Search strategy
We did a literature search (up to January 2012) of Medline and Embase for prospective cohort studies examining the association between rice intake and risk of type 2 diabetes. The search terms were (“Diabetes Mellitus”[Mesh] OR “diabetes”[All Fields]) AND (“Oryza sativa”[Mesh] OR “rice” OR “grain”) for Medline and (‘diabetes’/exp OR diabetes) AND (‘rice’/exp OR rice OR ‘grain’/exp OR grain) for Embase. We supplemented this search with a manual search of references cited by selected articles. One investigator (QS) did this literature search.
Study selection
We applied the following inclusion criteria: prospective cohort study, patients with self reported prevalent diabetes excluded at baseline, and point estimates of relative risk with 95% confidence intervals or standard errors available or derivable. We excluded animal studies, clinical trials, cross sectional studies, case-control studies, reviews, commentaries, letters, and studies that examined other associations. We also applied a criterion of study quality that all included cohort studies should have a loss to follow-up rate below 20%. Two investigators (EAH and AP) independently screened all studies by title or abstract and then by a full text evaluation. Any discrepancy between the two authors was solved by discussion with the senior investigator (QS).
Data extraction
We extracted the following information from each study: study’s characteristics (study name, authors, year of publication, journal, study location, duration of follow-up, person time, and number of participants and incident cases), participants’ characteristics (age and sex), exposure (rice intake levels for each category) and dietary assessment method, reproducibility and validity of assessment method, outcome (type 2 diabetes) ascertainment, and analysis strategy (statistical models, covariates included in the models, and risk estimates in each category). For studies that reported rice intake as servings per week or day, we converted it to grams per day by assuming that each serving was equivalent to 158 g of cooked rice. To convert raw rice intake levels to cooked rice consumption levels for Villegas et al’s study, we multiplied raw intake levels by a factor of 2.5.7
Two investigators (EAH and AP) extracted data independently, and any discrepancies were resolved by discussion. For two studies that expressed data separately for men and women or included data from multiple cohorts,5 13 we considered the analysis for each sex or cohort as an independent report and extracted data separately. As per our request, Hodge et al provided the number of cases and person years for each rice intake category in their study and confirmed that white rice accounted for the vast majority (approximately 95%) of total rice consumption in their study (Allison Hodge, email communication).12 To assess the quality of included studies, we derived a score that summarised 15 aspects of each study, including study design, response rate, follow-up rate, follow-up time, exposure and outcome measurements, and statistical analysis (supplementary table A).
Statistical analysis
Hazard ratios (equivalent to relative risks in cohort studies) were used as the common measure of association in all studies except for those by Nanri et al and Hodge et al, which used logistic regression to model the association of interest. Because of the very low incidence of type 2 diabetes in these two studies (five year risk was 1.9% for Nanri et al’s study and 1.2% for Hodge et al’s study), we considered odds ratios to be relatively accurate estimates of the true relative risks. We further derived estimates of person time for Nanri et al’s study by multiplying the number of participants by the follow-up time for each category of rice intake. We pooled all relative risks by using a random effects model comparing extreme categories of intake and set study weights to be equal to the inverse variance of each study’s effect estimate. We produced forest plots to assess the multivariate adjusted relative risks and corresponding 95% confidence intervals visually across studies. We evaluated heterogeneity of relative risks across studies by using the Cochrane Q statistic (we considered P<0.05 to be indicative of statistically significant heterogeneity) and the I2 statistic. We did stratified analyses according to ethnicity (Asian versus Western). To evaluate potential interactions between rice intake and ethnicity, we used Altman et al’s method to evaluate whether the pooled relative risks differed between different ethnic groups.14 We also fitted a fixed effects model to examine the between group heterogeneity and used the P value for heterogeneity to evaluate the interactions. We used Begg funnel plots and Egger’s tests to assess potential publication bias.15
To evaluate a potential non-linear dose-response relation, we first used a restricted cubic spline regression model (Stata RC_SPLINE command) with three knots to create spline variables that we subsequently used in our analysis to derive the generalised least squares trend estimation of pooled dose-response data (Stata GLST command).16 17 We then fitted another regression model without the spline terms. Lastly, we used the likelihood ratio test to examine the significance of any non-linearity by comparing the model with the linear term only and the model with both the linear and the cubic spline terms. This analysis used data from the relative risks and 95% confidence intervals, number of cases and person years, and median/mean of rice intake levels for each comparison group. Median intake levels for each rice consumption category were available in Nanri et al’s study.5 For all other studies that did not provide such data, we calculated the average of the two extreme values of each comparison group to determine mean rice intake levels. For the highest consumption category, we assumed that the average consumption level was the cut-off point plus a 25% increment, which was largely consistent with Nanri et al’s study. Because of unstable estimators in a random effects cubic spline model as a result of lack of power, we used a fixed effects model to evaluate the dose-response relation. We calculated absolute risk differences as background incidence rate×(relative risk−1).
We used Stata statistical software version 11.0 for all analyses. P values were two sided with a significance level of 0.05.
Results
Literature search
Figure 1⇓ shows results from the literature search and study selection process. We identified 825 articles from the Medline database and 2453 articles from the Embase database. After exclusion of duplicate records and studies that did not meet our inclusion criteria, 36 articles remained, and we further evaluated the full texts of these publications. Of these, we excluded five studies because they did not separate rice from other carbohydrate sources, 16 because no original data could be extracted (reviews, letters, or cross sectional studies), 10 that we deemed irrelevant, and one that had a loss to follow-up rate of 31.7%.18 Finally, four studies met the inclusion criteria and were included in the meta-analysis. A manual search of references cited by these studies did not yield new eligible articles. Among these four studies, Nanri et al’s study examined men and women separately and Sun et al’s report included data from three independent cohorts. Therefore, we included seven comparisons in the meta-analysis.
Study characteristics
Tables 1⇓ and 2⇓ show the characteristics of the included studies. All four studies were prospective cohort studies in participants who were free of self reported diabetes at baseline (total n=352 384). Hodge et al’s study further excluded any patients with diabetes whose date of diagnosis was before the study baseline even though they did not report having a diagnosis of diabetes at baseline interview. Among the participants, 13 284 incident cases of diabetes occurred during follow-up periods ranging from four to 22 years. Two studies were done in Asian populations (China and Japan) and the other two studies in Western populations (the United States and Australia). Average rice intake levels varied dramatically across studies. For example, in the Chinese study the mean intake level of cooked rice was approximately four servings (625 g) per day, whereas in the two studies done in the United States and Australia most (98% for the US study and 71% for the Australian study) participants consumed less than five servings a week. In all studies, dietary intake was assessed by food frequency questionnaires, which were validated against multiple day diet records or 24 hour recalls. Moderate correlation coefficients of dietary intake of rice have been found (ranged from 0.53 to 0.66), supporting reasonably good validity of rice intake assessment in these studies. Results from our assessment of study quality showed that most studies achieved a score of 12 or above (the maximum score was 15) except for Hodge et al’s study, which achieved a score of 7 (supplementary table A).
White rice intake and risk of type 2 diabetes
Figure 2⇓ summarises the comparisons of the highest and lowest categories of white rice intake levels. Overall, the random effects model summarising all seven comparisons suggested a positive association; the pooled relative risk was 1.27 (95% confidence interval 1.04 to 1.54), although significant heterogeneity was detected (I2=72.2%; Cochrane Q test P=0.001). With stratification by ethnicity, the association was stronger for Asian populations (pooled relative risk 1.55, 1.20 to 2.01) than for Western populations (1.12, 0.94 to 1.33). In both strata, the P for heterogeneity was not significant (P=0.17 and P=0.13). The difference in the pooled relative risks between these two groups reached statistical significance (P=0.038), suggesting an interaction between rice intake and ethnicity. Similarly, when we used a fixed effects model to examine between ethnic group heterogeneity, we found a significant P for heterogeneity (P<0.001). The funnel plot and Egger’s test (P=0.30) did not suggest evidence of publication bias (supplementary figure A).
Dose-response relation between white rice intake and diabetes risk
In a fixed effects cubic spline model that included all studies, we did not find evidence suggesting any non-linear relation between white rice consumption and risk of diabetes (P for non-linearity=0.51) (fig 3⇓). For each serving per day increment of white rice consumption, the relative risk was 1.11 (1.08 to 1.14; P for linear trend<0.001). Using the incidence rate of diabetes in the middle aged US population (15.2 cases/1000 population aged 45-64 years),19 we estimated that 167 cases of diabetes per 100 000 middle aged people would occur each year for each serving per day increase in consumption of white rice. (These estimates may underestimate the risk difference for Asian populations that are experiencing an accelerated incidence rate of diabetes.20) To illustrate this, we further plotted the incidence rate of diabetes by intake levels for each comparison category in Asian and Western populations (supplementary figure B).
Secondary analysis
We did several secondary analyses to examine the robustness of the primary results. Firstly, we evaluated a potential interaction by sex. Because Hodge et al’s did not separate men from women in their analysis, we excluded this study in this secondary analysis. The association was more pronounced among women (pooled relative risk 1.46, 1.16 to 1.83) than men (1.08, 0.87 to 1.34) (supplementary figure C). Secondly, we excluded Hodge et al’s study, which had a lower quality score than other studies and did not separate white rice from brown rice. This analysis yielded a pooled relative risk of 1.17 (0.97 to 1.42) comparing the highest and lowest categories in Western populations. The relative risk was 1.33 (1.09 to 1.63) for all populations, and a significant P for heterogeneity (P=0.006) was still present in this analysis. Lastly, when we included Yu et al’s study among Chinese people living in Hong Kong,18 which was excluded from the primary analysis because of a high rate of loss to follow-up, we found similar associations; comparing the highest and the lowest categories, the pooled relative risk was 1.45 (1.11 to 1.89) for Asian populations and 1.24 (1.03 to 1.50) for total populations.
Discussion
In this meta-analysis of prospective cohort studies, we found that higher white rice consumption was associated with a significantly elevated risk of type 2 diabetes. This association seems to be stronger for Asians than for Western populations. A dose-response analysis showed that each serving per day of white rice consumption was associated with an 11% increase in risk of diabetes in the overall population.
Strengths and limitations
Several caveats of this meta-analysis are worth discussing. Firstly, although the ethnicity stratified analysis did not show significant heterogeneity within each group, the limited number of studies may lead to diminished statistical power for detecting heterogeneity within each stratum. Secondly, although we included the results from only the fully adjusted models, because all individual studies were observational in nature the results of these studies may still be subject to residual confounding or other biases. Confounding by socioeconomic status is of particular concern because this is both a risk factor for type 2 diabetes and a predictor of rice consumption in Asian and Western populations.5 7 21 22 23 24 However, the US studies consisted of participants from the same professional background, so confounding by socioeconomic status was likely to be small. In addition, other studies controlled for indicators of socioeconomic status such as income and education. Nevertheless, residual confounding by socioeconomic status cannot be completely ruled out in these studies. Depending on the nature of uncontrolled or residual confounding, the associations seen in these individual studies and our meta-analysis could be biased in either direction. Large scale pooling projects, in which covariate adjustments and statistical analysis can be standardised, are needed to confirm the findings of this meta-analysis. Likewise, the dose-response relation could be more precisely modelled in such pooling projects.
Thirdly, all studies used food frequency questionnaires to assess levels of white rice intake. Although validation studies showed reasonable reproducibility and validity of self reported rice intake, measurement error is inevitable. Measurement error in assessment of exposure may lead to attenuation of true associations in a prospective study, especially when the exposure was assessed before disease assessment. Fourthly, although all studies excluded cases of self-reported diabetes at baseline, some undiagnosed cases may still be included in the analysis. However, the effect of such a bias is likely to be small. In the US studies, self reported diagnosis of diabetes was highly accurate; the Australian study further excluded any cases with a diagnosis date before baseline, even if they did not report diabetes at baseline interview; and in Asian studies, because rice is a staple food, substantial reduction of rice consumption after diagnosis of diabetes is unlikely. Lastly, we were unable to include brown rice in this meta-analysis or to evaluate the effects of substituting brown rice for white rice, because the association between brown rice and risk of diabetes was examined only in Sun et al’s study.13
The strengths of this meta-analysis include the large sample size and long duration of follow-up of the included studies. In addition, most established risk factors for type 2 diabetes were adjusted for in the fully adjusted models in these studies. Moreover, inclusion of studies in both Asian and Western countries allowed us to investigate the dose-response relation on the basis of a wide spectrum of white rice intake levels.
Results in relation to other studies
Several potential mechanisms could explain the association between white rice consumption and risk of type 2 diabetes. Among Asian populations, which consume white rice as a staple food, white rice is the predominant contributor to dietary glycaemic load. For example, in women living in Shanghai, white rice accounted for 73.9% of dietary glycaemic load7; in Japanese women, white rice explained 58.5% of dietary glycaemic load.25 In a meta-analysis that pooled data from cohort studies primarily done in Western populations, dietary glycaemic load was consistently associated with increased risk of developing type 2 diabetes.26 Similarly, recent investigations in Chinese and Japanese populations also support the hypothesis that high dietary glycaemic load is associated with increased risk of diabetes.7 25 27 The relatively weaker association for Western populations seen in this meta-analysis may be due to the fact that white rice intake was much lower than in Asians and, therefore, was only a minor contributor to dietary glycaemic load. In addition, the glycaemic index values of various white rice varieties depend on several factors including amylose content, other botanical structures, and processing methods.8 28 29 30 The contribution of white rice to dietary glycaemic load may vary substantially, especially when consumption levels are low. Nonetheless, high intake of white rice may also lead to increased risk of diabetes through mechanisms other than its contribution to dietary glycaemic load. Compared with minimally processed whole grains such as brown rice, white rice has a lower content of many nutrients including insoluble fibre, magnesium, vitamins, lignans, phytoestrogens, and phytic acid, which are lost during the refining process.31 Some of these nutrients, especially insoluble fibre and magnesium, have been associated with lower risk of type 2 diabetes in prospective cohort studies.9 10 32 33 34 35 36 Thus, a high consumption of white rice may lead to increased risk of diabetes because of the low intake of beneficial nutrients, in addition to its higher glycaemic load. Meanwhile, more data are needed to shed light on whether the interaction by ethnicity is due simply to substantially different white rice intake levels or to other mechanisms.
Data on the association between brown rice intake and type 2 diabetes are limited. In Sun et al’s work in Western populations, brown rice intake was associated with a modestly decreased risk of type 2 diabetes, and the substitution of brown rice or other whole grains for white rice was associated with a significantly lower risk of diabetes.13 Because Asian populations consume white rice almost exclusively, no data on the relation between brown rice and risk of diabetes are available in these populations. Nevertheless, a 16 week clinical trial in 76 Korean men showed that isocaloric replacement of white rice with whole grains and legume powder (composed of 66.6% whole grains, 22.2% legumes, 5.6% seeds, and 5.6% vegetables) led to significant reductions in serum glucose and insulin concentrations, whereas body weight remained unchanged.37 However, a recent study in Shanghai found that substituting brown rice for white rice for 16 weeks did not substantially affect metabolic markers in middle aged men and women, although high density lipoprotein cholesterol and diastolic blood pressure were significantly improved among people with diabetes through the brown rice intervention.38 More studies with larger sample sizes and longer durations of follow-up are warranted to examine the effects of substituting brown rice for white rice on risk of diabetes.
Conclusions
In summary, this meta-analysis suggests that higher white rice intake is associated with a significantly elevated risk of type 2 diabetes, especially among Asian populations. The recent transition in nutrition characterised by dramatically decreased physical activity levels and much improved security and variety of food has led to increased prevalence of obesity and insulin resistance in Asian countries.39 Although rice has been a staple food in Asian populations for thousands of years, this transition may render Asian populations more susceptible to the adverse effects of high intakes of white rice, as well as other sources of refined carbohydrates such as pastries, white bread, and sugar sweetened beverages. In addition, the dose-response relations indicate that even for Western populations with typically low intake levels, relatively high white rice consumption may still modestly increase risk of diabetes.
What is already known on this topic
The association between white rice consumption and risk of type 2 diabetes has been examined in both Asian and Western populations, but the findings are not entirely consistent across studies
The substantial difference in baseline rice intake levels between Asians and other populations for whom rice is not a staple food may contribute to the inconsistency of existing results
However, a systematic and quantitative summary of published studies to date is not available
What this study adds
Pooled data suggest that higher white rice consumption is associated with increased risk of developing type 2 diabetes in comparison with lower intake levels
This association is stronger for Asian (Chinese and Japanese) populations than for Western populations
Overall, there was a dose-response relation between higher intake of white rice and increasing risk of diabetes
Notes
Cite this as: BMJ 2012;344:e1454
Footnotes
We thank Allison M Hodge and Ruby Yu for providing data for the meta-analysis.
Contributors: EAH, AP, and QS searched the literature and extracted data. QS had the idea for the analysis. VM, AP, and QS provided statistical expertise. EAH and AP analysed the data. EAH wrote the first draft of the manuscript. All authors contributed to the interpretation of the results and critical revision of the manuscript for important intellectual content and approved the final version of the manuscript. QS is the guarantor.
Funding: QS was supported by career development award K99HL098459 from the US National Heart, Lung, and Blood Institute. The funding sources had no role in study design; in the collection, analysis, and interpretation of data; in the writing of the report; or in the decision to submit the article for publication. The authors are not affiliated with the funding institutions.
Competing interests: All authors have completed the Unified Competing Interest form at www.icmje.org/coi_disclosure.pdf (available on request from the corresponding author) and declare: no support from any organisation for the submitted work; no financial relationships with any organisations that might have an interest in the submitted work in the previous three years; no other relationships or activities that could appear to have influenced the submitted work.
Ethical approval: Not needed.
Data sharing: No additional data available.
This is an open-access article distributed under the terms of the Creative Commons Attribution Non-commercial License, which permits use, distribution, and reproduction in any medium, provided the original work is properly cited, the use is non commercial and is otherwise in compliance with the license. See: http://creativecommons.org/licenses/by-nc/2.0/ and http://creativecommons.org/licenses/by-nc/2.0/legalcode.