Next Article in Journal
Maternal Dietary Patterns and Gestational Diabetes Mellitus in a Multi-Ethnic Asian Cohort: The GUSTO Study
Previous Article in Journal
Antioxidant-Rich Extract from Plantaginis Semen Ameliorates Diabetic Retinal Injury in a Streptozotocin-Induced Diabetic Rat Model
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Nutrient Intake Is Associated with Longevity Characterization by Metabolites and Element Profiles of Healthy Centenarians

1
College of Light Industry and Food Engineering, Guangxi University, Nanning 530004, China
2
Institute of Agricultural Quality Standards and Testing Technology Research, Shandong Academy of Agricultural Sciences, Jinan 250100, China
3
College of Life Sciences, Tianjin Normal University, Tianjin 300387, China
4
Institute of Agro-Products Quality Testing, Jinan Academy of Agricultural Sciences, Jinan 250316, China
5
Food and Biological Technology Department, Guangxi Polytechnic College, Nanning 530226, China
6
Veterans Affairs Boston Healthcare System, Department of Neurology, Harvard Medical School, West Roxbury, MA 02132, USA
*
Authors to whom correspondence should be addressed.
Submission received: 15 July 2016 / Revised: 3 September 2016 / Accepted: 7 September 2016 / Published: 19 September 2016

Abstract

:
The relationships between diet and metabolites as well as element profiles in healthy centenarians are important but remain inconclusive. Therefore, to test the interesting hypothesis that there would be distinctive features of metabolites and element profiles in healthy centenarians, and that these would be associated with nutrient intake; the short chain fatty acids (SCFAs), total bile acids and ammonia in feces, phenol, p-cresol, uric acid, urea, creatinine and ammonia in urine, and element profiles in fingernails were determined in 90 healthy elderly people, including centenarians from Bama county (China)—a famous longevous region—and elderly people aged 80–99 from the longevous region and a non-longevous region. The partial least squares-discriminant analysis was used for pattern recognition. As a result, the centenarians showed a distinct metabolic pattern. Seven characteristic components closely related to the centenarians were identified, including acetic acid, total SCFA, Mn, Co, propionic acid, butyric acid and valeric acid. Their concentrations were significantly higher in the centenarians group (p < 0.05). Additionally, the dietary fiber intake was positively associated with butyric acid contents in feces (r = 0.896, p < 0.01), and negatively associated with phenol in urine (r = −0.326, p < 0.01). The results suggest that the specific metabolic pattern of centenarians may have an important and positive influence on the formation of the longevity phenomenon. Elevated dietary fiber intake should be a path toward health and longevity.

Graphical Abstract

1. Introduction

Health and longevity are the eternal pursuit of human beings. Diet is an important factor affecting health and longevity [1,2]. Appropriate dietary alterations, such as dietary restriction, can promote health and increase the life span of laboratory model organisms [3], but malnutrition, nutritional excess or imbalance accelerates the aging process and causes a variety of geriatric diseases in humans [4]. Besides, dietary compositions also influence the metabolism of an individual and contribute to a variation in metabolites.
Metabolite levels are important indicators of human health and longevity. Many diseases and inflammations are associated with metabolic imbalances [5]. Studies show that short chain fatty acids (SCFAs) (such as butyric acid) reduce the risk of gastrointestinal diseases and cancer [6], and also exert important immunomodulatory roles [7]. In addition, certain essential elements, such as Mn, Se and Zn, are critical to the human metabolism and aging [8,9]. Lv et al. reported that the higher Mn levels in the hair of centenarians were conducive to their long lifespan [10]. Besides, Se is associated with the prevention of some diseases such as cancer and Alzheimer's disease, and also has important impacts on longevity [11]. As for elderly people, appropriate levels of metabolites and elements in elderly people play important roles in the prevention of geriatric diseases and the maintenance of normal immune functions [12]. Therefore, monitoring the contents of relevant metabolites and elements is of great significance for human health and longevity. As centenarians represent the model of healthy and successful aging well [13], it is necessary to explore the features of metabolites and element profiles in healthy centenarians. Nevertheless, the relevant studies on the metabolites and element profiles in healthy centenarians are still limited. Furthermore, previous studies also have not specifically assessed the correlation between metabolites and nutrient intake of healthy elderly people over 80 years of age.
It is noteworthy that human health and longevity are associated with overall metabolites and element levels rather than individual ones. To comprehensively analyze the overall influence of various indexes on health and longevity, it is necessary to use a multivariate statistical analysis method to perform pattern recognition. Partial least squares-discriminant analysis (PLS-DA) can be applied to construct the classification model and identify characteristic components. More importantly, PLS-DA has advantages in solving the classification problem of a very small number of samples [14]. Thus, considering the limited sample size of this study, PLS-DA is used to analyze the comprehensive impact of metabolite and element levels on health and longevity. So far, the studies using PLS-DA for pattern recognition mainly focus on investigating the relationship between specific diseases and levels of metabolites or elements in the human body [15,16]. However, no work has been reported to use PLS-DA for comprehensive evaluation of the characteristics of metabolites and elements in healthy centenarians.
There is a marvelous phenomenon in Bama county, Guangxi province, China, where the ratio of centenarians was the highest in China. According to the Population Census of China in 2010 [17], there were 80 centenarians in the population of 224,637, reaching to a ratio of 35.6 centenarians per 1 × 105 persons [18]. The ratio of centenarians in this region far exceeds the world longevity county standard specified by the United Nations (7.5/100,000) [19]. Nevertheless, few researches have been focused on the metabolites and element profiles of elderly people living in Bama longevous region (LR) because of their traditional and conservative lifestyle.
In this study, based on the rare phenomenon of healthy aging and extreme longevity, the hypothesis that there would be remarkable and distinctive features of metabolites and element profiles in healthy centenarians from Bama LR, and that these would be associated with nutrient intake, should be of interest. To test this hypothesis, we enrolled 49 healthy centenarians in the LR (LRC group), 59 healthy elderly people aged 80–99 in the LR (LRE group) and 61 healthy elderly people aged 80–99 from a non-longevous region (NLR) (NLRE group). The levels of metabolites and elements in the three groups were analyzed. PLS-DA was used to further explore the overall impact of the metabolites and element levels on health and longevity, and to identify the characteristic components closely related to the centenarians. Moreover, we assessed the relationships between metabolites or element concentrations and dietary nutrient intake. This will be very important for a better understanding of the relationship between diet and longevity from the perspective of metabolism, which will provide a theoretical basis for a healthy diet.

2. Materials and Methods

2.1. Participants

The study was conducted in the Bama district and Xixiangtang district of the Guangxi province, China, during the period from 2013 to 2015. The Xixiangtang district was selected as the non-longevous region (NLR), where the ratio of centenarians was relatively lower [18] (only three centenarians in the population of 1,156,173 [17], a ratio of 0.26 centenarians per 1 × 105 persons), and the climatic environment was similar with the Bama district. By means of a distribution of flyers and information sessions in the residential community, and thorough community screenings according to the population census data, we enrolled healthy centenarians in the LR (LRC group), and elderly people aged 80–99 in the LR (LRE group) and the NLR (NLRE group). The age of participants was validated by the following steps. First, the names, dates of birth, and addresses of participants were obtained by the population census data. Second, the participants were required to provide a valid identification card or household register booklet, which lists the name, date of birth, and other important demographic information. Third, the Twelve Animal Year of birth according to the traditional Chinese calendar should match the official year of birth [20]. Fourth, the age of participants was also checked by using the generation of family and the familial reconstitution method. Finally, neighbors of participants and villagers were also consulted to confirm the age of participants. To guarantee that each participant was a relatively healthy elderly person, the following screening criteria were adopted. Ineligibility criteria included uncontrolled hypertension, unstable type 1 diabetes, major cardiovascular diseases, any form of cancer, severe dementia or dysphasia, hepatitis, renal impairment, endocrine disease, inflammatory diseases, compromised immune status (AIDS, multiple myeloma, chronic lymphatic leukaemias), a physical disability that limits walking or dependency on others for activities of daily living, cognitive impairment, gastrointestinal disorders and/or food allergies. In addition, no antibiotics or drugs known to affect experimental endpoints, or nutritional supplements were taken during the course of this study. During recruitment and screening, each participant was required to complete a questionnaire on their medical history to provide health information.
Forty-nine centenarians were enrolled in the LRC group, and 59 elderly people from the LR and 61 elderly people from the NLR were enrolled in the LRE group and NLRE group, respectively. Among these volunteers, 37 centenarians in the LRC group, 43 elderly people in the LRE group and 40 elderly people in the NLRE group met the screening criteria and started the study. Of these, 30 centenarians in the LRC group (11 male/19 female, age: 103 ± 3 years), 30 elderly people in the LRE group (12 male/18 female, age: 87 ± 5 years), and 30 elderly people in the NLRE group (13 male/17 female, age: 88 ± 4 years) completed the study protocol (Figure 1). The characteristics of the participants are shown in Table 1. There is no significant difference in age between the LRE group and NLRE group (p > 0.05), as well as in sex ratio among the three groups (p > 0.05).
All study procedures were reviewed and approved by the Institutional Ethics Review Board of Guangxi University prior to initiation of this study. Written informed consent was obtained from each participant.

