Research Paper Volume 16, Issue 3 pp 2638—2656

Serum metabolome and gut microbiome alterations are associated with low handgrip strength in older adults

Yan Guo1,2, *, , Qin Wang1,3, *, , Yifan Lv1, , Fan Xia1, , Xin Chen1, , Yan Mao1, , Xiaodong Wang4, , Guoxian Ding1, , Jing Yu1, ,

  • 1 Department of Geriatrics, Division of Geriatric Endocrinology, The First Affiliated Hospital of Nanjing Medical University, Nanjing, P.R. China
  • 2 Department of Neurology, Yancheng City No. 1 People’s Hospital, Yancheng, P.R. China
  • 3 Department of Geriatrics, Drum Tower Hospital Affiliated to Nanjing University Medical School, Nanjing, P.R. China
  • 4 Department of Endocrinology, The First Affiliated Hospital of Nanjing Medical University, Nanjing, P.R. China
* Equal contribution

Received: September 4, 2023       Accepted: January 3, 2024       Published: February 1, 2024      

https://doi.org/10.18632/aging.205501
How to Cite

Copyright: © 2024 Guo et al. This is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Abstract

Handgrip strength (HGS), which represents global muscle strength, is a powerful indicator of disability and mortality in older adults; it is also used for the diagnosis of possible- or probable- sarcopenia and physical frailty. This study aimed to explore the metabolic mechanisms and potential biomarkers associated with declining HGS among older adults. We recruited 15 age- and environment-matched inpatients (age, 77–90 years) with low or normal HGS. Liquid chromatography-mass spectrometry (LC-MS) and 16S ribosomal DNA (rDNA) gene sequencing were performed to analyze the metabolome of serum and stool samples and the gut microbiome composition of stool samples. Spearman’s correlation analysis was used to identify the potential serum and fecal metabolites associated with HGS. We assessed the levels of serum and fecal metabolites belonging to the class of cinnamic acids and derivatives and reported that the levels of carboxylic acids and their derivatives decreased in the low-HGS group. Serum levels of microbial metabolites, including cinnamoylglycine, 4-methoxycinnamic acid, and (e)-3,4,5-trimethoxycinnamic acid, were positively correlated with HGS. We found that gut microbial α-diversity was significantly higher in the low-HGS group, whereas higher β-diversity was observed in the normal group. The relative abundances of the genera Parabacteroides and Intestinibacter increased significantly in the low-HGS group and were negatively correlated with the serum levels of cinnamoylglycine. The identified metabolites whose levels were markedly altered, and intestinal flora associated with these metabolites suggest the potential metabolic underpinnings for HGS and provide a basis for the further identification of biomarkers of muscle strength decline in older adults.

Introduction

Age-related skeletal muscle strength decline is an independent risk factor for high mortality in older adults [1]. Because handgrip strength (HGS) is easy to measure and can represent global muscle strength in older people [2], it is not only a powerful indicator of physical activity, nutritional status, and even disability [3], but can also be used for the diagnosis of sarcopenia [4, 5] and physical frailty [6]. The inverse association between HGS assessed at a single point in time and subsequent mortality (higher values are associated with lower mortality risk) has been established among older people [7, 8]. Moreover, the clinical value of HGS measured longitudinally to assess the risk of mortality among older people has also been validated [9]. Hence, understanding the biological underpinnings and identifying biomarkers for declining HGS that accompany aging is imperative.

Previous studies have found that inflammation-induced loss of muscle mass, elevated oxidative stress, mitochondrial dysfunction, and accumulation of advanced glycation end-products are involved in the mechanisms of, and serve as biomarkers for, muscle dysfunction-associated conditions, including muscle mass and strength loss. In fact, age-related decline in muscle strength does not parallel the decline in muscle mass [10]. Biomarkers reflecting muscle strength decline in the early stage have only been scarcely investigated, especially in older adults with normal or low HGS. In recent years, studies have shown that all the pathological bases of aging can cause metabolic reactions and that there is a “metabolic clock” that controls aging [11]. Thus, as products of metabolic reactions, metabolites play an important role in physiological and pathological aging [12, 13], making metabolomics an important tool for the non-invasive identification and quantification of biomarkers in biological matrices. However, studies using metabolomics to investigate the metabolites associated with HGS decline in older adults are limited.

The metabolic capacity of the human body is influenced by the gut microbiota and its interactions with host cells [1416]. Gut microbiota dysbiosis, described as a profitless gut microbiota composition and diversity, has been proposed to occur during aging. Therefore, it has been proven that various age-related conditions, such as frailty and sarcopenia, are associated with the gut microbiome [17, 18]. The literature supports the possible presence of a “gut-muscle axis,” whereby gut microbial metabolism influences the functionality of muscle cells by producing mediators that drive the systemic effects of the gut microbiota [19, 20]. However, the potential relationship between the gut microbiome and metabolite profiles in older adults with low muscle strength is not well understood.

In this study, liquid chromatography-mass spectrometry (LC-MS) and 16S ribosomal DNA (rDNA) gene sequencing were performed to analyze the metabolome of serum and stool samples and the gut microbiome composition of stool samples from older adults with normal and low HGS. Correlation analysis was used to identify potential serum and fecal metabolites associated with the decline in HGS in older people. We further applied an integrated analysis of the microbiome and metabolome profiling to explore the correlations between altered gut microbiota and metabolites. The workflow is illustrated in Figure 1. The identified between-group differences for metabolites and their associated intestinal flora could indicate interactions between the host and gut microbiome in older adults with declining muscle strength.

Overview of the study.

Figure 1. Overview of the study.

Materials and Methods

Clinical assessment

Individuals included in the study were recruited from the Division of Geriatric Endocrinology of the First Affiliated Hospital of Nanjing Medical University from 2019 to 2020 [21]. Patients with malignant tumors in the acute phase of the disease, those with severe cognitive impairment, and those taking antibiotics within the past month were excluded. Participants were given a form with questions concerning several lifestyle variables (smoking, alcohol consumption, comorbidities, and medications). Body mass index (BMI) was calculated based on height and weight measurements.

Handgrip strength measurement

