Introduction

Japan has the world’s longest life expectancy (81.1 years for men and 87.1 years for women) [1]. This may be partly attributable to the Japanese dietary pattern. In fact, it has been reported that the Japanese dietary pattern was associated with a lower risk for all-cause and cardiovascular disease (CVD) mortality [2, 3]. However, evidence on the association between the Japanese diet and mortality risk is still lacking.

The 8-item Japanese Diet Index (JDI8) has been used to assess the degree of adherence to the Japanese diet [2, 4]. The JDI8 consists of the following eight components: high intake of rice, miso soup, seaweeds, pickles, green and yellow vegetables, fish, and green tea; and low intake of beef and pork. Studies have shown that a higher JDI8 score is associated with health benefits such as a lower risk of mortality [2] and disability [4]. However, these previous studies were conducted in a specific rural area of Japan. The Japan Public Health Center-based Prospective (JPHC) Study conducted a baseline survey of 140,420 registered residents aged 40–69 years within 11 public health center areas nationwide from 1990–1994 [5].

The aim of the present study was to examine the association between adherence to the Japanese diet assessed by JDI8 and the subsequent risk for all-cause and cause-specific mortality using data from a large-scale cohort from settings all over Japan.

Materials and methods

Study cohort

The JPHC Study was launched in 1990 and 1993 for cohorts I and II, respectively [5]. Participants in cohort I included residents aged 40–59 years from five Japanese public health center areas (Iwate, Akita, Nagano, Okinawa-Chubu, and Tokyo), and those in cohort II included residents aged 40–69 years from six other Japanese public health center areas (Ibaraki, Niigata, Kochi, Nagasaki, Okinawa-Miyako, and Osaka). A baseline survey questionnaire was distributed to 140,420 registered residents mostly by hand. Participants were informed of the objectives of the study, and that completion of the survey questionnaire was regarded as providing consent to participate. The 5- and 10-year follow-up surveys (second and third surveys, respectively) were conducted to update information on lifestyle habits and health conditions in 1995–1998 and 2000–2003, respectively. The present study used the second survey as the baseline.

A total of 136,163 subjects were eligible for participation in this study. Of those, 102,341 participants (75.2%) who provided valid responses formed the study cohort. We excluded 5116 participants who reported extreme total energy intake (sex-specific values outside of ± 2.5%) and 4256 participants who reported a history of disease, including cancer, stroke, myocardial infarction, and chronic liver disease. Thus 92,969 participants (42,700 men and 50,269 women) were analyzed (Fig. 1).

Fig. 1
figure 1

Flowchart of study participants

Exposure (8-item Japanese Diet Index score)

A semi-quantitative food frequency questionnaire (FFQ) was used to assess the average intake of 147 food and beverage items in the second survey [6, 7]. For most food items, participants were asked about consumption frequency and usual portion size. The validity and reproducibility of the FFQ have been established in previous studies [8,9,10].

Based on previous studies [2, 4], we used the JDI8 to assess the degree of adherence to the Japanese diet. The JDI8 consists of eight components: rice, miso soup, seaweeds, pickles, green and yellow vegetables (green vegetables, carrot, pumpkin, and tomato), fish (raw fish, salty fish, dried fish, seafood, canned tuna, and fish products), green tea, and beef and pork (beef, pork, and processed meat). The first seven components represent adherence to the Japanese diet, and participants received one point if their intake was more than or equal to the sex-specific median. The 8th component represents non-adherence to the Japanese diet, and participants received one point if their intake was less than the sex-specific median. Supplemental Table 1 shows the sex-specific median of intake of the JDI8 components. The JDI8 score ranged from 0 to 8, with higher scores indicating greater dietary conformity. Incidentally, the JDI8 was modified to exclude coffee from the original JDI, because lower coffee consumption was significantly associated with an increased risk of all-cause mortality [2]. Additionally, previous studies using a meta-analysis approach reported that coffee consumption was inversely associated with all-cause and CVD mortality [11].

Endpoints