2.2. Assessment of Dietary Nutrition Status

The assessment of dietary nutrition status was carried out by means of four season consecutive 7-day weighed dietary records (28-day WDRs) [21], which were conducted in January, April, July, and October. During the dietary survey, all of the participants were required to adopt the mode of individual dining, and were advised to maintain usual eating habits. All food and drinks consumed by the participants were weighed and recorded by the trained investigators using electronic food scales, measuring cups, and spoons. For unmeasurable food, standard portion size was used. For the purposes of quality control, ten percent of the dietary records were reviewed and checked by the trained dietitians on a dietary survey site. The dietary data was transformed to dietary nutrient intake based on the Chinese food composition tables [22]. Average daily nutrient intake was calculated by multiplying the quantities of food consumed (in g) or portion size by the nutrient contents per 100 g of food listed in the Chinese food composition tables [23].

2.3. Sample Collection and Preparation

The fecal, urine and fingernail samples were collected from the participants in January, April, July, and October. Fresh fecal samples and morning urine samples were maintained at 4 °C for ≤5 h before processing, and then stored at −80 °C until analysis. Fingernail samples were stored in clean separate bags.

2.4. Analysis of Metabolites in Feces

SCFAs including acetic acid, propionic acid, butyric acid, isobutyric acid, valeric acid and isovaleric acid in feces were analyzed as described previously [24]. After thawing at room temperature, samples were diluted, vortexed and then centrifuged. Supernatant fractions were then acidified with orthophosphoric acid. Individual SCFAs were separated and quantified by capillary GC (GC-2010 Plus, Shimadzu, Kyoto, Japan). The sum of the concentrations of acetic acid, propionic acid, butyric acid, isobutyric acid, valeric acid and isovaleric acid was calculated as the total SCFA concentration.
The contents of total bile acids in feces were measured by the enzymatic cycling method using a commercial kit (Total Bile Acids Kit; Mindray Bio-Medical Corp., Shenzhen, China) following the manufacturer’s instructions.
Ammonia contents in feces were determined using the indophenol blue procedure as described previously [25]. Fecal specimens were diluted to the appropriate concentration and centrifuged at 13,000× g for 20 min. The supernatant fractions were reacted with aqueous phenol plus sodium nitroprusside solution and alkaline sodium hypochlorite. The mixtures were heated for 15 min at 60 °C in a shaking water-bath. After cooling quickly to room temperature, the optical density (625 nm) of the endproduct (indophenol) formed was measured by colorimetry.

2.5. Analysis of Metabolites in Urine

Phenol and p-cresol in urine were analyzed as described previously [26]. After thawing at room temperature, the samples were centrifuged (9000× g, 20 min, at 4 °C). Supernatant fractions were acidified with hydrochloric acid, and then vortexed and heated for 60 min at 90 °C in a water-bath. After cooling quickly to room temperature, ethyl acetate was added for extraction. Phenol and p-cresol in the upper ethyl acetate phase were quantified by capillary GC.
The contents of uric acid in urine were analyzed by the uricase-peroxidase colorimetric method using a commercial kit (Uric Acid Kit; Mindray Bio-Medical Corp.) following the manufacturer’s instructions. Urea contents in urine were assayed by the urease-glutamic dehydrogenase colorimetric method using a commercial kit (Urea Kit; Mindray Bio-Medical Corp.) following the manufacturer’s instructions. Creatinine contents in urine were analyzed by the Jaffe colorimetric method using a commercial kit (Creatinine Kit; Mindray Bio-Medical Corp.) following the manufacturer’s instructions.
Ammonia contents in urine samples were also analyzed similarly using the abovementioned indophenol blue procedure.

2.6. Analysis of Element Profiles in Nails

The concentrations of elements including Na, Mg, K, Ca, Mn, Fe, Cu, Zn, As, Sn, Sb, Pb, Cr, Co, Ni, Se, Sr, and Ba in fingernails were determined by inductively coupled plasma mass spectrometry (ICP-MS) as described previously [27]. Nitric acid (65%) and hydrogen peroxide (30%) were added to fingernail samples, and then the microwave digestion was carried out. The element concentrations were determined by CRC-ICP-MS (Agilent 7700e ICP-MS with a collision/reaction cell) under optimized conditions.

2.7. Statistical Analysis

The results were presented as mean values with their standard deviation. The nonparametric Kruskal-Wallis test in SPSS 19.0 software (International Business Machines Corp., Armonk, NY, USA) was performed to analyze the statistical differences among the three groups. In addition, the SCFAs, total bile acids, phenol, p-cresol, uric acid, urea, creatinine, ammonia, and element contents were summarized by multivariate statistical analysis. PLS-DA was used to build the classification model among the three groups. PLS-DA was performed with the SIMCA-P 11.5 software package (Umetrics, Umea, Sweden). To minimize the influence of the heteroscedasticity resulting from the concentration variations of different components on the multivariate statistical model, the data set was mean-centered and scaled to unit variance prior to PLS-DA [28,29]. The values of the variable importance in the projection (VIP) obtained from the PLS-DA model were used to identify the differential components contributing to the variations among the three groups. In the VIP plot, each variable is provided with a VIP value and a 95% confidence interval (CI) derived from jack-knifing [30,31,32]. Variables, for which the VIP ± 95% CI exceeds 1, are designated as significant differential components in this analysis [31,32]. Consequently, the characteristic components closely related to centenarians from the LR were determined. For PLS-DA modeling, a typical sevenfold (Leave-1/7th Samples-Out) cross-validation procedure was performed to avoid model over-fitting. The critical p value was set as 0.05.
A Spearman correlation test was applied to evaluate correlations between the nutrient intake and metabolites as well as element levels. The statistical significance was set at p < 0.05 and p < 0.01.

3. Results

3.1. Nutrient Intake

The nutrient intake levels in the three groups are presented in Table 2. The intake of energy, protein, fat, cholesterol, riboflavin, nicotinic acid and sodium in the LRC group and LRE group were significantly lower than those in the NLRE group (p < 0.05). Of these nutrients, energy, protein, fat and cholesterol intake decreased by 20%, 27%, 38% and 54% respectively for the LRC group, and decreased by 19%, 31%, 42% and 48% respectively for the LRE group, compared with those for the NLRE group. However, it is noteworthy that the intake of dietary fiber and vitamin A increased 1.7-fold and 1.4-fold for the LRC group compared with those for the NLRE group (p < 0.05). In addition, the intake of vitamin B6, folic acid and magnesium in the LRC group was significantly higher than those in the LRE group (p < 0.05). Therefore, the results suggest that a lower intake of energy, fat, cholesterol and a higher intake of dietary fiber and vitamin A may have a positive effect on health and longevity. Interestingly, as for most of the nutrients, the intakes in the LRE group were more similar to those in the LRC group than to those in the NLRE group, probably due to the similar dietary habits and dietary pattern of elderly people in the LR. This suggests a major influence of geographical origin on the nutrient intakes of healthy elderly people.

3.2. Metabolites in Feces

The contents of acetic acid and total SCFA increased 1.4-fold and 1.3-fold for the LRC group compared with those for the LRE group (p < 0.05), and increased 1.8-fold and 1.9-fold for the LRE group compared with those for the NLRE group (p < 0.05) (Table 3). In addition, the contents of butyric acid and valeric acid in the LRC group were significantly higher than in the other two groups (p < 0.05). Butyric acid concentrations increased 2.1-fold and 3.4-fold for the LRC group, and valeric acid concentrations increased 1.4-fold and 1.8-fold for the LRC group, compared with those for the LRE group and NLRE group, respectively. Besides, the contents of propionic acid, isobutyric acid and isovaleric acid in the LRC group and LRE group were significantly higher than those in the NLRE group (p < 0.05). Additionally, the content of total bile acids in the LRC group was significantly higher than in the other two groups (p < 0.05). Based on the above mentioned indexes, the results showed that the contents of SCFAs and total bile acids in feces were significantly higher in the LRC group, indicating that the increase of fecal SCFAs and total bile acids has a positive impact on the health and longevity of the centenarians from the LR.

3.3. Metabolites in Urine