HGS of the dominant hand was measured using a portable hydraulic dynamometer (Jamar 5030J1, Jamar Technologies, Horsham, PA, USA). The participants adopted a seated upright position on a height-adjustable chair [22]. The test arm was positioned on a table to support the dynamometer’s weight. The participants were instructed and verbally encouraged to squeeze their handgrip as hard as possible. The best performance of the three trials was used in the analysis. Low HGS diagnostic cutoffs are <28.0 kg for men and <18.0 kg for women, according to the standard of the Asian Working Group for Sarcopenia (AWGS 2019) [4]. All measurements were performed by the same staff. Other methods for muscle mass measurement and physical performance, including the usual 4-m gait speed, 5-time chair stand test, and the Short Physical Performance Battery (SPPB), are specified in the Supplementary Materials and Methods.

Blood and stool sample collection

Blood samples were collected from all participants by nurses on the second day of overnight fasting. Samples were centrifuged at 3,000 rpm for 10 min at 4°C to obtain the serum, and then aliquoted and stored at −80°C until further analyses. Stool samples were collected within three days of admission using a Commode specimen collection system (Thermo Fisher Scientific, Waltham, MA, USA), aliquoted, frozen immediately in liquid nitrogen, and then stored at −80°C before further processing.

Untargeted metabolomics analysis of serum and stool samples

All the samples were thawed on ice before extraction using a solvent with a vortex or grinding step before centrifugation. After centrifugation, the reconstituted solution was added to the supernatant with a 60-s vortexing step before centrifugation for reconstitution. The supernatants were collected for metabolomic profiling using LC-MS analysis. A quality control (QC) sample was prepared by mixing equal volumes of the supernatant of each sample to assess the analytical variability. Untargeted metabolomics LC-MS analysis was performed, and data were collected in both the positive- and negative-ion modes to improve the coverage of metabolites. Refer to the Supplementary Materials and Methods for further details.

Raw LC-MS data were extracted and processed using multivariate statistical analysis (Partial Least Squares Method-Discriminant Analysis, PLS-DA) to establish a relationship model between the metabolite and sample groups. Univariate methods (Wilcoxon test and two-tailed Student’s t-test) were used to detect the metabolites whose levels were significantly altered; then, correction was performed by calculating the false discovery rate (FDR) to ensure that the metabolite peaks were reproducibly detected. Metabolites responsible for the difference in the metabolic profile scan between the groups were obtained based on a variable importance in the projection (VIP) threshold of 1 from the 7-fold cross-validated PLS-DA model. By combining the univariate and multivariate statistical analyses, metabolites whose levels were significantly altered between groups were acquired under the following conditions: p-value < 0.05, q-value < 0.05, fold change < 0.8 or > 1.2, and VIP > 1.