The participants’ residency and vital status were followed up from baseline to December 31, 2016, using the residential registry. Causes of death were confirmed by death certificates and were defined according to International Classification of Diseases, 10th revision (ICD-10). The major endpoints of this study were mortality from all-cause, cancer (ICD-10: C00–C97) and CVD (ICD-10: I00–I99). CVD mortality was subdivided into mortality from heart disease (ICD-10: I20–I52) and cerebrovascular disease (ICD-10: I60–I69).

Ethical issues

This study was approved by the Institutional Review Board of the National Cancer Center of Japan (approval code; 2015-085) and the Ethics Committee of Tohoku University Graduate School of Medicine (approval code; 2018-1-321).

Statistical analysis

We counted the person-years of follow-up for each subject from the date of response to the second survey questionnaire until the date of death or the end of the study period (December 31, 2016), whichever occurred first.

The adjusted Cox proportional hazards model was used to calculate the hazard ratios (HRs) and 95% confidence intervals (CIs) for all-cause and cause-specific mortality according to the groups of the JDI8 score (G1, score of 0–2; G2, score of 3; G3, score of 4–5; and G4, score of 6–8). The lowest group was used as the reference category. Multivariable-adjusted models were adjusted for the following variables. Model 1 was adjusted for age (45–49, 50–54, 55–59, 60–64, 65–69, or ≥ 70 years), sex, and study area (11 areas). To examine whether the association between adherence to the Japanese diet and mortality was attributable to physical health status or other lifestyle factors, Model 2 was further adjusted for body mass index (< 18.5, 18.5–24.9, 25–29.9, ≥ 30 kg/m2, or missing), smoking status (current, former, never, or missing), alcohol drinking (< 1 time/month, 1–3 times/month, 1–2, 3–4, 5–6 times/week, every day, or missing), total physical activity (quartile of metabolic equivalent task-hours/day), medication (antihypertensive, cholesterol-lowering, or hypoglycemic agents [yes or no for each item]), and occupation (agriculture, forestry, fishery, office work, self-employed, specialty work, housewife, unemployed, or other [yes or no for each item]). In addition to these adjustments, Model 3 was further adjusted for total energy intake (in kcal/d; sex-specific quartile categories). To test for linear trends, categories indicating an increase in the JDI8 score (scored as 1 for G1, 2 for G2, 3 for G3, and 4 for G4), and variables using the JDI8 score as a continuous variable were entered in the corresponding Cox model. The proportional hazards assumption was checked based on the Schoenfeld residual test using Stata command "estat phtest". In this analysis, the proportional hazards assumption was not rejected with the Schoenfeld residual test (P = 0.275).

We also conducted several sensitivity analyses to test the robustness of our findings. First, considering possible reverse causality, we conducted an analysis after exclusion of deaths occurring in the first 3 years of follow-up (n = 91,589). Second, we conducted analyses stratified by sex to assess whether effects of adherence to the Japanese diet on mortality varied by sex (male = 42,700, female = 50,269). Third, to minimize potential bias attributed to total energy intake, we calculated the JDI8 score using the data after adjustment by the residual method (n = 92,969). Fourth, for the same reason (i.e., to minimize potential bias attributed to total energy intake), we also analyzed the association including participants who reported extreme total energy intake (n = 97,802). Fifth, we performed an analysis that included sodium and soybean products (tofu, fermented soybean [natto], and soymilk, etc.) intakes in the covariates because it was reported that the Japanese diet was associated with high sodium intake [12, 13], and soybean products are characteristic of the Japanese diet [14] (n = 92,969). In addition, we investigated the relationship between each component of the JDI8 and all-cause and cause-specific mortality. Here, we conducted an analysis for each item by including it into the Cox model with full adjustment, respectively.

The Cox proportional hazards model was performed using SAS version 9.4 (SAS Inc., Cary, NC, USA). The proportional hazards assumption was performed by Stata/MP15 (StataCorp, Texas, USA). All statistical tests described here were two-sided, and differences at P < 0.05 were considered significant.

Results

Among the 92,969 participants, the proportion of men was 45.9%, the mean age (standard deviation [SD]) was 56.5 (7.8) years, median follow-up time was 18.9 years, follow-up rate was 99.7%, and follow-up for cause of death (in decedents) was over 98.0%.