The contents of phenol, p-cresol and urinary ammonia in the LRC group were significantly lower than in the other two groups (p < 0.05) (Table 4). Phenol, p-cresol and urinary ammonia concentrations decreased by 32%, 33% and 40% respectively for the LRC group compared with those for the LRE group, and decreased by 48%, 40% and 42% respectively for the LRC group compared with those for the NLRE group. Additionally, the contents of uric acid in the LRC group and LRE group were significantly lower than that in the NLRE group (p < 0.05). The contents of creatinine in the LRC group and LRE group were significantly higher than that in the NLRE group (p < 0.05), and the creatinine contents in the three groups were within the reference range provided in the manufacturer’s instructions. Consequently, the results indicate that lower contents of phenol, p-cresol, ammonia and uric acid have a positive influence on the health and longevity of the centenarians from the LR.

3.4. Element Levels in Fingernails

The contents of 18 chemical elements in fingernails were determined using ICP-MS, as shown in Table 5. The contents of Mn, Fe, Cu and Co in the LRC group and LRE group were significantly higher than those in the NLRE group (p < 0.05). Mn, Fe, Cu and Co contents in the LRC group increased 8.6, 2.4, 1.2 and 4.3-folds respectively, compared with those for the NLRE group (p < 0.05). Additionally, Zn content in the LRC group was significantly higher than that in the NLRE group (p < 0.05), and the content increased 1.2-fold. Besides, Se content increased 1.8-fold for the LRC group compared with that for the LRE group (p < 0.05). In addition, Pb content in the LRC group was the lowest among three groups (p < 0.05). Therefore, based on the element profiles in fingernails, the results revealed that higher levels of Mn, Fe, Cu, Co, Zn and Se, as well as a lower level of Pb have a positive influence on the health and longevity of the centenarians from the LR.

3.5. Pattern Recognition Analysis

In order to comprehensively analyze the overall impact of metabolite and element levels on health and longevity, pattern recognition technology was used to build the classification model and to identify the characteristic components closely related to the centenarians from the LR. The PLS-DA scores plot revealed a clear separation among the three groups, as shown in Figure 2. This indicated that the three groups had distinct differences in the metabolites and element profiles, and the intergroup differences were more significant, though there was some intrinsic physiological variability in each person. Therefore, the results showed that there were distinctive features of metabolites and element profiles in the LRC group.
In addition, the bi-plot of PLS-DA model showed that the contents of some beneficial components, such as SCFAs, Mn, Co, Zn and Se, in the LRC group were higher than those in the other two groups (Figure 3). Meanwhile, the contents of some harmful components, such as p-cresol, phenol and ammonia, in the NLRE group were higher than those in the other two groups.
To identify which components were accountable for the separation, the VIP plot of the PLS-DA model was used to select the differential components among the three groups (Figure 4). VIP statistics ranked the overall contribution of each component to the PLS-DA model. The VIP analysis identified seven components, including acetic acid (VIP = 1.914), total SCFA (VIP = 1.848), Mn (VIP = 1.753), Co (VIP = 1.587), propionic acid (VIP = 1.447), butyric acid (VIP = 1.418) and valeric acid (VIP = 1.377), as the significant and reliable differential components contributing to the variations among the three groups.
Besides, the differential components were further verified by the V plot of the PLS-DA model, as shown in Figure 5. In the V plot, each triangle denoted an individual component. The triangles far away from the origin represented the components responsible for the differences among the three groups. The abovementioned differential components, including acetic acid, total SCFA, Mn, Co, propionic acid, butyric acid and valeric acid, were the farthest away from the origin in the V plot, which further indicated that they were the differential components obtained from the PLS-DA model.
The univariate statistical analysis (Kruskal-Wallis test) was used to validate the statistically significant difference of the differential components derived from multivariate statistical analysis. There were significant differences in acetic acid (p = 0.000), total SCFA (p = 0.000), Mn (p = 0.000), Co (p = 0.000), propionic acid (p = 0.000), butyric acid (p = 0.000) and valeric acid (p = 0.000) among the three groups; the concentrations of these differential components in the LRC group were relatively higher. Therefore, the seven significant differential components could serve as the characteristic components closely related to centenarians from the LR.

3.6. Correlation Analysis

In order to assess the relationships between the nutrient intake and metabolites, as well as element levels, the Spearman correlation test was performed (Figure 6). There was a significant positive correlation between dietary fiber intake and butyric acid contents in feces, and the correlation coefficient was highest (r = 0.896, p < 0.01). The other short chain fatty acids were also positively associated with dietary fiber intake (p < 0.05). In addition, the contents of the SCFAs were also significantly associated with the intake of energy, protein, vitamin A and nicotinic acid (p < 0.05). Besides, the contents of phenol were negatively associated with the intake of dietary fiber (r = −0.326, p < 0.01). Additionally, Mn contents in fingernails were significantly associated with the intake of energy, protein, fat, dietary fiber, cholesterol, vitamin A, thiamine, nicotinic acid, sodium, magnesium, zinc and selenium (p < 0.05). Co contents in fingernails were significantly associated with the intake of energy, dietary fiber, cholesterol, vitamin A, thiamine, nicotinic acid, sodium and magnesium (p < 0.05).

4. Discussion

