- Open Access
Gut microbiota of the young ameliorates physical fitness of the aged in mice
Microbiome volume 10, Article number: 238 (2022)
Aging is a natural process that an organism gradually loses its physical fitness and functionality. Great efforts have been made to understand and intervene in this deteriorating process. The gut microbiota affects host physiology, and dysbiosis of the microbial community often underlies the pathogenesis of host disorders. The commensal microbiota also changes with aging; however, the interplay between the microbiota and host aging remains largely unexplored. Here, we systematically examined the ameliorating effects of the gut microbiota derived from the young on the physiology and phenotypes of the aged.
As the fecal microbiota was transplanted from young mice at 5 weeks after birth into 12-month-old ones, the thickness of the muscle fiber and grip strength were increased, and the water retention ability of the skin was enhanced with thickened stratum corneum. Muscle thickness was also marginally increased in 25-month-old mice after transferring the gut microbiota from the young. Bacteria enriched in 12-month-old mice that received the young-derived microbiota significantly correlated with the improved host fitness and altered gene expression. In the dermis of these mice, transcription of Dbn1 was most upregulated and DBN1-expressing cells increased twice. Dbn1-heterozygous mice exhibited impaired skin barrier function and hydration.
We revealed that the young-derived gut microbiota rejuvenates the physical fitness of the aged by altering the microbial composition of the gut and gene expression in muscle and skin. Dbn1, for the first time, was found to be induced by the young microbiota and to modulate skin hydration. Our results provide solid evidence that the gut microbiota from the young improves the vitality of the aged.
Aging is an inevitable deterioration of a structure and integrity of an organism, accompanying the functional decays in every aspect of physiology [1,2,3,4]. Understanding and regulating the aging process have been of great interest throughout the history [5, 6]. Manifestations of senescence occur not only at the cellular level such as shortened telomere and cell cycle arrest, but at the organ and organism levels—impaired cognitive function, dry eye, chronic gut inflammation, decreased skeletal muscle mass, skin wrinkling, et cetera .
Trillions of microorganisms reside in the mucosal surface of the host intestine, where they dynamically interact with a myriad of factors such as foods, host-derived molecules, and other microbes [7, 8]. The gut microbiota is considered a virtual organ of the host that plays a key role in maintaining immune and metabolic homeostasis [9, 10]; dysbiosis—disruption of microbiota homeostasis—has been repeatedly implicated in various diseases including inflammatory bowel disease, obesity, metabolic liver diseases, type 2 diabetes, cardiovascular diseases, cancers, and even neurological disorders [10, 11]. Accumulating evidence suggests that the gut microbiota changes with host aging . Moreover, the gut microbiota is emerging as a key player underlying several aging-related disorders [13, 14].
Recently, several studies have been conducted to demonstrate the causality of the gut microbiota on host health by transferring them to recipient groups of different ages. Some of them found that the gut microbiota derived from the aged promoted inflammation and in two cases cognitive decline as well [15,16,17]. Another study showed that the old-derived microbiota conferred obesogenic characteristics to the young germ-free mice . In the opposite direction, the gut microbiota from young African turquoise killifish was, for the first time in vertebrates, revealed to extend the life span of old killifish with improved muscle function and locomotor activity . Transcriptome analysis also indicated that gut aging is associated with increased inflammation and reduced proliferation. The killifish is a good model for aging-related studies of vertebrates as it only lives a few months. However, not only it is distantly related to mammals such as humans, but its gut microbiota differs from that of terrestrial vertebrates [20,21,22]. Further, underlying mechanisms for the counteracting effects of the young-derived microbiota on physical fitness including muscle and skin properties have yet to be explored in higher vertebrates.
The aim of this study is to systematically determine the revitalizing effects of the gut microbiota derived from the young on the physiology of the aged and to investigate the microbial and host genetic factors that play pivotal roles in the process. Thus, we transplanted the fecal microbiota from young (5 weeks after birth) mice to old (12 months) and very old (25 months) ones. We then monitored the gut microbiota profiles and examined the organ mass, serology, immune cellularity, gene expression profile, and aging manifestations in 2 months. We further analyzed microbial taxa that correlate with changes in the phenotypes and gene expression in mice. Mice heterozygous for Dbn1, which was the most upregulated gene in old mice after receiving the young-derived gut microbiota, were examined for selected phenotypes.
Animal experiments were approved by the Institutional Animal Care and Use Committee of Gachon University (IACUC number: GIACUC-R2019030) and the Institutional Animal Care and Use Committee of Yonsei University (IACUC number: 2018-0145), and the experiments were compliant with the Guide for the Care and Use of Laboratory Animals.
Reaching puberty, mammalian females are considered sexually mature after menarche, which occurs as early as 4 weeks after birth in mice [23, 24]. Reproductively senescent perimenopause females (between 9 and 12 months in mice) undergo a drastic decrease in sex hormone levels, which significantly contributes to the initiation of programmed aging [25,26,27]; menopause marks the final and irreversible termination of the female reproductive life span . Survivorship drops off markedly for mice over 24 months of age and can be considered very old . Based on these dispositions, female mice at 5 weeks were defined as young and two aged ones at 12 and 25 months as old and very old, respectively, in this study.
In late July 2019, C57BL6 female mice at 3 weeks, 12 months, and 25 months were obtained from the Laboratory Animal Resource Center, Daejeon, Korea, maintained in a specific pathogen-free (SPF) facility under a 12-h light/dark cycle. Mice were acclimated to the laboratory animal facility at Gachon University for 2 weeks before starting the experiment. Mice were fed with a normal Chow diet. Each individual mouse was tail-tagged for identification and randomly assigned to the treatment groups, which allows for matching different types of data including microbial structure, RNA sequencing, and aging-associated phenotypes.
Dbn1Het C57BL/6 mice (16-week-old) used for validation of Dbn1 for its involvement in skin properties were available from Korea Mouse Phenotyping Center (KMPC). Dbn1WT and Dbn1Het mice were maintained by crossbreeding Dbn1Het; complete loss of the Dbn1 gene (Dbn1-/-) resulted in embryonic lethality. Four weeks after birth, offspring were separated from the cage housing the mother and genotyped using the following primers: 5′–GTCCTTCCTCTTGGTCATTCCC and 5′–TGGAGAAACCAGGGAGATGTTG (for wild type), 5′–GCTACCATTACCAGTTGGTCTGGTGTC and 5′–CAGAGCCCAAGACTAATACACCC (for knockout).
Fecal microbiota transplantation
To prepare donor samples for fecal microbiota transplantation (FMT), we collected fresh feces after 2 weeks of acclimation by gently rubbing the belly from young (Y; 5 weeks: adolescence), old (O; 12 months: perimenopause), and very old (VO; 25 months: senescence) C57BL6 female mice, respectively, and maintained them under the strict anaerobic condition throughout the process. The feces were pooled by groups and mixed with PBS supplemented with 10% glycerol (30 mg feces/mL PBS; ca. 5 × 108 colony-forming units/mL) and thoroughly homogenized via vortexing. After sedimentation for 10 min, the supernatant was filtered through a 40-μm cell strainer and the homogenous aliquots were stored in a refrigerator at −80 °C.
Approximately 100 μL of the microbial suspension was transferred into each C57BL6 female mouse by oral gavage twice a week for a total of 8 weeks using autoclaved feeding needles. Superscript was used to describe the FMT donor. Thus, if feces from the young mice (Y) are given to 12-month-old or 25-month-old groups, they are designated as OldY or vOldY, respectively. As a control (C), 100 μL of PBS with 10% glycerol was administrated to old (OldC) or very old (vOldC) mice. In the same manner, feces from old or very old mice (O or VO) were respectively given to 12- or 25-month-old groups, resulting in OldO or vOldVO mice. All the recipient mice were housed under the SPF condition and weighed twice a week during the FMT timeline. Examination of passive avoidance, grip strength, skin properties, and visual abilities was performed 3 days after the last FMT. Then, all mice were euthanized and sacrificed for the measurement of physiological changes.
To generate microbiota-depleted old (xOld) mice by clearing out pre-existing gut microbiota, drinking water containing antibiotics (1 g/L ampicillin, 0.5 g/L vancomycin, 1 g/L metronidazole, 1 g/L neomycin, and 2.5 mg/L amphotericin B) was given to mice ad libitum during 7 days and withdrawn 24 h before FMT. When we tried the bacterial clearance in very old mice to make xvOld, they failed to tolerate the antibiotic treatment and died.
Fecal microbiota analysis
DNA was extracted using a Fast DNA SPIN kit for fecal samples (MP Biomedicals, Irvine, CA, USA), according to the manufacturer’s instructions. The V3-V4 region of the 16S ribosomal RNA (rRNA) gene was sequenced on an Illumina MiSeq platform by Chun Laboratory (Seoul, Korea) using sequence-specific primers (337F: CCTACGGGA(N)GGCWGCAG, 806R: GACTACHVGGGTM(A)TCTAAT). Raw sequencing data were analyzed using the QIIME2 pipeline (version 2017.10). The q2-demux plugin with DADA was used to de-multiplex and trim the adaptor sequences, and sequence variants were clustered into amplicon sequence variants (ASVs). The phylogenetic tree was constructed using q2-alignment and q2-phylogeny plugins. Sequence variants were taxonomically assigned using the SILVA reference database . The feature table and phylogenetic distances were imported into R for downstream analysis.
Alpha diversity indices including the number of observed ASVs, Chao1, and Shannon’s diversity were calculated using phyloseq and vegan R packages. Linear discriminant analysis effect size (LEfSe) was used to identify microbial biomarkers using the factorial Kruskal-Wallis test (P < 0.05) ; the threshold of logarithmic LDA score was 2.0. Spearman’s correlation coefficient was calculated for bacterial species that have more than 0.5% relative abundance in 12-month-old recipient mice using the psych R package.
Passive avoidance test
To examine the cognitive response of mice, a passive avoidance test was performed. A mouse was placed in a light chamber for 30 s and allowed to move to a dark chamber upon gate opening. The gate was immediately closed after the mouse moved to the dark chamber, and an electric foot shock (0.5 mA) was administered for 2 s. The time between the gate opening and entering of the mouse into the dark chamber was recorded as a training time. Then, the mouse was returned to the home cage. After 24 h, the mouse was returned to the light chamber for 30 s, and then, the gate of the dark chamber was opened. The time for the mouse to enter the dark chamber after the gate opening was recorded as a retention time.
Measurement of grip strength
Forelimb grip strength test was performed using GSM Grip-Strength Meter (Ugo Basile, 47200). The mouse was placed on an acrylic panel and the peak grip force of the mouse on the metal wire of the apparatus was measured in grams. The peak force was averaged from three independent trials.
Measurement of skin properties
The day before measuring skin hydration phenotypes, mouse dorsal hair was carefully shaved. Skin hydration was measured using MoistureMeterSC (Delfin) which registers a probe to the skin surface until skin hydration is stabilized. The capacitance is expressed in arbitrary units (AU). Transepidermal water loss (TEWL) was measured using VapoMeter (Delfin) that registers a probe to the skin surface until TEWL is stabilized (less than 1 min). The result is expressed in g/hm2.
Spatial frequency thresholds (i.e., visual acuity) were assessed by optokinetic nystagmus (OKN) using a virtual optokinetic system (OptoMotry, Cerebral Mechanics, Medicine Hat, Alberta, Canada). A video camera at the ceiling of the device records and transfers images to the connected computer. The clockwise movement of the drift grid tracked the movement of the left eye, and the counterclockwise movement tracked the reaction of the right eye of the mouse. The experimenter judged whether the head and body of the mouse were tracking the direction of the drift grid rotation. If tracking was unclear or absent, the process was repeated. The maximum spatial frequency capable of driving head tracking was determined.
For the intraocular pressure (IOP) test, mice were anesthetized using an intraperitoneal injection of xylazine (10 mg/kg; Rompun®, Bayer Animal Health) and zolazepam and tiletamine (30 mg/kg; Zoletil 50®, Virbac, Carros, France). IOP was measured using a rebound tonometer (Icare® TONOLAB tonometer, Colonial Medical Supply, Franconia, NH, USA), according to the manufacturer’s instructions. One trial result was obtained after six consecutive measurements, and the mean of consecutive trials was used for analyses.
Electroretinogram (ERG) analysis was performed using Micron Ganzfeld ERG (Phoenix Research Labs, Pleasanton, USA). The mouse was dark-adapted at least 12 h before the experiment for scotopic testing (rod cell response). After anesthesia, the pupil dilated as previously described. Once the pupil was fully dilated, we applied hypromellose 2.5% (Goniovisc®) and inserted the electrodes. ERG was recorded using Micron Ganzfeld ERG, according to the manufacturer’s instructions. Scotopic ERG was obtained with increasing flash intensity at a range of −1.7 log cd/s/m2 to 1.9 log cd/s/m2. Mice were light-adapted for 15 min prior to the cone cell response experiment. Photopic ERG was performed with increasing flash intensity at a range of −0.5 log cd/s/m2 to 4.1 log cd/s/m2. Values were based on the average of 10 responses to light stimuli. The implicit times of rod and cone cell response were determined.
For blood sample collection, the three groups of FMT recipients were euthanized using CO2 gas, and blood was directly obtained from the heart using 26G 1 mL syringe. To isolate serum from the total blood, blood samples were placed in the MiniCollect® tube (Greiner Bio-One, Kremsmünster, Austria) and centrifuged at 1 × 104 rpm for 5 min at 25 °C. The levels of aspartate aminotransferase (GOT), alanine aminotransferase (GPT), creatinine phosphokinase (CPK), alkaline phosphatase (ALP), blood urea nitrogen (BUN), amylase (AMYL), creatinine (CRE), albumin (ALB), triglyceride (TG), total bilirubin (TBIL), total cholesterol (TCHO), r-glutamyltransferase (GGT), lactate dehydrogenase (LDH), high-density lipoprotein-cholesterol (HDLC), and creatine kinase-mb (CKMB) in the serum samples (10 μL each) were measured using DRI-chem 4000i (Fuji, Minato, Japan) automated clinical chemistry analyzer.
The mouse tissues were fixed with ice-cold 4% paraformaldehyde and embedded in paraffin. The tissues were cut into 5-μm thick sections; deparaffinized with xylene three times for 20 min each, 100% EtOH three times for 10 min each, 90% EtOH twice for 10 min each, and 75% EtOH for 10 min; and stained with hematoxylin and eosin (H&E). Slides containing the stained tissue were dehydrated and mounted on Shandon Synthetic Mount (Thermo, Inc., Waltham, MA, USA). The diameter of the myofiber in the hindlimb skeletal muscle was determined using digital images of the sections using the ToupView program (http://www.touptek.com).
Immunohistochemical and immunofluorescence staining
For immunohistochemistry, the tissues were fixed with ice-cold 4% paraformaldehyde in PBS and mounted on paraffin blocks. Samples were cut into 3 μm sections, deparaffinized, and rehydrated in PBS. Antigens were then retrieved for 15 min at high pressure in Target Retrieval Solution (Dako, Santa Clara, CA, USA). Then, the specimens were chilled on ice for 1 h, washed with PBS three times for 5 min each and blocked with 3% H2O2 in PBS for 30 min to eliminate the endogenous peroxidase. The slides were washed with PBS, blocked for 2 h at 25 °C with Serum-Free Protein Block (Dako), probed at 4 °C overnight with the primary antibodies (1/1000 dilution), stained for 30 min with horseradish peroxidase-conjugated anti-mouse (Dako) or anti-rabbit IgG (Dako) secondary antibody, and developed with Liquid DAB+ Substrate Chromogen System (Dako). Finally, the specimens were counterstained with Mayer’s hematoxylin (Dako) and mounted on Shandon Synthetic Mount (Thermo).
For immunofluorescence, tissue sample preparation and blocking were performed in the same manner as described in immunohistochemistry. Then, the slides were probed at 4 °C overnight with the primary antibodies (1/1000 dilution), labeled for 2 h with Cy3-conjugated goat anti-rabbit IgG (Thermo) or A488 donkey anti-mouse IgG (Thermo), and stained with DAPI (Sigma, Burlington, MA, USA) for 15 min. Finally, the slides were mounted on ProLong Gold antifade reagent (Thermo), imaged on an Axio Imager M2 fluorescence microscope (Zeiss, Oberkochen, Land Baden-Württemberg, Germany), and analyzed using Zeiss Zen Blue Edition. The following primary antibodies were used for staining: anti-Ki67 (Abcam, Cambridge, UK), anti-DBN1 (Novus, Denver, CO, USA), anti-ITGB4 (Abcam), and anti-KRT10 (Abcam).
Multiplex cytokine and chemokine bead array assays
Multiplex bead array assays for serum cytokine quantification were performed using (#740446; BioLegend, San Diego, CA, USA) according to the manufacturer’s instructions. Serum samples were analyzed via a CytoFLEX platform (Beckman Coulter, Brea, CA, USA).
Immune cell preparation
Single-cell suspensions of the spleen, lymph nodes, and liver were prepared by mechanical homogenization of the organs and passing through 40-μm or 70-μm cell strainer (BD Biosciences, San Jose, CA, USA). For the spleen, erythrocytes were removed using ACK lysing buffer (Thermo). Liver cells were purified by density gradient centrifugation on a 67% and 44% Percoll gradient (GE Healthcare, Chicago, IL, USA).
Cells were seeded in a 96-well plate and stained with fluorochrome-conjugated antibodies for 20 min at 4 °C in the dark. Cells were subsequently washed with PBS containing 2% fetal bovine serum. The following antibodies were used: CD8 (53-6.7; BD Biosciences), CD62L (MEL-14; BD Biosciences), CCR7 (4B12; eBioscience, San Diego, CA, USA), CD103 (2E7; eBioscience), CD25 (PC61; BioLegend), CD45.2 (104; BioLegend), CD44 (IM7; BioLegend), CD69 (H1.2F3; BioLegend), and CD4 (RM4-5; BioLegend). Cells were stained with LIVE/DEAD® Fixable Blue Dead Cell Stain (Invitrogen, Carlsbad, CA, USA) to discriminate dead cells. All stained samples were acquired using CytoFLEX LX instrument (Beckman Coulter), and the data were analyzed using the FlowJo software (BD Biosciences).
Total RNA was isolated from the distal colon, quadriceps femoris muscle (quad muscle), and back skin, respectively, using Trizol reagent (Invitrogen). RNA quality was assessed by Agilent 2100 bioanalyzer using the RNA 6000 Nano Chip (Agilent Technologies, Amstelveen, The Netherlands), and RNA quantification was performed using ND-2000 Spectrophotometer (Thermo). For messenger RNA sequencing, a complementary DNA library was constructed using QuantSeq 3’ mRNA-Seq Library Prep Kit (Lexogen, Inc., Austria), according to the manufacturer’s instructions. Briefly, 500 ng of the total RNA was prepared and an oligo-dT primer containing an Illumina-compatible sequence at its 5′ end was hybridized to the RNA and reverse transcribed. After degradation of the RNA template, second strand synthesis was initiated using a random primer containing an Illumina-compatible linker sequence at its 5′ end. The double-stranded DNA library was purified using magnetic beads to remove all reaction components. The library was amplified to obtain the complete adapter sequences required for cluster generation. The finished library was purified to remove PCR components. High-throughput sequencing was performed using single-end 75-bp reads using NextSeq 500 (Illumina, Inc., USA).
QuantSeq 3′ mRNA-Seq reads were aligned using Bowtie2. Bowtie2 indices were either generated from the genome assembly sequence or representative transcript sequences for aligning to the genome and transcriptome. The alignment file was used for assembling transcripts, estimating their abundances, and detecting differential expression of genes. Raw read counts were normalized and differentially expressed genes (DEGs) were identified using the edgeR package (P < 0.01). The DEGs were determined based on read counts from unique and multiple alignments using coverage in Bedtools. Functional annotation of DEGs was performed using DAVID (http://david.abcc.ncifcrf.gov/) and MEDLINE (http://www.ncbi.nlm.nih.gov/) databases.
Complete hierarchical clustering based on Euclidean distance was performed to classify co-occurred species into clusters using the hclust function in the R package, and four clusters were obtained by cutting off the dendrogram tree at height 6. Spearman’s correlation coefficient was calculated to examine the relationship between DEGs and microbial composition. P values were corrected for multiple comparison problems using Benjamini-Hochberg method with a 0.05 cutoff.
Statistical analyses of host phenotype data were performed using Prism software (GraphPad). Raw data were subjected to one-way ANOVA to evaluate statistical significance between at least three groups, and pairwise comparison was conducted using Student’s t test. Means were considered significantly different at P < 0.05. Statistical details, including sample sizes for each experiment, are provided in the relevant figure legends.
Microbiota of the young improves the muscle and skin properties of the aged
To investigate whether the phenotypes of aged mice could be improved by the gut microbiota of the young, we first designated three different age groups of female C57BL/6 mice and compared the baseline microbiota via rRNA gene sequencing. The three groups from which fresh fecal samples were collected are as follows: 5 weeks after birth for young (represents adolescence), 12 months for old (perimenopause), and 25 months for very old (senescence). Differentially abundant microbes included Prevotellaceae (belongs to Bacteroidetes) and Ruminococcaceae (Oscillospiraceae; Firmicutes) in the young mice and Erysipelotrichaceae (Firmicutes) and Clostridaceae 1 (Firmicutes) in the two old groups; however, the abundance of Muribaculaceae (Bacteroidetes) was comparable between the groups (Fig. 1A).
Fecal contents were repeatedly transferred to mice at 12-month-old mice (Old) by oral gavage for 8 weeks (twice a week) to ensure the establishment of donor microbiota: vehicle control (C), 5 weeks (Y), and 12 months (O) (Fig. 1B). The fecal microbiota of the recipient mice was longitudinally collected at weeks 0, 4, and 8. The abundance of bacterial families including Muribaculaceae, Lachnospiraceae, Lactobacillaceae, Bacteroidaceae, Ruminococcaceae, and Rikenellaceae differently changed over time depending on the source of microbiota (Fig. 1C). To screen for the microbial biomarkers that are specific for each group, LEfSe was applied to the 8-week microbiota data. Results showed that species belonging to Muribaculaceae, Bacteroidaceae, and Prevotellaceae are most discriminative in old mice that received microbiota derived from the young (OldY), whereas species of Lachnospiraceae, Desulfovibrionaceae, Ruminococcaceae, and Tannerellaceae were discriminative in old mice that received microbiota from the old (OldO) (LDA effect size ≥ 2.0; Fig. 1D). There was no significant difference in observed ASVs, Chao1, and Shannon’s diversity index across the recipient groups (Fig. S1A).
During the experiment, the recipient mice appeared healthy and maintained their body and organ weights, implying that the repeated oral gavages were not excessively stressful for the recipients (Fig. S1B, C). Data from the blood chemistry analysis showed that FMT did not change the levels of most systemic metabolic indicators except a slight increase in amylase level in OldY (Fig. S1D). In addition, the exogenous microbiota did not alter two well-known aging-associated phenotypes, which are avoidance-learning memory and visual abilities such as visual acuity, IOP, and response time of rod and cone cells (Fig. S1E, F). However, when we performed forelimb grip strength test that is a common method to measure neuromuscular function of rodents, OldY mice had 30–50% stronger grip strength than OldC and OldO (PYvsO = 0.000079; Fig. 1E). In line with the result, the diameter of the muscle fiber was significantly increased in OldY, suggesting the improved muscle function by the young-derived microbiota (PYvsO = 0.0042; Fig. 1F). In skin, the number of stratum corneum layer (outermost barrier of the epidermis), was highly increased in OldY (5.3±2.23 in OldY, 2.83±1.36 in OldO; Fig. 1G). Consequently, the ability to retain moisture was dramatically enhanced as shown in the skin hydration assessment (7.47±1.29 AU in OldY, 3.66±0.38 AU in OldO) and TEWL measurement (3.71±0.52 g/hm2 in OldY, 6.24±1.34 g/hm2 in OldO; Fig. 1H). The epithelium of OldY had more Ki-67-positive proliferating cells, and in situ expression of keratin 10 (KRT10) was also elevated in OldY compared to OldO, but integrin beta 4 (ITGB4) was comparable (Fig. 1I–K).
To investigate the effects of the young-derived microbiota on the very old mice, we also performed young-to-old microbiota transfer with 25-month-old mice (vOld) using three donor groups: vehicle control, 5 weeks (Y), and 25 months (VO) (Fig. S2A). During 0, 4, and 8 weeks, most of the microbial families that are abundant in Old were present in vOld as well, with the exception of Bifidobacteriaceae in vOld (Fig. S2B). LEfSe indicated that biomarkers belonging to Prevotellaceae and Muribaculaceae that were retrieved in OldY were replicated in vOldY at week 8, as well as Roseburia (Lachnospiraceae) in vOldVO, but a species of Eubacteriaceae was newly represented in vOldVO (Fig. S2C). There was no significant difference in the number of observed ASVs, Chao1, and Shannon’s diversity index across the recipients (Fig. S2D). No significant changes were observed in body weight, organ mass, and the number of stratum corneum layers; only the relative brain weight (PYvsVO = 0.004) and the thickness of muscle fiber (PYvsVO = 0.085) were increased in vOldY (Fig. S2E–H).
Further, we designed an experiment in which mice were pre-treated with an antibiotic cocktail (Abx) before FMT to remove the resident gut microbiota and to exclude the involvement of pre-existing microbes (Fig. 2A). By 4 weeks after the first FMT to Abx-treated 12-month-old mice (xOld), the taxonomic structure was restored harboring high proportion of Muribaculaceae, Lachnospiraceae, Lactobacillaceae, and Bacteroidaceae, which was similar to those of the gut microbiome of SPF mice described above and the common murine gut microbiome (Fig. 2B). Consistent with the results of SPF mice, Alloprevotella (Prevotellaceae) was found to be most discriminative in Abx-treated old mice that received microbiota from the young (xOldY), and Akkermansia that is predominant in microbiota derived from the young was also observed (Fig. 2C). Body weight of the recipients was stable during FMT, and most of the phenotypic changes were reproduced without pre-existing microbiota. There were no significant differences in organ weights except brain in xOldY mice (2.04±0.11 % organ weight in xOldY, 1.91±0.07 % in xOldO), and the levels of metabolic indicators except amylase and aversive memory were not changed (Fig. 2D–G). Most importantly, muscle fitness and skin properties were consistently improved (Fig. 2H–K). Very old (25 months) mice that also received the same dose of Abx for microbiota clearance failed to tolerate the toxicity of Abx and most died; results from xvOld were excluded in further analysis.
Although FMT has been reportedly implicated in the intestinal immune responses , our study showed that there was no significant alteration in the serum cytokine profile or T cell population in the spleen, liver, and mesenteric or inguinal lymph nodes (Fig. S3). These results highlight that transplanting the young-derived microbiota can rejuvenate the aging-associated phenotypes of the skin and muscle without eliciting severe inflammation.
Young microbiota alters host gene expression in the colon, muscle, and skin
To characterize the changes in host gene expression after FMT, we performed RNA sequencing (RNA-seq) of three tissue homogenates of the colon, quad muscle, and back skin samples at week 8. Analysis of the RNA-seq data from the colon showed 89 and 91 genes expressed more in either OldY or OldO, respectively (Fig. 3A). More genes annotated to be involved in cell differentiation or K+ ion transport were upregulated in OldY, and more death-associated genes were expressed higher in OldO (Fig. 3B and Table S1). Furthermore, gene set enrichment analysis (GSEA) based on the KEGG database showed that receptor-signaling pathways including those related to a serotonin receptor and a G-protein-coupled neurotransmitter receptor had more genes significantly upregulated in OldY than OldO; however, no inflammation-related signaling pathways were detected (Fig. 3C).
In the muscle, we found that 62 and 109 genes, including Tfap2b and Serpinb6c, respectively, were expressed more in either OldY or OldO (P < 0.05; Fig. 3D). More genes involved in cell differentiation, proliferation, fatty acid synthesis, or neutral lipid metabolism pathways showed higher expression in OldY (Fig. 3E and Table S2), and glutathione (the most abundant intracellular anti-oxidant thiol) metabolism was overrepresented in OldY according to GSEA (Fig. 3F). In skin, 47 and 73 genes, including Dbn1 and Serpinb6c, respectively, were expressed more in OldY or OldO (Fig. 3G). While more aging- or death-associated genes were expressed higher in OldO, more genes involved in cell differentiation and migration were highly expressed in OldY (Fig. 3H and Table S3). In addition, GSEA showed that more genes related to cell junction maintenance and epithelial structure maintenance were expressed in OldY than in OldO (Fig. 3I).
Certain microbial clusters link to rejuvenation
To define a consortium of core microbial taxa that are involved in aging-associated phenotypes, we clustered the entire families of the gut microbiota at week 8 according to microbial abundance correlation (Fig. 4A). Among the four clusters retrieved by hierarchical clustering, Lachnospiraceae-dominant cluster 1 was prevalent in OldC and OldO groups, whereas cluster 2 containing Muribaculaceae, Bacteroidaceae, Lactobacillaceae, and Prevotellaceae was abundant in OldY. While cluster 3 that distinctly harbored Rikenellaceae, Ruminococcaceae, and Erysipelotrichaceae was enriched in OldY, it was also significantly enriched in OldO compared to non-FMT control. In addition, Desulfovibrionaceae-containing cluster 4 had no preference for any treatments, implying a minor impact of its members on the host (Fig. 4A).
To better understand the association of microbial clusters with host physiology, we calculated Spearman’s correlation coefficient between the cumulative abundance of clusters and the indicators of host fitness (Fig. 4B). Prevotellaceae and Muribaculaceae in clusters 2 and 4 showed a positive correlation with both grip strength and skin moisture, the most improved phenotypes in OldY. Bacteroidaceae in cluster 2 positively correlated with skin moisture and brain size, and Erysipelotrichaceae in cluster 3 and Akkermansiaceae in cluster 4 with the inverse of TEWL (TEWL-1). In contrast, Ruminococcaceae in cluster 1 negatively correlated with grip strength and TEWL-1, Tannerellaceae in cluster 4 with skin moisture, TEWL-1, and Ki-67+, and Ruminococcaceae in cluster 4 with skin moisture and brain size (Fig. 4B). Thus, Bacteroidetes families in clusters 2 and 3 are strongly associated with restored phenotypes of the aged hosts, while Firmicutes families mainly present in clusters 1 and 4 may promote functional deterioration.
Next, we explored the correlation between tissue-specific DEGs in OldY and OldO and clustered taxa. Muribaculaceae, Prevotellaceae, and Bacteroidaceae (all in the phylum Bacteroidetes) in cluster 2, as well as Ruminococcaceae and Erysipelotrichaceae in cluster 3, had a strong association with upregulated genes in the muscle and skin of OldY (Fig. 4C). Conversely, upregulated genes in OldO negatively correlated with cluster 2 families, but positively with Lachnospiraceae in clusters 1 and 4 and Ruminococcaceae in cluster 1. Similar patterns were observed in the colon tissues (Fig. S4A). Correlation coefficients between the top DEGs and bacterial clusters were compared in each tissue. Cluster 2 showed the highest positive correlation with Tfap2b and Dbn1, which are most upregulated in the muscle and skin of OldY, respectively. Serpinb6c and Rreb1, most downregulated in each tissue, respectively, negatively correlated with cluster 2 (Fig. 4D). Similar results were produced for Trpm3 and Wdr54 in the colon (Fig. S4B). Consistently, members of Muribaculaceae, Prevotellaceae, and Bacteroidaceae were associated with increased gene expression and improved host phenotypes in OldY, suggesting that they may have an ability to rejuvenate host physiology.
Dbn1 is involved in water retention in the skin
Notably, among the sets of DEGs, Dbn1 in the skin of OldY was most upregulated across the three tissues (171.1 fold increase; P = 0.001) (Table S3). Thus, we examined the in situ expression of DBN1 via immunofluorescence and found that the dermis of OldY harbored about two times more DBN1-positive (DBN1+) cells than those of OldC and OldO (Fig. 5A, B). To validate the functionality of Dbn1, we measured the water retention capacity of the wild type and Dbn1-heterozygous (Dbn1Het) C57BL/6 littermates (8 weeks old). While the levels of skin hydration and TEWL were comparable between young Dbn1WT and OldY, the halved gene dosage of Dbn1 resulted in impairment of skin hydration and barrier function as much as those observed in OldO (Fig. 5C). These data indicate that the gut microbiota derived from the young improves water-holding capacity in perimenopause female mice by modulating host gene expression like upregulation of Dbn1.
Perspectives on understanding aging are important for intervention and management of the process. The concept of “programmed aging” that considers aging as a designed biological program allows, at least partially, the rejuvenation of old body through the modulation of effector(s) . Since Smith and colleagues first reported extension of life span and delayed behavioral decline in short-lived killifish by the young-derived microbiota , a few groups have added evidence on the counteracting effects of fecal microbiota transfer on hallmarks of aging gut, brain, and retina using the mouse model [34,35,36]. Our results highlight that aging phenotypes of physical fitness in aged mice can be ameliorated by transplanting the young-derived microbiota.
Indicators of muscle and skin properties including forelimb grip strength, the thickness of muscle fiber, number of stratum corneum layers, skin water retention, and skin cell markers (Ki-67 and KRT10) were all improved in old mice at 12 months with or without Abx pre-treatment. Further, a statistically significant increase in muscle fiber and brain size was observed in very old mice at 25 months. A couple of studies reported that gut microbiota is implicated in the skeletal muscle function. Germ-free mice suffer from reduced muscle mass and locomotor activity as compared to mice raised under the SPF condition , and xenogeneic transplantation of the pig gut microbiota to germ-free mice synchronized lipid metabolisms in muscle, collectively supporting the effect of the gut microbiota in the muscle [27, 37].
Among all organs measured, only the brain, usually weighing 2% of the body weight, showed a moderate change by fecal microbiota transfer. Old mice at 12 months receiving the gut microbiota from the young maintained their relative brain mass regardless of Abx pre-treatment, but those receiving the old-derived microbiota showed a slight decrease. Despite the change, avoidance-learning memory was not different among them, implying that the behavioral relevance of the reduced brain size in xOldO is trivial. Of note, there was a significant increase in brain size for vOldY as compared to vOldVO and vOldC; however, the physical activity of the very old mice was too low for the passive avoidance test as well as the grip strength test. Given the pre-documented association between the gut microbiota and brain function [35, 36], these results may serve as preliminary data for a further study on the physical and functional fitness of the aging brain. The two studies revealed the beneficial effects of the young microbiota, altering the immune profiles in the gut and brain. The results are in line with the “inflammaging”| theory; however, it is hard to distinguish the true anti-inflammatory activity of the young-derived microbiota from the absence of the old-derived immunogenic molecules. Transplantation of the fecal microbiota from old mice did not induce significant changes in body weight, organ mass, the level of biomarkers for tissue damage, pro-inflammatory cytokines, and the number of T cells in our analyses (data not shown), suggesting that ours are not mediated by the lack of the old-derived immunogenicity.
We identified that a group of microorganisms consisting of Muribaculaceae, Prevotellaceae, and Bacteroidaceae (all Bacteroidetes) represents the young-derived microbiota that are discriminative in old mice at 12 months after fecal transfer and strongly correlates with enhanced gene expression of the recipients’ muscle and skin. This is in line with our big data analysis of the human gut microbiota with conserved enrichment of Bacteroidaceae in the young population (unpublished data) and a previous report on high frailty of old people having reduced Bacteroides and Prevotella . Moreover, these three families co-occurring in a Bacteroidetes-rich cluster based on hierarchical clustering are significantly associated with higher grip strength and improved skin water retention. Microbial taxa more abundant in aged mice and likely to contribute to aging phenotypes were also defined. Some species of Lachnospiraceae and Ruminococcaceae in the phylum Firmicutes were consistently over-represented in old mice that received the old-derived microbiota, and correlated with upregulated genes in the muscle and skin. In addition, a few members of Ruminococcaceae are negatively associated with grip strength, skin water retention, and brain size. Our results provide a list of causative bacteria that may induce host gene expression and subsequent phenotypic changes with regards to revitalizing the aged.
Finally, we elucidated that Dbn1, highly induced in the skin of OldY mice, exerts a protective role in skin hydration. Although DBN1, a cytoplasmic actin-binding protein, is known to be expressed ubiquitously, its function has been documented mainly in neurobiology such as neuronal stem cell differentiation and cognitive function in Alzheimer’s disease [38, 39]. Our results are the first evidence that relates Dbn1 to skin hydration; however, considering the predefined role of Dbn1, it is plausible that DBN1 may facilitate the proliferation or differentiation of epidermal stem cells, leading to the improved skin phenotypes, which awaits further investigation.
The interaction between the gut microbiome and the host has been shown to modulate host physiology, which inspired people to examine the revitalizing potential of the gut microbiota from the young. While impaired physical fitness such as sarcopenia and skin barrier dysfunction in older people has corroborated the “damage theory” of aging, our study for the first time demonstrates that microbial intervention strengthens the aged muscle and rejuvenates the skin through changes in the microbial structure and host gene expression. We proposed Dbn1 as an important player in maintaining skin integrity and further identified microbial members related to pro- or anti-aging phenotypes. Our study provides comprehensive evidence of the revitalization of multiple organs in the host and points to host and microbial targets for anti-aging treatment.
Availability of data and materials
The sequencing reads of the 16S rRNA gene for the mouse fecal microbiota have been deposited in GenBank under the BioProject number PRJNA789433 (https://dataview.ncbi.nlm.nih.gov/object/PRJNA789433?reviewer=p6bokr5gnifd4atmgjj9f6dik2).
Amplicon sequence variant
Blood urea nitrogen
Differentially expressed gene
Fecal microbiota transplantation
Gene set enrichment analysis
Hematoxylin and eosin
Kyoto Encyclopedia of Genes and Genomes
Linear discriminant analysis
Linear discriminant analysis effect size
Operational taxonomic unit
Principal coordinate analysis
Transepidermal water loss
López-Otín C, Blasco MA, Partridge L, Serrano M, Kroemer G. The hallmarks of aging. Cell. 2013;153:1194–217.
McLean AJ, Le Couteur DG. Aging biology and geriatric clinical pharmacology. Pharmacol Rev. 2004;56:163–84.
Cruz-Jentoft AJ, Bahat G, Bauer J, Boirie Y, Bruyere O, et al. Sarcopenia: revised European consensus on definition and diagnosis. Age Ageing. 2019;48:16–31.
Marengoni A, Angleman S, Melis R, Mangialasche F, Karp A, et al. Aging with multimorbidity: a systematic review of the literature. Ageing Res Rev. 2011;10:430–9.
Depp CA, Jeste DV. Definitions and predictors of successful aging: a comprehensive review of larger quantitative studies. Am J Geriatr Psychiatry. 2006;14:6–20.
Cruz-Jentoft AJ, Landi F, Schneider SM, Zúñiga C, Arai H, et al. Prevalence of and interventions for sarcopenia in ageing adults: a systematic review. Report of the International Sarcopenia Initiative (EWGSOP and IWGS). Age Ageing. 2014;43:748–59.
Gilbert JA, Blaser MJ, Caporaso JG, Jansson JK, Lynch SV, et al. Current understanding of the human microbiome. Nat Med. 2018;24:392–400.
Sommer F, Bäckhed F. The gut microbiota--masters of host development and physiology. Nat Rev Microbiol. 2013;11:227–38.
Cani PD. Human gut microbiome: hopes, threats and promises. Gut. 2018;67:1716–25.
Cho I, Blaser MJ. The human microbiome: at the interface of health and disease. Nat Rev Genet. 2012;13:260–70.
Fan Y, Pedersen O. Gut microbiota in human metabolic health and disease. Nat Rev Microbiol. 2021;19:55–71.
Nagpal R, Mainali R, Ahmadi S, Wang S, Singh R, et al. Gut microbiome and aging: Physiological and mechanistic insights. Nutr Healthy Aging. 2018;4:267–85.
Wilmanski T, Diener C, Rappaport N, Patwardhan S, Wiedrick J, et al. Gut microbiome pattern reflects healthy ageing and predicts survival in humans. Nat Metab. 2021;3:274–86.
Bana B, Cabreiro F. The microbiome and aging. Annu Rev Genet. 2019;53:239–61.
Fransen F, van Beek AA, Borghuis T, Aidy SE, Hugenholtz F, et al. Aged gut microbiota contributes to systemical inflammaging after transfer to germ-free mice. Front Immunol. 2017;8:1385.
Li Y, Ning L, Yin Y, Wang R, Zhang Z, et al. Age-related shifts in gut microbiota contribute to cognitive decline in aged rats. Aging. 2020;12:7801–17.
D'Amato A, Di Cesare ML, Lucarini E, Man AL, Le Gall G, et al. Faecal microbiota transplant from aged donor mice affects spatial learning and memory via modulating hippocampal synaptic plasticity- and neurotransmission-related proteins in young recipients. Microbiome. 2020;8:140.
Binyamin D, Werbner N, Nuriel-Ohayon M, Uzan A, Mor H, et al. The aging mouse microbiome has obesogenic characteristics. Genome Med. 2020;12:87.
Smith P, Willemsen D, Popkes M, Metge F, Gandiwa E, et al. Regulation of life span by the gut microbiota in the short-lived African turquoise killifish. Elife. 2017;6:e27014.
Douglas AE. Simple animal models for microbiome research. Nat Rev Microbiol. 2019;17:764–75.
Lee WJ, Brey PT. How microbiomes influence metazoan development: insights from history and Drosophila modeling of gut-microbe interactions. Annu Rev Cell Dev Biol. 2013;29:571–92.
Kostic AD, Howitt MR, Garrett WS. Exploring host-microbiota interactions in animal models and humans. Genes Dev. 2013;27:701–18.
Silver LM. Mouse genetics – concepts and applications. 1st ed: Oxford University Press; 1995.
Dutta S, Sengupta P. Men and mice: relating their ages. Life Sci. 2016;152:244–8.
Felicio LS, Nelson JF, Finch CE. Longitudinal studies of estrous cyclicity in aging C57BL/6J mice: II. Cessation of cyclicity and the duration of persistent vaginal cornification. Biol Reprod. 1984;31:446–53.
Cruz G, Fernandois D, Paredes AH. Ovarian function and reproductive senescence in the rat: role of ovarian sympathetic innervation. Reproduction. 2017;153:R59–68.
Lahiri S, Kim H, Garcia-Perez I, Reza MM, Martin KA, et al. The gut microbiota influences skeletal muscle mass and function in mice. Sci Transl Med. 2019:11:eaan5662.
Gosden RG. Menopause. In: Birren JE, editor. Encyclopedia of gerontology. 2nd ed. New York: Elsevier; 2007. p. 151–9.
Flurkey K, Currer JM, Harrison DE. Mouse models in aging research. In: Fox JG, editor. The mouse in biomedical research. 2nd ed: American College of Laboratory Animal Medicine; 2007. p. 637–72.
Quast C, Pruesse E, Yilmaz P, Gerken J, Schweer T, et al. The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res. 2013;41:D590–6.
Segata N, Izard J, Waldron L, Gevers D, Miropolsky L, et al. Metagenomic biomarker discovery and explanation. Genome Biol. 2011;12:R60.
Quraishi MN, Shaheen W, Oo YH, Iqbal TH. Immunological mechanisms underpinning faecal microbiota transplantation for the treatment of inflammatory bowel disease. Clin Exp Immunol. 2020;199:24–38.
Jin K. Modern biological theories of aging. Aging Dis. 2010;1:72–4.
Stebegg M, Silva-Cayetano A, Innocentin S, Jenkins TP, Cantacessi C, et al. Heterochronic faecal transplantation boosts gut germinal centres in aged mice. Nat Commun. 2019;10:2443.
Boehme M, Guzzetta KE, Bastiaanssen TFS, van de Wouw M, Moloney GM, et al. Microbiota from young mice counteracts selective age-associated behavioral deficits. Nat Aging. 2021;1:666–76.
Parker A, Romano S, Ansorge R, Aboelnour A, Le Gall G, et al. Fecal microbiota transfer between young and aged mice reverses hallmarks of the aging gut, eye, and brain. Microbiome. 2022;10:68.
Yan H, Diao H, Xiao Y, Li W, Yu B, et al. Gut microbiota can transfer fiber characteristics and lipid metabolic profiles of skeletal muscle from pigs to germ-free mice. Sci Rep. 2016;6:31786.
Ao X, Liu Y, Qin M, Li C, Chen X, et al. Expression of Dbn1 during mouse brain development and neural stem cell differentiation. Biochem Biophys Res Commun. 2014;449:81–7.
Liu Y, Xu YF, Zhang L, Huang L, Yu P, et al. Effective expression of Drebrin in hippocampus improves cognitive function and alleviates lesions of Alzheimer's disease in APP (swe)/PS1 (ΔE9) mice. CNS Neurosci Ther. 2017;23:590–604.
We thank In Kwon Chung and Young Jun Oh for helping with the discussions, Boyoung Lee and Lae-Guen Jang for their involvement in the project, Sung-Rae Cho for supplying the Grip-Strength Meter, Dae Youn Hwang for assisting with the passive avoidance test, Irin M. Kim for providing mouse illustrations, and ebiogen for the next-generation sequencing service.
This study was financially supported by grants from the National Research Foundation of Korea (NRF-2021M3H9A2098025, NRF-2021M3A9I4021432, and NRF-2014M3C9A33068822) and the Yonsei Signature Research Cluster Program to J.F.K., and the NRF of Korea (NRF-2022R1A2C3007850 and NRF-2022M3A9F3016364) and the Korea Mouse Phenotyping Project (NRF-2016M3A9D5A01952416) to K.T.N. Research and publication were also supported in part by the Brain Korea 21 project, and Y.C. and J.K.Y. are fellowship awardees of this program.
Ethics approval and consent to participate
Consent for publication
All authors read and approved this version of the manuscript for publication.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file 1: Figure S1.
Physiological characteristics of old mice exposed to the three treatments. Figure S2. Physiological characteristics of very old mice transplanted with microbiota from control, very old, and young mice. Figure S3. Immune profile of FMT mice. Figure S4. Correlation analysis between microbial clusters and DEGs in the colon tissue.
Additional file 2: Table S1.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Kim, K.H., Chung, Y., Huh, JW. et al. Gut microbiota of the young ameliorates physical fitness of the aged in mice. Microbiome 10, 238 (2022). https://doi.org/10.1186/s40168-022-01386-w
- Programmed aging