The online HMDB database (http://www.hmdb.ca) was used to annotate metabolites by matching the exact molecular mass data (m/z) of the samples with those from the database [23]. The KEGG database (https://www.genome.jp/kegg/) was used to understand the functional characteristics of the differential metabolites and to determine the main biochemical metabolic pathways and signal transduction associated with the metabolites. A hypergeometric test was used to identify the significantly enriched pathway entries. Metabolic pathways with a p-value < 0.05 were deemed to be significantly enriched by differential metabolites.

16S rDNA microbiome analysis of stool samples

DNA was extracted from all fecal samples and the 16S rDNA was PCR-amplified and sequenced on the MiSeq system (Illumina, San Diego, CA, USA). The sequence reads with a similarity greater than 97% were identified for further analysis. Operational taxonomic unit (OTU) representative sequences were taxonomically classified using the Ribosomal Database Project (RDP) Classifier, and community composition was analyzed for each taxonomic rank: domain, kingdom, phylum, class, order, family, genus, and species. The complete protocol can be found in the Supplementary Materials and Methods.

Spearman’s multiomics correlation analysis

Spearman correlation analysis based on the R package (v3.4.1) was used to analyze the correlation between the metabolites and microbiota. Spearman’s correlation coefficients were also computed for the relationships between the HGS values and individual metabolomic features. Differences were considered significant when the p-value was < 0.05 and |r| was > 0.5. If r < 0, there was a negative correlation; otherwise, there was a positive one.

Data availability

Raw LC-MS data are available in the MetaboLights repository (accession no. MTBLS4367 and MTBLS4372). The sequencing data of the 16S rDNA were deposited in the BioProject database (accession no. PRJNA787524).

Results

Clinical data of study participants

Eight inpatients with low HGS and seven age- and environment-matched inpatients with normal HGS (age, 77–90 years) were recruited. HGS is an important component of sarcopenia diagnosis and a representative indicator of muscle strength. Thus, we also collected information regarding other components of sarcopenia, including appendicular skeletal muscle (ASM), 4-m gait speed, 5-time chair-stand test, and SPPB (Table 1). The clinical characteristics of the study participants are shown in Table 2. There were no statistically significant differences between the two groups with regard to age, BMI, health behaviors, ASM, 4-m gait speed, or SPPB. The total number of comorbidities and medications, especially, diabetes and antidiabetic medications, was higher in the low-HGS group than in the normal group. Similar to the HGS, a 5-time chair stand test was another evaluation method used to determine “Possible sarcopenia”. The time for the test increased in the low-HGS group.

Table 1. Sarcopenia-associated parameters in accordance with handgrip strength in 15 participants.

Normal HGS
GenderAge (years)Muscle strengthMuscle massPhysical performance
HGS (kg)RASM (kg/m2)4-m gait speed (m/s)5-time chair stand test (s)SPPB
M8229.805.800.9917.018
M8129.906.061.2015.4010
M7929.937.071.2611.1812
M8331.106.411.729.3912
M8734.305.560.80-5
F8220.005.351.2911.559
F7720.104.751.15-8
Avg ± std81.57 ± 3.1527.88 ± 5.575.86 ± 0.751.20 ± 0.2812.91 ± 3.179 ± 2
Low HGS (Male: < 28 kg, Female: < 18 kg)
GenderAge (years)Muscle strengthMuscle massPhysical performance
HGS (kg)RASM (kg/m2)4-m gait speed (m/s)5-time chair stand test (s)SPPB
M8127.104.661.0418.309
M9018.805.45--0
M8415.005.95--0
M8026.675.380.7218.976
M9023.375.191.1615.3410
M8524.405.980.9824.509
F*8015.905.981.2612.0711
F*8710.375.820.5723.174
Avg ± std84.63 ± 4.1420.20 ± 6.115.55 ± 0.470.95 ± 0.2618.73 ± 4.676 ± 4
P-value0.13640.02520.35590.13640.04270.1326
Abbreviations: M: male; F: female; HGS: handgrip strength; RASM: relative appendicular skeletal muscle mass; SPPB: Short Physical Performance Battery; -: Participants cannot complete the test. *: Possible sarcopenia, : Sarcopenia, : Severe sarcopenia (according to the diagnosis criteria of Asian Working Group for Sarcopenia, 2019).

Table 2. Clinical characteristics of the study participants.

Handgrip strengthP-value
Normal (n = 7)Low (n = 8)
Age, years81.57 ± 3.1584.63 ± 4.140.1364
Sex (male/female)5/26/2N/A
BMI, kg/m224.00 ± 1.6123.26 ± 2.430.5062
Health behavior (n (%))
 Smoking2 (29%)3 (38%)0.7144
 Alcohol drinking1 (14%)2 (25%)0.6048
Total number of comorbidities (n (%))0.96 ± 0.112.75 ± 0.370.0024
 Hypertension2 (28.57%)5 (62.50%)0.1888
 Diabetes mellitus0 (0%)5 (62.50%)0.0104
 Stroke0 (0%)3 (37.50%)0.0701
 Coronary artery disease0 (0%)1 (12.50%)0.3329
Total number of medications (n (%))1.43 ± 0.535.25 ± 0.840.0026
 Psycho-drugs0 (0%)0 (0%)N/A
 Antidiabetic agent0 (0%)5 (62.50%)0.0104
 Statin3 (42.86%)1 (12.50%)0.1847
 β-blocker0 (0%)1 (12.50%)0.3329
 Antiplatelet2 (28.57%)1 (12.50%)0.4376
 ACEi and/or ARB1 (14.29%)3 (37.50%)0.3104
 β2 sympathomimetic inhaler0 (0%)0 (0%)N/A

Serum metabolomics analysis

Metabolomics analysis yielded 386 metabolites in the negative-ion mode and 861 metabolites in the positive-ion mode. Supervised PLS-DA revealed that the metabolomic profiles differed between the normal and low HGS groups (Figure 2A). In total, the abundance of 67 annotated metabolites was found to be significantly different between the two groups, with 25 metabolites having higher concentrations and 42 metabolites having lower concentrations in the serum of older adults with low HGS (Supplementary Table 1). Class information was available for only 47 of these 67 metabolites, according to the Human Metabolome Database (HMDB version 5.0). Analysis of these metabolites showed that compared with the normal group, the shifts in the low-HGS group mainly included benzene and their (substituted) derivatives, cinnamic acids and their derivatives, and fatty acids (Figure 2B). In particular, the levels of metabolites belonging to the class of cinnamic acids and their derivatives were lower in the serum of older adults with low HGS, while the levels of metabolites belonging to the class of fatty acids were higher in the normal group (Supplementary Table 1). KEGG metabolic pathway analysis revealed a significant enrichment of 10 pathways, including linoleic acid metabolism, biosynthesis of unsaturated fatty acids, and ABC transporters (Figure 2C, Supplementary Table 2). The results revealed alterations in serum metabolic characteristics in subjects with decreased HGS.

Serum metabolomics analysis. (A) PLS-DA analysis of the grouped discrimination by the first two principal components (PCs) in negative and positive ion modes. (B) Pie graph of the class composition according to the number of altered metabolites. (C) Bubble chart of the pathway enrichment analysis of differential metabolites analyzed in the negative- and positive-ion modes. RichFactor was the number of differential metabolites divided by the identified metabolites annotated to the pathway.

Figure 2. Serum metabolomics analysis. (A) PLS-DA analysis of the grouped discrimination by the first two principal components (PCs) in negative and positive ion modes. (B) Pie graph of the class composition according to the number of altered metabolites. (C) Bubble chart of the pathway enrichment analysis of differential metabolites analyzed in the negative- and positive-ion modes. RichFactor was the number of differential metabolites divided by the identified metabolites annotated to the pathway.

To evaluate the relationship between the HGS value and the 67 annotated metabolites, Spearman’s correlation analysis was performed. Six metabolites were positively correlated with the HGS values: cinnamoylglycine, 4-methoxycinnamic acid, (e)-3,4,5-trimethoxycinnamic acid, dimethyl ((e)-(1-methoxy-2-oxo-1,2-dihydro-3h-indol-3-ylidene)methyl)carbonodithioimidate, dalpanin, and Af4878000 (Figure 3). Among these metabolites, the appearance of cinnamoylglycine, 4-methoxycinnamic acid, and (e)-3,4,5-trimethoxycinnamic acid suggested that cinnamic acid metabolism may play a role in the decline of HGS in older adults.

Correlation analysis of the HGS values and differential metabolites in serum samples.

Figure 3. Correlation analysis of the HGS values and differential metabolites in serum samples.

Fecal metabolomics analysis

In total, 1268 metabolites were identified in the negative-ion mode and 2255 metabolites were identified in the positive-ion mode. PLS-DA showed a clear separation trend between the normal and low HGS groups (Figure 4A). All fecal metabolites with annotations that were higher (n = 78) or lower (n = 99) in the low-HGS group are shown in Supplementary Table 3. Analysis of 96 metabolites for which class information was available revealed that the levels of fecal metabolites were altered in the low-HGS group; these metabolites mainly included benzene and its (substituted) derivatives, carboxylic acids and their derivatives, flavonoids, and phenol ethers (Figure 4B). Moreover, the relative concentration of metabolites belonging to the class of cinnamic acids and their derivatives also decreased in the feces of older subjects with low HGS (Supplementary Table 3). KEGG metabolic pathway mapping showed that the metabolites were significantly enriched in seven pathways, including cholesterol metabolism, primary bile acid biosynthesis, and bile secretion (Figure 4C, Supplementary Table 4).

Fecal metabolomics analysis. (A) PLS-DA analysis of the grouped discrimination by the first two principal components (PCs) in the negative- and positive-ion modes. (B) Pie graph of the class composition according to the number of altered metabolites. (C) Bubble chart of pathway enrichment analysis of differential metabolites analyzed in the negative- and positive-ion modes. RichFactor was the number of differential metabolites divided by the identified metabolites annotated to the pathway.

Figure 4. Fecal metabolomics analysis. (A) PLS-DA analysis of the grouped discrimination by the first two principal components (PCs) in the negative- and positive-ion modes. (B) Pie graph of the class composition according to the number of altered metabolites. (C) Bubble chart of pathway enrichment analysis of differential metabolites analyzed in the negative- and positive-ion modes. RichFactor was the number of differential metabolites divided by the identified metabolites annotated to the pathway.

Similarly, to determine the differences that could be associated with HGS, we examined the correlations between the relative concentrations of the annotated metabolites and HGS values. Thirteen metabolites were found to be associated with changes in HGS values, of which four were positively correlated with HGS and nine were negatively correlated with HGS (Supplementary Figure 1). Most of the metabolites were drugs and their metabolic derivatives, such as sitagliptin, pioglitazone, emedastine, and venlafaxine. Gamma-glu-gln and histidylphenylalanine, which were included in the class of carboxylic acids and their derivatives, showed a positive correlation with HGS. Indole-2-carboxylic acid was positively correlated with HGS. Importantly, one metabolite from the class of cinnamic acids and their derivatives, (2e)-5-hydroxyferulic acid, was also positively correlated with HGS. These results suggest that the fecal metabolite profiles between the normal and low-HGS groups differed.

Microbiome analysis

Previous studies demonstrated a surprisingly notable effect of the gut microbiome on mammalian blood metabolites [24] and approximately 58% of metabolites detected in human body originate from microbiota [25]. To identify the changes in the gut microbiota of older adults with low HGS, we performed 16S rDNA amplicon sequencing of fecal samples from 15 participants. In total, 543 OTUs were detected, and 349 OTUs overlapped between the normal and low-HGS groups (Figure 5A). The number of unique OTUs in the low-HGS group was slightly higher than that in the normal group, indicating a more abundant microbiota in individuals with low HGS. According to the species associated OTUs and sequence number, rarefaction curves of the samples were calculated, and a flat trend indicated that the sampling size was reasonable (Figure 5B). Gut microbial α-diversity was lower in the low-HGS group than in the normal group, as calculated using the Simpson indices (Figure 5C). For β-diversity, UPGMA cluster analysis based on weighted UniFrac analysis was performed, and the phylogenetic distance between the samples was calculated (Figure 5D). We observed a higher β-diversity in the gut microbiota of older adults with normal HGS, indicating a more heterogeneous community structure among participants with normal HGS than among those with low HGS (Figure 5E).

Microbiome analysis. (A) The Core-Pan graph of OTU distribution between the two groups. (B) The rarefaction curve of random sequences per sample and their corresponding number of observed species. (C) Species diversity differences estimated on the basis of the observed Sobs, Chao, Ace, Shannon, Simpson, and Coverage indices. (D) UPGMA cluster analysis of 15 samples at the genus level. A1–A7 represent the normal handgrip strength group; B8–B15 represents the low handgrip strength group. (E) β-diversity box-plot constructed on the basis of Weighted UniFrac analysis between the two groups. (F) The percentages of gut microbiota diversity at the phylum level. (G) The percentages of gut microbiota diversity at the genus level. (H) Linear discriminant analysis (LDA) integrated with effect size (LEfSe). Left: the phylogenetic distribution of microbiota in the Cladogram. Right: the differences in the abundance of microbiota.

Figure 5. Microbiome analysis. (A) The Core-Pan graph of OTU distribution between the two groups. (B) The rarefaction curve of random sequences per sample and their corresponding number of observed species. (C) Species diversity differences estimated on the basis of the observed Sobs, Chao, Ace, Shannon, Simpson, and Coverage indices. (D) UPGMA cluster analysis of 15 samples at the genus level. A1–A7 represent the normal handgrip strength group; B8–B15 represents the low handgrip strength group. (E) β-diversity box-plot constructed on the basis of Weighted UniFrac analysis between the two groups. (F) The percentages of gut microbiota diversity at the phylum level. (G) The percentages of gut microbiota diversity at the genus level. (H) Linear discriminant analysis (LDA) integrated with effect size (LEfSe). Left: the phylogenetic distribution of microbiota in the Cladogram. Right: the differences in the abundance of microbiota.

Next, the gut microbial composition abundances in the stool samples were analyzed at the phylum and genus levels (Figure 5F, 5G). At the phylum level, Bacteroidetes and Firmicutes were the predominant phyla. There were no significant differences in microbiome composition abundance at the phylum level between the normal and low-HGS groups. At the genus level, the relative abundance of the genera Parabacteroides and Intestinibacter increased in the low-HGS group (Table 3). Furthermore, we used linear discriminant analysis effect size (LEfSe) to generate a cladogram to identify specific bacteria associated with decreased HGS. The results showed that members from Bacteroidaceae, Bacteroides, Selenomonadales, Negativicutes, Porphyromonadaceae, Parabacteroides, Intestinibacter, Enterococcus, and Enterococcaceae were the most abundant microbiota in the low-HGS group (all | Linear discriminant analysis (LDA) scores (log10) | > 3.6). However, participants with normal HGS were mainly characterized by higher abundances of Sporobacter and Terrisporobacter (all | LDA scores (log10) | > 3.6) (Figure 5H). The data revealed pronounced differences in microbiota between the two groups.

Table 3. Relative abundance of fecal microbiota at the genus level.

Bacteria speciesNormal HGSLow HGSP-valueFDR
Relative abundance (%)
Parabacteroides0.4000736.8898530.0065360.466388
Intestinibacter0.0044980.0278540.0089690.466388

Correlations between gut microbiota and serum/fecal metabolites

To explore the potential dependencies between microbiome composition, host metabolism, and the metabolome, we examined the correlations between the two datasets. Spearman’s correlation coefficient was computed between the relative abundance of HGS-associated species at the genus level and the different serum and fecal metabolites. Enrichment of members from the genus Parabacteroides in stool samples from older adults with low HGS was negatively correlated with the six serum metabolites whose levels increased with the HGS value (Figure 6A). The correlation coefficient between cinnamoylglycine and Parabacteroides was the most significant (|r| = 0.82). Moreover, Intestinibacter, whose abundance also increased in older adults with low HGS, correlated negatively with cinnamoylglycine. For fecal metabolites associated with low HGS, the abundance of Parabacteroides was also negatively correlated with (2e)-5-hydroxyferulic acid, a cinnamic acid derivative (Figure 6B). The abundances of the genera Parabacteroides and Intestinibacter were negatively correlated with indole-2-carboxylic acid. Correlations between serum and fecal metabolites altered by HGS were also analyzed. Six serum metabolites showed positive correlations with HGS and were positively correlated with fecal indole-2-carboxylic acid and gamma-glu-gln levels. Serum cinnamoylglycine was positively correlated with fecal (2e)-5-hydroxyferulic acid, which also belongs to the class cinnamic acids and their derivatives (Figure 6C). These results indicated that alterations in the levels of serum metabolites and their related fecal metabolites in the low-HGS group may be associated with the gut microbiota.

Correlations between gut microbiota and serum/fecal metabolites. (A) Spearman’s correlation coefficient between the genera Parabacteroides and Intestinibacter and serum metabolites. (B) Spearman’s correlation coefficient between the genera Parabacteroides and Intestinibacter and fecal metabolites. (C) Spearman’s correlation coefficient between the serum and fecal metabolites. Differences were considered significant when the p-value was  0.5. If r

Figure 6. Correlations between gut microbiota and serum/fecal metabolites. (A) Spearman’s correlation coefficient between the genera Parabacteroides and Intestinibacter and serum metabolites. (B) Spearman’s correlation coefficient between the genera Parabacteroides and Intestinibacter and fecal metabolites. (C) Spearman’s correlation coefficient between the serum and fecal metabolites. Differences were considered significant when the p-value was < 0.05 and |r| was > 0.5. If r < 0, there was a negative correlation.

Discussion

Declining muscle strength with aging is almost inevitable; however, the metabolic mechanisms underlying this phenomenon remain poorly understood. Using an LC-MS-based metabolomics approach, we investigated the serum and fecal metabolites whose levels were altered in older adults with normal and low HGS. In addition to the metabolites belonging to the class of benzene and its (substituted) derivatives, the shifts in the low-HGS group mainly included cinnamic acids and derivatives. 4-methoxycinnamic acid, (e)-3,4,5-trimethoxycinnamic acid, sinapinic acid, and (2e)-3-(3,4-dimethoxyphenyl) acrylic acid were the four metabolites belonging to the class of cinnamic acids and their derivatives, and their relative concentrations decreased in the serum of older adults with low HGS. In particular, 4-methoxycinnamic acid and (e)-3,4,5-trimethoxycinnamic acid levels were positively correlated with HGS. Fecal metabolites of this class also decreased in the low-HGS group. Cinnamic acids and their derivatives, both synthetic and derived from natural sources, have been reported to exert a wide range of biological activities, including anticancer, hepatoprotective, neuroprotective, cardioprotective, and antidiabetic effects, and potent functions in muscle cell proliferation, differentiation, and development [2630]. Sinapinic acid is one of the most common hydroxycinnamic acids; it is widespread in the plant kingdom [31]. It has also been reported to be a major active component of traditional Chinese remedies [32]. The antioxidant and anti-inflammatory activities of sinapinic acid have been reported to be highly significant. Sinapinic acid scavenges native peroxynitrite (ONOO) and is more efficient than its alkyl esters [33, 34]. Lee et al. [35] showed that sinapinic acid modulates inflammation by suppressing NOD-like receptor pyrin domain-containing 3 (NLRP3) inflammasome activation. In addition, supplementation with sinapinic acid affected the intestinal microbiome by improving the proportion of the butyrate acid producers Blautia and Dorea and inhibiting the growth of bacterial species associated with diseases and inflammation, such as Bacteroides [36]. Although further studies on other cinnamic acids and their derivatives are needed, the above findings, including those of our study, fully suggest the potential role of this class of metabolites in aging and related diseases.

Cinnamoylglycine, another metabolite involved in the metabolism of cinnamic acid, is a glycine conjugate of cinnamic acid; it belongs to the class of carboxylic acids and their derivatives. In the present study, the relative concentration of cinnamoylglycine decreased with the decrease in the HGS. We also found that the levels of hippurate, a glycine conjugate of benzoic acid, decreased in the serum of older adults with declining HGS. These two metabolites are derived from catabolism of dietary polyphenols found in plant-based foods, which is performed by the intestinal microflora [24, 37]. Hippuric acid in plasma has been recognized as a plausible hallmark of frailty and geriatric syndromes [38]. Lustgarten et al. found that cinnamoylglycine, which is related to gut bacterial metabolism, was negatively associated with muscle quality [39] and physical function in older adults [40]. Cinnamoylglycine also predicted a higher gut microbiome diversity and was linked to a lower incidence of type 2 diabetes [41]. Although the current findings regarding the role of cinnamoylglycine in muscle strength and body function are inconsistent, cinnamoylglycine has the potential to serve as a biomarker of age-related changes in muscle strength; this warrants further studies.

A few integrated studies regarding the relationship between serum, fecal metabolomics, and gut microbiome and muscle dysfunction have been published [42, 43]. The microbial species and their associated metabolites involved in HGS-associated alterations were identified in the present study. We found that the abundance of two genera, Parabacteroides and Intestinibacter, which belong to the phyla Bacteroidetes and Firmicutes, respectively, increased in the low-HGS group. The two gut commensal bacteria have been reported to produce short-chain fatty acids [44, 45]. Parabacteroides produces acetate to reduce neutrophil infiltration [45] and can be affected by prebiotic supplementation in frail older subjects [46]. The metabolic benefits of Parabacteroides distasonis in decreasing weight gain, hyperglycemia, and hepatic steatosis have been previously studied [47]. A higher abundance of Intestinibacter, which can produce butyrate, is associated with the reduced incidence of type 2 diabetes [44] and its abundance decreases immediately after metformin treatment [48, 49]. A recent study on coronavirus disease 2019 (COVID-19) found that Intestinibacter bartlettii was positively correlated with anorexia and fatigue in survivors of COVID-19 after discharge from the hospital [50]. We found that the members of Parabacteroides, which were enriched in stool samples of low-HGS older adults, were negatively related to six serum metabolites that were associated with HGS, including cinnamoylglycine, a cinnamic acid derivative. The abundance of Parabacteroides was also negatively correlated with fecal (2e)-5-hydroxyferulic acid, which is a cinnamic acid derivative. The abundances of members from Parabacteroides and Intestinibacter were negatively correlated with the fecal metabolite indole-2-carboxylic acid. The discovery of these microbial metabolites and their correlated microbiota in older adults with low HGS indicates that microbial metabolites, along with other traditional markers, can be used as potential biomarkers for muscle strength decline.

This study has limitations including small sample size. Although environmental-matched controls were recruited and participated, an important limitation was the lack of dietary control among participants. Dietary patterns consistently correlate with groups of bacteria with shared functional roles in both health and disease [51]. The influence of dietary protein on the gut microbiome and its impact on sarcopenia also has been explored [52]. We will expand the sample size and consider various influencing factors, including diet, in future studies.

In summary, this work demonstrates that the levels of serum and fecal metabolites, especially, gut microbiota-related metabolites, are notably altered in older people with low HGS. The findings of the current study suggest potential metabolic underpinnings for HGS and provide a basis for the further identification of biomarkers of muscle strength decline in older adults.

Author Contributions

All authors participated in the design, interpretation of the studies and analysis of the data and review of the manuscript; JY, GD and XW designed the experiments, YG, QW and JY wrote and revised the manuscript. YG, QW, YL and FX collected the samples and analyzed the data, and XC and YM performed the geriatric assessments.

Conflicts of Interest

The authors declare no conflicts of interest related to this study.

Ethical Statement and Consent

The Ethics Committee of the First Affiliated Hospital of Nanjing Medical University approved this study (approval number 2019-NT-48), which was performed in accordance with the guidelines for the protection of human subjects. Informed consent was obtained from all the participants included in the study.

Funding

This work was supported by the grants from the National Natural Science Foundation of China (81870613 to G.D.) and the Cadre Health Research Project of Jiangsu Province (BJ19013 to J.Y.).

References

  • 1. Arango-Lopera VE, Arroyo P, Gutiérrez-Robledo LM, Pérez-Zepeda MU, Cesari M. Mortality as an adverse outcome of sarcopenia. J Nutr Health Aging. 2013; 17:259–62. https://doi.org/10.1007/s12603-012-0434-0 [PubMed]
  • 2. Porto JM, Nakaishi APM, Cangussu-Oliveira LM, Freire Júnior RC, Spilla SB, Abreu DCC. Relationship between grip strength and global muscle strength in community-dwelling older people. Arch Gerontol Geriatr. 2019; 82:273–8. https://doi.org/10.1016/j.archger.2019.03.005 [PubMed]
  • 3. Syddall H, Cooper C, Martin F, Briggs R, Aihie Sayer A. Is grip strength a useful single marker of frailty? Age Ageing. 2003; 32:650–6. https://doi.org/10.1093/ageing/afg111 [PubMed]
  • 4. Chen LK, Woo J, Assantachai P, Auyeung TW, Chou MY, Iijima K, Jang HC, Kang L, Kim M, Kim S, Kojima T, Kuzuya M, Lee JSW, et al. Asian Working Group for Sarcopenia: 2019 Consensus Update on Sarcopenia Diagnosis and Treatment. J Am Med Dir Assoc. 2020; 21:300–7.e2. https://doi.org/10.1016/j.jamda.2019.12.012 [PubMed]
  • 5. Cruz-Jentoft AJ, Bahat G, Bauer J, Boirie Y, Bruyère O, Cederholm T, Cooper C, Landi F, Rolland Y, Sayer AA, Schneider SM, Sieber CC, Topinkova E, et al, and Writing Group for the European Working Group on Sarcopenia in Older People 2 (EWGSOP2), and the Extended Group for EWGSOP2. Sarcopenia: revised European consensus on definition and diagnosis. Age Ageing. 2019; 48:16–31. https://doi.org/10.1093/ageing/afy169 [PubMed]
  • 6. Dudzińska-Griszek J, Szuster K, Szewieczek J. Grip strength as a frailty diagnostic component in geriatric inpatients. Clin Interv Aging. 2017; 12:1151–7. https://doi.org/10.2147/CIA.S140192 [PubMed]
  • 7. Gale CR, Martyn CN, Cooper C, Sayer AA. Grip strength, body composition, and mortality. Int J Epidemiol. 2007; 36:228–35. https://doi.org/10.1093/ije/dyl224 [PubMed]
  • 8. Rantanen T, Harris T, Leveille SG, Visser M, Foley D, Masaki K, Guralnik JM. Muscle strength and body mass index as long-term predictors of mortality in initially healthy men. J Gerontol A Biol Sci Med Sci. 2000; 55:M168–73. https://doi.org/10.1093/gerona/55.3.m168 [PubMed]
  • 9. Malhotra R, Tareque MI, Tan NC, Ma S. Association of baseline hand grip strength and annual change in hand grip strength with mortality among older people. Arch Gerontol Geriatr. 2020; 86:103961. https://doi.org/10.1016/j.archger.2019.103961 [PubMed]
  • 10. Delmonico MJ, Harris TB, Visser M, Park SW, Conroy MB, Velasquez-Mieyer P, Boudreau R, Manini TM, Nevitt M, Newman AB, Goodpaster BH, and Health, Aging, and Body. Longitudinal study of muscle strength, quality, and adipose tissue infiltration. Am J Clin Nutr. 2009; 90:1579–85. https://doi.org/10.3945/ajcn.2009.28047 [PubMed]
  • 11. López-Otín C, Galluzzi L, Freije JMP, Madeo F, Kroemer G. Metabolic Control of Longevity. Cell. 2016; 166:802–21. https://doi.org/10.1016/j.cell.2016.07.031 [PubMed]
  • 12. Hertel J, Friedrich N, Wittfeld K, Pietzner M, Budde K, Van der Auwera S, Lohmann T, Teumer A, Völzke H, Nauck M, Grabe HJ. Measuring Biological Age via Metabonomics: The Metabolic Age Score. J Proteome Res. 2016; 15:400–10. https://doi.org/10.1021/acs.jproteome.5b00561 [PubMed]
  • 13. Ma S, Yim SH, Lee SG, Kim EB, Lee SR, Chang KT, Buffenstein R, Lewis KN, Park TJ, Miller RA, Clish CB, Gladyshev VN. Organization of the Mammalian Metabolome according to Organ Function, Lineage Specialization, and Longevity. Cell Metab. 2015; 22:332–43. https://doi.org/10.1016/j.cmet.2015.07.005 [PubMed]
  • 14. Krautkramer KA, Fan J, Bäckhed F. Gut microbial metabolites as multi-kingdom intermediates. Nat Rev Microbiol. 2021; 19:77–94. https://doi.org/10.1038/s41579-020-0438-4 [PubMed]
  • 15. Nicholson JK, Holmes E, Kinross J, Burcelin R, Gibson G, Jia W, Pettersson S. Host-gut microbiota metabolic interactions. Science. 2012; 336:1262–7. https://doi.org/10.1126/science.1223813 [PubMed]
  • 16. Holmes E, Li JV, Marchesi JR, Nicholson JK. Gut microbiota composition and activity in relation to host metabolic phenotype and disease risk. Cell Metab. 2012; 16:559–64. https://doi.org/10.1016/j.cmet.2012.10.007 [PubMed]
  • 17. Ticinesi A, Nouvenne A, Cerundolo N, Catania P, Prati B, Tana C, Meschi T. Gut Microbiota, Muscle Mass and Function in Aging: A Focus on Physical Frailty and Sarcopenia. Nutrients. 2019; 11:1633. https://doi.org/10.3390/nu11071633 [PubMed]
  • 18. Casati M, Ferri E, Azzolino D, Cesari M, Arosio B. Gut microbiota and physical frailty through the mediation of sarcopenia. Exp Gerontol. 2019; 124:110639. https://doi.org/10.1016/j.exger.2019.110639 [PubMed]
  • 19. Prokopidis K, Chambers E, Ni Lochlainn M, Witard OC. Mechanisms Linking the Gut-Muscle Axis With Muscle Protein Metabolism and Anabolic Resistance: Implications for Older Adults at Risk of Sarcopenia. Front Physiol. 2021; 12:770455. https://doi.org/10.3389/fphys.2021.770455 [PubMed]
  • 20. Ticinesi A, Lauretani F, Milani C, Nouvenne A, Tana C, Del Rio D, Maggio M, Ventura M, Meschi T. Aging Gut Microbiota at the Cross-Road between Nutrition, Physical Frailty, and Sarcopenia: Is There a Gut-Muscle Axis? Nutrients. 2017; 9:1303. https://doi.org/10.3390/nu9121303 [PubMed]
  • 21. Guo Y, Zhu G, Wang F, Zhang H, Chen X, Mao Y, Lv Y, Xia F, Jin Y, Ding G, Yu J. Distinct Serum and Fecal Metabolite Profiles Linking With Gut Microbiome in Older Adults With Frailty. Front Med (Lausanne). 2022; 9:827174. https://doi.org/10.3389/fmed.2022.827174 [PubMed]
  • 22. Roberts HC, Denison HJ, Martin HJ, Patel HP, Syddall H, Cooper C, Sayer AA. A review of the measurement of grip strength in clinical and epidemiological studies: towards a standardised approach. Age Ageing. 2011; 40:423–9. https://doi.org/10.1093/ageing/afr051 [PubMed]
  • 23. Feng Q, Liu Z, Zhong S, Li R, Xia H, Jie Z, Wen B, Chen X, Yan W, Fan Y, Guo Z, Meng N, Chen J, et al. Integrated metabolomics and metagenomics analysis of plasma and urine identified microbial metabolites associated with coronary heart disease. Sci Rep. 2016; 6:22525. https://doi.org/10.1038/srep22525 [PubMed]
  • 24. Wikoff WR, Anfora AT, Liu J, Schultz PG, Lesley SA, Peters EC, Siuzdak G. Metabolomics analysis reveals large effects of gut microflora on mammalian blood metabolites. Proc Natl Acad Sci U S A. 2009; 106:3698–703. https://doi.org/10.1073/pnas.0812874106 [PubMed]
  • 25. Dekkers KF, Sayols-Baixeras S, Baldanzi G, Nowak C, Hammar U, Nguyen D, Varotsis G, Brunkwall L, Nielsen N, Eklund AC, Bak Holm J, Nielsen HB, Ottosson F, et al. An online atlas of human plasma metabolite signatures of gut microbiome composition. Nat Commun. 2022; 13:5370. https://doi.org/10.1038/s41467-022-33050-0 [PubMed]
  • 26. Adisakwattana S. Cinnamic Acid and Its Derivatives: Mechanisms for Prevention and Management of Diabetes and Its Complications. Nutrients. 2017; 9:163. https://doi.org/10.3390/nu9020163 [PubMed]
  • 27. Kuppusamy P, Soundharrajan I, Kim DH, Hwang I, Choi KC. 4-hydroxy-3-methoxy cinnamic acid accelerate myoblasts differentiation on C2C12 mouse skeletal muscle cells via AKT and ERK 1/2 activation. Phytomedicine. 2019; 60:152873. https://doi.org/10.1016/j.phymed.2019.152873 [PubMed]
  • 28. Pandi A, Kalappan VM. Pharmacological and therapeutic applications of Sinapic acid-an updated review. Mol Biol Rep. 2021; 48:3733–45. https://doi.org/10.1007/s11033-021-06367-0 [PubMed]
  • 29. Song F, Li H, Sun J, Wang S. Protective effects of cinnamic acid and cinnamic aldehyde on isoproterenol-induced acute myocardial ischemia in rats. J Ethnopharmacol. 2013; 150:125–30. https://doi.org/10.1016/j.jep.2013.08.019 [PubMed]
  • 30. Zhang WX, Wang H, Cui HR, Guo WB, Zhou F, Cai DS, Xu B, Jia XH, Huang XM, Yang YQ, Chen HS, Qi JC, Wang PL, Lei HM. Design, synthesis and biological evaluation of cinnamic acid derivatives with synergetic neuroprotection and angiogenesis effect. Eur J Med Chem. 2019; 183:111695. https://doi.org/10.1016/j.ejmech.2019.111695 [PubMed]
  • 31. Herrmann K. Occurrence and content of hydroxycinnamic and hydroxybenzoic acid compounds in foods. Crit Rev Food Sci Nutr. 1989; 28:315–47. https://doi.org/10.1080/10408398909527504 [PubMed]
  • 32. Chen C. Sinapic Acid and Its Derivatives as Medicine in Oxidative Stress-Induced Diseases and Aging. Oxid Med Cell Longev. 2016; 2016:3571614. https://doi.org/10.1155/2016/3571614 [PubMed]
  • 33. Gaspar A, Martins M, Silva P, Garrido EM, Garrido J, Firuzi O, Miri R, Saso L, Borges F. Dietary phenolic acids and derivatives. Evaluation of the antioxidant activity of sinapic acid and its alkyl esters. J Agric Food Chem. 2010; 58:11273–80. https://doi.org/10.1021/jf103075r [PubMed]
  • 34. Zou Y, Kim AR, Kim JE, Choi JS, Chung HY. Peroxynitrite scavenging activity of sinapic acid (3,5-dimethoxy-4-hydroxycinnamic acid) isolated from Brassica juncea. J Agric Food Chem. 2002; 50:5884–90. https://doi.org/10.1021/jf020496z [PubMed]
  • 35. Lee EH, Shin JH, Kim SS, Seo SR. Sinapic Acid Controls Inflammation by Suppressing NLRP3 Inflammasome Activation. Cells. 2021; 10:2327. https://doi.org/10.3390/cells10092327 [PubMed]
  • 36. Yang C, Deng Q, Xu J, Wang X, Hu C, Tang H, Huang F. Sinapic acid and resveratrol alleviate oxidative stress with modulation of gut microbiota in high-fat diet-fed rats. Food Res Int. 2019; 116:1202–11. https://doi.org/10.1016/j.foodres.2018.10.003 [PubMed]
  • 37. Lees HJ, Swann JR, Wilson ID, Nicholson JK, Holmes E. Hippurate: the natural history of a mammalian-microbial cometabolite. J Proteome Res. 2013; 12:1527–46. https://doi.org/10.1021/pr300900b [PubMed]
  • 38. De Simone G, Balducci C, Forloni G, Pastorelli R, Brunelli L. Hippuric acid: Could became a barometer for frailty and geriatric syndromes? Ageing Res Rev. 2021; 72:101466. https://doi.org/10.1016/j.arr.2021.101466 [PubMed]
  • 39. Lustgarten MS, Price LL, Fielding RA. Analytes and Metabolites Associated with Muscle Quality in Young, Healthy Adults. Med Sci Sports Exerc. 2015; 47:1659–64. https://doi.org/10.1249/MSS.0000000000000578 [PubMed]
  • 40. Lustgarten MS, Price LL, Chalé A, Fielding RA. Metabolites related to gut bacterial metabolism, peroxisome proliferator-activated receptor-alpha activation, and insulin sensitivity are associated with physical function in functionally-limited older adults. Aging Cell. 2014; 13:918–25. https://doi.org/10.1111/acel.12251 [PubMed]
  • 41. Menni C, Zhu J, Le Roy CI, Mompeo O, Young K, Rebholz CM, Selvin E, North KE, Mohney RP, Bell JT, Boerwinkle E, Spector TD, Mangino M, et al. Serum metabolites reflecting gut microbiome alpha diversity predict type 2 diabetes. Gut Microbes. 2020; 11:1632–42. https://doi.org/10.1080/19490976.2020.1778261 [PubMed]
  • 42. Ma L, Keng J, Cheng M, Pan H, Feng B, Hu Y, Feng T, Yang F. Gut Microbiome and Serum Metabolome Alterations Associated with Isolated Dystonia. mSphere. 2021; 6:e0028321. https://doi.org/10.1128/mSphere.00283-21 [PubMed]
  • 43. Clos-Garcia M, Andrés-Marin N, Fernández-Eulate G, Abecia L, Lavín JL, van Liempd S, Cabrera D, Royo F, Valero A, Errazquin N, Vega MCG, Govillard L, Tackett MR, et al. Gut microbiome and serum metabolome analyses identify molecular biomarkers and altered glutamate metabolism in fibromyalgia. EBioMedicine. 2019; 46:499–511. https://doi.org/10.1016/j.ebiom.2019.07.031 [PubMed]
  • 44. Chen Z, Radjabzadeh D, Chen L, Kurilshikov A, Kavousi M, Ahmadizar F, Ikram MA, Uitterlinden AG, Zhernakova A, Fu J, Kraaij R, Voortman T. Association of Insulin Resistance and Type 2 Diabetes With Gut Microbial Diversity: A Microbiome-Wide Analysis From Population Studies. JAMA Netw Open. 2021; 4:e2118811. https://doi.org/10.1001/jamanetworkopen.2021.18811 [PubMed]
  • 45. Lei Y, Tang L, Liu S, Hu S, Wu L, Liu Y, Yang M, Huang S, Tang X, Tang T, Zhao X, Vlodavsky I, Zeng S, et al. Parabacteroides produces acetate to alleviate heparanase-exacerbated acute pancreatitis through reducing neutrophil infiltration. Microbiome. 2021; 9:115. https://doi.org/10.1186/s40168-021-01065-2 [PubMed]
  • 46. Tran TTT, Cousin FJ, Lynch DB, Menon R, Brulc J, Brown JR, O'Herlihy E, Butto LF, Power K, Jeffery IB, O'Connor EM, O'Toole PW. Prebiotic supplementation in frail older people affects specific gut microbiota taxa but not global diversity. Microbiome. 2019; 7:39. https://doi.org/10.1186/s40168-019-0654-1 [PubMed]
  • 47. Wang K, Liao M, Zhou N, Bao L, Ma K, Zheng Z, Wang Y, Liu C, Wang W, Wang J, Liu SJ, Liu H. Parabacteroides distasonis Alleviates Obesity and Metabolic Dysfunctions via Production of Succinate and Secondary Bile Acids. Cell Rep. 2019; 26:222–35.e5. https://doi.org/10.1016/j.celrep.2018.12.028 [PubMed]
  • 48. Bryrup T, Thomsen CW, Kern T, Allin KH, Brandslund I, Jørgensen NR, Vestergaard H, Hansen T, Hansen TH, Pedersen O, Nielsen T. Metformin-induced changes of the gut microbiota in healthy young men: results of a non-blinded, one-armed intervention study. Diabetologia. 2019; 62:1024–35. https://doi.org/10.1007/s00125-019-4848-7 [PubMed]
  • 49. Rodriguez J, Hiel S, Delzenne NM. Metformin: old friend, new ways of action-implication of the gut microbiome? Curr Opin Clin Nutr Metab Care. 2018; 21:294–301. https://doi.org/10.1097/MCO.0000000000000468 [PubMed]
  • 50. Zhou Y, Zhang J, Zhang D, Ma WL, Wang X. Linking the gut microbiota to persistent symptoms in survivors of COVID-19 after discharge. J Microbiol. 2021; 59:941–8. https://doi.org/10.1007/s12275-021-1206-5 [PubMed]
  • 51. Bolte LA, Vich Vila A, Imhann F, Collij V, Gacesa R, Peters V, Wijmenga C, Kurilshikov A, Campmans-Kuijpers MJE, Fu J, Dijkstra G, Zhernakova A, Weersma RK. Long-term dietary patterns are associated with pro-inflammatory and anti-inflammatory features of the gut microbiome. Gut. 2021; 70:1287–98. https://doi.org/10.1136/gutjnl-2020-322670 [PubMed]
  • 52. Prokopidis K, Cervo MM, Gandham A, Scott D. Impact of Protein Intake in Older Adults with Sarcopenia and Obesity: A Gut Microbiota Perspective. Nutrients. 2020; 12:2285. https://doi.org/10.3390/nu12082285 [PubMed]