Skip to main content

Population-based studies of relationships between dietary acidity load, insulin resistance and incident diabetes in Danes

Abstract

Background

It has been suggested that the acidity of the diet may be related to increased risk of type 2 diabetes. To investigate this hypothesis, we tested if the acidity of the diet, measured as the Potential Renal Acid Load (PRAL) score, was associated with incident diabetes and diabetes-related intermediary traits.

Methods

A total of 54,651 individuals from the Danish Diet, Cancer and Health (DCH) cohort were included in the prospective cox regression analyses of incident diabetes over a 15 years follow-up period. Moreover, 5724 Danish individuals with baseline data from the Inter99 cohort were included in the cross sectional, multivariate and logistic regression analyses of measures of insulin sensitivity, insulin release and glucose tolerance status derived from an oral glucose tolerance test (OGTT).

Results

In the DCH cohort a trend analysis showed that quintiles of PRAL score were, after multifactorial adjustment, associated with a higher incidence of diabetes (ptrend = 6 × 10− 7). HR for incident diabetes was 1.24 (1.14; 1.35) (p = 7 × 10− 7) between first and fifth PRAL score quintile.

In Inter99 higher PRAL score associated with insulin resistance as estimated by lower BIGTT-Si (an OGTT-derived index of insulin sensitivity) (p = 4 × 10− 7) and Matsuda index of insulin sensitivity (p = 2 × 10− 5) as well as higher HOMA-IR (p = 0.001). No association was observed for measures of insulin release, but higher PRAL score was associated with lower OGTT-based disposition index.

Conclusions

A high dietary acidity load is associated with a higher risk of diabetes among middle-aged Danes. Although adjustment for BMI attenuated the effect sizes the association remained significant. The increased risk of diabetes may be related to our finding that a high dietary acidity load associates with impaired insulin sensitivity.

Peer Review reports

Background

Accumulating evidence suggests that a high dietary acidity load results in chronic tissue metabolic acidosis which, in turn, may contribute to the development of insulin resistance and type 2 diabetes (T2D) [1,2,3,4,5].

Observational studies report that a high dietary acidity load associates with the risk of developing T2D: an epidemiological study of ~ 66,000 middle-aged French women, including 1372 incident T2D cases, reported a higher incidence of T2D during 14 years of follow-up in those study participants with a high dietary acidity load [6]. Also, in an analysis combining data from three observational studies of American health professionals with a total of 15,305 cases of T2D in 4,025,131 person years of follow-up, the authors reported an increased risk of T2D with a higher PRAL score [7]. Yet, this finding could not be reproduced in a Swedish study of 911 elderly men with 115 cases through 18 years of follow-up [8]. Gender discrepancies have been suggested since a study of 1191 incident cases of T2D among ~ 65,000 Japanese showed an association with T2D in men only during 5 years of follow-up [9]. Finally, higher dietary acidity load was reported to associate with increased insulin resistance in a study of 1732 Japanese (> 90% men) [10].

A widely used approach to estimate dietary acidity load is the Potential Renal Acid Load (PRAL) score which is a validated proxy for renal net acid excretion [11]. Another estimation for the acidity of the diet is the NEAP score (Net Endogenous Acid Production). Both estimates seem to reflect a similar risk of incident diabetes [6, 7, 9, 10].

The PRAL score is based on dietary intake of protein and various micronutrients, phosphorus, potassium, calcium and magnesium, and takes into account the absorption rate of the nutrients in the gut, unlike the NEAP score, which only operates with potassium and protein intake [6]. A negative PRAL score reflects an alkalizing potential of the diet whereas a positive PRAL score reflects an acidifying potential of the diet.

The aim of the present study was to substantiate and elaborate previous findings and test if PRAL score associates with impaired glucose tolerance and incident diabetes in our study sample of middle-aged people from the general Danish population. Furthermore, in a cross sectional study of middle-aged individuals from the Danish general population, we aimed to test the hypothesis that a higher PRAL score associates with diabetes-related intermediary traits, including impaired beta-cell function and insulin resistance, derived from Oral Glucose Tolerance Tests (OGTT).

Methods