Table 1 shows characteristics of study participants according to the JDI8 score groups. Participants with a higher JDI8 score were older, more likely to be male, to drink alcohol every day, to engage in more physical activity, to have greater total energy intake, and less likely to have high BMI (25.0 to < 30.0 and ≥ 30 kg/m2). In contrast, participants with a lower JDI8 score were younger and current smokers.

Table 1 Characteristics according to groups of the JDI8 score (n = 92,969)

During 1,635,302 person-years of follow-up, we identified 20,596 all-cause deaths, 7148 cancer deaths, 4990 CVD deaths, including 2600 heart disease deaths and 1950 cerebrovascular disease deaths. Table 2 shows the association between groups by JDI8 score and all-cause, cancer, CVD, heart disease, and cerebrovascular disease mortality, along with HRs and 95% CIs. The JDI8 score was inversely associated with risk of all-cause, CVD, and heart disease mortality. The multivariable-adjusted HR (95% CI) of all-cause, CVD, and heart disease mortality for the highest group of the JDI8 score was 0.86 (0.81–0.90, P trend < 0.001), 0.89 (0.80–0.99, P trend = 0.007) and 0.89 (0.77–1.03, P trend = 0.037), respectively compared to the lowest group. In addition, each one-point increase in the JDI8 score was linearly and inversely associated with mortality risk. The multivariable-adjusted HR (95% CI) of all-cause, CVD, and heart disease mortality for each one-point increase of the JDI8 score was 0.97 (0.96–0.98, P trend < 0.001), 0.97 (0.95–0.99, P trend = 0.003) and 0.96 (0.94–0.99, P trend = 0.008), respectively. The risk of cancer and cerebrovascular disease mortality decreased with higher JDI8 score after adjustment for age, sex, and area. However, these associations were attenuated after adjustment for covariates. The multivariable-adjusted HR (95% CI) of cancer and cerebrovascular disease mortality for the highest JDI8 score group was 0.94 (0.86–1.03, P trend = 0.175) and 0.89 (0.75–1.05, P trend = 0.237), respectively compared to the lowest JDI8 score group.

Table 2 Relationship between the JDI8 score and mortality (n = 92,969)

These results did not change substantially after excluding 1,380 participants who died in the first 3 years of follow up (Supplemental Table 2). In addition, these inverse associations did not differ between the sexes (Supplemental Table 3). To minimize potential bias attributed to total energy intake, even when we used the data after adjustment by the residual method (Supplemental Table 4) and also when including participants who reported extreme total energy intake (Supplemental Table 5), the JDI8 score was inversely associated with the risk of all-cause, CVD, and heart disease mortality.

Table 3 shows the associations between each component of the JDI8 and all-cause, cancer, CVD, heart disease, and cerebrovascular disease mortality. Except for rice and miso soup, consuming the five other adherence components resulted in a significant inverse association with all-cause mortality; the multivariable-adjusted HRs (95% CI) were 0.94 (0.92–0.97) for seaweeds, 0.95 (0.92–0.98) for pickles, 0.94 (0.91–0.96) for green and yellow vegetables, 0.97 (0.94–0.997) for fish, and 0.89 (0.86–0.91) for green tea. No significant association was observed between the non-adherence component and all-cause mortality. The multivariable-adjusted HR (95% CI) was 1.02 (0.99–1.05) for beef and pork.

Table 3 Relationship between the JDI8 score components and mortality (n = 92,969)

In addition, when we performed the analysis that included sodium and soybean products intake in the covariates, we observed an inverse association between the JDI8 score and the risk of all-cause, CVD, and heart disease mortality as well as cancer mortality risk (Supplemental Table 6).

Discussion

We investigated the association between adherence to the Japanese diet and all-cause and cause-specific mortality in this cohort study covering multiple areas of Japan. We observed that a higher JDI8 score was significantly associated with lower risks of all-cause, CVD, and heart disease mortality.