In this study, we captured the characteristics of metabolites in feces and urine and element profiles in the fingernails of healthy centenarians from the LR using pattern recognition analysis combined with a univariate statistical test. As a result, the LRC group showed a distinct metabolic pattern. Seven characteristic components closely related to the centenarians from the LR were identified. The concentrations of the characteristic components, including acetic acid, total SCFA, Mn, Co, propionic acid, butyric acid and valeric acid, were significantly higher in the LRC group (p < 0.05). Additionally, the further correlation analysis revealed that there were some significant correlations between nutrient intake and metabolites. Therefore, these results confirm the original research hypothesis that there would be distinctive features of metabolites and element profiles in healthy centenarians from the LR, and that these would be associated with nutrient intake.
SCFAs in the gut perform various physiological functions [33,34], which are recognized as essential energy sources and act as signal transduction molecules via G-protein coupled receptors [35]. Butyric acid and propionic acid can activate intestinal gluconeogenesis, which promote metabolic benefits via gut-brain neural circuits [36]. Besides, SCFAs have anti-inflammatory, antitumorigenic, and antimicrobial effects [33], and are also associated with a lower risk for some diseases, such as irritable bowel syndrome, inflammatory bowel diseases, cardiovascular diseases and cancer [33,37]. By contrast, phenol, p-cresol and ammonia are potential cytotoxins and carcinogens, and are related to a greater risk of cancers [38]. However, the relevant reports on the metabolites including SCFAs, phenol, p-cresol, and ammonia of elderly people are rather limited, especially for centenarians. As for this study, we discovered that the levels of acetic acid, propionic acid, butyric acid, valeric acid and total SCFA in the LRC group were significantly higher (p < 0.05), which were identified as characteristic components closely related to the centenarians, whereas the contents of phenol, p-cresol and urinary ammonia in the LRC group were significantly lower than in the other two groups (p < 0.05). Together, the results indicate that for elderly people, the relatively higher levels of SCFAs and lower levels of phenol, p-cresol and ammonia are conducive to the longevity of the centenarians from the LR. More interestingly, as for most of the SCFAs in feces, the contents in the LRE group were more similar to those in the LRC group than to those in the NLRE group, which could be attributed, firstly, to the similar dietary habits, dietary pattern and lifestyle of elderly people in the LR; secondly, to the same geographical environment in the LR. This suggests an important influence of geographical origin on the metabolites of healthy elderly people. It has been reported that diet, lifestyle and other environmental factors impact on the metabolite levels [39,40] and give rise to regional metabolomic phenotypes [41].
As for the elements in fingernails, Mn is an essential cofactor of mitochondrial superoxide dismutase (SOD) which is a key antioxidant enzyme [42,43]. Moreover, Mn is known as an important anti-aging element [19]. Co plays a biologically essential role as a cofactor in a number of proteins, and also a component of the vitamin B12 complex [44]. In this study, we discovered that the contents of Mn and Co in the LRC groups were significantly higher (p < 0.05), which were identified as characteristic components closely related to the centenarians. Therefore, the results suggested that the two characteristic elements also provided an essential material guarantee for health and longevity of the centenarians from the LR.
More importantly, we comprehensively analyzed the overall effect of metabolites and elements using the PLS-DA method. This pattern recognition method can distinguish the ‘among-groups’ and ‘within-groups’ variability, and thereby extract the most relevant information in differences among groups to build the classification model [14,45]. Therefore, the PLS-DA model can maximize systematic variations among groups, and thus achieve global separation profile among groups [46]. Additionally, PLS-DA plays a unique role in the identification of the most important variables and handling of highly collinear data [47,48]. Moreover, the method has distinct advantages in constructing the classification model for a small number of samples, which can overcome the limitation of a small sample size [49,50,51]. Relevant studies have shown that PLS-DA is much more robust for a small sample size [45]. For example, the Irs1−/− model group (n = 5) and WT control group (n = 5) were clearly distinguished [51]. Therefore, based on the advantages of PLS-DA, the method was used for pattern recognition in the present study and the results of the analysis are reliable and efficient. PLS-DA results revealed a unique metabolic pattern of the centenarians from the LR, and seven characteristic components were identified that are closely related to the centenarians. In consideration of the important physiological functions of these characteristic components, the specific metabolic pattern of the LRC group should have an important and positive influence on health and longevity. Additionally, the results suggested that PLS-DA could be used to manifest the different metabolic patterns of different populations. In a previous study on urine metabolic profiling of northern Italian centenarians, Collino et al. found that compared with the elderly group, the contents of phenylacetylglutamine, p-cresol-sulfate, and 2-hydroxybenzoate in centenarians were significantly higher [13]. There are some differences from our results, indicating that the centenarians from Bama LR have a unique metabolic pattern. However, there are very few researches on using PLS-DA to depict the comprehensive feature of metabolites and elements in healthy centenarians.
More remarkably, we discovered that dietary fiber intake was positively associated with SCFAs (p < 0.05), especially butyric acid (p < 0.01), and negatively associated with phenol concentrations (p < 0.01). It is therefore concluded that an appropriate increase of dietary fibers in a daily diet is a path toward the formation of the longevity phenomenon. Nevertheless, to date, we have not yet found published data on associating diet with metabolites as well as elements in the fingernails of healthy elderly people over 80 years of age. A previous study reported that SCFAs are produced during bacterial fermentation of undigested carbohydrates, such as non-starch polysaccharides, dietary fibers and resistant starch [33]. β-Glucan, fructo-oligosaccharide, and a mixture of dietary fibers can be converted to butyric acid in vivo [52]. Besides, a study showed that the high-fat diet reduced formation of butyric acid in rats, and the dietary fiber counteracted the harmful effect and stimulated the butyric acid formation [53]. To date, most of the relevant studies focus on investigating the relationships between a specific diet and metabolites by means of dietary intervention [53,54,55]. However, the study on the relations between metabolites and nutrient intakes of healthy elderly people maintaining their usual eating patterns is rather limited. McOrist et al. discovered that the total SCFA and butyric acid concentrations in feces from eight healthy volunteers (age: 31–59 years, n = 8) were weakly associated with dietary fat, sugar and carbohydrate intake that were obtained by self-report [56]. There are some differences from our results, which reveals distinct characteristics of metabolites and nutrient intake of the healthy people over 80 years of age. In addition, the differences of the dietary assessment method also could lead to the differences in results. Among available nutrition assessment approaches, the WDR method is the most precise and is accepted as the gold standard [57], though it is time-consuming, expensive, and generally requires considerable commitment on the part of the volunteers [23]. In this study, four season consecutive 7-day WDRs were carried out to assess the habitual intake of nutrients, which represents real life intake and minimizes variances in dietary intake according to seasons and days.
The present study has several strengths, including the use of the four season consecutive 7-day WDRs method, which provides a more accurate and robust report of elderly people’s usual intake. Furthermore, we pooled data of metabolites and elements obtained in January, April, July, and October to minimize intra-individual variation according to seasons and days. Additionally, the three kinds of different biological samples, including feces, urine, and fingernails, were analyzed to comprehensively assess the features of metabolites and element profiles in healthy centenarians from different perspective. Lastly, seven characteristic components closely related to centenarians from the LR were discovered using the PLS-DA method which is robust and reliable for a small sample size.
However, some limitations of this preliminary exploratory study need to be noted. First, health and longevity also depend on many other factors, such as genes and environment factors [58], while our study is limited by its cross-sectional design, therefore a causal relation between metabolites or elements and longevity cannot be concluded. Nevertheless, this study discovered the unique metabolic pattern of the centenarians from the LR, which provided new clues for further research on the relationship between longevity and metabolism. Future studies are needed to explore the mechanisms of metabolites or elements impacting healthy aging and the important metabolic pathways affecting aging. Second, the detailed data regarding the smoking history and amount were unavailable, and hence it is necessary to exclude the confounding influence of smoking in future studies. Third, the rather limited number of healthy centenarians and elderly people aged 80–99 years, as well as the strict screening criteria, led to a relatively small sample size. Besides, because of the tedious process of the 28-day weighed dietary records method and the complex protocols of sample collection, many participants withdrew from the study. However, we applied the PLS-DA method to perform pattern recognition, which overcomes the limitation of a small sample size. Further studies including populations from other different longevous regions are required to confirm these findings.

5. Conclusions

In summary, to our knowledge, this is the only study on the relationships between the nutrient intake and metabolites of healthy elderly people over 80 years of age, and the only study to use the PLS-DA method to assess the overall effect of metabolites and elements. The results indicated that the unique metabolic pattern of the LRC group may have an important and positive influence on the formation of the longevity phenomenon. Elevated dietary fiber intake should be a path toward health and longevity. These findings will provide new clues for future research on nutritional metabonomics and element metabolism in healthy elderly people, which will provide precious information for future human metabolism and longevity studies.

Acknowledgments

This study was supported by the grants from the National Natural Science Foundation of China (No. 31371762) and the “Bagui Scholars Distinguished Professor” Special Project. We are deeply grateful to all the participants and their offspring for their collaboration in 28-day WDRs and sample collection.

Author Contributions

