Introduction

Aging traditionally involves a decline in executive functions, which in some cases progresses towards cognitive impairment. Given that the world’s population is getting older, the prevalence of subjective cognitive decline and mild cognitive impairment (MCI) is on the rise1,2. Therefore, the development of strategies to maintain normal cognitive function and delay cognitive decline is relevant.

Cardiovascular risk factors like type 2 diabetes (T2DM) or obesity, have been consistently associated with cognitive decline and MCI2,3, as well as the progression to dementia4. Thus, it has been postulated that therapeutic management of these factors can be useful to maintain cognitive function and even slow down cognitive decline2. Therapies based on lifestyle interventions, specifically those designed to prevent cardiovascular disease, could be a safe approach to maintain cognitive function and delay aging-associated decline. Although approaches using dietary supplements did not report significant benefits5,6, studies with traditional Mediterranean Diet (MedDiet) have shown promising results7,8,9,10. MedDiet has also been reported to be associated with better self-reported quality of life11.

Physical activity (PA) is another modifiable lifestyle factor with potentially large beneficial effects on cognitive functioning. It has been postulated that exercise promotes healthy aging and has beneficial effects on cognition in the elderly12. Several cross-sectional analyses and interventional studies have found positive associations between PA, especially aerobic exercise, and cognitive function in the elderly13,14,15,16. These studies strongly suggest that the PA association with cognitive function is relevant enough to be taken into consideration when designing lifestyle interventions addressed to improve health in the elderly. However, reported associations between PA and cognitive health are often small to moderate, depending on the cognitive domain, and inconsistent across studies15,16,17,18,19,20,21,22,23. PA has also been associated with a better subjective quality of life24,25.

Although there is compelling evidence pointing towards an association between PA and cognitive function, the heterogeneity of the findings calls for additional research to unravel those conditions where the effect is beneficial. These inconsistencies can be due to the different baseline health status of participants, environmental factors, tools used to assess cognition and PA, different PA programs, or statistical models. In this regard, physical fitness (PF) needs to be considered. PF is a measure of physical condition and well-being, and it has been reported that PA with an effective improvement in cardiovascular fitness has a greater effect on cognitive function than PA focused on strength26. We aimed to investigate the relationship between PA and PF with cognitive function and quality of life, independently of the impact of the adherence to an energy-restricted MedDiet, using baseline data from the PREDIMED-Plus trial27.

Results

Association between PF and cognitive function

Participants in the higher physical PF quartile (Q4) were mainly younger males (Table 1). Higher PF was associated with lower BMI (p < 0.001), longer sleeping time at weekends (p = 0.027), lengthier schooling time (p < 0.001) and higher education (p < 0.001). Former smokers were also more frequently classified in the Q4 (p < 0.001) as well as married participants (p < 0.001) and active workers (p < 0.001) (Table 1). Moreover, participants in the highest PF quartile had a lower prevalence of diabetes (p < 0.001), depression (p < 0.001) and hypertension (p = 0.002) (Table 1).

Table 1 Baseline characteristics of the participants according to baseline fitness and PA levels.

We found an association between the language domain of the MMSE test and the total MMSE score and PF (p = 0.024 in both cases) (Table 2). Both were higher in the upper PF quartile, although the increase in the score was below 1%. However, when we considered the PF score as a continuous variable, we lost significance after complete adjustment (Supplementary Table S1). We also found a higher calculation score associated with higher PF, although significance was lost after complete adjustment. Similar results were found with the PF score as a continuous variable (Supplementary Table S1). We found that better PF was associated with a higher score in phonemic (4.5%, 5.57% and 7.54% higher score for Q2, Q3 and Q4 vs. Q1, respectively) and semantic (3.24%, 5.86% and 6.53% higher score for Q2, Q3 and Q4 vs. Q1, respectively) verbal fluency tests (p < 0.001) (Table 2). Results were confirmed when we treated PF as a continuous variable (β = 0.0434 (0.0228–0.064) p = 0.001 for phonemic and 0.0428 (0.0203–0.0653) p = 0.003 for semantic fluency tests respectively) (Supplementary Table S1). TMT test time performance was also lower in higher PF quartiles. Specifically, TMT A time was 5.94%, 9.25%, and 10.32% lower and TMT B time was 6.14%, 10.04% and 11.05% lower in Q2, Q3 and Q4, respectively (p < 0.001 in both cases). When we analyzed PF as a continuous variable, we confirmed the lower TMT A time with higher PF score (β = −0.23 (−0.36–0.094), p = 0.0012), but we did not confirm the association between PF score and TMT B time (Supplementary Table S1). As sensitivity analyses, we excluded those with very high values of chair test (>2 SD, n = 237) as may represent implausible values/outliers and we lost the significant association with the language domain but not with total MMSE test score (p = 0.006), while association with TMT A and phonemic and semantic fluency tests remained significant, either analyzing PF quartiles or continuous scores (data not shown).