The present study is based on two Danish population-based cohorts: the Diet, Cancer and Health cohort (DCH) and the Inter99 cohort (ClinicalTrials.gov ID no. NCT00289237).

Diet, Cancer and health cohort

During 1993 to 1997, 160,725 Danish men and women were invited to participate in the DCH cohort; the inclusion criteria being age 50–64 years, living in the greater Copenhagen and Aarhus areas, born in Denmark and not registered with a previous cancer diagnosis in the Danish Cancer Registry. In total, 27,178 men and 29,875 women participated. However, 574 individuals were excluded due to cancer diagnosis before baseline, leaving 56,479 participants available for analysis.

The present study is based on data from 25,808 men and 28,843 women after exclusion of patients with known diabetes at baseline (n = 1371), participants with incomplete dietary registration (n = 53), participants with extreme values of self-reported energy intake (< 1000 KJ/day (n = 0) or > 20,000 KJ/day (n = 197)) and participants with missing values for BMI and lifestyle characteristics (diet, smoking and physical activity) (n = 207). Information on incident diabetes during the 15 years of follow-up was obtained from The National Diabetes Register [12] and dates of death were obtained from the Danish Civil Registration System. All other information was collected at baseline. A flowchart is given in Additional file 1: FigureS1.

In the DCH cohort, diet was monitored at recruitment by a 192-item FFQ that each participant received by mail before their visit to the study centre. The FFQ was designed specifically for this study population, aiming to capture the average intake of different food and beverage items over the past 12 months before study inclusion. Daily intakes of foods and nutrients were calculated for each participant by the software programme FoodCalc (www.ibt.ku.dk/jesper/foodcalc/). A description of the development and validation of this FFQ, and a detailed description of the estimation of the dietary intake of the population have been published [13,14,15,16]. Information on smoking habits and physical activity was obtained from questionnaires.

Inter99 cohort

The Inter99 cohort is a non-pharmacological intervention study for the prevention of ischaemic heart disease [17]. Detailed description of the Inter99 study is given in the Additional file 1.

The present study is based on data from 2843 men and 2881 women after exclusion of participants with incomplete dietary registration (n = 150), participants with missing data from the OGTT (n = 359) and participants with extreme values of self-reported energy intake (< 1000 KJ/day (n = 6) or > 20,000 KJ/day (n = 93)). Additionally, 19 individuals had fasting serum C-peptide levels lower than 150 pmol/l and were, due to suspicion of type 1 diabetes, excluded from further analyses. Four hundred thirty-three individuals had missing information on smoking, physical activity, dietary intake or body mass index (BMI) and were thus excluded. A flowchart is given in Additional file 1: FigureS2.

Based on OGTT derived data, participants were characterised as having normal glucose tolerance (NGT) (n = 4288), impaired fasting glucose (IFG) (n = 474), impaired glucose tolerance (IGT) (n = 655) or screen detected, treatment-naive T2D (n = 214) according to the 1999 WHO criteria [18]. Additionally, a self-reported diabetes diagnosis was reported for 93 participants. In the present analytical protocol, individuals with combined IFG and IGT are presented together in the IGT group.

In the multivariate analyses of diabetes-related intermediary traits, participants with self-reported diabetes at baseline were excluded (n = 93) leaving 5631 participants eligible for analysis.

To assess beta cell function we used insulinogenic index and corrected insulin response as well as disposition index. To review insulin sensitivity we used Homeostatic Model Assessment of Insulin Resistance (HOMA-IR), Matsuda index of insulin sensitivity (ISIMatsuda) and BIGTT-Si. Calculations of these indices are given in Additional file 1.

The Inter99 study participants completed, at recruitment, a validated and self-administered food frequency questionnaire (FFQ) during their visit to the Research Centre [19]. They were asked to report their dietary intake during the month before examination. The FFQ included 198 questions on food items and beverages with additional questions regarding portion sizes of selected food items. All food items in the FFQ were linked to a food item in the Danish Food Composition Databank [20]. A detailed description of the questionnaire and estimation of the dietary intake of the population has been published [19]. Smoking status and physical activities were obtained from validated questionnaires as reported [17].

The potential renal acid load (PRAL) score

PRAL score was derived based on the estimated intake of several nutrients calculated from the FFQ used in Inter99 and DCH [6]:

$$ PRAL\ \left(\frac{mEq}{day}\right)=0.49\times protein\left(\frac{g}{day}\right)+0.037\times phosphorus\ \left(\frac{mg}{day}\right)-0.021\times potassium\ \left(\frac{mg}{day}\right)-0.013\times calcium\left(\frac{mg}{day}\right)-0.026\times magnesium\left(\frac{mg}{day}\right) $$

Statistical analyses

All statistical tests were performed using the R statistical package (http://cran.r-project.org/, version 3.1.3). A p-value of < 0.05 was considered statistically significant.

In the DCH cohort, Cox proportional hazards regression models with age as the time scale were used to estimate hazard ratios (HRs) of incident diabetes. Participants were censored at their date of death or emigration. PRAL score was analysed as categorised into quintiles, with the lowest category as the reference group. Two multivariate models were used: model 1 was adjusted for age (as time scale), smoking status (never-smoker/ex-smoker/current smoker), physical activity (metabolic equivalent of task (MET)) [21] and total fat, carbohydrate and energy intake, and model 2 was further adjusted for BMI. We tested the assumption of proportional hazards by visual inspection. Due to an uneven distribution of men and women across the quintiles and to address the previous diverging findings on men and women we conducted the analyses in combined datasets and stratified by sex.

In Inter99 logistic regression models were used to assess the association between PRAL score and IFG, IGT, and diabetes. Additionally, in Inter99 linear regression models were used to test whether PRAL score (as a continuous variable) was associated with diabetes-related intermediary traits in non-diabetic individuals only (n = 5631). Two multivariate models were used: model 1 was adjusted for age, sex, total fat, carbohydrate and energy intake, smoking status (never-smoker/ex-smoker/current smoker) and level of physical activity (0–2 h/week, 2–4 h/week, 4–7 h/week, 7–12 h/week) [22] and model 2 was further adjusted for BMI. To test whether the assumptions for linear regression analysis were fulfilled, all variables were inspected with regards to linearity as well as homogeneity of variance and normality of the residuals. Natural logarithmic transformation was used to approximate a normal distribution when needed. Test of effect modification was performed by introducing the interaction term in the linear models.

Results

Baseline characteristics of included participants in the two cohorts are shown in Table 1 as stratified by quintiles of PRAL score. Individuals with a higher PRAL score had a higher intake of total and saturated fat, lower intake of fruits and vegetables and higher daily energy intake (see Table 1). In the subsequent analyses, we adjusted for these differences.

Table 1 Baseline characteristics of study participants in the Inter99 study and in the DCH cohort by quintiles of PRAL score

Outcome from analyses of the DCH cohort

During a mean follow-up period of 15 years 7201 incident cases of diabetes occurred in the total study population. Multifactor adjusted HRs (model 1) of diabetes were 1.06 (95% CI 0.98; 1.15) (p = 0.12), 1.10 (95% CI 1.02; 1.19) (p = 0.02), 1.13 (95% CI 1.04; 1.22) (p = 0.003) and 1.24 (95% CI 1.14; 1.35) (p = 7 × 10− 7), respectively, for the second, third, fourth and fifth quintile versus the first quintile of PRAL score (ptrend = 6 × 10− 7) (Fig. 1). Further adjustment for BMI attenuated the HRs as HR for the 5th vs. the 1st quintile was 1.10 (95% CI 1.01; 1.20) (p = 0.03) (Fig. 1, model 2), but the same trend was observed across the quintiles (ptrend = 0.04). In the sex stratified analyses a similar association was observed for both men and women (Fig. 2). However, when BMI was included in the model the association with incident diabetes disappeared for men, but remained significant for women (ptrend = 0.02).

Fig. 1
figure 1

Hazard ratio (HR) of incident diabetes according to PRAL score quintiles in the full DCH cohort. Age is used as the time scale in the cox model. Model 1 is adjusted for fat, energy and carbohydrate intake, smoking and physical activity, while model 2 is adjusted for BMI, fat, energy and carbohydrate intake, smoking and physical activity. Squares indicate HR for model 1, and circles indicate HR for model 2. The p-value shown indicates significance of a trend test

Fig. 2
figure 2

Hazard ratio (HR) of incident diabetes according to PRAL score quintiles in the DCH cohort stratified by sex. Age is used as the time scale in the cox model. Model 1 is adjusted for fat, energy and carbohydrate intake, smoking and physical activity, while model 2 is adjusted for BMI, fat, energy and carbohydrate intake, smoking and physical activity. Squares indicate HR for model 1, and circles indicate HR for model 2. The p-value shown indicates significance of a trend test

Outcome from analyses of the Inter99 study

The multifactor-adjusted odds ratio (OR) for a one standard deviation higher PRAL score was 1.11 (95% CI 1.01; 1.22) (p = 0.03) for IGT (Fig. 3). The OR diminished with further adjustment for BMI (p = 0.11). No significant association was observed between PRAL score and T2D (Fig. 3) and PRAL score and IFG. However, when adjusting for BMI we found an association between lowered PRAL score and IFG with an OR of 0.89 (95% CI 0.80; 0.99) (p = 0.03).

Fig. 3
figure 3

PRAL score and glucose tolerance status among participants of the Inter99 study. Odds ratios (OR) of impaired fasting glucose (IFG), impaired glucose tolerance (IGT) and type 2 diabetes (T2D) according to PRAL score as shown for one standard deviation of PRAL score in the total population (16.36 mEq/day). Model 1 is adjusted for age, sex, smoking, physical activity and fat, energy and carbohydrate intake, while model 2 is adjusted for age, sex, BMI, smoking, physical activity and fat, energy and carbohydrate intake. NGT = Normal Glucose Tolerance. Participants having combined IFG and IGT were classified as IGT. Squares indicate OR for model 1, and circles indicate OR for model 2

A higher PRAL score was associated with lower OGTT-based measures of insulin resistance (model 1) as expressed by decreased BIGTT-Si (p = 4 × 10− 7) and ISIMatsuda (p = 2 × 10− 5) as well as increased HOMA-IR (p = 0.001) (Table 2). Accordingly, a higher PRAL score was associated with higher serum insulin levels at fasting and 120 min during the OGTT (p = 2 × 10− 4 and p = 8 × 10− 16, respectively), and with higher plasma glucose levels at 120 min (p = 4 × 10− 10) (Table 2). No associations were observed between PRAL score and corrected insulin response (p = 0.2) or insulinogenic index (p = 0.3), both of which are indices of beta cell function. In addition, a higher PRAL score was associated with a lower OGTT-based disposition index (p = 0.004). All findings remained significant after further adjustment for BMI (model 2) and with similar effect sizes.

Table 2 Associations between PRAL score and diabetes-related intermediary traits in the baseline part of Inter99 study as well as means and medians for the individual quintiles

When testing for effect modification we introduced the interaction term PRAL×sex in the linear regression models, but the interaction terms were not significant (data not shown). However; to address diverging previous findings in men and women we further stratified by sex. We found that a higher PRAL score was associated with indices of insulin sensitivity in both women and men (Table 3), but the association with increased HOMA-IR with increased PRAL score was only apparent in women. Furthermore, an association between higher PRAL score and lower HbA1c was seen in men, but only after adjustment for BMI. No associations were found between PRAL score and indices of beta cell function in the sex stratified analyses, but a higher PRAL score was associated with lower disposition index in women (Table 3).

Table 3 Associations in the Inter99 study between PRAL score and diabetes-related intermediary traits stratified by sex as well as means and medians for the individual quintiles

Sensitivity analysis

In the Inter99 cohort we further did a sensitivity analysis where we adjusted for family history of diabetes, hypertension and dietary patterns (see Additional file 1: TableS1). This information was not available to us in the DCH cohort, and so the sensitivity analysis was only conducted in the Inter99 cohort. The sensitivity analysis showed comparable associations as demonstrated in the primary analyses except the associations between PRAL score and HbA1c and disposition index. In the sensitivity analysis the two latter associations became non-significant (see Additional file 1: TableS1).

Discussion

In the DCH cohort of more than 54,000 individuals with ~ 7200 incident cases of diabetes after 15 years of follow-up, we demonstrated a positive relationship between PRAL score and development of incident diabetes, substantiating previous findings [6,7,8,9,10]. Moreover, in baseline data from the Inter99 study of 5631 non-diabetic individuals with OGTT-derived data, we found that a higher PRAL score was associated with reduced insulin sensitivity, but not with changes in proxies of beta cell function. In addition, higher PRAL score was associated with a lower OGTT-based disposition index. We did not observe any major diverging results between the sexes, and so we cannot confirm that any specific differences exist in regards to PRAL score and sex. When adjusting for BMI in the DCH cohort, only the increased HR between the first and the fifth quintiles remained significant. Furthermore, when adjusting for BMI, the significant trend seen in males disappeared, but remained significant in women. In the Inter99 cohort, adjustment for BMI had no effect on the analyses. Since the individuals with a higher PRAL score generally had a higher overall energy intake, the increased BMI may be part of the explanation for the increased incidence of diabetes. However, since the associations persist even after adjustment for BMI, it shows that there is indeed an association between PRAL score and incident diabetes as well as diabetes intermediary traits which is independent of BMI.

Despite a lack of association with indices of beta-cell function in the Inter99 study, we did find an association with lower disposition index as an indication of a decreased beta cell secretion of insulin at the concomitant level of insulin resistance. Hence, we cannot completely discard an association between PRAL score and beta-cell function. Higher PRAL score is consistent with a generally unhealthy diet, which is reflected in our finding that individuals with a higher PRAL score also had a higher intake of total and saturated fat, lower intake of fruits and vegetables and higher daily energy intake. We have accounted for these confounders by adjusting for dietary fat, carbohydrate and total energy intake in our analyses.

The different nutrients, especially protein, embedded in the calculation of the PRAL score may come from various sources. Since the amino acid composition of plant proteins is different from the composition of animal proteins, it is possible that ingestion of plant proteins will have a different effect on the tissue acidity than animal proteins. Indeed, a previous study showed an association between incident T2D and animal protein intake, but not with plant protein intake [23]. One could therefore speculate that a type of measure where proteins were reported as of their origin or corresponding amino acids would provide useful information about the underlying theories of dietary acidity load and diabetes.

The amount of dietary fiber ingested are surprisingly similar across the quintiles. However, the food items from where the dietary fiber is obtained may vary considerably and may therefore have different impacts on the PRAL score. For instance, dietary fibers originating from vegetables and fruits will have an alkalizing effect and thus decrease the PRAL score, whereas dietary fiber originating from whole grain may have a different effect on the PRAL score [1]. As appears from the baseline table (Table 1) of the DCH cohort, the dietary fiber ingested in the lower quintiles must come from fruits and vegetables, whereas the dietary fiber in the higher quintiles originate from other food items, although we cannot say for certain which ones. Additionally, it is noteworthy that the quintiles with the highest PRAL scores are also the ones with the highest number of current smokers. This is the case in both the DCH cohort and the Inter99 cohort. It cannot be excluded that smoking status might influence food preferences. Indeed, a study by Endoh et al. shows, that the dietary patterns in Japanese men and women differ between smokers and non-smokers [24].

An intervention study conducted in American adults showed that individuals would be able to lower their PRAL score by 13 units when following a vegan diet for two to three days a week, while a reduction of almost 30 units was seen in individuals following the vegan diet every day in a week [25]. If we examine our data from the DCH cohort, we find that a reduction of 30 units would be sufficient to move an individual from Q4 to Q1. With a reduction of just 12 units an individual would step down an entire quintile. Despite the modest increased risk of diabetes between the lowest and the highest quintile, diabetes is still major burden in society today, and even slight reductions in this risk through a diet change might have considerable impact on public health.

Our study has limitations. First, in the present epidemiological studies, we did not use the hyperinsulinemic euglycemic clamp to measure insulin sensitivity, as this would not be feasible in such large cohorts; still the surrogate measures applied in the current study, HOMA-IR, ISIMatsuda and BIGTT-Si, have all been validated in previous studies and provide physiologically relevant proxies for insulin sensitivity [26,27,28]. Additionally, to overcome any uncertainties in the individual surrogate measures, we have used several different proxies of insulin sensitivity to measure the same outcome. Second, and most important, in our study of incident diabetes, information was collected at baseline only. Consequently, all conclusions based on the DCH cohort rely on the assumption that the study participants have not changed their diet or lifestyle substantially over the following 15 years of follow-up. Third, the diagnosis of diabetes was based on The Danish National Diabetes Register which does not clearly distinguish between patients with type 1 diabetes or T2D [12]. Yet, individuals in our study were at least 50 years of age at baseline, and since newly diagnosed type 1 diabetes is quite rare in Denmark above this age we do not expect this lack of diagnostic accuracy to bias our analytical outcomes. The PRAL score in our two cohorts ranged from around − 100 to 100 and with a median range of − 20.6 to 18.6 in the Inter99 cohort and − 16.3 to 18.0 in the DCH cohort, respectively. These median ranges are roughly the same as in the two studies by Kiefte-de Jong et al. (− 17.4 to 21.3) and Fagherazzi et al. (− 23.0 to 14.3). However, the PRAL score range in the study by Xu et al. was from − 40 to 323. For histograms of the PRAL scores in the two cohorts, please see Additional file 1: FigureS3 and FigureS4. We believe that these between study differences may be related to differences in dietary patterns and the demographics of the study populations. Furthermore, the variations in PRAL score may occur due to differences in how the questionnaires and food composition tables are structured between studies.

Conclusions

We confirm the association of a high acidity load with incident T2D and suggest that the risk of T2D might be mediated partially through a decrease in insulin sensitivity. Although the present study and previous reports suggest a link between dietary acidity load and incident diabetes among middle-aged people, carefully designed and conducted dietary interventions are needed to elucidate whether a causal link exists between the two and which mechanisms might be involved.

Abbreviations

BIGTT:

Beta-cell function; Insulin sensitivity; Glucose; Tolerance; Test

BIGTT-Si:

An OGTT-derived index of insulin sensitivity (Si)

DCH:

The Diet, Cancer and Health cohort

FFQ:

Food frequency questionnaire

IFG:

Impaired fasting glucose

IGT:

Impaired glucose tolerance

ISIMatsuda :

Matsuda index of insulin sensitivity

MET:

Metabolic equivalent of task

NEAP:

Net Endogenous Acid Production

OGTT:

Oral glucose tolerance test

PRAL:

Potential Renal Acid Load

T2D:

Type 2 diabetes

References

  1. Adeva MM, Souto G. Diet-induced metabolic acidosis. Clin Nutr. 2011;30:416–21.

    Article  CAS  Google Scholar 

  2. Engberink MF, Bakker SJ, Brink EJ, van Baak MA, van Rooij FJ, Hofman A, et al. Dietary acid load and risk of hypertension: the Rotterdam study. Am J Clin Nutr. 2012;95:1438–44.

    Article  CAS  Google Scholar 

  3. Souto G, Donapetry C, Calvino J, Adeva MM. Metabolic acidosis-induced insulin resistance and cardiovascular risk. Metab Syndr Relat Disord. 2011;9:247–53.

    Article  CAS  Google Scholar 

  4. Williams RS, Kozan P, Samocha-Bonet D. The role of dietary acid load and mild metabolic acidosis in insulin resistance in humans. Biochimie. 2016;124:171–7.

    Article  CAS  Google Scholar 

  5. Della, Guardia L, Thomas MA, Cena H. Insulin sensitivity and glucose homeostasis can be influenced by metabolic acid load. Nutrients. 2018;10(5):618.

    Article  Google Scholar 

  6. Fagherazzi G, Vilier A, Bonnet F, Lajous M, Balkau B, Boutron-Rualt MC, et al. Dietary acid load and risk of type 2 diabetes: the E3N-EPIC cohort study. Diabetologia. 2014;57:313–20.

    Article  CAS  Google Scholar 

  7. Kiefte-de Jong JC, Li Y, Chen M, Curhan GC, Mattei J, Malik VS, et al. Diet-dependent acid load and type 2 diabetes: pooled results from three prospective cohort studies. Diabetologia. 2017;60:270–9.

    Article  CAS  Google Scholar 

  8. Xu H, Jia T, Huang X, Riserus U, Cederholm T, Arnlov J, et al. Dietary acid load, insulin sensitivity and risk of type 2 diabetes in community-dwelling older men. Diabetologia. 2014;57:1561–8.

    Article  CAS  Google Scholar 

  9. Akter S, Kurotani K, Kashino I, Goto A, Mizoue T, Noda M, et al. High dietary acid load score is associated with increased risk of type 2 diabetes in Japanese men: the Japan public health center-based prospective study. J Nutr. 2016;146:1076–83.

    Article  CAS  Google Scholar 

  10. Akter S, Eguchi M, Kuwahara K, Kochi T, Ito R, Kurotani K, et al. High dietary acid load is associated with insulin resistance: the Furukawa nutrition and health study. Clin Nutr. 2016;35:453–9.

    Article  CAS  Google Scholar 

  11. Remer T, Manz F. Estimation of the renal net acid excretion by adults consuming diets containing variable amounts of protein. Am J Clin Nutr. 1994;59:1356–61.

    Article  CAS  Google Scholar 

  12. Carstensen B, Kristensen JK, Marcussen MM, Borch-Johnsen K. The National Diabetes Register. Scand J Public Health. 2011;39:58–61.

    Article  Google Scholar 

  13. Overvad K, Tjonneland A, Haraldsdottir J, Ewertz M, Jensen OM. Development of a semiquantitative food frequency questionnaire to assess food, energy and nutrient intake in Denmark. Int J Epidemiol. 1991;20:900–5.

    Article  CAS  Google Scholar 

  14. Tjonneland A, Overvad K, Haraldsdottir J, Bang S, Ewertz M, Jensen OM. Validation of a semiquantitative food frequency questionnaire developed in Denmark. Int J Epidemiol. 1991;20:906–12.

    Article  CAS  Google Scholar 

  15. Tjonneland A, Haraldsdottir J, Overvad K, Stripp C, Ewertz M, Jensen OM. Influence of individually estimated portion size data on the validity of a semiquantitative food frequency questionnaire. Int J Epidemiol. 1992;21:770–7.

    Article  CAS  Google Scholar 

  16. Haraldsdottir J, Tjonneland A, Overvad K. Validity of individual portion size estimates in a food frequency questionnaire. Int J Epidemiol. 1994;23:786–96.

    CAS  PubMed  Google Scholar 

  17. Jorgensen T, Borch-Johnsen K, Thomsen TF, Ibsen H, Glumer C, Pisinger C. A randomized non-pharmacological intervention study for prevention of ischaemic heart disease: baseline results Inter99. Eur J Cardiovasc Prev Rehabil. 2003;10:377–86.

    Article  Google Scholar 

  18. WHO. Definition, diagnosis and classification of diabetes mellitus and its complications. In: Part 1: diagnosis and classification of diabetes mellitus; 1999.

    Google Scholar 

  19. Toft U, Kristoffersen L, Ladelund S, Bysted A, Jakobsen J, Lau C, et al. Relative validity of a food frequency questionnaire used in the Inter99 study. Eur J Clin Nutr. 2008;62:1038–46.

    Article  CAS  Google Scholar 

  20. Møller A SE, Christensen A, Hartkopp H. Danish Food Composition Databank version 5.0 (Fødevaredatabanken V 5.0). Danish Institute for Food and Veterinary research, Ministry of Food, Fisheries, and AgricultureSøborg, Denmark 2002.

  21. Friedenreich C, Norat T, Steindorf K, Boutron-Ruault MC, Pischon T, Mazuir M, et al. Physical activity and risk of colon and rectal cancers: the European prospective investigation into cancer and nutrition. Cancer Epidemiol Biomark Prev. 2006;15:2398–407.

    Article  Google Scholar 

  22. von Huth Smith L, Borch-Johnsen K, Jorgensen T. Commuting physical activity is favourably associated with biological risk factors for cardiovascular disease. Eur J Epidemiol. 2007;22:771–9.

    Article  Google Scholar 

  23. van Nielen M, Feskens EJ, Mensink M, Sluijs I, Molina E, Amiano P, et al. Dietary protein intake and incidence of type 2 diabetes in Europe: the EPIC-InterAct case-cohort study. Diabetes Care. 2014;37:1854–62.

    Article  CAS  Google Scholar 

  24. Endoh K, Kuriki K, Kasezawa N, Tohyama K, Goda T. Association between smoking status and food and nutrient consumption in Japanese: a large-scale cross-sectional study. Asian Pac J Cancer Prev. 2015;16:6527–34.

    Article  Google Scholar 

  25. Cosgrove K, Johnston CS. Examining the impact of adherence to a vegan diet on Acid-Base balance in healthy adults. Plant Foods Hum Nutr. 2017;72(3):308–13.

    Article  CAS  Google Scholar 

  26. Gutch M, Kumar S, Razi SM, Gupta KK, Gupta A. Assessment of insulin sensitivity/resistance. Indian J Endocrinol Metab. 2015;19:160–4.

    Article  CAS  Google Scholar 

  27. Hansen T, Drivsholm T, Urhammer SA, Palacios RT, Volund A, Borch-Johnsen K, et al. The BIGTT test: a novel test for simultaneous measurement of pancreatic beta-cell function, insulin sensitivity, and glucose tolerance. Diabetes Care. 2007;30:257–62.

    Article  CAS  Google Scholar 

  28. Matthews DR, Hosker JP, Rudenski AS, Naylor BA, Treacher DF, Turner RC. Homeostasis model assessment: insulin resistance and beta-cell function from fasting plasma glucose and insulin concentrations in man. Diabetologia. 1985;28:412–9.

    Article  CAS  Google Scholar 

  29. Tjonneland A, Olsen A, Boll K, Stripp C, Christensen J, Engholm G, et al. Study design, exposure variables, and socioeconomic determinants of participation in diet, Cancer and health: a population-based prospective cohort study of 57,053 men and women in Denmark. Scand J Public Health. 2007;35:432–41.

    Article  Google Scholar 

Download references

Acknowledgements

We are indebted to the staff and participants of the Inter99 and DCH study for their important contributions. The authors wish to thank A. Forman, T. Lorentzen, B. Andreasen, G.J. Klavsen and M.J. Nielsen for technical assistance and T.F. Toldsted and G. Lademann for management assistance. The DCH study was funded by the Danish Cancer Society.

Funding

The Novo Nordisk Foundation Center for Basic Metabolic Research is an independent Research Center at the University of Copenhagen partially funded by an unrestricted donation from the Novo Nordisk Foundation (www.metabol.ku.dk).

Author information

Authors and Affiliations

Authors

Contributions

JG: compiled and interpreted all data analyses, made tables and figures, wrote the first draft of the manuscript and delivered subsequently substantial contributions to the manuscript. TN: supported and supervised JG in analysis and interpretation of the data and delivered substantial contributions to the manuscript. MM: interpreted all data analyses and delivered substantial contributions to the manuscript. UT: took part in the execution and epidemiological analysis of the Inter99 Study. TJ: took part in the execution and epidemiological analysis of the Inter99 Study. KO: took part in the execution and epidemiological analysis of the Diet, Cancer and Health cohort. AT: took part in the execution and epidemiological analysis of the Diet, Cancer and Health cohort. TH: took part in the execution and epidemiological analysis of the Inter99 Study. KA: supported and supervised JG in analysis and interpretation of the data and delivered substantial contributions to the manuscript. OP: conceived and designed and supervised the present study and provided substantial contributions to the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Oluf Pedersen.

Ethics declarations

Ethics approval and consent to participate

The DCH study was approved by the regional Ethical Committees (jr.nr (KF) 01–0345/93) and was conducted in accordance with the principles of the Declaration of Helsinki [29].The Inter99 study was approved by the Ethical Committee of Copenhagen County (KA 98155) and was conducted in accordance with the principles of the Declaration of Helsinki.

Consent for publication

Not applicable

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional files

Additional file 1:

Supplements. (DOCX 85 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Gæde, J., Nielsen, T., Madsen, M.L. et al. Population-based studies of relationships between dietary acidity load, insulin resistance and incident diabetes in Danes. Nutr J 17, 91 (2018). https://doi.org/10.1186/s12937-018-0395-1

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12937-018-0395-1

Keywords