Q.L., M.Z. and X.D. conceived and designed the experiments; D.C., S.Z., F.C. and X.T. performed the experiments; D.C., Q.L., D.L., G.H. and Z.Z. recruited the subjects; D.C., Q.L., D.L., G.H., Z.Z. and D.L. performed the dietary survey; D.C. analyzed the data; D.C. wrote the paper; Q.L. and S.L. assisted with manuscript preparation.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Péter, S.; Saris, W.H.M.; Mathers, J.C.; Feskens, E.; Schols, A.; Navis, G.; Kuipers, F.; Weber, P.; Eggersdorfer, M. Nutrient status assessment in individuals and populations for healthy aging—Statement from an expert workshop. Nutrients 2015, 7, 10491–10500. [Google Scholar] [CrossRef] [PubMed]
  2. Shi, Z.; Zhang, T.; Byles, J.; Martin, S.; Avery, J.C.; Taylor, A.W. Food habits, lifestyle factors and mortality among oldest old Chinese: The Chinese longitudinal healthy longevity survey (CLHLS). Nutrients 2015, 7, 7562–7579. [Google Scholar] [CrossRef] [PubMed]
  3. Fontana, L.; Partridge, L.; Longo, V.D. Extending healthy life span-from yeast to humans. Science 2010, 328, 321–326. [Google Scholar] [CrossRef] [PubMed]
  4. Brunner, E.J.; Mosdol, A.; Witte, D.R.; Martikainen, P.; Stafford, M.; Shipley, M.J.; Marmot, M.G. Dietary patterns and 15-year risks of major coronary events, diabetes, and mortality. Am. J. Clin. Nutr. 2008, 87, 1414–1421. [Google Scholar] [PubMed]
  5. McNiven, E.M.; German, J.B.; Slupsky, C.M. Analytical metabolomics: Nutritional opportunities for personalized health. J. Nutr. Biochem. 2011, 22, 995–1002. [Google Scholar] [CrossRef] [PubMed]
  6. Wong, J.M.W.; de Souza, R.; Kendall, C.W.C.; Emam, A.; Jenkins, D.J.A. Colonic health: Fermentation and short chain fatty acids. J. Clin. Gastroenterol. 2006, 40, 235–243. [Google Scholar] [CrossRef] [PubMed]
  7. Nastasi, C.; Candela, M.; Bonefeld, C.M.; Geisler, C.; Hansen, M.; Krejsgaard, T.; Biagi, E.; Andersen, M.H.; Brigidi, P.; Ødum, N.; et al. The effect of short-chain fatty acids on human monocyte-derived dendritic cells. Sci. Rep. 2015, 5, 16148. [Google Scholar] [CrossRef] [PubMed]
  8. Hambidge, M. Biomarkers of trace mineral intake and status. J. Nutr. 2003, 133, 948S–955S. [Google Scholar] [PubMed]
  9. Rodrigues, J.L.; Batista, B.L.; Fillion, M.; Passos, C.J.; Mergler, D.; Barbosa, F., Jr. Trace element levels in whole blood of riparian villagers of the Brazilian Amazon. Sci. Total Environ. 2009, 407, 4168–4173. [Google Scholar] [CrossRef] [PubMed]
  10. Lv, J.; Wang, W.; Zhang, F.; Krafft, T.; Yuan, F.; Li, Y. Identification of human age using trace element concentrations in hair and the support vector machine method. Biol. Trace Elem. Res. 2011, 143, 1441–1450. [Google Scholar] [CrossRef] [PubMed]
  11. Papp, L.V.; Holmgren, A.; Khanna, K.K. Selenium and selenoproteins in health and disease. Antioxid. Redox Signal. 2010, 12, 793–795. [Google Scholar] [CrossRef] [PubMed]
  12. Savarino, L.; Granchi, D.; Ciapetti, G.; Cenni, E.; Ravaglia, G.; Forti, P.; Maioli, F.; Mattioli, R. Serum concentrations of zinc and selenium in elderly people: Results in healthy nonagenarians/centenarians. Exp. Gerontol. 2001, 36, 327–339. [Google Scholar] [CrossRef]
  13. Collino, S.; Montoliu, I.; Martin, F.P.; Scherer, M.; Mari, D.; Salvioli, S.; Bucci, L.; Ostan, R.; Monti, D.; Biagi, E.; et al. Metabolic signatures of extreme longevity in northern Italian centenarians reveal a complex remodeling of lipids, amino acids, and gut microbiota metabolism. PLoS ONE 2013, 8, e56564. [Google Scholar] [CrossRef]
  14. Boulesteix, A.L.; Strimmer, K. Partial least squares: A versatile tool for the analysis of high-dimensional genomic data. Brief. Bioinform. 2006, 8, 32–44. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Tan, B.; Qiu, Y.; Zou, X.; Chen, T.; Xie, G.; Cheng, Y.; Dong, T.; Zhao, L.; Feng, B.; Hu, X.; et al. Metabonomics identifies serum metabolite markers of colorectal cancer. J. Proteome Res. 2013, 12, 3000–3009. [Google Scholar] [CrossRef] [PubMed]
  16. Zhao, T.; Chen, T.; Qiu, Y.; Zou, X.; Li, X.; Su, M.; Yan, C.; Zhao, A.; Jia, W. Trace element profiling using inductively coupled plasma mass spectrometry and its application in an osteoarthritis study. Anal. Chem. 2009, 81, 3683–3692. [Google Scholar] [CrossRef] [PubMed]
  17. Shi, S.; Liu, S.; Huang, Q.; Huang, F.; Wei, Y. Guangxi Zhuang Autonomous Region 2010 Population Census; China Statistics Press: Beijing, China, 2012. [Google Scholar]
  18. Chu, H. The Spatial-Temporal Evolution and Influence Factors Analysis of Longevity Level in Southern China. Master’s Thesis, Central China Normal University, Wuhan, China, May 2015. [Google Scholar]
  19. Lv, J.; Wang, W.; Krafft, T.; Li, Y.; Zhang, F.; Yuan, F. Effects of several environmental factors on longevity and health of the human population of Zhongxiang, Hubei, China. Biol. Trace Elem. Res. 2011, 143, 702–716. [Google Scholar] [CrossRef] [PubMed]
  20. Wang, Z.; Zeng, Y.; Jeune, B.; Vaupel, J.W. Age validation of Han Chinese centenarians. Genus 1998, 54, 123–141. [Google Scholar] [PubMed]
  21. Tokudome, S.; Imaeda, N.; Tokudome, Y.; Fujiwara, N.; Nagaya, T.; Sato, J.; Kuriki, K.; Ikeda, M.; Maki, S. Relative validity of a semi-quantitative food frequency questionnaire versus 28 day weighed diet records in Japanese female dietitians. Eur. J. Clin. Nutr. 2001, 55, 735–742. [Google Scholar] [CrossRef] [PubMed]
  22. Yang, Y.; Wang, G.; Pan, X. China Food Composition; The Peking University Medical Press: Beijing, China, 2009. [Google Scholar]
  23. Zhuang, M.; Yuan, Z.; Lin, L.; Hu, B.; Wang, X.; Yang, Y.; Chen, X.; Jin, L.; Lu, M.; Ye, W. Reproducibility and relative validity of a food frequency questionnaire developed for adults in Taizhou, China. PLoS ONE 2012, 7, e48341. [Google Scholar] [CrossRef] [PubMed]
  24. Brinkworth, G.D.; Noakes, M.; Clifton, P.M.; Bird, A.R. Comparative effects of very low-carbohydrate, high-fat and high-carbohydrate, low-fat weight-loss diets on bowel habit and faecal short-chain fatty acids and bacterial populations. Br. J. Nutr. 2009, 101, 1493–1502. [Google Scholar] [CrossRef] [PubMed]
  25. Bird, A.R.; Vuaran, M.S.; King, R.A.; Noakes, M.; Keogh, J.; Morell, M.K.; Topping, D.L. Wholegrain foods made from a novel high-amylose barley variety (Himalaya 292) improve indices of bowel health in human subjects. Br. J. Nutr. 2008, 99, 1032–1040. [Google Scholar] [CrossRef] [PubMed]
  26. Ministry of Public Health of China. Determination of Phenol in Urine by Gas Chromatographic Method. WS/T 50-1996; China Standards Press: Beijing, China, 1996. [Google Scholar]
  27. Samanta, G.; Sharma, R.; Roychowdhury, T.; Chakraborti, D. Arsenic and other elements in hair, nails, and skin-scales of arsenic victims in West Bengal, India. Sci. Total Environ. 2004, 326, 33–47. [Google Scholar] [CrossRef] [PubMed]
  28. Trygg, J.; Holmes, E.; Lundstedt, T. Chemometrics in metabonomics. J. Proteome Res. 2007, 6, 469–479. [Google Scholar] [CrossRef] [PubMed]
  29. Wold, S.; Trygg, J.; Berglund, A.; Antti, H. Some recent developments in PLS modeling. Chemometr. Intell. Lab. Syst. 2001, 58, 131–150. [Google Scholar] [CrossRef]
  30. Wiklund, S.; Johansson, E.; Sjostrom, L.; Mellerowicz, E.J.; Edlund, U.; Shockcor, J.P.; Gottfries, J.; Moritz, T.; Trygg, J. Visualization of GC/TOF-MS-based metabolomics data for identification of biochemically interesting compounds using OPLS class models. Anal. Chem. 2008, 80, 115–122. [Google Scholar] [CrossRef] [PubMed]
  31. Rios-Avila, L.; Coats, B.; Chi, Y.Y.; Midttun, O.; Ueland, P.M.; Stacpoole, P.W.; Gregory, J.F. Metabolite profile analysis reveals association of vitamin B-6 with metabolites related to one-carbon metabolism and tryptophan catabolism but not with biomarkers of inflammation in oral contraceptive users and reveals the effects of oral contraceptives on these processes. J. Nutr. 2015, 145, 87–95. [Google Scholar] [PubMed]
  32. Da Silva, V.R.; Rios-Avila, L.; Lamers, Y.; Ralat, M.A.; Midttun, O.; Quinlivan, E.P.; Garrett, T.J.; Coats, B.; Shankar, M.N.; Percival, S.S.; et al. Metabolite profile analysis reveals functional effects of 28-day vitamin B-6 restriction on one-carbon metabolism and tryptophan catabolic pathways in healthy men and women. J. Nutr. 2013, 143, 1719–1727. [Google Scholar] [CrossRef] [PubMed]
  33. Tan, J.; McKenzie, C.; Potamitis, M.; Thorburn, A.N.; Mackay, C.R.; Macia, L. The role of short-chain fatty acids in health and disease. Adv. Immunol. 2014, 121, 91–119. [Google Scholar] [PubMed]
  34. Boets, E.; Deroover, L.; Houben, E.; Vermeulen, K.; Gomand, S.V.; Delcour, J.A.; Verbeke, K. Quantification of in vivo colonic short chain fatty acid production from inulin. Nutrients 2015, 7, 8916–8929. [Google Scholar] [CrossRef] [PubMed]
  35. Kasubuchi, M.; Hasegawa, S.; Hiramatsu, T.; Ichimura, A.; Kimura, I. Dietary gut microbial metabolites, short-chain fatty acids, and host metabolic regulation. Nutrients 2015, 7, 2839–2849. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Vadder, F.D.; Kovatcheva-Datchary, P.; Goncalves, D.; Vinera, J.; Zitoun, C.; Duchampt, A.; Backhed, F.; Mithieux, G. Microbiota-generated metabolites promote metabolic benefits via gut-brain neural circuits. Cell 2014, 156, 84–96. [Google Scholar] [CrossRef] [PubMed]
  37. Thibault, R.; Blachier, F.; Darcy-Vrillon, B.; de Coppet, P.; Bourreille, A.; Segain, J.P. Butyrate utilization by the colonic mucosa in inflammatory bowel diseases: A transport deficiency. Inflamm. Bowel Dis. 2010, 16, 684–695. [Google Scholar] [CrossRef] [PubMed]
  38. Bingham, S.A. Diet and colorectal cancer prevention. Biochem. Soc. Trans. 2000, 28, 12–16. [Google Scholar] [CrossRef] [PubMed]
  39. Xie, G.; Li, X.; Li, H.; Jia, W. Toward personalized nutrition: comprehensive phytoprofiling and metabotyping. J. Proteome Res. 2013, 12, 1547–1559. [Google Scholar] [CrossRef] [PubMed]
  40. Walsh, M.C.; Brennan, L.; Pujos-Guillot, E.; Sébédio, J.L.; Scalbert, A.; Fagan, A.; Higgins, D.G.; Gibney, M.J. Influence of acute phytochemical intake on human urinary metabolomic profiles. Am. J. Clin. Nutr. 2007, 86, 1687–1693. [Google Scholar] [PubMed]
  41. Holmes, E.; Loo, R.L.; Stamler, J.; Bictash, M.; Yap, I.K.; Chan, Q.; Ebbels, T.; de Iorio, M.; Brown, I.J.; Veselkov, K.A.; et al. Human metabolic phenotype diversity and its association with diet and blood pressure. Nature 2008, 453, 396–400. [Google Scholar] [CrossRef] [PubMed]
  42. Miller, A.F. Superoxide dismutases: Ancient enzymes and new insights. FEBS Lett. 2012, 586, 585–595. [Google Scholar] [CrossRef] [PubMed]
  43. Carneiro, M.F.H.; Grotto, D.; Batista, B.L.; Rhoden, C.R.; Barbosa, F., Jr. Background values for essential and toxic elements in children’s nails and correlation with hair levels. Biol. Trace Elem. Res. 2011, 144, 339–350. [Google Scholar] [CrossRef] [PubMed]
  44. Okamoto, S.; Eltis, L.D. The biological occurrence and trafficking of cobalt. Metallomics 2011, 3, 963–970. [Google Scholar] [CrossRef] [PubMed]
  45. Vignon, M. Inference in morphological taxonomy using collinear data and small sample sizes: Monogenean sclerites (Platyhelminthes) as a case study. Zool. Scr. 2011, 40, 306–316. [Google Scholar] [CrossRef]
  46. Su, M.; Zheng, X.Y.; Zhang, T.; Pei, L.; Wang, F.; Zheng, X.; Gu, X.; Song, X.; Lu, X.; Chen, G.; et al. Integrated profiling of metabolites and trace elements reveals a multifaceted malnutrition in pregnant women from a region with a high prevalence of congenital malformations. Metabolomics 2012, 8, 831–844. [Google Scholar] [CrossRef]
  47. Brereton, R.G. Pattern recognition in chemometrics. Chemometr. Intell. Lab. Syst. 2015, 149, 90–96. [Google Scholar] [CrossRef]
  48. Gromski, P.S.; Muhamadali, H.; Ellis, D.I.; Xu, Y.; Correa, E.; Turner, M.L.; Goodacre, R. A tutorial review: Metabolomics and partial least squares-discriminant analysis—A marriage of convenience or a shotgun wedding. Anal. Chim. Acta 2015, 879, 10–23. [Google Scholar] [CrossRef] [PubMed]
  49. Gou, X.; Tao, Q.; Feng, Q.; Peng, J.; Sun, S.; Cao, H.; Zheng, N.; Zhang, Y.; Hu, Y.; Liu, P. Urinary metabonomics characterization of liver fibrosis induced by CCl4 in rats and intervention effects of Xia Yu Xue Decoction. J. Pharmaceut. Biomed. Anal. 2013, 74, 62–65. [Google Scholar] [CrossRef] [PubMed]
  50. Selman, C.; Kerrison, N.D.; Cooray, A.; Piper, M.D.W.; Lingard, S.J.; Barton, R.H.; Schuster, E.F.; Blanc, E.; Gems, D.; Nicholson, J.K.; et al. Coordinated multitissue transcriptional and plasma metabonomic profiles following acute caloric restriction in mice. Physiol. Genom. 2006, 27, 187–200. [Google Scholar] [CrossRef] [PubMed]
  51. Wijeyesekera, A.; Selman, C.; Barton, R.H.; Holmes, E.; Nicholson, J.K.; Withers, D.J. Metabotyping of long-lived mice using 1H NMR spectroscopy. J. Proteome Res. 2012, 11, 2224–2235. [Google Scholar] [CrossRef] [PubMed]
  52. Nilsson, U.; Nyman, M. Short-chain fatty acid formation in the hindgut of rats fed oligosaccharides varying in monomeric composition, degree of polymerisation and solubility. Br. J. Nutr. 2005, 94, 705–713. [Google Scholar] [CrossRef] [PubMed]
  53. Jakobsdottir, G.; Xu, J.; Molin, G.; Ahrne, S.; Nyman, M. High-fat diet reduces the formation of butyrate, but increases succinate, inflammation, liver fat and cholesterol in rats, while dietary fibre counteracts these effects. PLoS ONE 2013, 8, e80476. [Google Scholar] [CrossRef] [PubMed]
  54. Bird, A.R.; Flory, C.; Davies, D.A.; Usher, S.; Topping, D.L. A novel barley cultivar (Himalaya 292) with a specific gene mutation in starch synthase lla raises large bowel starch and short-chain fatty acids in rats. J. Nutr. 2004, 134, 831–835. [Google Scholar] [PubMed]
  55. Bird, A.R.; Hayakawa, T.; Marsono, Y.; Gooden, J.M.; Record, I.R.; Correll, R.L.; Topping, D.L. Coarse brown rice increases fecal and large bowel short-chain fatty acids and starch but lowers calcium in the large bowel of pigs. J. Nutr. 2000, 130, 1780–1787. [Google Scholar] [PubMed]
  56. McOrist, A.L.; Abell, G.C.; Cooke, C.; Nyland, K. Bacterial population dynamics and faecal short-chain fatty acid (SCFA) concentrations in healthy humans. Br. J. Nutr. 2008, 100, 138–146. [Google Scholar] [CrossRef] [PubMed]
  57. Tokudome, Y.; Imaeda, N.; Nagaya, T.; Ikeda, M.; Fujiwara, N.; Sato, J.; Kuriki, K.; Kikuchi, S.; Maki, S.; Tokudome, S. Daily, weekly, seasonal, within- and between-individual variation in nutrient intake according to four season consecutive 7 day weighed diet records in Japanese female dietitians. J. Epidemiol. 2002, 12, 85–92. [Google Scholar] [CrossRef] [PubMed]
  58. Davinelli, S.; Willcox, D.C.; Scapagnini, G. Extending healthy ageing: nutrient sensitive pathway and centenarian population. Immun. Ageing 2012, 9, 9. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Flow diagram of participants during the study.