Table 2 Comparison of the cognitive assessment variables of the participants according to their physical fitness (quartiles).

Association between PA levels and energy expenditure and cognitive function

Participants with the highest level of PA were mainly men (p < 0.001), former smokers (p < 0.001), with high adherence to MedDiet (p < 0.001), slept longer time in weekdays (p = 0.014) and showed lower weight (p < 0.001) and BMI (p < 0.001) (Table 1). Married participants were more frequently classified as high PA level, whereas single, divorced or widow participants were more frequently classified as light PA level (p = 0.018). Additionally, current workers showed a lower PA level than retired participants (p < 0.001). Finally, the prevalence of diabetes and depression was lower among participants with a high PA level (p = 0.017 and p = 003, respectively) (Table 1). We found a positive association between higher PA and higher scores in phonemic (2.39% and 5.45% for moderate and high PA level vs. low PA level, p < 0.001) and semantic fluency tests (0.25% and 3.88% for moderate and high PA level vs. low PA level, p = 0.003) under the minimally adjusted model (Table 3). However, the significance was lost after complete adjustment. We did not find any association between the level of PA and the score in any of the neurocognitive tests in the complete adjustment model (Table 3). Similarly, we found higher scores in the direct numeric series associated with higher PA level under the minimally adjusted model (1.73% and 3.46% for moderate and high PA levels compared with low PA level, p = 0.032), but significance was lost after complete adjustment (Table 3). Similar results were observed when analyzing the association between scoring in neurocognitive tests and energy expenditure with total PA (Supplementary Table S1).

Table 3 Comparison of the cognitive assessment variables of the participants according to their PA level.

Association between both, PF and PA and health-related quality of life

Higher levels of PA and higher energy expenditure with total PA were associated with higher scores in all domains of the SF36-HRQL questionnaire (Table 4 and Supplementary Table S2). Accordingly, individuals in the upper quartile of PF or with higher Chair Stand repetitions showed higher scores in all domains of the HRQL-SF36 questionnaire, including the aggregated physical and mental health domains (Table 5 and Supplementary Table S2). We repeated the analysis excluding those with very high values of chair test (>2 SD, n = 237) and results did not change (data not shown).

Table 4 Adjusted means for HRQL dimensions according to baseline PA levels.
Table 5 Adjusted means for HRQL dimensions according to the baseline level of PF (quartiles).

Discussion

Our study showed that PF was positively associated with a better performance in the TMT, measuring executive function, and the COWAT, measuring language ability and executive control. These results are in accordance with those reported by Alosco et al. in 200 adult patients of heart failure using a similar battery of cognitive and PF tests28. The authors suggested that the negative effect of obesity on cognitive function was likely a result of the combination of clinical factors (obesity, diabetes, and hypertension) and PF28. Because we included BMI, diabetes, and hypertension as adjusting co-variables, our results suggest that PF influences cognitive function independently of these other factors.

We did not find any association between the level of PA and cognitive function. The association between PA and cognitive function is still controversial. Some previous reports have suggested a positive association. For instance, a study conducted in 2,315 cognitively healthy elder participants from CFAS-Wales study showed that PA was positively associated with cognitive function measured by the Cambridge Cognitive Examination battery of neuropsychological tests29 and a larger observational study carried out in 93,082 respondents ≥45 years old from 21 US states reported that people with self-reported cognitive decline were more frequently inactive than those without this subjective cognitive decline status30. Moreover, the study conducted by Xu et al. in 27,651 participants (50–85 years old) showed a dose-response relationship between PA and cognitive function14. Other population studies and meta-analyses support an association between PA and cognitive function20,22,31. Conversely, Kooistra et al. did not find any association between PA and cognitive performance or cognitive decline in a baseline analysis of 1,232 participants of the SMART-RM study16. These inconsistencies could be related to differences in populations and tools used to measure cognitive function and PA. The intensity of the PA could also be a determinant factor. The OSHPE study found that moderate activity, but not light activity, was associated with hippocampal volume but not with memory15 and a recent systematic review concluded that moderate and vigorous leisure-time PA was associated with benefits in different cognitive domains in older adults17. Therefore, it is possible that an association between PA and cognitive function was significant only in individuals at a high PA level. Indeed, our results show a trend towards a better scoring in TMT and COWAT test in participants categorized as high-level PA.

