Overweight and obesity are the fifth leading causes of death globally( 1 ). The effectiveness of diets low in carbohydrate and high in protein and fat to promote weight loss and reduce cardiometabolic risk has been shown in several systematic review and meta-analysis studies( Reference Santos, Esteves and da Costa Pereira 2 – Reference Clifton, Condo and Keogh 5 ). However, the safety of consuming low-carbohydrate diets (LCD) remains controversial. Although some recent cohort studies have supported the long-term beneficial effects of plant-based LCD on cardiovascular morbidity or mortality( Reference Halton, Willett and Liu 6 ), other studies have suggested that LCD may increase the risk of CVD morbidity and mortality( Reference Fung, van Dam and Hankinson 7 – Reference Trichopoulou, Psaltopoulou and Orfanos 10 ). Differential effects of animal-based and plant-based LCD on health outcomes have also been indicated( Reference Halton, Willett and Liu 6 , Reference Fung, van Dam and Hankinson 7 , Reference Willett 11 ).
Very recently, a meta-analysis of observational studies on the association of LCD with total mortality has been carried out by Noto et al. ( Reference Noto, Goto and Tsujimoto 12 ). This study showed LCD to be associated with a significantly higher risk of total mortality without any changes in CVD mortality and incidence. Surprisingly, however, these authors analysed data of individuals mostly living in European countries or in the USA; no data of those living in Japan were analysed. Compared with Western populations, Japanese people living in Japan consume a higher amount of carbohydrate and a lower amount of saturated fat. Fish intake in Japan is very high compared with that in Western countries( Reference Nakamura, Ueshima and Okamura 13 – Reference Zhou, Stamler and Dennis 16 ). Thus, it is quite possible that LCD may have very different effects on health outcomes in Japanese people living in Japan.
The aim of the present study was to assess the association of usual, animal-based, and plant–fish-based LCD with CVD and total mortality in a large cohort of participants randomly recruited from selected districts in Japan.
Participants and methods
Participants
Cohort studies of the National Survey on Circulatory Disorders, Japan, are known as NIPPON DATA (National Integrated Project for Prospective Observation of Non-communicable Disease and Its Trends in the Aged). The present study analysed data obtained from NIPPON DATA80, in which baseline surveys were carried out in 1980. Details of this cohort have been reported elsewhere( Reference Okamura, Kadowaki and Hayakawa 17 ).
A total of 300 districts were randomly selected throughout Japan. The selection process is outlined below. For a comprehensive survey of living conditions, 1800 residential areas were sampled from over 900 thousand Population Census areas of the year 1975 by stratified sampling. Each area was divided into two districts. From 3600 districts, 300 were chosen by stratified sampling for the present study. A total of 13 771 individuals aged ≥ 30 years resided in the participating districts. All of them were invited to participate in the study. Of them, 10 546 community-based individuals agreed to participate in the study. The participation rate was 76·6 % (10 546 of 13 771) before exclusion for reasons mentioned below. The survey consisted of history-taking, physical examinations, blood tests and a self-administered questionnaire on lifestyle. For the present study, the participants were followed up to 2009 (NIPPON DATA80, 1980–2009).
The participants were excluded from follow-up because of missing baseline data (n 84), a past history of CHD or stroke at baseline (n 153), or loss to follow-up (n 1109). The latter group was excluded because of the absence of a permanent address that was needed to link to vital statistical records. The final sample comprised 9200 participants (5160 (56 %) women and 4040 men). There were no significant differences between participants who were lost to follow-up and those who were included in the present study in terms of several risk factors.
Biochemical and baseline examinations
The baseline surveys were conducted at public health centres according to a standardised manual. Blood pressure (BP) was measured by trained research nurses using a standard mercury sphygmomanometer by attaching it to the right arm of seated participants after at least 5 min of rest. Hypertension was defined as systolic BP ≥ 140 mmHg or diastolic BP ≥ 90 mmHg or use of antihypertensive agents. Height and weight were measured in stocking feet and light clothing. BMI was calculated as weight (kg) divided by the square of height (m2).
A lifestyle survey was also carried out using a self-administered questionnaire. The participants were also asked about their employee class (executive, professional and others (including homemaker)). Reported information was confirmed by public health nurses through face-to-face interviews.
Casual blood samples were drawn and centrifuged within 60 min of collection and stored at − 70°C until analyses as described previously( Reference Okamura, Kadowaki and Hayakawa 17 , Reference Nakamura, Okamura and Higashiyama 18 ). Serum concentrations of glucose were measured using the cupric-neocuproine method, and the value was converted so as to better correspond with those obtained from the more widely used hexokinase method( Reference Iso, Imano and Kitamura 19 ).
Nutritional survey
We used the data of the National Nutritional Survey in Japan (NNSJ). In the NNSJ, a food intake survey using weighed food records over three consecutive representative weekdays was conducted by specially trained dietary interviewers. Dietary interviewers visited participants at their homes at least once during the survey. We used the Standard Tables for Food Composition in Japan, 4th edition, with matched fatty acid and micronutrient intake values to calculate Japanese nutrient intake values. We estimated the nutrient intake of each household member by dividing household intake data of the NNSJ80 conducted in 1980 proportionally using average intakes by sex and age groups calculated for the NNSJ conducted in 1995, when collection of information on dietary intakes of individual household members started. These methods have been described in detail elsewhere( Reference Okuda, Miura and Yoshita 20 ). For each individual, means of the estimated intake of individual nutrients from the 3 d records were used in the analyses. The validation of this method was ascertained in a previous study( Reference Iwaoka, Yoshiike and Date 21 ). Although alcohol intake is an important dietary factor, it was not included in the integrated datasets. Intakes of alcoholic beverages were recorded in the NNSJ, and the alcohol intake per household was included in the NNSJ datasets. However, we thought that calculation of the alcohol intake of individuals from the intake per household would not be appropriate for the investigation of the alcohol–disease association, due to the large inter-individual differences in alcohol consumption within a household( Reference Okuda, Miura and Yoshita 20 ). Therefore, we used the values of alcohol consumption frequency recorded in NIPPON DATA. We present data as the contribution to total energy intake (%kJ) from carbohydrate, total fat, SFA, MUFA+PUFA, total protein, and vegetable protein; green and yellow vegetables, fruits, fish, meats, and eggs (g/d); total dietary fibre (g/4186 kJ (1000 kcal)); and Na and K (mg/4186 kJ (1000 kcal)).
Calculation of the low-carbohydrate diet score
The LCD score was calculated by modifying the methods of Halton et al. ( Reference Halton, Willett and Liu 6 ). Instead of calculating scores based on animal or vegetable fat, we used more specifically SFA or MUFA+PUFA values, because the harmful effects of SFA and the beneficial effects of MUFA+PUFA are known. Data are reported as a percentage of energy. Women and men were analysed separately, by dividing them into eleven strata each for fat, protein and carbohydrate intake. For fat and protein, the points ranged from 10, given to participants in the highest stratum, to 0, given to participants in the lowest stratum. For carbohydrate, the order of the strata was reversed. The points for each of the three macronutrients were then summed to obtain the overall diet score, which ranged from 0 to 30. Herein, this score is called the usual LCD score. Because previous studies have shown differential effects of animal-based and plant-based LCD on health outcomes, we also created two additional LCD scores. One was calculated according to the percentage of energy as carbohydrate, animal protein and SFA and the other according to the percentage of energy as carbohydrate, vegetable protein and MUFA+PUFA. Herein, the former is called the animal-based LCD score and the latter is called the plant–fish-based LCD score.
Endpoint determination
To determine the cause of death after 29 years of follow-up, we used the National Vital Statistics Database of Japan with permission from the Management and Coordination Agency, Government of Japan. The underlying causes of death were coded according to the 9th International Classification of Disease (ICD-9) through the end of 1994 and ICD-10 from the beginning of 1995. The details of the classification used in the present study have been described elsewhere( Reference Okamura, Hayakawa and Kadowaki 22 ). The causes of CVD deaths (ICD-9: 390–459 and ICD-10: I00–I99) were identified. Approval for the study was obtained from the Institutional Review Board of Shiga University of Medical Science (no. 12–18, 2000).
Statistical analyses
SAS version 9.4 for Windows (SAS Institute, Inc.) was used for the analyses. We divided the participants into ten deciles according to their LCD score. Age data are reported as means and standard deviations. For other continuous variables, we obtained age-adjusted data by ANCOVA, and these data are reported as means with their standard errors. Prevalence is reported as a percentage. Variables were compared across the ten deciles according to the LCD score. The χ2 test was used to compare dichotomous variables, and when statistically significant, followed by a post hoc logistic analysis adjusted for age, taking decile 1 as a reference. A one-way ANOVA was used to compare means among the groups, and when significant, followed by a post hoc Dunnett's test, taking decile 1 as a reference.
To assess the association between the LCD score and CVD and total mortality, multivariable-adjusted hazard ratios (HR) and P values were calculated using a Cox proportional-hazards model. Women and men were analysed separately and combined. Variables in model 1 were age and deciles of diet scores. Model 2 comprised model 1 variables and BMI (five categories divided at 18·5, 23, 25 and 30 kg/m2; 18·5–23 kg/m2: a reference), hypertension, cigarette smoking (never-smokers, ex-smokers, and three categories of current smokers divided at twenty and forty cigarettes per d; never-smokers: a reference), and alcohol drinking (ex-drinkers, current drinkers and never-drinkers; never-drinkers: a reference). Model 3 comprised model 2 variables and serum total cholesterol and blood glucose concentrations (standardised to have the mean at 0 (sd 1)) and serum creatinine concentrations (divided at the 75th percentile, 10 mg/l). Model 4 comprised model 3 variables and total dietary fibre intake, Na:K ratio and employee classes (executive, professional and others; others: a reference). Model 3 and 4 in women and men combined comprised model 3 and 4 variables and sex. To determine whether the food records obtained in 1980 were representative of diet in the subsequent 29 years, model 4 analyses of usual LCD score for CVD mortality in women at 14 and 19 years of follow-up were carried out. All the variables in model 4 were retained using the PROC GLMSELECT option with the LASSO selection method. Thus, collinearity was not a problem in the analysis.
Results
Descriptive statistics
During follow-up for 29 years (224 610 person-years), 1171 CVD deaths (52 % in women) and 3443 total deaths (48 % in women) were confirmed.
Criteria for determining the LCD score for women are given in Table 1 and those for men are given in Table S2 (available online). In both women and men, about 60 % of energy was from carbohydrate, about 15 % was from protein and about 20 % was from fat. For female participants in the lowest stratum of carbohydrate intake, the intake values ranged from 17·3 to 53·5 %, and for male participants, the values ranged from 18·8 to 51·6 %. For female participants in the highest stratum of carbohydrate intake, the intake values ranged from 71·3 to 83·3 %, and for male participants, the values ranged from 68·3 to 81·0 %.
* The LCD score for women was calculated by modifying the methods of Halton et al. ( Reference Halton, Willett and Liu 6 ). We divided the study participants into eleven strata each for fat, protein and carbohydrate intake, expressed as a percentage of energy. For fat and protein, the points ranged from 10, given to participants in the highest stratum, to 0, given to those in the lowest stratum. For carbohydrate, the order of the strata was reversed. The points for each of the three macronutrients were then summed to obtain the overall diet score, which ranged from 0 to 30 (usual LCD score). We also created two additional LCD scores. One score was calculated according to the percentage of energy as carbohydrate, the percentage of energy as animal protein, and the percentage of energy as SFA (animal-based LCD score) and the other score according to the percentage of energy as carbohydrate, the percentage of energy as vegetable protein, and the percentage of energy as MUFA+PUFA (plant–fish-based LCD score).
The characteristics of participants in 1980 according to the LCD score for women are given in Table 2. The mean intakes of rice, fish and shellfish, and meats were about 200, 100 and 50 g/d, respectively. For women in the higher-LCD score groups of all the three LCD scores, mean age was lower, age-adjusted rice intake was lower; age-adjusted intakes of fruits, vegetables, fish and shellfish, meats, and eggs were higher; and the prevalence of professional employee class was higher. For women in the higher-LCD score groups of two of the LCD scores, except the animal-based LCD score, age-adjusted total fibre, Na and K intakes were higher. Women in the higher-animal-based LCD score groups consumed a higher amount of age-adjusted K. The intakes of flour products were higher in the higher-LCD score groups. Age-adjusted mean serum total cholesterol concentrations were significantly higher in the higher-LCD score groups. There were no differences in mean BMI, blood glucose concentrations and total energy intake among the groups.
TCH, total cholesterol.
Mean values were significantly different from those of the decile 1 group: * P< 0·05; ** P< 0·01. A one-way ANOVA was used to compare means among the groups, followed by a post hoc Dunnett's test, taking decile 1 as a reference, when there was a significant difference in F value (P< 0·05).
† Continuous variables, except age, are age-adjusted data by ANCOVA.
‡ The χ2 test was used to compare dichotomous variables, followed by a post hoc logistic analysis adjusted for age.
The characteristics of participants in 1980 according to the LCD score for men are given in Table S3 (available online). The mean intakes of rice, fish and shellfish, and meats were about 300, 120 and 65 g/d, respectively. The characteristics of men were similar to those of women.
Associations of low-carbohydrate diets with CVD and total mortality
The results of Cox analyses on the associations of usual LCD with CVD and total mortality in women and in women and men combined are given in Table 3. The usual LCD score for women was significantly inversely associated with CVD mortality in age-adjusted model 1 (P trend= 0·030). Adjustments for confounding risk factors in women resulted in a more significant trend in model 2 (P trend= 0·024). Further adjustments for possible mediators, serum glucose and blood sugar, and a confounder, serum creatinine, resulted in a lower P trend value in model 3 (HR in decile 10 = 0·59, 95 % CI 0·38, 0·92, HR per decile = 0·97, P trend= 0·019). Finally, with further adjustment for total fibre intake, Na:K ratio and employee class included in model 4, the inverse association remained significant (HR in decile 10 = 0·60, 95 % CI 0·38, 0·94, HR per decile = 0·97, P trend= 0·021). In model 3 analysis in women and men combined with adjustment for sex also, a significant inverse association with CVD mortality was observed (HR in decile 10 = 0·74, 95 % CI 0·55, 0·99, HR per decile = 0·98, P trend= 0·033). In model 4 analysis in women and men combined, the inverse association with CVD mortality was attenuated (HR in decile 10 = 0·78, 95 % CI 0·458, 1·05, HR per decile = 0·98, P trend= 0·079).
* We calculated HR using a Cox proportional-hazards model.
† Variables in model 1 were age and deciles of diet scores.
‡ Model 1 variables and BMI (five categories divided at 18·5, 23, 25 and 30 kg/m2; 18·5–23 kg/m2: a reference), hypertension, cigarette smoking (never-smokers, ex-smokers, and three categories of current smokers divided at twenty and forty cigarettes per d; never-smokers: a reference) and alcohol consumption (ex-drinkers, current drinkers, and never-drinkers; never-drinkers: a reference).
§ Model 2 variables and serum total cholesterol and blood glucose concentrations (standardised to have the mean at 0 (sd 1)) and serum creatinine concentrations (divided at the 75th percentile, 10 mg/l).
∥ Model 3 variables and total dietary fibre intake, Na:K ratio and employee classes (executive, professional and others; others: a reference).
¶ Model 3 and 4 in women and men combined: model 3 and 4 variables and sex.
During follow-up for 14 and 19 years, 242 and 355 CVD deaths, respectively, in women were confirmed. In model 4 analyses at 14 years of follow-up, the inverse association of usual LCD score with CVD mortality was not significant (HR per decile = 0·96, P trend= 0·127), but it was statistically significant at 19 years of follow-up (HR per decile = 0·96, P trend= 0·024).
The animal-based LCD score for women was significantly inversely associated with CVD mortality in models 2 and 3 (P trend= 0·046 and 0·028), but not in model 1 (P trend= 0·064). In model 4 analysis in women, the inverse association with CVD mortality remained significant (HR in decile 10 = 0·73, 95 % CI 0·48, 1·10, HR per decile = 0·97, P trend= 0·038). In model 3 analysis in women and men combined with adjustment for sex also, a significant inverse association with CVD mortality was observed (HR in decile 10 = 0·79, 95 % CI 0·60, 1·05, HR per decile = 0·98, P trend= 0·046). In model 4 analysis in women and men combined, the inverse association with CVD mortality was attenuated (HR in decile 10 = 0·83, 95 % CI 0·63, 1·11, HR per decile = 0·98, P trend= 0·140) (Table S1, available online).
The plant–fish-based LCD score for women was marginally inversely associated with CVD mortality in models 1–3 (P trend= 0·051, 0·050, 0·050). In model 4 analysis in women, the inverse association became significant (HR in decile 10 = 0·73, 95 % CI 0·49, 1·09, HR per decile = 0·97, P trend= 0·046). In model 3 analysis in women and men combined with adjustment for sex, a significant inverse association with CVD mortality was observed (HR in decile 10 = 0·83, 95 % CI 0·63, 1·08, HR per decile = 0·98, P trend= 0·016). In model 4 analysis in women and men combined, the inverse association with CVD mortality remained significant (HR in decile 10 = 0·84, 95 % CI 0·64, 1·11, HR per decile = 0·98, P trend= 0·026) (Table S1, available online)
The usual LCD score for women was significantly inversely associated with total mortality in age-adjusted model 1 (P trend= 0·018). Adjustments for confounding risk factors in women resulted in a more significant trend in model 2 (P trend= 0·016). Further adjustments for possible mediators, serum glucose and blood sugar, and a confounder, serum creatinine, resulted in a lower P trend value in model 3 (HR in decile 10 = 0·73, 95 % CI 0·57, 0·93, HR per decile = 0·98, P trend= 0·020). In model 4 analysis in women, the inverse association remained significant (HR in decile 10 = 0·74, 95 % CI 0·57, 0·95, HR per decile = 0·98, P trend= 0·029). In model 3 analysis in women and men combined with adjustment for sex also, a significant inverse association with total mortality was observed (HR in decile 10 = 0·84, 95 % CI 0·72, 0·99, P trend= 0·030). However, in model 4 analysis in women and men combined, the inverse association with total mortality was attenuated (HR in decile 10 = 0·87, 95 % CI 0·74, 1·02, HR per decile = 0·99, P trend= 0·090) (Table 3).
The results of Cox analyses on the associations of LCD with CVD and total mortality in men are given in Table S4 (available online). Although all the three types of LCD tended to be inversely associated with CVD or total mortality, none of the associations was statistically significant.
Discussion
We found that a usual LCD was significantly inversely associated with CVD and total mortality in women and in women and men combined, independent of confounders in our large cohort. We did not observe any differential effects of animal-based and plant-based LCD. Inclusion of total fibre intake, Na:K ratio and employee class, a marker of socio-economic status (SES), in models resulted in the attenuation of statistical significance in some, but not in all, analyses.
This is the first long-term study to show a beneficial effect of a moderate usual LCD on CVD and total mortality. Despite several studies showing that LCD are effective at promoting weight loss and inducing a better cardiovascular risk profile( Reference Santos, Esteves and da Costa Pereira 2 – Reference Clifton, Condo and Keogh 5 ), previous long-term studies in people mostly living in Western countries have shown either neutral or adverse effects of LCD on total mortality, CVD morbidity and/or mortality( Reference Fung, van Dam and Hankinson 7 – Reference Trichopoulou, Psaltopoulou and Orfanos 10 ). A recent meta-analysis carried out by Noto et al. ( Reference Noto, Goto and Tsujimoto 12 ) has found that LCD are associated with a significantly higher risk of total mortality without any changes in CVD morbidity or mortality. Fung et al. ( Reference Fung, van Dam and Hankinson 7 ) showed the animal-based LCD score to be associated with higher CVD and total mortality and cancer mortality in the Nurses' Health Study and Health Professionals Follow-up Study. In contrast, Fung et al. ( Reference Fung, van Dam and Hankinson 7 ) also showed a higher vegetable-based LCD score to be associated with lower CVD and total mortality. Lagiou et al. ( Reference Lagiou, Sandin and Lof 9 ) followed up 43 396 middle-aged participants in Sweden and found LCD to be associated with an increased risk of CVD.
Exploring the reasons for differences in the results of the present study and those of previous studies carried out in Western countries is of interest. Compared with Western populations, Japanese people living in Japan consume a higher amount of carbohydrate. For both women and men, about 60 % of energy was from carbohydrate, which was about 10 points higher than that reported in previous studies in Western populations( Reference Halton, Willett and Liu 6 , Reference Zhou, Stamler and Dennis 16 , 23 ). In the present study, even for women in the lowest stratum of carbohydrate intake, the intake values ranged from 17·3 to 53·5 %, and for men, the values ranged from 18·8 to 51·6 %. Even for women and men in decile 10 of LCD score, the mean carbohydrate intake was about 50 % of total energy. Because very few participants were obese and mean BMI was not different across the deciles of LCD score, they appeared to have not followed strict carbohydrate restrictions as recommended by Atkins( Reference Atkins 24 ). In fact, this diet was not known in 1980 in Japan. For women and men in the higher-LCD score groups of all the three LCD types, age-adjusted intakes of fruits, vegetables, fish and shellfish, meats, and eggs were higher. It is understandable that the former three dietary tendencies resulted in better cardiometabolic risk profiles and outcomes judging from the results of previous studies( Reference Kelemen, Kushi and Jacobs 25 , Reference Hu and Willett 26 ). However, the latter two dietary tendencies could have resulted in worse cardiometabolic risk profiles and outcomes as reported in previous studies( Reference De Stefani, Boffetta and Ronco 27 , Reference Lam, Cross and Consonni 28 ).
A large-scale INTERMAP study has shown that dietary animal protein intake is not associated with higher BP after adjustment for height and weight, while vegetable protein intake is inversely associated with BP( Reference Elliott, Stamler and Dyer 29 ). Interestingly, several longitudinal studies have suggested a beneficial effect of meat intake on cardiovascular outcomes. Using the Hiroshima/Nagasaki Life Span Study cohort data, Sauvaget et al. ( Reference Sauvaget, Nagano and Allen 30 , Reference Sauvaget, Nagano and Hayashi 31 ) showed that the intake of animal products has protective effects against intracerebral haemorrhage and cerebral infarction. A pooled analysis of Asian prospective cohort studies has shown that red meat intake is inversely associated with CVD mortality in men( Reference Lee, McLerran and Rolland 32 ). Recently, the intake of animal protein and fat has increased significantly in Japan; however, the current mean intake of these nutrients in Japan is still low compared with that in Western countries( Reference Iso, Stampfer and Manson 33 , Reference Ueshima 34 ). It is possible that even if a very high meat intake is associated with negative outcomes, a higher intake in a lower range may have beneficial effects.
In the present study, we did not observe any differential effects of animal-based and plant-based LCD that have been indicated in studies carried out mostly in the USA. In a study carried out in Sweden, Lagiou et al. ( Reference Lagiou, Sandin and Lof 9 ) also failed to find such differential effects. This may have resulted from the above-mentioned reasons.
Although all the three types of LCD tended to be inversely associated with CVD or total mortality, none of the associations was statistically significant in men, unlike in women. Variability in dietary estimates caused by chance eating out may have been higher in men than in women because the proportion of men in full-time employment was higher than that of women; a higher prevalence of current smoking and alcohol drinking in men than in women resulted in some dilution of the effects of LCD on CHD mortality in men.
We observed an influence of SES on LCD scores: percentages of participants with a higher SES were higher in the higher-LCD score groups. A higher intake of protein and a lower intake of carbohydrates in participants in higher-SES groups have been shown in a previous study( Reference Stallone, Brunner and Bingham 35 ). However, inclusion of total dietary fibre intake, Na:K ratio and employee class, a marker of SES, in models in the present study resulted in the attenuation of statistical significance in some, but in not all, the analyses.
We recruited a large cohort of participants from randomly selected health districts in Japan. The participants of the present study were followed up for 29 years, which is a long period, and this increases the value of the present study.
Study limitations
The present study has several limitations. First, the follow-up period was long and dietary exposure may have changed over time( Reference Stallone, Brunner and Bingham 35 ). To avoid misclassification of dietary exposure over time, repeated dietary intake assessments would have been better. However, this possible variability generally tends to result in underestimated risk. The results of model 4 analyses of usual LCD score for CVD mortality in women at 14 and 19 years, in addition to 29 years, of follow-up did not indicate any significant dilution of effects over long follow-up duration. Second, information on fish oil supplementation was not available in the baseline surveys, although supplement use was not common among Japanese in 1980. Third, residual confounding from unmeasured or unknown factors cannot be excluded. For example, physical activity and amount of alcohol consumption (although we used values of alcohol consumption frequency) were not considered in the present study. Lastly, data on CVD incidence were not available; therefore, we could not examine the relationships with CVD incidence.
Conclusions
Moderate diets lower in carbohydrate and higher in protein and fat are significantly inversely associated with CVD and total mortality in women.
Supplementary material
To view supplementary material for this article, please visit http://dx.doi.org/10.1017/S0007114514001627
Acknowledgements
The authors thank all members of Japanese Association of Public Health Center Directors and all staff members of the public health centres that contributed to the present study. Details regarding the investigators and members of the research group are given in the appendix of Nakamura et al. ( Reference Nakamura, Okamura and Higashiyama 18 ).
The present study was supported in part by grants from the Ministry of Health, Labor and Welfare under the auspices of Japanese Association for Cerebro-cardiovascular Disease Control, a Research Grant for Cardiovascular Diseases (7A-2) from the Ministry of Health, Labor and Welfare, and a Health and Labor Sciences Research Grant, Japan (Comprehensive Research on Aging and Health (H11-Chouju-046, H14-Chouju-003, H17-Chouju-012, H19-Chouju-Ippan-014) and Comprehensive Research on Life-Style Related Diseases including Cardiovascular Diseases and Diabetes Mellitus (H22-Jyunkankitou-Seisyu-Sitei-017, H25-Jyunkankitou-Seisyu-Sitei)).
The authors' contributions are as follows: Y. N. and H. U. designed the study; Y. N., H. U., N. O., T. O., A. K., K. M., N. M., T. H., Y. K., A. F., M. N., N. T., T. O. and A. O. conducted the study; Y. N. and H. U. analysed the data; Y. N., H. U. and K. M. wrote the article; Y. N. had primary responsibility for the final content. All authors read and approved the final manuscript.
None of the authors has any conflicts of interest to declare.