Figure 1. Flow diagram of participants during the study.
Nutrients 08 00564 g001
Figure 2. Scores plot from the partial least squares-discriminant analysis of overall indexes including the metabolites and elements in the three groups. Red symbols represent the participants in the LRC group. Green symbols represent the participants in the LRE group. Blue symbols represent the participants in the NLRE group. Their spatial distribution reveals the variations of the metabolites and element profiles among the three groups. The three groups exhibit a clear separation, indicating that there are differences in the metabolites and element profiles among the three groups. LRC, centenarians from longevous region; LRE, elderly people aged 80–99 years from longevous region; NLRE, elderly people aged 80–99 years from a non-longevous region.
Figure 2. Scores plot from the partial least squares-discriminant analysis of overall indexes including the metabolites and elements in the three groups. Red symbols represent the participants in the LRC group. Green symbols represent the participants in the LRE group. Blue symbols represent the participants in the NLRE group. Their spatial distribution reveals the variations of the metabolites and element profiles among the three groups. The three groups exhibit a clear separation, indicating that there are differences in the metabolites and element profiles among the three groups. LRC, centenarians from longevous region; LRE, elderly people aged 80–99 years from longevous region; NLRE, elderly people aged 80–99 years from a non-longevous region.
Nutrients 08 00564 g002
Figure 3. Bi-plot from the partial least squares-discriminant analysis of overall indexes including the metabolites and elements in the three groups. Red triangles represent the participants in the LRC group. Green squares represent the participants in the LRE group. Blue rhombuses represent the participants in the NLRE group. Their spatial distribution reveals the variations of the metabolites and element profiles among the three groups. The three groups exhibit a clear separation, indicating that there are differences in the metabolites and element profiles among the three groups. Moreover, the bi-plot showed that the contents of some beneficial components, such as short chain fatty acids (SCFAs), Mn, Co, Zn and Se, in the LRC group were higher than those in the other two groups. Meanwhile, the contents of some harmful components, such as p-cresol, phenol and ammonia, in the NLRE group were higher. LRC, centenarians from longevous region; LRE, elderly people aged 80–99 years from longevous region; NLRE, elderly people aged 80–99 years from a non-longevous region; Acet, acetic acid; Prop, propionic acid; Buty, butyric acid; Isobu, isobutyric acid; Vale, valeric acid; Isova, isovaleric acid; T-SCFA, total SCFA; TBA, total bile acids; Phen, phenol; p-Cre, p-cresol; UA, uric acid; CREA, creatinine; Fec Amm, fecal ammonia; Uri Amm, urinary ammonia.
Figure 3. Bi-plot from the partial least squares-discriminant analysis of overall indexes including the metabolites and elements in the three groups. Red triangles represent the participants in the LRC group. Green squares represent the participants in the LRE group. Blue rhombuses represent the participants in the NLRE group. Their spatial distribution reveals the variations of the metabolites and element profiles among the three groups. The three groups exhibit a clear separation, indicating that there are differences in the metabolites and element profiles among the three groups. Moreover, the bi-plot showed that the contents of some beneficial components, such as short chain fatty acids (SCFAs), Mn, Co, Zn and Se, in the LRC group were higher than those in the other two groups. Meanwhile, the contents of some harmful components, such as p-cresol, phenol and ammonia, in the NLRE group were higher. LRC, centenarians from longevous region; LRE, elderly people aged 80–99 years from longevous region; NLRE, elderly people aged 80–99 years from a non-longevous region; Acet, acetic acid; Prop, propionic acid; Buty, butyric acid; Isobu, isobutyric acid; Vale, valeric acid; Isova, isovaleric acid; T-SCFA, total SCFA; TBA, total bile acids; Phen, phenol; p-Cre, p-cresol; UA, uric acid; CREA, creatinine; Fec Amm, fecal ammonia; Uri Amm, urinary ammonia.
Nutrients 08 00564 g003
Figure 4. Variable importance in the projection (VIP) plot from the partial least squares-discriminant analysis of overall indexes including the metabolites and elements in the three groups. Variables for which the VIP ± 95% confidence interval (CI) exceeds 1 are designated as significant differential components, including acetic acid, total SCFA, Mn, Co, propionic acid, butyric acid and valeric acid. Acet, acetic acid; Prop, propionic acid; Buty, butyric acid; Isobu, isobutyric acid; Vale, valeric acid; Isova, isovaleric acid; T-SCFA, total SCFA; TBA, total bile acids; Phen, phenol; p-Cre, p-cresol; UA, uric acid; CREA, creatinine; Fec Amm, fecal ammonia; Uri Amm, urinary ammonia.
Figure 4. Variable importance in the projection (VIP) plot from the partial least squares-discriminant analysis of overall indexes including the metabolites and elements in the three groups. Variables for which the VIP ± 95% confidence interval (CI) exceeds 1 are designated as significant differential components, including acetic acid, total SCFA, Mn, Co, propionic acid, butyric acid and valeric acid. Acet, acetic acid; Prop, propionic acid; Buty, butyric acid; Isobu, isobutyric acid; Vale, valeric acid; Isova, isovaleric acid; T-SCFA, total SCFA; TBA, total bile acids; Phen, phenol; p-Cre, p-cresol; UA, uric acid; CREA, creatinine; Fec Amm, fecal ammonia; Uri Amm, urinary ammonia.
Nutrients 08 00564 g004
Figure 5. V plot of the partial least squares-discriminant analysis of overall indexes including the metabolites and elements in the three groups. Each triangle denotes an individual component. The triangles far away from the origin represent the components responsible for the differences among the three groups. The acetic acid, total SCFA, Mn, Co, propionic acid, butyric acid and valeric acid are the farthest away from the origin in the V plot, which further indicates that they are the differential components obtained from the PLS-DA model. Acet, acetic acid; Prop, propionic acid; Buty, butyric acid; Isobu, isobutyric acid; Vale, valeric acid; Isova, isovaleric acid; T-SCFA, total SCFA; TBA, total bile acids; Phen, phenol; p-Cre, p-cresol; UA, uric acid; CREA, creatinine; Fec Amm, fecal ammonia; Uri Amm, urinary ammonia.
Figure 5. V plot of the partial least squares-discriminant analysis of overall indexes including the metabolites and elements in the three groups. Each triangle denotes an individual component. The triangles far away from the origin represent the components responsible for the differences among the three groups. The acetic acid, total SCFA, Mn, Co, propionic acid, butyric acid and valeric acid are the farthest away from the origin in the V plot, which further indicates that they are the differential components obtained from the PLS-DA model. Acet, acetic acid; Prop, propionic acid; Buty, butyric acid; Isobu, isobutyric acid; Vale, valeric acid; Isova, isovaleric acid; T-SCFA, total SCFA; TBA, total bile acids; Phen, phenol; p-Cre, p-cresol; UA, uric acid; CREA, creatinine; Fec Amm, fecal ammonia; Uri Amm, urinary ammonia.
Nutrients 08 00564 g005
Figure 6. Correlations between the nutrient intake and metabolites as well as element levels. Red denotes positive correlation, blue denotes negative correlation, and yellow denotes no correlation. There was a significant positive correlation between dietary fiber intake and butyric acid contents in feces, and the correlation coefficient was highest (r = 0.896, p < 0.01). The other short chain fatty acids also were positively associated with dietary fiber intake (p < 0.05). Besides, the dietary fiber intake was negatively associated with the contents of phenol in urine (r = −0.326, p < 0.01).
Figure 6. Correlations between the nutrient intake and metabolites as well as element levels. Red denotes positive correlation, blue denotes negative correlation, and yellow denotes no correlation. There was a significant positive correlation between dietary fiber intake and butyric acid contents in feces, and the correlation coefficient was highest (r = 0.896, p < 0.01). The other short chain fatty acids also were positively associated with dietary fiber intake (p < 0.05). Besides, the dietary fiber intake was negatively associated with the contents of phenol in urine (r = −0.326, p < 0.01).
Nutrients 08 00564 g006
Table 1. Participant characteristics 1.
Table 1. Participant characteristics 1.
CharacteristicsLRC GroupLRE GroupNLRE Group
Age (year)103 ± 387 ± 588 ± 4
Sex (M/F)11/1912/1813/17
Height (cm)145.9 ± 10.6150.8 ± 7.7157.0 ± 12.0
Weight (kg)43.1 ± 10.045.6 ± 6.559.1 ± 8.3
Body mass index (kg/m2)20.0 ± 2.820.1 ± 3.123.9 ± 1.4
1 Values are means ± standard deviation (SD). LRC, centenarians from longevous region; LRE, elderly people aged 80–99 years from longevous region; NLRE, elderly people aged 80–99 years from a non-longevous region; M, male; F, female.
Table 2. Nutrient intakes in the LRC group, LRE group and NLRE group 1.
Table 2. Nutrient intakes in the LRC group, LRE group and NLRE group 1.
NutrientLRC GroupLRE GroupNLRE Groupp
Energy (Kcal)1220.30 ± 134.60 a1237.20 ± 154.45 a1520.10 ± 215.62 b0.000
Protein (g)38.90 ± 7.39 a36.83 ± 8.61 a53.48 ± 14.84 b0.000
Fat (g)42.24 ± 15.78 a39.16 ± 13.02 a67.93 ± 25.00 b0.000
Carbohydrate (g)172.56 ± 20.54 a180.12 ± 24.91 a167.80 ± 32.05 a0.343
Dietary fiber (g)23.48 ± 8.26 a13.90 ± 6.21 b13.77 ± 5.86 b0.000
Cholesterol (mg)110.73 ± 71.64 a124.43 ± 121.97 a238.67 ± 130.94 b0.001
Vitamin A (μgRE)1308.37 ± 439.39 a1181.73 ± 370.05 a,b956.47 ± 496.79 b0.001
Thiamine (mg)0.48 ± 0.16 a0.57 ± 0.14 a,b0.60 ± 0.11 b0.002
Riboflavin (mg)0.61 ± 0.11 a0.62 ± 0.17 a0.84 ± 0.30 b0.007
Vitamin B6 (mg)0.18 ± 0.08 a0.13 ± 0.06 b0.17 ± 0.10 a,b0.040
Folic acid (μg)67.36 ± 41.75 a37.81 ± 30.09 b57.48 ± 23.88 a0.002
Nicotinic acid (mg)7.93 ± 2.30 a7.65 ± 1.86 a11.88 ± 3.38 b0.000
Vitamin C (mg)61.45 ± 20.02 a,b51.64 ± 19.06 a68.49 ± 23.20 b0.014
Vitamin E (mg)8.16 ± 3.65 a,b6.28 ± 3.41 a9.15 ± 3.95 b0.014
Calcium (mg)481.90 ± 87.48 a,b421.93 ± 127.94 a511.83 ± 158.33 b0.003
Phosphorus (mg)602.77 ± 74.36 a,b569.13 ± 113.71 a779.93 ± 223.26 b0.001
Potassium (mg)1433.00 ± 203.42 a,b1269.13 ± 197.30 a1546.33 ± 252.80 b0.000
Sodium (mg)1817.67 ± 222.17 a1854.26 ± 474.14 a2144.78 ± 471.59 b0.000
Magnesium (mg)354.73 ± 71.19 a283.73 ± 78.90 b276.95 ± 80.48 b0.000
Iron (mg)14.64 ± 4.40 a13.85 ± 3.69 a15.26 ± 3.11 a0.058
Zinc (mg)5.40 ± 1.35 a6.10 ± 1.74 a6.35 ± 1.96 a0.062
Selenium (μg)13.86 ± 5.36 a15.36 ± 4.82 a,b17.68 ± 4.10 b0.001
Copper (mg)2.08 ± 2.11 a3.09 ± 2.12 a2.06 ± 1.89 a0.067
Manganese (mg)3.19 ± 0.83 a3.60 ± 0.70 a3.68 ± 1.16 a0.145
1 Values are means ± SD. Values without a common superscript letter in a row were significantly different, p < 0.05. The p values are obtained by Kruskal-Wallis test. LRC, centenarians from longevous region; LRE, elderly people aged 80–99 years from longevous region; NLRE, elderly people aged 80–99 years from a non-longevous region.
Table 3. Metabolites in feces in the LRC group, LRE group and NLRE group 1.
Table 3. Metabolites in feces in the LRC group, LRE group and NLRE group 1.
MetaboliteLRC GroupLRE GroupNLRE Groupp
Acetic acid (µg/g)2539.47 ± 875.80 a1825.13 ± 527.79 b1016.17 ± 644.02 c0.000
Propionic acid (µg/g)875.53 ± 363.69 a830.80 ± 506.01 a326.67 ± 214.75 b0.000
Isobutyric acid (µg/g)195.03 ± 75.92 a186.18 ± 113.04 a109.51 ± 67.54 b0.000
Butyric acid (µg/g)780.61 ± 587.01 a365.33 ± 291.05 b226.99 ± 153.17 b0.000
Isovaleric acid (µg/g)358.19 ± 184.83 a388.27 ± 254.08 a185.15 ± 138.60 b0.000
Valeric acid (µg/g)223.48 ± 76.80 a157.93 ± 93.42 b121.77 ± 49.31 b0.000
Total SCFA (µg/g)4972.31 ± 1773.99 a3753.63 ± 1355.86 b1986.27 ± 1175.16 c0.000
Total bile acids (μmol/g)0.15 ± 0.06 a0.10 ± 0.04 b0.08 ± 0.07 b0.000
Fecal ammonia (mg/g)0.68 ± 0.28 a0.59 ± 0.22 a0.56 ± 0.33 a0.126
1 Values are means ± SD. Values without a common superscript letter in a row were significantly different, p < 0.05. The p values are obtained by Kruskal-Wallis test. LRC, centenarians from longevous region; LRE, elderly people aged 80–99 years from longevous region; NLRE, elderly people aged 80–99 years from a non-longevous region.
Table 4. Metabolites in urine in the LRC group, LRE group and NLRE group 1.
Table 4. Metabolites in urine in the LRC group, LRE group and NLRE group 1.
MetaboliteLRC GroupLRE GroupNLRE Groupp
Phenol (mg/L)13.32 ± 16.02 a19.60 ± 8.26 b25.79 ± 29.46 b0.003
p-Cresol (mg/L)65.52 ± 25.08 a98.26 ± 31.50 b109.52 ± 40.77 b0.000
Uric acid (μmol/L)974.63 ± 525.26 a830.97 ± 251.23 a1147.03 ± 191.99 b0.000
Urea (mmol/L)285.02 ± 132.62 a371.80 ± 106.73 b325.96 ± 151.43 ab0.027
Creatinine (μmol/L)7672.50 ± 5840.73 a5776.94 ± 1258.51 a3791.72 ± 1119.87 b0.000
Urinary ammonia (μg/μL)0.29 ± 0.15 a0.48 ± 0.17 b0.50 ± 0.35 b0.001
1 Values are means ± SD. Values without a common superscript letter in a row were significantly different, p < 0.05. The p values are obtained by Kruskal-Wallis test. LRC, centenarians from longevous region; LRE, elderly people aged 80–99 years from longevous region; NLRE, elderly people aged 80–99 years from a non-longevous region.
Table 5. Element concentrations in fingernails in the LRC group, LRE group and NLRE group 1.
Table 5. Element concentrations in fingernails in the LRC group, LRE group and NLRE group 1.
ElementLRC GroupLRE GroupNLRE Groupp
Na (μg/L)107.836 ± 106.206 a136.393 ± 122.729 a164.905 ± 127.408 a0.208
Mg (μg/L)61.143 ± 30.351 a62.514 ± 18.946 a53.848 ± 26.503 a0.332
K (μg/L)5.658 ± 8.982 a4.759 ± 7.180 a18.521 ± 48.124 a0.355
Ca (μg/L)869.077 ± 452.755 a850.038 ± 488.529 a650.867 ± 407.325 a0.170
Mn (μg/L)3.051 ± 2.294 a1.499 ± 0.822 a0.353 ± 0.368 b0.000
Fe (μg/L)38.153 ± 15.967 a40.478 ± 25.144 a16.162 ± 15.351 b0.000
Cu (μg/L)4.579 ± 1.409 a5.547 ± 2.592 a3.936 ± 3.153 b0.000
Zn (μg/L)137.406 ± 43.331 a118.979 ± 24.022 a,b112.877 ± 40.170 b0.012
As (μg/L)0.138 ± 0.091 a0.241 ± 0.315 a0.289 ± 0.521 a0.061
Sn (μg/L)1.029 ± 2.994 a0.567 ± 1.172 a0.866 ± 1.679 a0.730
Sb (μg/L)0.126 ± 0.167 a0.059 ± 0.071 a0.111 ± 0.131 a0.306
Pb (μg/L)0.140 ± 0.153 a0.311 ± 0.244 b0.279 ± 0.181 b0.000
Cr (μg/L)0.957 ± 1.058 a0.763 ± 0.824 a0.743 ± 1.268 a0.071
Co (μg/L)0.030 ± 0.015 a0.024 ± 0.017 a0.007 ± 0.010 b0.000
Ni (μg/L)2.130 ± 4.765 a0.410 ± 0.301 b0.785 ± 0.689 a,b0.009
Se (μg/L)0.476 ± 0.293 a0.267 ± 0.255 b0.315 ± 0.314 a,b0.020
Sr (μg/L)0.561 ± 0.391 a0.449 ± 0.325 a0.329 ± 0.246 a0.058
Ba (μg/L)0.879 ± 0.783 a0.745 ± 0.603 a0.636 ± 0.473 a0.605
1 Values are means ± SD. Values without a common superscript letter in a row were significantly different, p < 0.05. The p values are obtained by Kruskal-Wallis test. LRC, centenarians from longevous region; LRE, elderly people aged 80–99 years from longevous region; NLRE, elderly people aged 80–99 years from a non-longevous region.