Our findings of inverse associations between adherence to the Japanese diet and all-cause and CVD mortality are consistent with previous studies, which reported that both sexes who adhered to the Japanese diet had a 9% lower risk of all-cause mortality [2] and a 26% lower risk of CVD mortality [3]. Therefore, the results of our large-scale cohort from settings all over Japan support the findings in these previous studies.

In our analysis using each component of the JDI8 as an exposure variable, we found that higher consumption of seaweeds, pickles, green and yellow vegetables, fish, and green tea were associated with a decreased risk of all-cause, CVD, and heart disease mortality. These foods contain many beneficial nutrients, for instance, vegetables contains a myriad of nutrients and phytochemicals, including fiber, vitamin C, potassium, carotenoids, antioxidants, and flavonoids [15]; fish contains very-long-chain fatty acids, EPA, and DHA [16]; green tea contains polyphenol catechins with antioxidative activity [17, 18]; and seaweeds contains dietary fiber [19].

Although the impact of individual components of the JDI8 score upon mortality risk was small, the impact of the JDI8 score as a whole (i.e., the Japanese dietary pattern) was stronger. This finding may be explained by the cumulative effects of individual components of the Japanese diet. A previous study suggested that owing to the complex biological interactions between different components of the diet, the use of a whole diet approach rather than individual nutrients or food groups might help to understand the role of diet in health outcomes [20].

In this study, participants with higher JDI8 score consumed sodium, which increases the risk of gastric cancer [21, 22] and hypertension [23, 24], and rice that raises glycemic indices, which increases the risk of colorectal cancer [21] and cardiometabolic complications [25]. On the other hand, participants with higher JDI8 score may also consume such beneficial nutrients as protein, fiber, vitamins A, C, and E, calcium, iron, potassium, and magnesium which have been associated with lower risk for CVD, including heart disease and stroke [26,27,28,29,30,31,32,33,34], because their intake of all adhering components of JDI8 are high.

Additionally, the present study indicated that intake of potassium as well as sodium increased among the subjects with higher JDI8 score, and the sodium–potassium ratio was almost the same among all the groups. It was reported that sodium–potassium ratio predicted the risk for cardiovascular events better than sodium intake alone [35]. Evidence also indicates that high potassium intake can mitigate the risk for high blood pressure despite a high-sodium diet [36, 37]. Therefore, the inverse association of the Japanese diet assessed by the JDI8 with mortality may be partly explained by the combined effect of these nutrients.

On the other hand, we did not observe a significant association between the JDI8 score and cancer mortality. According to the World Cancer Research Fund [21] and the National Cancer Center of Japan [22], vegetable intake was associated with a lower risk for stomach and esophagus cancer; fish intake was associated with a lower risk for cervical cancer; pickles intake was associated with a higher risk for stomach cancer; red meat and processed meat intake were associated with a higher risk for colorectal cancer; and dietary fiber and calcium intake were associated with a lower risk for colorectal cancer. These findings of cancer incidence imply that the association between foods or nutrients and cancer differ by site of cancer. However, less is known about the association between diet and cancer mortality. Other factors, such as stage at diagnosis or treatment, are important contributors to cancer mortality outcomes and were not adjusted for in these analyses. Thus, further investigation of the association between the Japanese diet and cause-specific cancer is required.

Our study had several strengths: (1) it was a large population-based cohort study including 92,969 persons; (2) the response rate was relatively high (79.2%); (3) the follow-up period was long enough for us to eliminate possible biases by conducting analysis excluding participants who died in the first 3 years of follow-up; (4) few participants were lost during follow-up (0.3%) and an almost complete follow-up for cause of death (in decedents) was possible among all subjects (over 98%); (5) we used the validated semi-quantitative FFQ; and (6) many confounding factors were taken into account.

However, there were also several limitations. First, because dietary intake was assessed at only one time point, changes in dietary habits were not considered. Second, although we adjusted for many potential confounders, residual confounding such as socioeconomic status [38,39,40] could not be completely ruled out.

In conclusion, the present study indicated that the Japanese diet assessed by the JDI8 is associated with a decreased risk of all-cause, CVD, and heart disease mortality among Japanese.