Study design and participants
This cross-sectional study was based on data from a baseline survey of Hitachi Health Study II. The Hitachi Health Study II is an ongoing prospective study conducted on current and retired employees and their spouses at Hitachi, Ltd., Japan, a manufacturing company in the Ibaraki Prefecture. Data collection for the baseline survey was conducted as part of the health check-up from April 2017 to March 2020. Of those who underwent health check-ups, individuals aged 60 years and over (as of 31st of March each year) and screened for their cognitive function were instructed to participate in the baseline survey. We further instructed participants aged 60, 63, 66, or 69 years to fill out two questionnaires for overall health-related lifestyle and dietary habits on the day of the health check-up. Because of the aim of the present study, we considered those who underwent cognitive function screening and questionnaire surveys as the target of the present analysis (n = 1,581). No women were included because of the small proportion of the original sample (10.4%). For the present analysis, after excluding those aged 70 years and over (n = 7), we included participants who completed both questionnaire surveys (n = 1,551). We further excluded those with a medical history of stroke (n = 25) and those with no information on the variables of interest (n = 32). In total, 1,494 men were included in the final analysis (Figure 1).
This study was conducted in accordance with the guidelines of the Declaration of Helsinki. All procedures involving human subjects were approved by the Ethics Committee of the National Center for Global Health and Medicine (approval number: NCGM-G-002208) and Hitachi Health Care Center. Written informed consent was obtained from all participants before participating in the present study.
Estimation of dietary hardness
Dietary hardness in the present study was defined as an estimate of the masticatory muscle activity required for the consumption of the habitual diet, according to previous studies [10–13]. Habitual dietary intake was assessed using a previously validated, brief-type, self-administered diet history questionnaire (BDHQ) [14,15]. The BDHQ estimates daily intake of 58 items (consisting of 38 foods, 12 beverages, and 8 seasonings) during the preceding month, which were commonly consumed in Japan [14]. The BDHQ is a fixed-portion size questionnaire that asks about the frequency of selected food consumption, but not portion size. Intakes of energy and selected nutrients were calculated using an ad hoc computer algorithm for the BDHQ [15], according to the Standard Tables of Food Composition in Japan, 2010 [16]. The BDHQ has a satisfactory ranking ability of energy-adjusted dietary intake using the density method in Japanese men aged 32–76 years [14,15]. Median Spearman’s correlation coefficient with 16-day dietary records was 0.48 (interquartile range, 0.33–0.56) for food groups [14], while Pearson’s correlation coefficient was 0.56 (interquartile range, 0.41–0.63) for nutrients [15].
Subsequently, the dietary hardness of each participant’s habitual dietary intake (mV·s/day) was calculated by summing the products of the hardness of each food item in the BDHQ (mV·s/cm3) and its volume consumed (cm3/day) [10,11]. Briefly, to estimate the hardness of each food item (mV·s/cm3), 34 out of 38 food items in the BDHQ were directly matched to an equivalent food item for which information on masticatory muscle activity (mV·s/2.197 cm3) was available from Yanagisawa et al. [13] and then divided by 2.197 [10]. The corresponding values for similar food items were used as proxies for the remaining four food items, whereas the hardness of beverages and seasonings was not estimated [10]. Given the great influence of cooking methods on the hardness of vegetables [13], the observed ratio of consumption of raw/cooked form (S. Sasaki, unpublished observations, 2006) was considered as much as possible in the matching procedure [10,11]. Food volume consumed (cm3/day) was estimated based on the weight in grams (g/day) assessed using the BDHQ by assuming that the density for all foods was 1 (g/cm3) [10,11].
To consider differences in dietary intake due to varying body sizes and energy requirements and to attenuate the influence of misreporting and a high correlation between energy intake and the crude estimate of dietary hardness, dietary hardness and dietary intake were energy-adjusted using the density methods [17]. Because of the higher proportion of alcohol consumers in the present study (75.2%) than in a previous study (37.0%) [10] and the non-contribution of fluids, including alcoholic beverages, to the estimation of dietary hardness, energy-adjusted dietary hardness was provided as the value per 1000 kcal of energy intake from solid foods (i.e. foods and seasonings). Energy-adjusted nutrient intake is presented as units/1000 kcal of total energy intake.
Assessment of cognitive dysfunction
Cognitive dysfunction was assessed using a computerised test battery for screening individuals at risk of Alzheimer’s disease (MSP-1100, Nihon Kohden Corporation, Tokyo, Japan) [18]. Briefly, the test battery was developed based on the revised version of Hasegawa’s Dementia Scale [19]. It consisted of 4 tasks for examining temporal memory (3 items), temporal orientation (4 items), three-dimensional visual-spatial perception (2 items), and short-term memory (3 items) of participants [18]. Each item was scored as 1 (for the former 3 tasks) or 2 (for the latter 1 task) points for each correct response. The score ranged from 0 to 15 points, with a higher score indicating cognitive improvement. According to previous studies [18,20], participants who scored ≤13 points on the test battery were defined as having cognitive dysfunction.
Assessment of covariates
For the assessment of covariates, we referred to the participants’ health check-up data, including anthropometric and biochemical measurements and information on the history of diseases, and the overall health-related lifestyle questionnaire. Body height and weight were measured to the nearest 0.1 kg and 0.1 cm, respectively, while the participants wore light clothes and no shoes. Body mass index (BMI, kg/m2) was calculated as body weight divided by height squared. Blood pressure was measured using an automatic sphygmomanometer. Fasting plasma glucose (FPG) level was measured using the glucose oxidase enzyme electrode method (A&T, Tokyo, Japan), and haemoglobin A1c (HbA1c) level was measured using high-performance liquid chromatography (HLC723-G9, TOSOH, Tokyo, Japan). Alcohol consumption was calculated based on information on the frequency and amount of alcohol consumption collected during the health check-up. Hypertension was defined as present when participants had systolic blood pressure ≥ 140 mmHg, diastolic blood pressure ≥ 90 mmHg, and/or self-reported medication use for hypertension [21]. Diabetes was defined as present when participants had FPG level ≥ 126 mg/dL, HbA1c level ≥ 6.5%, and/or self-reported medication use for diabetes [22]. Depressive symptoms were measured using the Japanese version of the short version of the Center for Epidemiologic Studies Depression Scale, which consists of 11 of the original 20 items [23–25]. In line with a previous study that used arithmetic conversion to define a cut-off score, those who scored ≥ 9 points were defined as having depressive symptoms [26]. Dietary reporting status was evaluated based on the ratio of the reported energy intake to basal metabolic rate (BMR), using the Goldberg cut-off for 30-day (i.e. 1 month) dietary data and a physical activity level for sedentary lifestyle (i.e. 1.55) [27]. BMR was estimated using sex- and age-specific equations developed for the Japanese population [28,29].
Statistical analyses
Descriptive data are presented as means and standard deviations (SDs) for continuous variables or numbers and percentages of participants for categorical variables. Energy-adjusted dietary hardness (mV·s/1000 kcal) was categorised into tertiles and used to compare the selected characteristics of the participants.
Odds ratios (ORs) and 95% confidence intervals (CIs) for cognitive dysfunction were estimated for each tertile of dietary hardness by logistic regression analysis, using the lowest category as the reference. Three models were considered in the analysis. Model 1 was adjusted for age (years, continuous variables). In Model 2, we further adjusted for the following potential confounding factors: education (< 10, 10–12, or ≥ 13 years), current employment (yes or no), living alone (yes or no), smoking status (current or past/non-smoking), alcohol consumption status (none, > 0 to < 46, or ≥ 46 g/day), habitual exercise (yes or no), BMI (kg/m2, continuous), hypertension (yes or no), diabetes (yes or no), depressive symptoms (yes or no), dietary counselling from a doctor or dietitian (yes or no), dietary reporting status (under-, plausible, or overreporting), and energy intake (kcal, continuous) [4,5]. In Model 3, we further adjusted for intake of nutrients, including n-3 PUFA; vitamins A, D, E, B6, B12,and C; and folate (unit/1000 kcal, continuous) [5,30], which may reduce the risk of cognitive impairment, to consider whether the observed association would be independent of nutrient intake. We tested linear trends using dietary hardness as a continuous variable.
In the sensitivity analysis, we classified participants into halves or quartiles according to dietary hardness or used ≤ 12 points on the MSP-1100 score as a definition of cognitive dysfunction based on instructions. We repeated the same analysis based on energy-adjusted dietary hardness using total energy intake. Furthermore, we investigated the contribution of each component score of the MSP-1100 (i.e. temporal memory, temporal orientation memory, three-dimensional visual-spatial perception, and short-term memory) based on ordered logistic regression analysis using the highest scores of each component as a reference. The results are presented in terms of beta coefficients and 95% CI. All statistical analyses were performed using SAS version 9.4 (SAS Institute Inc., Cary, NC, USA). All reported P-values were two-tailed, and statistical significance was set at P < 0.05.