Multidimensional associations between nutrient intake and healthy ageing in humans
BMC Biology volume 20, Article number: 196 (2022)
Little is known about how normal variation in dietary patterns in humans affects the ageing process. To date, most analyses of the problem have used a unidimensional paradigm, being concerned with the effects of a single nutrient on a single outcome. Perhaps then, our ability to understand the problem has been complicated by the fact that both nutrition and the physiology of ageing are highly complex and multidimensional, involving a high number of functional interactions. Here we apply the multidimensional geometric framework for nutrition to data on biological ageing from 1560 older adults followed over four years to assess on a large-scale how nutrient intake associates with the ageing process.
Ageing and age-related loss of homeostasis (physiological dysregulation) were quantified via the integration of blood biomarkers. The effects of diet were modelled using the geometric framework for nutrition, applied to macronutrients and 19 micronutrients/nutrient subclasses. We observed four broad patterns: (1) The optimal level of nutrient intake was dependent on the ageing metric used. Elevated protein intake improved/depressed some ageing parameters, whereas elevated carbohydrate levels improved/depressed others; (2) There were non-linearities where intermediate levels of nutrients performed well for many outcomes (i.e. arguing against a simple more/less is better perspective); (3) There is broad tolerance for nutrient intake patterns that don’t deviate too much from norms (‘homeostatic plateaus’). (4) Optimal levels of one nutrient often depend on levels of another (e.g. vitamin E and vitamin C). Simpler linear/univariate analytical approaches are insufficient to capture such associations. We present an interactive tool to explore the results in the high-dimensional nutritional space.
Using multidimensional modelling techniques to test the effects of nutrient intake on physiological dysregulation in an aged population, we identified key patterns of specific nutrients associated with minimal biological ageing. Our approach presents a roadmap for future studies to explore the full complexity of the nutrition-ageing landscape.
How does what we eat affect our healthspan and longevity? The answer to this relatively concise question is unavoidably complex. Conventional approaches to understanding the effects of diet on health and ageing, particularly in human nutrition, have usually focussed on single nutrients or a handful of dietary attributes/patterns [1,2,3]. Yet, nutrients have both individual and interactive effects. For example, at the most macro-level, protein, carbohydrate and fat energy sources interact to determine metabolic, physiological and cognitive functioning (e.g. [4,5,6]). Similarly, the numerous phenotypic changes that occur with age are increasingly recognised as interconnected and multidimensional [7,8,9,10]. Thus seemingly distinct physiological components of ageing likely reflect a broader loss of homeostasis in a complex dynamic system rather than independent processes . Such interdependencies mean that the atomised interpretation of the effects of a single nutrient, diet, molecular mechanism or biomarker is likely to be context-dependent [12,13,14,15]; a consequence being that the results of univariate studies are spurious and/or more difficult to reproduce, leading to inconsistent conclusions between studies.
The Geometric Framework for Nutrition (GFN) is a state-space approach to nutrition that deals with dietary complexity by considering multiple dimensions of nutrient intake simultaneously (Fig. 1) [16, 17]. Using fully factorial experiments, the GFN has shown that components of biological ageing are affected by the ratio of dietary macronutrients independently of the effects of net energy intake. These effects have been observed across taxa [18,19,20]. For example, mice subjected to a lifetime of low-protein (within the boundaries of what can support growth and development) high-carbohydrate intake display improved cardiometabolic health and increased median age at death relative to animals with higher protein or fat intakes [4, 21]. Any such benefits seemingly disappear in old age, though, when higher protein reduces late-life mortality . These experimental findings in mice are consistent with epidemiological studies in humans that emphasise the importance of dietary protein for the elderly (reviewed in [22, 23]). Ecological analyses of national-level data in humans have also found that even at this highest order level, the macronutrient composition of the food supply is a powerful predictor of international variation in patterns of mortality . Collectively, these studies emphasise the importance of multi-dimensional thinking in nutrition.
The biological ageing process is no more tractable than nutrition. There is no clear consensus as to what ageing is , though most researchers now agree it is multi-factorial [26, 27]. Different methods for quantifying ageing correlate poorly with one another after chronological age adjustment , implying that ageing is a compound process. An emerging approach is to measure the effects of ageing via the breakdown in homeostatic regulation (i.e. dysregulation) across physiological systems [10, 29], an approach complementary to ‘biological age’. A statistical distance can be used to quantify how abnormal, or ‘dysregulated’, an individual’s biomarker profile is, either globally or within specific systems. In this way, it is possible to generate dysregulation scores that are predictive of a wide array of health outcomes during ageing, including mortality, frailty, and chronic diseases [10, 30,31,32,33]. Dysregulation scores thus simultaneously provide a proxy for the ageing process and general health state. In this sense, dysregulation may be a better metric of health than measures of the ageing process per se, if such a thing actually exists.
Here, we apply the GFN to model the effects of nutrient intake on dysregulation and ageing scores in community-dwelling older adults (>67 years old). We use data from the Quebec Longitudinal Study on Nutrition and Successful Aging (NuAge) , which is extensive enough to permit the application of the GFN in an epidemiological context. We show that combining these two methods provides a means to integrate the complexity inherent to nutrition and ageing physiology. We begin by testing the hypothesis, suggested by the mouse data, that higher protein intakes during old age are associated with markers of improved health. Simultaneously, we also test whether protein interacts with the intake of the other macronutrients, notably whether an increase in the ratio of carbohydrate to protein reduces markers of ageing. We then show how the same approach can be used inductively to detect micronutrient interactions that have systemic effects.
Briefly, NuAge participants were community-dwelling men and women, aged 67-84 years in the Montreal, Laval, or Sherbrooke areas in Quebec (Canada). They were selected randomly from the Quebec Medicare database (n=36,183), after stratification for age and sex. Individuals with good general health were recruited (n=1793) between November 2003 and June 2005 (T1) . Participants were re-examined annually for 3 years (T2, T3 and T4). Dietary intake data were collected annually using 3 non-consecutive 24-h diet recalls. Of the original recruits, 1754 (98%) provided consent for the integration of their data and biological specimens into the NuAge Database and Biobank for future studies. Measured from serum/plasma samples, 30 biomarkers were used to calculate dysregulation globally and for five systems that a previous study validated as being largely independent : (1) oxygen transport; (2) liver/kidney function; (3) leukopoiesis; (4) micronutrients; and (5) lipids (see Table 1 for biomarkers in each score). We also calculated two other integrative, clinical-biomarker-based measures of biological ageing: phenotypic age (PhenoAge) and Klemera-Doubal biological age [35,36,37].
We assessed the effects of intakes of macronutrients and micronutrients, as well as their interactions, on measures of dysregulation and ageing. Our primary tool was the generalised additive model (GAM), a form of multiple regression. GAMs test for non-linear multidimensional effects using ‘smooth’ terms, which can revert back to simple linear terms (identical to those in linear regression) where the simpler effect gives the best fit to the data [38, 39]. GAMs are particularly useful in nutrition research where there is now abundant evidence that nutrient intakes can have non-linear effects on health outcomes (e.g. [4, 16, 40]). We explored three-dimensional effects of nutrient intake. Because GAMs can estimate non-linear effects, qualitative interpretation of the sign of estimated effects comes through visualisation (rather than a single regression coefficient). Here effects were visualised using the nutrient intake surfaces common to the GFN (Fig. 1). Importantly, GAMs can be used to correct for factors (e.g. sociodemographic status) that might be expected to confound relationships, in the same way as standard linear regression can be used in epidemiology. For each outcome, we fitted a series of eight models exploring different nutritional predictors and making statistical corrections for different factors. Factors explored were income, education level, age, physical activity, number of comorbidities, sex and current smoking status. Models 1 through 4 explored the effects of macronutrient intake with differing degrees of correction. Models 5 and 6 explored the 3-way effects of different combinations of micronutrients, while models 7 and 8 contained terms for both macronutrients and micro-nutrients simultaneously (see Additional file 1, Text S1 for a full description of all models) [41,42,43,44,45,46].
The main text contains a complete case analysis comprising 3569 observations from 1560 people. In Additional file 1 (Texts S2 and S3) we report sensitivity analyses, where we have imputed missing income data, and also analysed a more exclusive subset of the complete cases dataset (see Additional file 1: Table S1 for population summaries; the exclusive dataset excluded diabetics, individuals on prescribed diets, BMI <22 or >29.9, or with substantial weight fluctuation). These sensitivity analyses estimated similar effects to those in the main text. However, for the exclusive dataset, in many places, these effects (although qualitatively similar) do not meet the criteria for statistical significance. We interpret this latter point as evidence that the effects in the two datasets are similar, but that the power of the complete case analysis is required to detect statistical significance.
Our first model, model 1, tested for the effects of macronutrient intake (in kJ/day) on outcomes without any statistical corrections (see Additional file 1: Text S1 for parameters of all models). Because this model does not make correction for any confounders its output shows the unadjusted association between macronutrient intake and dysregulation/ageing scores within the data. We detected statistically significant effects of macronutrient intake on liver/kidney and micronutrient system dysregulation scores, as well as biological age (Fig. 2A–C; see Additional file 1: Table S2 for all statistical model output). Model 1 predicted a relatively minor effect of protein on liver/kidney function dysregulation, and non-linear effects of all carbohydrates and fats on liver/kidney function dysregulation. Individuals who consumed high (> 6000 kJ/day) or low (< 3000 kJ/day) levels of carbohydrates typically had slightly elevated (around 0.25 SD above the population mean) dysregulation scores (Fig. 2A). Very high intakes of lipids (> 4000 kJ/day; note this is within 2SD above the mean lipid intake) had the highest liver/kidney function dysregulation (~0.4 SD above the mean; Fig. 2A). With regard to micronutrient dysregulation scores, individuals with moderately high carbohydrate (5000 to 6000 kJ/day) and low lipid (< 2000 kJ/day) intake had minimal dysregulation. Again, high lipid intake was associated with maximal dysregulation (Fig. 2B). Subjects with protein intake of around 2000 kJ/day were predicted by model 1 to have low dysregulation (Fig. 2B). Biological age was predicted to be the lowest for individuals with high levels of intake for all three macronutrients (Fig. 2C); note that in this analysis nutrient intakes are not expressed relative to any measure of individual requirements.
Because model 1 shows unadjusted associations that may be confounded, in model 2 we considered dysregulation score as a function of macronutrient intake relative to what we define as the typical energy intake given height, weight, age, sex and physical activity (Additional file 1: Text S4). We again detected effects of macronutrient intake on liver/kidney function and micronutrient dysregulation and biological age score, but in addition, effects on PhenoAge were detected (Fig. 3A–D; Additional file 1: Table S3). Again, this analysis indicated high protein intake relative to what is typical (100% above average) had low liver/kidney function dysregulation scores (Fig. 3A). Interestingly, both PhenoAge and biological age were predicted to be low at elevated carbohydrate levels and typical levels of lipid and protein (Fig. 3C, D). These effects remained after making corrections for potential confounders, both excluding and including the number of comorbidities (models 3 and 4; Additional file 1: Tables S4 and S5), and showed a similar trend within the exclusive dataset (see Additional File 1: Text S3 and Figure S3).
Dietary micronutrients and nutrient subclasses
Intakes of a number of the micronutrient and nutrient subclasses (hereafter referred to collectively as micronutrients) explored were strongly correlated (Fig. 4). Hierarchical clustering based on correlation distance suggested seven clusters of micronutrients with very highly correlated intakes (r > 0.65). The first principal component (PC1) of intakes within each of the seven clusters explained between 90 and 100% of the variation of intake of micronutrients within that cluster (Additional file 1: Table S6). PC1 estimates uniformly displayed positive correlations with intakes of micronutrients within the clusters. For subsequent analyses, cluster-specific PC1 values were used as measures of intakes of micronutrients within clusters.
After using PCA to reduce dimensionality, we were left with 19 minimally correlated micronutrient/subclass variables (PC1 of clusters or individual micronutrients), making it feasible to run GAMs for all 969 3-way micronutrient combinations for each dysregulation score (micronutrient-specific models; see Additional file 1: Text S1). For all scores except lipid dysregulation, there were a greater number of significant micronutrient smooth terms than would be expected under the null hypothesis (Fig. 5A–H). After correction for the false discovery rate (FDR), there were no significant effects of micronutrient intake on oxygen transport or lipid dysregulation scores. There were 905 models that detected significant effects for at least one score, although 363 of these were solely related to ageing scores (i.e. no significant effect on any dysregulation score; Additional file 1: Table S7). There were 17 combinations for which we detected significant effects on all scores (except oxygen transport and lipid dysregulation). Interestingly all 17 models contained α-tocopherol (vitamin E) as one of the three micronutrients. Arguably any one of these micronutrient combinations may be of interest and could warrant further investigation in response to a priori hypotheses (results from all models can be found at https://cohenaginglab.github.io/micronuage/). Here however we focus on the effects of α-tocopherol, vitamin C and trans-fatty acids because this combination had the highest mean per cent deviance explained across all scores.
In models 5 and 6 we tested for effects of these three micronutrients with correction for confounders, excluding and including comorbidities respectively. In these models, we detected the effects of micronutrients on leukopoiesis, liver/kidney function, micronutrients and global dysregulation (Additional file 1: Tables S8 and S9). The effects varied slightly across the four systems. (Fig. 6A–D for model 6). For example, elevated intakes of trans-fatty acids are predicted to be detrimental for liver/kidney dysregulation, but have a minor beneficial effect on leukopoiesis dysregulation. Nevertheless, consuming around 2SD of α-tocopherol above the population mean, while consuming vitamin C and trans-fatty acids at the population mean results in low dysregulation across the four scores, suggesting a systemic benefit of high, but not excessive, α-tocopherol intake (Fig. 6). To ensure that increased α-tocopherol does not have detrimental interactions with other micronutrients, we screened the complete list of significant micronutrient combinations (after FDR adjustment) for those that included α-tocopherol intake; 153 combinations were identified for at least one outcome score. Figure 5I shows the distribution of dysregulation scores (all traits) predicted for a 2SD increase in α-tocopherol for any identified effects of micronutrients (after inclusion for potential confounders). For all scores, intakes of α-tocopherol 2SD above our sample mean (mean ± SD = 4.75 ± 2.73) is predicted to lead to reductions in dysregulation or ageing Z-scores (below the population average) or changes close to 0 (the population average).
Macro, micro or both
We found statistical support for the effects of both macro and micronutrients on liver/kidney function and micronutrient dysregulation. There can be strong covariances amongst macro and micronutrient intake owing to their co-occurrences in foods (see Additional file 1: Text S5). For example, vitamin C intake is often considered to be a marker of fruit and vegetable consumption , and unsurprisingly in this dataset high vitamin C intake coincides with a high carbohydrate diet (Additional file 1: Text S5). One may thus question which of these two nutritional levels is a better predictor of dysregulation, or if both must be considered. We evaluated this by fitting the two three-dimensional effects for macro and micronutrients simultaneously (models 7 and 8, excluding and including comorbidities respectively) and assessing model fit relative to corresponding previous models using the Akaike information criterion (AIC; models with minimal AIC are favoured). For liver/kidney function dysregulation a model including macronutrients is favoured by AIC (Table 2). In contrast, for the micronutrient dysregulation score, a model including only micronutrients is favoured by AIC (Table 2). We also note here that even the most complex models fitted, which include both macro and micronutrients, as well as other potential predictors of health status, explain less than 5% of the deviance in these dysregulation scores for this population (Table 2).
Here we used multidimensional modelling techniques to test for associations between nutrient intake and physiological dysregulation in an aged population. We find that consuming above-average protein and α-tocopherol intake, which is the most active form of ‘vitamin E’, is associated with lower levels of physiological dysregulation. In contrast, individuals with above-average carbohydrates (given what is typical for height, weight, age sex and physical activity) coupled with typical protein and lipid intakes had minimal measures of biological ageing.
Analysis of the effects of dietary macronutrients on age-specific mortality in mice has shown that protein restriction when coupled with increased carbohydrate intake, extends median lifespan through reduced mortality in middle life, but that higher protein in late life may reduce mortality . Given that our study population consisted entirely of older adults, our findings with regard to physiological dysregulation are consistent with this experimental literature. Our results are also concordant with numerous studies highlighting the need for increased protein intake in older people, in particular, to offset sarcopenia and decreased physical performance associated with ageing (reviewed in [22, 23]). It is interesting to speculate, though, whether reapplying our methods to middle-aged cohorts would detect similar benefits of elevated protein intake; the experimental literature, and some data in humans, suggests perhaps not [20, 21, 48].
Numerous micronutrient interactions have been demonstrated in the experimental absorption literature (e.g. calcium-iron and vitamin C-iron ). Our analysis indicates that such interactions may have strong enough effects on health in old age to be detected at the level of the population. For example, elevated vitamin E intake can produce either low or high micronutrient dysregulation scores, depending upon the level of vitamin C intake (Fig. 6C). In vitro and in vivo studies, including a supplementation trial in healthy adults, have detected vitamin C-E absorption interactions, including where simultaneous supplementation leads to elevated circulating levels, and where vitamin C can reduce oxidised vitamin-E recovering its role as an antioxidant [50,51,52]. Thus, there is a mechanistic basis for what we have detected, and it is certainly worth testing for interactions between these micronutrients in other epidemiological cohorts.
We found that consuming α-tocopherol at 2SD above the population average is associated with benefits; in the subject population, this level corresponds to 10.21mg/day of α-tocopherol. The World Health Organisation recommended intake of α-tocopherol for those aged 65+ is 7.5mg/day in females and 10mg/day in males, thus the value we highlight is not substantially beyond current guidelines . Highlighting the importance of considering non-linearity, more extreme intake patterns (e.g. >4 SDs above average) were associated with harm (Fig. 5). This finding accords with the results of RCTs suggesting excessive vitamin E supplementation may increase all-cause mortality . We also detected a non-linear effect of carbohydrate intake on dysregulation, which suggests individuals at the upper/lower extremes of the observed carbohydrate intakes suffer poor health. Epidemiological and meta-analytic study of the effects of carbohydrates on all-cause mortality in humans has found identical patterns , possibly moderated by carbohydrate quality. Generally speaking, our study therefore provides further support to the importance of looking beyond ‘single nutrient at a time’ and monotonic ‘more is better’ analyses [1, 55], to detect the interactive effects of nutrients. For space reasons, we have not presented and discussed the full range of results for all micronutrient intake interactions detected here. However, we have illustrated an approach that can be used to study the effects of multiple micronutrient intakes on health. The GAM-based approach can either be used for discovery (i.e. as means to detect micronutrient interactions) or as a targeted mechanism-driven paradigm to test specific hypotheses about micronutrient interactions derived from experimental biology or nutrition science [49, 56]. We encourage readers to explore http://cohenaginglab.github.io/micronuage where the effects of micronutrient intakes on dysregulation and ageing scores in this population can be visualised interactively, a priori hypotheses explored, and complete models downloaded.
Even in our most complex models, the deviance explained remains at around 5%. This is within the bounds of what analyses of other outcomes in this cohort have found (e.g. ). Nevertheless, the surfaces we present here show relatively large effects of nutrition on dysregulation levels, which is potentially important as nutrition is readily modifiable. The low deviance explained may be due to a number of factors including the reporting errors, and systematic reporting biases inherent in epidemiological studies of nutrition [58,59,60]; the fact we consider additive effects of macro and micronutrients rather than interactions (more complex models are theoretically possible but are limited by sample size and visualisation beyond three-dimensions); and the rarity of the extreme dietary profiles at the edges of the surfaces, where the strongest effects are found (the standard error of the surface is a proxy for sample density; e.g. Additional file 1: Fig. S6). An implication of this latter point is that some of the largest effects of nutrient intake we report will only apply to a small proportion of the population and that our physiology is often robust enough to tolerate relatively wide variation without much consequence. Similar patterns are observed when using the GFN to map evolutionary fitness to organisms that ecologists pre-define as dietary ‘generalists’ [16, 61,62,63]. This is consistent with an understanding of nutrition in which our ancestors evolved to tolerate an array of dietary patterns . Accordingly, homeostasis can be maintained across a wide array of nutritional states, with the caveat that when diet becomes too extreme dysregulation can increase very rapidly (‘falling off a dietary cliff’). The tolerance for different diets—the size of the plateau in this analogy—could of course vary as a function of genetic or environmental factors that predispose us to greater risk .
In contrast to dysregulation, we found that PhenoAge and biological age were minimised on relatively high carbohydrate and lower to moderate protein intakes. Similar results have been observed in numerous experimental studies on the effects of protein to carbohydrate ratio on ageing in model organisms . The somewhat discordant results between dysregulation measures and PhenoAge/Biological Age are not necessarily that surprising. Increasingly, the literature shows that ageing is multivariate and heterogeneous [10, 28]. In this context, our results imply that different dietary patterns come with their own benefits and drawbacks in their effects on the different facets of ageing. A practical consequence of this specificity is that dietary recommendations could be tailored to slow the most advanced ageing processes based on an individual’s biological profile. Given the strength of the apparent trade-off between different ageing/dysregulation scores, it seems unlikely that any one intake pattern can simultaneously minimise all scores. However, once a key target score of interest is identified, multi-criteria optimisation may help define a set of ‘non-dominated’ dietary intakes that maximally improve that score without degrading other outcomes unnecessarily. Regardless of the methods used to advance these findings, replication in other cohorts is needed to confirm the precise patterns we report, and even then, context-specific clinical recommendations will be essential.
An important next step for analyses of this nature will be to complement them with analyses of the effects of food/broad dietary patterns on physiological dysregulation . The intakes of individual nutrients are often markers for broader dietary patterns. High vitamin C/E intake can reflect a diet rich in fresh fruit and vegetables , whereas a high trans-fatty acid intake may reflect a diet comprised of processed foods. Such analyses may help to elucidate why, for example, we detect a minor beneficial effect of trans-fatty acids on leukopoiesis dysregulation, where other studies have detected the opposite ; it is possible that in this population high trans-fatty acid intakes represent a diet high in dairy, as opposed to processed foods.
Previously, studies have applied GFN thinking to experimental and observational contexts in humans. They have shown that the macronutrient composition of the diet is associated with biomarkers of cardiovascular health, energy intake, obesity and specific chronic diseases [68,69,70,71,72,73]. Important as they are, these analyses were concerned with unidimensional measures of health and healthspan. Here, for the first time, we have used the GFN to model holistic measures of systemic functioning during ageing. How results from experimental nutrition and geroscience map onto humans living in the community remains unclear. In part, our understanding has been hampered by a lack of techniques that cut across the complexity inherent to nutrition and biological ageing in real-world contexts; the approach we present here is a promising start. Future applications could include personalised approaches to aid in healthy ageing, and screening of at-risk older adults to ensure they do not fall off the ‘dietary cliff’. Finally, our results advocate against the popular practice of eating to maximise or minimise certain nutrients. The dose-response relationship is often U-shaped, and highly dependent on context (e.g. age, other aspects of diet); targeting in the absence of clear evidence is likely to do more harm than good.
Participant information and dietary intake data
The NuAge Database and Biobank and the present study have been approved by the Research Ethics Board (REB) of the CIUSSS-de-l’Estrie-CHUS (Quebec, Canada; projects 2019-2832 and 2015-868/14-141, respectively). The original sample includes 1587 individuals recruited between November 2003 and June 2005 (T1) to which 206 volunteers were added . They were community-dwelling men and women, aged 67–84 years, able to speak English or French, and in good general health at recruitment. Notably, they had to be free of disabilities in activities of daily living, not cognitively impaired and able to walk 300 metres or to climb 10 stairs without rest. A structured interview was conducted annually in the NuAge Study at baseline (T1) and for the next 3 years (T2, T3 and T4) to gather the following data. Sociodemographic (actual income, education) and lifestyle (smoking, alcohol) information was obtained using a general study questionnaire developed from standard health survey questions . The number of self-reported chronic health conditions (i.e. comorbidities) was computed from an adaptation of the Older Americans Resources and Services (OARS) Multidimensional Functional Assessment questionnaire . Chronic health conditions considered were self-reported cancer (within the last 5 years); hypertension; self-reported liver or gallbladder disease; self-reported surgery of the digestive system; self-reported heart trouble; self-reported circulation trouble in arms or legs; self-reported thrombosis, cerebral haemorrhage, or cerebrovascular accident; self-reported transient cognitive impairment; self-reported Parkinson's disease; self-reported diabetes; self-reported emphysema or chronic bronchitis; self-reported asthma; rheumatoid arthritis, arthritis, or rheumatism; self-reported osteoporosis; self-reported kidney diseases; self-reported thyroid and gland problems; self-reported surgery of the circulatory system; and mini-mental state examination (MMSE) total score between 22 and 26. Usual physical activity was assessed using the Physical Activity Scale for the Elderly (PASE) questionnaire, with higher scores corresponding to higher physical activity levels . Weight and height were measured, and body mass index (BMI = weight [kg]/height [m]2) calculated for each participant .
Dietary intake data were collected annually (T1 to T4) using 3 non-consecutive 24-h diet recalls [77,78,79]. Each set included 2 weekdays and 1 weekend day, with the first administered during the annual face-to-face interview and the others by telephone interviews without prior notice. Based on the USDA 5-step multiple-pass method , interviewers recorded a detailed description and portion sizes of all foods and beverages consumed by each participant the day before the interview. Only energy and nutrients coming from foods were considered (i.e. excluding supplements). All interviewers were trained registered dietitians. Energy and nutrient intake were computed from the 24-h diet recalls using the CANDAT-Nutrient Calculation System (version 10, ©Godin London Inc.) based on the 2007b version of the Canadian Nutrient File (CNF) from Health Canada and a database of 1200 additional foods that were developed on site . This software produced intakes of the three macronutrients, as well as a full suite of micronutrients and several nutrient subclasses (principally lipid subclasses). We assessed the effects of intakes of all nutrient dimensions available to us.
Biomarker data and dysregulation score
Biomarkers were chosen based on our previous study [10, 82], their clinical use, and according to their availability in NuAge. Biomarkers were not chosen for their a priori expected relationship to either diet or ageing; indeed, our previous work has shown that dysregulation measures are largely robust to the choice of component markers, reflecting instead an emergent property of the system as a whole. Accordingly, while future work will look at markers of specific interest such as FGF-21, this study does not use markers that are expected individually to integrate metabolism (fatty acids, hydroxybutyrate, insulin, adiponectin, etc.). In total, 30 biomarkers were used to calculate dysregulation globally and for five systems that a previous study validated as being largely independent : (1) oxygen transport; (2) liver/kidney function; (3) leukopoiesis; (4) micronutrients; and (5) lipids (Table 1).
Not all individuals had complete data on biomarkers and nutrition, so sample size depended on the physiological system and the analysis. For the oxygen transport and leukopoiesis systems, the sample size varied between 1224 and 1331. For the other systems, the sample size varied between 654 and 730. The difference in sample size is due to a more extensive serum biomarker analysis conducted in 2016 on a subsample of ~750 individuals that were selected at random amongst the 904 of the 1754 participants who met the following criteria: (1) the individual needed to have blood sampling conducted without any missing intermediate visits; (2) the individual needed at least two visits with blood samples; and (3) there had to be a sufficient number of stored aliquots at all visits to analyse the selected biomarkers.
Dysregulation was calculated using the Mahalanobis distance (DM), a statistical distance which measures how deviant or aberrant an individual’s biomarker profile is compared to a reference population. Using DM, a global dysregulation score (‘Global’) was calculated with all 30 biomarkers and by the physiological system for each participant at each visit, as previously described . Higher scores indicate higher dysregulation. Age at recruitment in NuAge is restricted from 67 to 84 years, thus, to construct our reference population we used data from the first visit, which represents a population in general good health, according to the inclusion and exclusion criteria (22). Some biomarkers were log- or square-root- transformed as necessary to approach normality. All biomarkers were then centred at the mean of the reference population and divided by the standard deviation of the reference population to standardise them before calculation of dysregulation scores (i.e. UV scaled).
Alongside the 6 measures of physiological dysregulation, we also included two other measures of biological ageing: phenotypic age (PhenoAge) and Klemera-Doubal biological age. PhenoAge measures the predicted age of the individual based on the individual’s mortality risk as assessed via 9 biomarkers and calibrated in NHANES IV. Biological age measures the age of an individual as predicted based on linear projections of a set of age-associated biomarkers. We calculated PhenoAge as described elsewhere [36, 37], although because CRP (C-reactive protein) is not available in the NuAge cohort, we used the mean CRP value from NHANES IV 1999–2010  for all individuals. Validation with individually-imputed values showed a very strong correlation (r=0.99) with PhenoAge scores using this mean value; we chose the mean value as the more conservative option and simpler model. PhenoAge showed a weaker correlation with chronological age in our population than previously reported (r=0.67 vs 0.96 ;). Exploratory analysis of other datasets in our possession indicated the discrepancy was likely due to the limited age range rather than the absence of CRP (data not shown). Calculation of biological age was based on previous work by Levine . We first searched for biomarkers that correlated with chronological age and found 13 with r > 0.1 (p ≤ 0.05). After removing three biomarkers with numerous missing values (non-HDL cholesterol, LDL, and estimated glomerular filtration rate), we were left with the following list: haemoglobin, haematocrit, red blood cell count, red cell distribution width, monocyte count, albumin, folate, creatinine, blood urea nitrogen, and lymphocyte percentage. Using these biomarkers, we calculated biological age as previously described [35, 83].
All analyses were performed in the statistical programming environment R . In all cases, dysregulation scores were first log-transformed and then Z-transformed (mean centred and divided by one SD) prior to model fitting. Analysis was performed using GAMs [38, 85]. GAMs are a form of the generalised linear (mixed) models (GLMs, and GLMMs) that allow the user to model complex non-linear effects through the inclusion of non-parametric ‘smooth’ terms (usually implemented as a spline). These terms can be specified alongside the conventional linear parametric terms in a linear model. Smooth terms can be additive in a single dimension, or multi-dimensional to test for synergistic interactions. The statistical significance of smooth terms can be interpreted via a p-value, but their effect must be interpreted visually (as opposed to numerically via a regression coefficient). GAMs were implemented using the ‘gam’ function in the R package mgcv, and terms estimated by restricted maximum likelihood [38, 39].
We first tested the hypothesis that increased protein is associated with improved health in old age. We estimated the effects of daily intake of macronutrients (protein, carbohydrate and fat, in kJ) on each dysregulation/ageing score using GAMs. In all models, dysregulation/ageing score at an observation was the outcome. In model 1, macronutrient intakes at those observations were fitted as a three-dimensional smooth term (thin-plate spline), and individual subject ID was fitted as a random effect (there are multiple observations per individual). However, the dataset is made up of individuals whose energy requirements likely vary due to height, weight, age, sex and physical activity level. Thus, we also estimated macronutrient intakes relative to what is typical for an individual in this population given the aforementioned variables (estimated from the residuals of a model of intake; see Additional file 1: Text S4). In model 2, we modelled dysregulation scores as a function of relative macronutrient intake using a GAM as above.
In model 3, we tested for the effects of confounders by refitting model 2 to include potential confounders as additive effects; sex (men/women), smoking status (current/not current), current income, number of years of education, alcohol intake (g/day), physical activity level (PASE) and age were included. We ran a separate model with comorbidities as an additional potential confounder (model 4). Numeric predictors were Z-transformed prior to fitting and were included in the model as smooth terms, and categorical predictors as parametric terms.
In the second part of our analyses, we considered the effects of those micronutrients for which intake data were available. Micronutrient intakes are likely to be highly correlated. Therefore, we performed hierarchical clustering on correlations (using 1—correlation coefficient) of micronutrient intakes based on correlation distance using the ‘hclust’ function in the stats package in base R. For any clusters of highly correlated micronutrients, we performed a principal component analysis (PCA; ‘prcomp’ function in base R) and used the first principal component (PC1) as our measure of intake for the nutrients within that cluster. The number of micronutrients adds complexity to understating their interactive effects because (1) the estimation of multidimensional smooth terms in GAMs becomes challenging as the number of dimensions grows and (2) the smooth terms in GAMs must be interpreted visually. For these reasons, we restricted ourselves to considering 3 micronutrient dimensions at a time. For each three-way combination of micronutrient intake (PC1 for clusters) we fitted a GAM with a three-dimensional smooth term for the three micronutrients and a random effect for the individual. We then identified those models with significant effects after correction for the false discovery rate (FDR, q < 0.05 ;). We note that correction by FDR assumes that p-values are independent, which is not the case here. However, overlooking this non-independence is conservative in that it will result in fewer significant effects rather than more, and thus we proceeded without correction. In the main text, we interpret specific cases of interest, however, all results and models can be accessed at http://cohenaginglab.github.io/micronuage. We also tested for effects of micronutrients alongside correction for the potential confounders discussed above (models 5 and 6), and by modelling micro and macronutrient intakes simultaneously (models 7 and 8). Note that including two three-dimensional smooth terms (one for micronutrients and one for macronutrients) is not the same as a six-dimensional smooth term; this approach is substantially less power-hungry but still allows us to adjust micronutrient models for macronutrient intake, and vice-versa. See Additional file 1 (Text S1) for a concise list of all models implemented.
Statistical significance was inferred for terms in GAMs based on α of 0.05. Where we found it necessary to compare amongst models containing different predictors we used the Akaike Information Criterion (AIC), where smaller values (beyond a margin of two points) indicate a better model fit . Smooth terms from GAMs were interpreted from model predictions. For three-dimensional effects of nutrient intake, we present two-dimensional surface plots, where intakes of nutrients are given on the x- and y-axes, while the intake of the third nutrient is held constant at a value stated on the figure panel. On all surfaces deep blue/red areas indicate low/high Z bio-scores (i.e. good/bad health), respectively.
The main text contains a complete cases analysis where we report results from analyses on all observations for which relevant predictor and dysregulation score data were available (i.e. inclusive dataset with missing income data excluded). Additional file 1, Supplementary Texts S2 and S3 contain the results of two sets of sensitivity analyses. In the first, we have imputed missing income data (31%) as the participant-specific mean value to increase our sample size. In the second we have analysed a subset of the data (<50% of the total data) where we have excluded any observation at which the subject was recorded as either; diabetic (type-1 or 2), reported as being on a medically prescribed diet, having a BMI outside the range of 22 to 29.9, or that came from a subject with a coefficient of variation of weight > 0.04 over the course of all the observations (i.e. for whom weight fluctuated substantially).
Key analyses were performed by two separate members of the team, using broadly the same approach but with no consultation on the detailed analytical decisions (i.e. what might be inferred from seeing the Results but not the Methods). Discrepancies leading to qualitative changes in the conclusions were flagged and resolved. VL performed this validation, with AMS as primary data analyst. Datasets containing all participants’ variables used in this study and those from dietary intakes estimated from 24-h dietary recalls were transmitted by the NuAge team in November–December 2018 and October 2019, respectively. Data for some biomarkers were obtained in a previous secondary project .
Availability of data and materials
All data generated during this study are included in this published article, its supplementary information files and publicly available repositories. Code is available at https://github.com/AlistairMcNairSenior/NuAGE_GFN (DOI: 10.5281/zenodo.6984539) All surfaces and associated models from which conclusions are drawn are available at http://cohenaginglab.github.io/micronuage. The raw data are subject to strict restrictions due to ethical concerns and their use must be assessed by the NuAge Research Committee. Requests for access to then raw data should be made to NuAgefirstname.lastname@example.org for review.
Raubenheimer D, Simpson SJ. Nutritional ecology and human health. Ann Rev Nutr. 2016;36:603–26.
Hoffmann I. Transcending reductionism in nutrition research. Am J Clin Nutr. 2003;78(3):514S–6S.
Samieri C, Sun Q, Townsend MK, Chiuve SE, Okereke OI, Willett WC, et al. The association between dietary patterns at midlife and health in aging. Ann Intern Med. 2013;159(9):584–91.
Solon-Biet SM, McMahon AC, Ballard JWO, Ruohonen K, Wu L, Cogger VC, et al. The ratio of macronutrients, not caloric intake, dictates cardiometabolic health, aging, and longevity in ad libitum-fed mice. Cell Metab. 2014;19(3):418–30.
Fischer K, Colombani PC, Langhans W, Wenk C. Cognitive performance and its relationship with postprandial metabolic changes after ingestion of different macronutrients in the morning. Brit J Nutr. 2001;85(3):393–405.
Simpson SJ, Le Couteur DG, Raubenheimer D. Putting the balance back in diet. Cell. 2015;161(1):18–23.
Cappola AR, Maggio M, Ferrucci L. Is research on hormones and aging finished? No! Just started! J Gerontol A Biol Sci Med Sci. 2008;63(7):696–7.
Fulop T, Larbi A, Dupuis G, Le Page A, Frost EH, Cohen AA, et al. Immunosenescence and inflamm-aging as two sides of the same coin: Friends or foes? Front Immunol. 2018;8:1960.
Fried LP, Xue Q-L, Cappola AR, Ferrucci L, Chaves P, Varadhan R, et al. Nonlinear multisystem physiological dysregulation associated with frailty in older women: Implications for etiology and treatment. J Gerontol A Biol Sci Med Sci. 2009;64(10):1049–57.
Li Q, Wang S, Milot E, Bergeron P, Ferrucci L, Fried LP, et al. Homeostatic dysregulation proceeds in parallel in multiple physiological systems. Aging Cell. 2015;14(6):1103–12.
Cohen AA. Complex systems dynamics in aging: New evidence, continuing questions. Biogerontology. 2016;17(1):205–20.
Cohen AA, Legault V, Fuellen G, Fülöp T, Fried LP, Ferrucci L. The risks of biomarker-based epidemiology: Associations of circulating calcium levels with age, mortality, and frailty vary substantially across populations. Exp Gerontol. 2018;107:11–7.
Gurven M, Kaplan H, Winking J, Eid Rodriguez D, Vasunilashorn S, Kim JK, et al. Inflammation and infection do not promote arterial aging and cardiovascular disease risk factors among lean horticulturalists. PLoS One. 2009;4(8):e6590.
Raubenheimer D, Simpson SJ, Le Couteur DG, Solon-Biet SM, Coogan SCP. Nutritional ecology and the evolution of aging. Exp Gerontol. 2016;86:50–61.
Piper Matthew DW, Partridge L, Raubenheimer D, Simpson SJ. Dietary restriction and aging: A unifying perspective. Cell Metab. 2011;14(2):154–60.
Simpson SJ, Raubenheimer D. The nature of nutrition: A unifying framework from animal adaptations to human obesity. Oxford: Princeton University Press; 2012.
Raubenheimer D, Simpson SJ. The geometry of compensatory feeding in the locust. Anim Behav. 1993;45(5):953–64.
Solon-Biet SM, Walters KA, Simanainen U, McMahon AC, Ruohonen K, Ballard JWO, et al. Macronutrient balance, reproductive function and lifespan in aging mice. Proc Natl Acad Sci U S A. 2015;112:3481–6.
Lee KP, Simpson SJ, Clissold FJ, Brooks R, Ballard JWO, Taylor PW, et al. Lifespan and reproduction in drosophila: New insights from nutritional geometry. Proc Natl Acad Sci U S A. 2008;105(7):2498–503.
Le Couteur DG, Solon-Biet S, Cogger VC, Mitchell SJ, Senior A, de Cabo R, et al. The impact of low-protein high-carbohydrate diets on aging and lifespan. Cell Mol Life Sci. 2016;73(6):1237–52.
Senior AM, Solon-Biet SM, Cogger VC, Le Couteur DG, Nakagawa S, Raubenheimer D, et al. Dietary macronutrient content, age-specific mortality and lifespan. Proc Roy Soc B. 1902;2019(286):20190393.
Landi F, Calvani R, Tosato M, Martone MA, Ortolani E, Savera G, et al. Protein intake and muscle health in old age: From biological plausibility to clinical evidence. Nutrients. 2016;8(5):295.
Bauer J, Biolo G, Cederholm T, Cesari M, Cruz-Jentoft AJ, Morley JE, et al. Evidence-based recommendations for optimal dietary protein intake in older people: A position paper from the prot-age study group. J Am Med Dir Assoc. 2013;14(8):542–59.
Senior AM, Nakagawa S, Raubenheimer D, Simpson SJ. Global associations between macronutrient supply and age-specific mortality. Proc Natl Acad Sci U S A. 2020;117(48):30824.
Cohen AA, Kennedy BK, Anglas U, Bronikowski AM, Deelen J, Dufour F, et al. Lack of consensus on an aging biology paradigm? A global survey reveals an agreement to disagree, and the need for an interdisciplinary framework. Mech Ageing Dev. 2020;191:111316.
López-Otín C, Blasco MA, Partridge L, Serrano M, Kroemer G. The hallmarks of aging. Cell. 2013;153(6):1194–217.
Kennedy BK, Berger SL, Brunet A, Campisi J, Cuervo AM, Epel ES, et al. Geroscience: Linking aging to chronic disease. Cell. 2014;159(4):709–13.
Belsky DW, Moffitt TE, Cohen AA, Corcoran DL, Levine ME, Prinz JA, et al. Eleven telomere, epigenetic clock, and biomarker-composite quantifications of biological aging: Do they measure the same thing? Am J Epidemiol. 2017;187(6):1220–30.
Cohen AA, Martin LB, Wingfield JC, McWilliams SR, Dunne JA. Physiological regulatory networks: Ecological roles and evolutionary constraints. Trends Ecol Evol. 2012;27(8):428–35.
Cohen AA, Milot E, Yong J, Seplaki CL, Fülöp T, Bandeen-Roche K, et al. A novel statistical approach shows evidence for multi-system physiological dysregulation during aging. Mech Ageing Dev. 2013;134(3):110–7.
Cohen AA, Li Q, Milot E, Leroux M, Faucher S, Morissette-Thomas V, et al. Statistical distance as a measure of physiological dysregulation is largely robust to variation in its biomarker composition. PLoS One. 2015;10(4):e0122541.
Cohen AA, Milot E, Li Q, Legault V, Fried LP, Ferrucci L. Cross-population validation of statistical distance as a measure of physiological dysregulation during aging. Exp Gerontol. 2014;57:203–10.
Milot E, Morissette-Thomas V, Li Q, Fried LP, Ferrucci L, Cohen AA. Trajectories of physiological dysregulation predicts mortality and health outcomes in a consistent manner across three populations. Mech Ageing Dev. 2014;141-142:56–63.
Gaudreau P, Morais JA, Shatenstein B, Gray-Donald K, Khalil A, Dionne I, et al. Nutrition as a determinant of successful aging: Description of the quebec longitudinal study nuage and results from cross-sectional pilot studies. Rejuvenation Res. 2007;10(3):377–86.
Levine ME. Modeling the rate of senescence: Can estimated biological age predict mortality more accurately than chronological age? J Gerontol Ser A Biol Sci Med Sci. 2013;68(6):667–74.
Levine ME, Lu AT, Quach A, Chen BH, Assimes TL, Bandinelli S, et al. An epigenetic biomarker of aging for lifespan and healthspan. Aging. 2018;10(4):573–91.
Liu Z, Kuo P-L, Horvath S, Crimmins E, Ferrucci L, Levine M. A new aging measure captures morbidity and mortality risk across diverse subpopulations from nhanes iv: A cohort study. PLoS Med. 2019;15(12):e1002718.
Wood SN. Generalized additive models: An introduction with r: Chapman and Hall/CRC; 2017.
Wood SN. Fast stable restricted maximum likelihood and marginal likelihood estimation of semiparametric generalized linear models. J R Stat Soc Series B. Stat Methodol. 2011;73:3–36.
Seidelmann SB, Claggett B, Cheng S, Henglin M, Shah A, Steffen LM, et al. Dietary carbohydrate intake and mortality: A prospective cohort study and meta-analysis. Lancet Public Health. 2018;3(9):e419–28.
Lawson J, Willden C. Mixture experiments in r using mixexp. J Stat Soft. 2016;72:1–20 Code snippet 2.
Saner C, Tassoni D, Harcourt BE, Kao K-T, Alexander EJ, McCallum Z, et al. Evidence for protein leverage in children and adolescents with obesity. Obesity. 2020;28(4):822–9.
Song M, Giovannucci E. Substitution analysis in nutritional epidemiology: Proceed with caution. Eur J Epidemiol. 2018;33(2):137–40.
Akaike H. Information theory as an extension of the maximum likelihood principle. In: Second International Symposium on Information Theory. Budapest: Akademiai Kiado; 1973.
Raubenheimer D. Toward a quantitative nutritional ecology: The right-angled mixture triangle. Ecol Monogr. 2011;81(3):407–27.
Nakagawa S, Cuthill IC. Effect size, confidence interval and statistical significance: A practical guide for biologists. Biol Revs. 2007;82(4):591–605.
Michels KB, Welch AA, Luben R, Bingham SA, Day NE. Measurement of fruit and vegetable consumption with diet questionnaires and implications for analyses and interpretation. Am J Epidemiol. 2005;161(10):987–94.
Levine Morgan E, Suarez Jorge A, Brandhorst S, Balasubramanian P, Cheng C-W, Madia F, et al. Low protein intake is associated with a major reduction in igf-1, cancer, and overall mortality in the 65 and younger but not older population. Cell Metab. 2014;19(3):407–17.
Sandström B. Micronutrient interactions: Effects on absorption and bioavailability. Brit J Nutr. 2001;85(S2):S181–5.
Hamilton IMJ, Gilmore WS, Benzie IFF, Mulholland CW, Strain JJ. Interactions between vitamins c and e in human subjects. Brit J Nutr. 2000;84(3):261–7.
Tanaka K, Hashimoto T, Tokumaru S, Iguchi H, Kojo S. Interactions between vitamin c and vitamin e are observed in tissues of inherently scorbutic rats. J Nutr. 1997;127(10):2060–4.
Traber MG: Vitamin e. In: Modern nutrition in health and disease. Ross AC, Caballero B, Cousins RJ, Tucker KL, Ziegler TR, 11th. Philadelphia: Lippincott Williams & Wilkins; 2012: 293-304.
Nations WHOaFaAOotU. Vitamin and mineral requirments in human nutrition. 2nd ed; 2004.
Miller ER, Pastor-Barriuso R, Dalal D, Riemersma RA, Appel LJ, Guallar E. Meta-analysis: High-dosage vitamin e supplementation may increase all-cause mortality. Ann Intern Med. 2005;142(1):37–46.
Fardet A, Rock E. Perspective: Reductionist nutrition research has meaning only within the framework of holistic and ethical thinking. Adv Nutr. 2018;9(6):655–70.
Lönnerdal B. Interactions between micronutrients: Synergies and antagonisms. In: Pettifor JM, Zlotkin S, editors. Micronutrient deficiencies during the weaning period and the first years of life. Basel: Karger; 2004. p. 67–81.
Mercille G, Richard L, Gauvin L, Kestens Y, Shatenstein B, Daniel M, et al. Associations between residential food environment and dietary patterns in urban-dwelling older adults: Results from the voisinuage study. Public Health Nutr. 2012;15(11):2026–39.
Hebert JR, Clemow L, Pbert L, Ockene IS, Ockene JK. Social desirability bias in dietary self-report may compromise the validity of dietary intake measures. Int J Epidemiol. 1995;24(2):389–98.
Dahm CC. Correcting measurement error in dietary exposure assessments: No piece of cake. Am J Clin Nutr. 2020;112(1):11–2.
Ferrari P, Slimani N, Ciampi A, Trichopoulou A, Naska A, Lauria C, et al. Evaluation of under- and overreporting of energy intake in the 24-hour diet recalls in the european prospective investigation into cancer and nutrition (epic). Public Health Nutr. 2002;5(6b):1329–45.
Machovsky-Capuska GE, Senior AM, Simpson SJ, Raubenheimer D. The multidimensional nutritional niche. Trends Ecol Evol. 2016;31:355–65.
Raubenheimer D, Simpson SJ. Nutrient balancing in grasshoppers: Behavioural and physiological correlates of dietary breadth. J Exp Biol. 2003;206(10):1669–81.
Warbrick-Smith J, Raubenheimer D, Simpson SJ, Behmer ST. Three hundred and fifty generations of extreme food specialisation: Testing predictions of nutritional ecology. Entomol Exp Appl. 2009;132(1):65–75.
Luca F, Perry GH, Di Rienzo A. Evolutionary adaptations to dietary changes. Ann Rev Nutr. 2010;30(1):291–314.
Nijhout HF, Sadre-Marandi F, Best J, Reed MC. Systems biology of phenotypic robustness and plasticity. Integr Comp Biol. 2017;57(2):171–84.
Tapsell LC, Neale EP, Satija A, Hu FB. Foods, nutrients, and dietary patterns: Interconnections and implications for dietary guidelines. Adv Nutr. 2016;7(3):445–54.
Mozaffarian D, Pischon T, Hankinson SE, Rifai N, Joshipura K, Willett WC, et al. Dietary intake of trans fatty acids and systemic inflammation in women. Am J Clin Nutr. 2004;79(4):606–12.
Martinez-Cordero C, Kuzawa CW, Sloboda DM, Stewart J, Simpson SJ, Raubenheimer D. Testing the protein leverage hypothesis in a free-living human population. Appetite. 2012;59(2):312–5.
Blumfield M, Hure A, MacDonald-Wicks L, Smith R, Simpson S, Raubenheimer D, et al. The association between the macronutrient content of maternal diet and the adequacy of micronutrients during pregnancy in the women and their children’s health (watch) study. Nutrients. 2012;4(12):1958–76.
Raubenheimer D, Machovsky-Capuska GE, Gosby AK, Simpson S. Nutritional ecology of obesity: From humans to companion animals. Brit J Nutr. 2015;113(S1):S26–39.
Gosby AK, Lau NS, Tam CS, Iglesias MA, Morrison CD, Caterson ID, et al. Raised fgf-21 and triglycerides accompany increased energy intake driven by protein leverage in lean, healthy individuals: A randomised trial. PLoS One. 2016;11(8):e0161003.
Simpson SJ, Raubenheimer D, Cogger VC, Macia L, Solon-Biet SM, Le Couteur DG, et al. The nutritional geometry of liver disease including non-alcoholic fatty liver disease. J Hepatol. 2018;68(2):316–25.
Saner C, Tassoni D, Harcourt BE, Kao KT, Alexander EJ, McCallum Z, et al. Evidence for protein leverage in children and adolescents with obesity. Obesity. 2020; In. Press.
Levasseur M. Sources et justifications des questions utilisées dans l’enquête santé québec. In: Cahier technique 87-03: Ministère de la Santé et des Services sociaux, Gouvernement du Québec; 1987.
Fillenbaum GG, Smyer MA. The development, validity, and reliability of the oars multidimensional functional assessment questionnaire1. J Gerontol. 1981;36(4):428–34.
Washburn RA, McAuley E, Katula J, Mihalko SL, Boileau RA. The physical activity scale for the elderly (pase): Evidence for validity. J Clin Epidemiol. 1999;52(7):643–51.
Presse N, Payette H, Shatenstein B, Greenwood CE, Kergoat MJ, Ferland G. A minimum of six days of diet recording is needed to assess usual vitamin k intake among older adults. J Nutr. 2010;141(2):341–6.
Allès B, Samieri C, Jutand M-A, Carmichael P-H, Shatenstein B, Gaudreau P, et al. Nutrient patterns, cognitive function, and decline in older persons: Results from the three-city and nuage studies. Nutrients. 2019;11(8):1808.
Shatenstein B, Gauvin L, Keller H, Richard L, Gaudreau P, Giroux F, et al. Individual and collective factors predicting change in diet quality over 3 years in a subset of older men and women from the nuage cohort. Eur J Nutr. 2016;55(4):1671–81.
Moshfegh M, Borrud L, Perloff B, LaComb R. Improved method for the 24-hour dietary recall for use in national surveys. FASEB J. 1999;13(4):A603.
Johnson-Down L, Ritter H, Starkey LJ, Gray-Donald K. Primary food sources of nutrients in the diet of canadian adults. Can J Diet Pract Res. 2006;67(1):7–13.
Ghachem A, Fried LP, Legault V, Bandeen-Roche K, Presse N, Gaudreau P, et al. Evidence from two cohorts for the frailty syndrome as an emergent state of parallel dysregulation in multiple physiological systems. Biogerontology. 2020.
Klemera P, Doubal S. A new approach to the concept and computation of biological age. Mech Ageing Dev. 2006;127(3):240–8.
R-Development-Core-Team: R: A language and environemnt for statistical computing. In., 4.0.3 edn. Available at http://www.r-project.org. 2020.
Hastie TJ, Tibshirani RJ. Generalized additive models. London: Chapman and Hall/CRC; 1990.
Benjamani Y, Hochberg Y. Controlling the false discovery rate: A practical and powerful approach to multiple testing. J Royl Stat Soc Ser B. 1995;57(1):289–300.
We would like to thank members of the ageing physiology research group at the University of Sherbrooke, as well as Dr Rosilene Ventura Ribeiro at the Charles Perkins Centre, for their thoughtful discussion about this work.
AMS was supported by a discovery early career researcher award from the Australian Research Council (ARC DECRA: DE180101520). The NuAge cohort was supported by the Canadian Institutes of Health Research (CIHR) grant #62842. The NuAge Database and Biobank are supported by the Fonds de recherche du Québec (FRQ) grant #2020-VICO-279753, the Quebec Network for Research on Aging, a thematic network funded by the FRQ - Santé (FRQ-S) and by the Merck-Frosst Chair funded by La Fondation de l’Université de Sherbrooke. This work was supported by the Canadian Institutes of Health Research (CIHR, grant #153011). AAC is supported by a CIHR New Investigator Salary Award and FRQ-S Senior Fellowship and is a member of the FRQ-S funded Centre de recherche du CHUS and Centre de recherche sur le vieillissement. NP is a Junior 1 Research Scholar of the FRQ-S.
Ethics approval and consent to participate
The NuAge Database and Biobank and the present study have been approved by the Research Ethics Board (REB) of the CIUSSS-de-l’Estrie-CHUS (Quebec, Canada; projects 2019-2832 and 2015-868/14-141, respectively). All participants gave written consent for their data and samples to be included in the NuAge Database and Biobank.
Consent for publication
This manuscript does not include details, images, or videos relating to an individual person.
AAC declares a conflict of interest as founder and CEO at Oken Health. The other authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Contains the following. Text S1. Model Descriptions. Text S2. Analysis of Inclusive Dataset with Imputation of Missing Income Data. Text S3. Analysis of Exclusive Dataset. Text S4. Estimation of Relative Intake. Text S5. Dietary Macronutrient Composition and Micronutrient Intake. Table S1. Variables and Summary Statistics. Tables S2 though S5, S8 & S9. GAM outputs. Table S6. Principal Component Analysis. Table S7. Models with Significant Effects of Micronutrients. Table S10 & S11. Model AICs. Figure S1. & S3. Effects of Relative Macronutrient Intake. Figure S2 & S4. Effects of Dietary Micronutrients. Figure S5. Effects of Diet Macronutrient Composition on Micronutrient intake. Figure S6. Standard Errors for Effects of Macronutrients.
About this article
Cite this article
Senior, A.M., Legault, V., Lavoie, F.B. et al. Multidimensional associations between nutrient intake and healthy ageing in humans. BMC Biol 20, 196 (2022). https://doi.org/10.1186/s12915-022-01395-z
- Geometric framework
- Systems epidemiology