The discordant results could also be attributed to the type of PA, aerobic, or resistance/strength exercises, with aerobic exercises likely having a greater effect. Protocols for each type of exercise are different in intensity and length and produce different physical adaptations, with aerobic exercise having a greater impact on %VO2max and energy expenditure32,33,34. We found that, although higher PA and PF levels were both associated with better self-reported quality of life, only PF was associated with cognitive function. This differential effect could be attributed to a different impact of aerobic vs. anaerobic exercise on cognitive function. Aerobic fitness has been positively correlated with better high-interference memory in older adults18, and aerobic exercises have been associated with better improvement in the score of neuropsychological tests in comparison with anaerobic exercise35. Aerobic exercise in the elderly has also been associated with better performance on tasks requiring executive control in comparison with anaerobic exercises36. Other reports also showed a specific association between aerobic exercise, but not anaerobic exercise, and executive control37,38.

Moreover, a meta-analysis published in 2003 found robust, but specific, benefits of aerobic fitness on cognitive function19. Another meta-analysis including 29 randomized interventional trials involving aerobic exercise training (n = 2049) found a moderate but consistent improvement in attention and processing speed, executive function and memory associated with aerobic exercise. However, it failed to find any association with working memory21. These studies strongly support the hypothesis that aerobic exercise, but not anaerobic one, has a specific impact on executive function. We did not discriminate between aerobic and anaerobic exercise, neither measured cardiorespiratory fitness; however, we found a significant association between TMT and COWAT tests and PF. TMT test measures the speed of information processing and the COWAT test measures verbal memory, attention, and inhibition. Together, they define the individual condition regarding executive function and cognitive control. Therefore, our results strongly suggest that good PF, achieved through an aerobic exercise program, could be a significant predictor of executive function in the elderly. The mechanisms behind these associations could be the increase in the hippocampal volume associated with aerobic fitness39. In this regard, a randomized clinical trial conducted in 120 adults without dementia showed an increase in the volume of the left and right hippocampus in the aerobic exercise group and a decrease in the stretching control group after one year34.

Interestingly, authors showed a selective impact of aerobic exercise on the volume of the anterior hippocampus, involved in spatial memory and submitted to age-related atrophy34. However, a Cochrane systematic review showed no effect of aerobic exercise, even if the intervention leads to an improvement in cardiorespiratory fitness, on cognitive domains23. More extensive interventional studies are, therefore, needed to clarify the impact of aerobic exercise leading to an improved PF on cognitive function.

We also found that PF was associated with a better quality of life, with positive associations in all domains, including the aggregated ones. PF has been reported to account for up to 22.5% of the variability in the quality of life scoring of obese individuals with binge eating disorder40, and cardiorespiratory fitness has been suggested to be a good predictor of physical and mental quality of life in individuals with bipolar disorder41. Our results reinforce the hypothesis that PF is associated with a better cognitive function and quality of life and encourage the development of PA programs addressed to increase PF in the late adulthood to prevent morbidity in the elderly. Although PA levels seemed not to be associated with cognitive function, they were associated with a higher score in all domains of the quality of life questionnaire. These results are consistent with the previously reported association between PA and quality of life in the elderly42,43.

In conclusion, our results showed that PF, but not PA level, was associated with cognitive function, specifically with those domains related to language ability, attention, and processing speed. Undoubtedly, PA, specifically aerobic exercises, contribute to PF and thus PA could be associated with cognitive function through an improvement in PF. Our results highlight the importance of maintaining a good physical condition to retain cognitive function. In this regard, the primary aim of PA programs addressed to maintain cognitive health should be the improvement of PF, and the type of exercise should be considered when designing PA programs.

This study has some limitations. First, this is an observational cross-sectional study that does not allow to infer causality. Second, this is an elderly population with obesity and metabolic syndrome with an average low level of PA and without cognitive impairment and this could be masking a significant relationship between PA and cognitive function. However, this population is representative of a high cardiovascualr disease (CVD) risk elderly population in a medium-high income Western country and this study has been carried out with a large sample and with a stringent adjustment for confounding factors. Third, HRQL and PA activities are self-reported; thus, we cannot rule out the overestimation or underestimation of those variables. Although we have adjusted for potential confounding factors (age, gender, recruitment center, adherence to an er-MedDiet, BMI, civil status, educational level, sleep hours, working status, TV hours, prevalence of T2DM, depression, hypertension and cancer) our population is heterogeneous, especially regarding sociodemographic parameters like civil status, education level and working status. Thus, we cannot rule out the possibility that this heterogeneity produced a moderation or exacerbation of the observed association between PF and cognitive function.

Finally, we should bear in mind that past experiences or the accumulation of other factors, not considered within this analysis, may also have an impact on cognitive function in the elderly.

Methods

Study design and participants

We used baseline data of the PREDIMED-Plus trial, a 6-year multicenter, randomized, parallel-group clinical trial conducted in Spain for primary prevention of CVD through lifestyle modification. The cohort and protocol, including sample size calculation and recruitment flowchart, have been described27, and the protocol can be found at the PREDIMED-Plus website44. This trial was retrospectively registered at the International Standard Randomized Controlled Trial (ISRCTN89898870, https://www.isrctn.com/ISRCTN89898870?q=ISRCTN89898870&filters=&sort=&offset=1&totalResults=1&page=1&pageSize=10&searchType=basic-search) on 07/24/2014. Participants were recruited by primary care health doctors and nurses associated with study centers who assessed the inclusion/exclusion criteria. Eligible participants were community-dwelling adults (55–75 years men and 60–75 years women) with overweight/obesity (27 ≤ BMI ≥ 40 Kg/m2) who met at least three criteria of metabolic syndrome (according to the International Diabetes Federation and the American Heart Association)45 without cardiovascular and neurodegenerative disease at baseline. We included baseline data from 6874 participants recruited in 23 Spanish study centers from October 2013 to December 2016. Data include cognitive function, PA, quality of life, anthropometrical and biochemical measures, PF conditions, history of cancer or lung, disease sociodemographic parameters, and depression status27.

The study protocol and procedures were approved according to the ethical standards of the Declaration of Helsinki by the Institutional Review Boards (IRBs) of all the participating institutions: CEI Provincial de Málaga, CEI de los Hospitales Universitarios Virgen Macarena y Virgen del Rocío, CEI de la Universidad de Navarra, CEI de las Illes Balears, CEIC del Hospital Clínic de Barcelona, CEIC del Parc de Salut Mar, CEIC del Hospital Universitari Sant Joan de Reus, CEI del Hospital Universitario San Cecilio, CEIC de la Fundación Jiménez Díaz, CEIC Euskadi, CEI en Humanos de la Universidad de Valencia, CEIC del Hospital Universitario de Gran Canaria Doctor Negrín, CEIC del Hospital Universitario de Bellvitge, CEIC de IMDEA Alimentación, CEIC del Hospital Clínico San Carlos, CEI Provincial de Málaga, CCEIBA de la Investigación Biomédica de Andalucía, CEIC del Hospital General Universitario de Elche, Comité de Ética del Hospital Universitario Reina Sofía and CEIC de León. All participants provided informed written consent.

Dietary and anthropometrical measures

Adherence to an er-MedDiet diet was measured through a 17-item questionnaire in which an adequate consumption of typical traditional Mediterranean foods adds add one point, and low consumption of foods that are not characteristic of the traditional Mediterranean diet also adds one point. This questionnaire has been specially developed for the PREDIMED-Plus trial11 and was applied by a trained nutritionist in a face-to-face interview. Weight and height were measured with light clothing and no shoes, with calibrated scales and a wall-mounted stadiometer, respectively. BMI was calculated as weight in kilograms divided by the square of height in meters.

Physical activity and fitness

PF was assessed using the 30 seconds Chair Stand test designed to test the functional fitness of senior population46. Participants were classified according to quartiles of test scoring: Q1 (0–11 repetitions), Q2 (12–13 repetitions), Q3 (14–16 repetitions) and Q4 (17–37 repetitions).

PA was assessed using the validated REGICOR questionnaire and the Rapid Assessment of Physical Activity Questionnaire (RAPA)47. REGICOR questionnaire provides information about PA intensity according to energy expenditure (MET/min/week) distinguishing light (<4 MET), moderate (4–5.5 MET) and vigorous (≥6 MET) PA and recording information about the PA at work and in everyday life and about sedentary behaviors42. Frequency and intensity of PA are quantified by requesting participants to provide information about the frequency (days per month) and the intensity (minutes per session) of physical activities like light walking, brisk walking, gardening, trail/hiking, climbing stairs, other indoor or outdoor sports activities. RAPA questionnaire allows an easy identification of PA level as low (light-moderate PA activities but not every week), moderate (moderate PA activities less than 150 minutes per week or vigorous PA activities less than 75 minutes per week) or high (more than 150 minutes per week of moderate PA activities or more than 75 minutes per week of vigorous activities). Trained interviewers administered PA questionnaires in individual face-to-face sessions.

Neuropsychological assessment

Executive function was measured with the Trail-making Test (TMT) A and B3,48. TMT-A provides information about cognitive processing skills (psychomotor speed, processing speed, and visuospatial skills), whereas TMT-B provides information about executive functioning49. In task A, participants are requested to connect circled numbers 1–25 in the correct order. In task B, numbers and letters have to be alternatively connected following an ascending order in each case (1-A-2-B-3-C-…).

Semantic and phonemic verbal ability and executive control were measured with the Controlled Oral Word Association Test (COWAT)3,50. Semantic fluency is evaluated by requesting participants to mention as many animals as possible in 1 minute. Phonemic fluency is evaluated by requesting participants to say as many words starting with “P” as they can in 1 minute50.

Short-term memory and attention and working memory were measured by the forward and backward Digit Span (DS) test of the Wechsler Adult Intelligence Scale-III battery3,51, respectively. In the DS-Forward recall, participants are requested to repeat a series of random single digits in the same order they heard them. In the DS-Backward recall they are requested to repeat another series of a random single digit in the inverse order they heard them.

Cognitive impairment was assessed with the Mini-mental State Examination test (MMSE) and the Clock-Drawing test (CDT)52,53,54. The MMSE test is a 30-point questionnaire that examines domains like time and spatial orientation, immediate and deferred recall, attention, calculation, and language. A score is calculated for each domain, and all are added to obtain a total MMSE score. A threshold of 24 is defined to identify cognitive impairment. The CDT is useful to examine visuospatial functioning and memory. Participants are required to draw a clock with a specific time, and a numeric score of 0 or 1 is given depending if they correctly or incorrectly draw the numbers, minute and seconds clock hands and the requested time49. All scores are sum to get a total CDT score.

All cognitive tests were applied by nutritionists who had been previously trained.

Health-related quality of life (HRQL)

We used an adapted version of previously published 36-items HRQL questionnaire (SF36-HRQL)55, which was validated for the Spanish population and extensively used to measure the subjective perception of health and capability or physical limitations to manage with daily tasks56, as described before11,57.

Statistical analyses

Baseline characteristics of the participants were described as means ± SD for quantitative variables and as proportions for qualitative variables. Chair Stand test scores and MET/Min/week of total PA were used to quantitatively measure PF and PA, respectively, and were analyzed as continuous variables. Additionally, participants were classified according to the quartile of PF measured by the Chair Stand test.: Q1 (0–11 repetitions), Q2 (12–13 repetitions), Q3 (14–16 repetitions) and Q4 (17–37 repetitions). PA activity level was also analyzed as a categorical variable: low, moderately, and highly active. The participant’s categorization has been described above.

Characteristics of the participants were compared among PF quartiles and PA levels with the χ2 test for categorical variables and ANOVA for continuous variables.

The score for each domain of the SF36-HRQL tests as well as for the aggregated dimensions was treated as a continuous variable and compared with PA and PF with ANCOVA. Neurocognitive scores were also treated as continuous variables and compared with PA levels and PF with ANCOVA.

SF36-HRQL and neurocognitive scores (independent variables) were analyzed for their association with Chair Stand test score and MET/min/week energy expenditure as continuous variables (dependent variables) with a linear regression model.

Two adjustment models were applied. Model 1 adjustment (minimal adjustment) included age, gender, and recruitment center as cofactors. Model 2 adjustment (complete adjustment) included, additionally, smoking status (current, never or former smoker), adherence to er-MedDiet11, BMI, civil status (married, single, divorced or widow), educational level (primary, secondary or college), sleep hours, working status (active, long-term inactive or retired), TV hours, prevalence of T2DM3, depression, hypertension and cancer. We included data from 6874 participants. Missing values for predictive and descriptive variables were imputed using Predictive Mean Matching. The specific sample size for outcome variables is detailed in corresponding tables. The percentage of missing values range from 0.029 to 9%, except for the DS test and the aggregated physical and mental dimensions of the SF36-HRQL test, for whose percentage of missing values ranged from 20.99% to 67.99%. To prevent from type-I error inflation, a multiple testing correction with the Holm’s method was applied58. Holm’s method is a sequential hypothesis rejection test that controls the family-wise error rate like Bonferroni’s one but with much increased power. As sensitivity analyses, we excluded participants with >2 SD values of the Chair Stand Test (n = 237) as may represent implausible values/outliers and repeated ANCOVA analyses. We considered the two-tailed significant level of 0.05 as a threshold for statistical significance. All analyses were conducted with data from database version 201903131142_PREDIMEDplus_2019-03-12 and with R.