Share and Cite

MDPI and ACS Style

Cai, D.; Zhao, S.; Li, D.; Chang, F.; Tian, X.; Huang, G.; Zhu, Z.; Liu, D.; Dou, X.; Li, S.; et al. Nutrient Intake Is Associated with Longevity Characterization by Metabolites and Element Profiles of Healthy Centenarians. Nutrients 2016, 8, 564. https://0-doi-org.brum.beds.ac.uk/10.3390/nu8090564

AMA Style

Cai D, Zhao S, Li D, Chang F, Tian X, Huang G, Zhu Z, Liu D, Dou X, Li S, et al. Nutrient Intake Is Associated with Longevity Characterization by Metabolites and Element Profiles of Healthy Centenarians. Nutrients. 2016; 8(9):564. https://0-doi-org.brum.beds.ac.uk/10.3390/nu8090564

Chicago/Turabian Style

Cai, Da, Shancang Zhao, Danlei Li, Fang Chang, Xiangxu Tian, Guohong Huang, Zhenjun Zhu, Dong Liu, Xiaowei Dou, Shubo Li, and et al. 2016. "Nutrient Intake Is Associated with Longevity Characterization by Metabolites and Element Profiles of Healthy Centenarians" Nutrients 8, no. 9: 564. https://0-doi-org.brum.beds.ac.uk/10.3390/nu8090564

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop