- Open Access
Microbial shifts in the aging mouse gut
Microbiome volume 2, Article number: 50 (2014)
The changes that occur in the microbiome of aging individuals are unclear, especially in light of the imperfect correlation of frailty with age. Studies in older human subjects have reported subtle effects, but these results may be confounded by other variables that often change with age such as diet and place of residence. To test these associations in a more controlled model system, we examined the relationship between age, frailty, and the gut microbiome of female C57BL/6 J mice.
The frailty index, which is based on the evaluation of 31 clinical signs of deterioration in mice, showed a near-perfect correlation with age. We observed a statistically significant relationship between age and the taxonomic composition of the corresponding microbiome. Consistent with previous human studies, the Rikenellaceae family, which includes the Alistipes genus, was the most significantly overrepresented taxon within middle-aged and older mice.
The functional profile of the mouse gut microbiome also varied with host age and frailty. Bacterial-encoded functions that were underrepresented in older mice included cobalamin (B12) and biotin (B7) biosynthesis, and bacterial SOS genes associated with DNA repair. Conversely, creatine degradation, associated with muscle wasting, was overrepresented within the gut microbiomes of the older mice, as were bacterial-encoded β-glucuronidases, which can influence drug-induced epithelial cell toxicity. Older mice also showed an overabundance of monosaccharide utilization genes relative to di-, oligo-, and polysaccharide utilization genes, which may have a substantial impact on gut homeostasis.
We have identified taxonomic and functional patterns that correlate with age and frailty in the mouse microbiome. Differences in functions related to host nutrition and drug pharmacology vary in an age-dependent manner, suggesting that the availability and timing of essential functions may differ significantly with age and frailty. Future work with larger cohorts of mice will aim to separate the effects of age and frailty, and other factors.
The human microbiome influences and is influenced by several aspects of the host’s health and development . Perturbations to human-associated bacterial communities are associated with many disorders such as colon cancer, autoimmune diseases, inflammatory bowel disease, and Clostridium difficile infection -. The influence of the gut microbiota on human health is driven by interactions between microbes and the host: for example, different groups of bacteria can synthesize energy sources such as butyrate and other short-chain fatty acids ,, stimulate the immune system ,, and provide protection from pathogens through competitive exclusion and the production of protective compounds such as bacteriocins ,.
The age of the host appears to be linked with the composition of the associated microbiome ,. Although many studies have focused on associations between the microbiome and early life stages ,,, relatively few studies have looked into the effect of the microbiome on aging and frailty in later life -. It is known that the microbiome changes drastically between infant and adult stages of life, with a shift from dominance by Bifidobacterium to genera within Bacteroidetes and Clostridia , that reflects a change from primarily lactate metabolism and plant polysaccharide breakdown to short-chain fatty acid (SCFA) production and vitamin (such as cobalamin/B12) and carbohydrate metabolism. Previous studies have observed changes to the bacterial communities between young/middle-aged adults and older subjects ,-. Centenarians were found to have a decrease in Clostridia and an increase in Proteobacteria and Bacilli in their gut microbiomes . Such shifts reduce the abundance of SCFA producers and increase the number of facultative anaerobes and opportunistic pathogens (primarily Proteobacteria) within the gut microbiome of the older population . Such shifts would likely lead to increased inflammation, which along with aging-associated changes is a major contributor to the overall frailty of an individual . Studies focusing on the frailty of the host as a factor in microbiome composition found that members of the Oscillibacter and Alistipes genera were in high abundance in the most frail individuals  whereas the abundance of Eubacteriaceae, Faecalibacterium, and Lactobacillus was reduced . Although these shifts in constituent microbes are likely associated with changes in microbial metabolism and interactions with the host, the consequences of such shifts are not yet understood. Beyond these studies, little is known to date about the link between the microbiome, aging and frailty in terms of perturbations to the microbial communities. Since different microbial taxa often carry out different molecular functions and take on different ecological roles in the gut, we may expect to see commensurate changes in the metabolic potential of the microbes and the manner in which they interact with the host.
Given the potential for the gut microbiome to change as an individual ages, and its possible role in the health of older individuals, our aim is to investigate linkages between aging, frailty, and the microbiome. There are many probable confounding factors in human studies, including change in diet, medications, and housing status (e.g., home residence versus long-term care facilities), making it difficult to identify direct effects of aging and frailty. Although not free of confounding factors themselves, mouse studies allow for better-controlled observations and experiments and have been shown to serve as good models of the human microbiome ,. A frailty model has recently been demonstrated to provide similar information about aging mice  as traditional frailty assessments do for older human patients . A previous study of the gut microbiome in older mice related dietary intervention to shifts in bacterial composition  but did not investigate changes in relation to age itself. Using metagenomic sequencing of fecal samples, we find shifts both in microbial composition and specific molecular functions that correlate with the age and frailty of the host. Our results corroborate previous findings of microbiome shifts in aging adults and identify changes in taxonomy and function with relevance to aging and frailty.
Results and discussion
Murine age groupings and frailty index scores
A total of 21 stool samples for metagenomic analysis were collected from ten different mice with varying ages and murine frailty index (FI) scores (Additional file 1). Age correlated with frailty (Spearman correlation = 0.86, p = 1.064 × 10−5), and fell into three natural groupings which we refer to as ‘young’ (age in days: mean 174 ± s.d. 15; FI: 0.024 ± 0.016; samples = 9; mice = 5), ‘middle’ (age in days: 589 ± 18; FI: 0.097 ± 0.030; samples = 6; mice = 2), and ‘old’ (age in days: 857 ± 16; FI: 0.302 ± 0.088, samples = 6; mice = 3) (Figure 1, Additional file 2). The strong correlation between FI and age reinforces that both are suitable measurements to compare with changes in the gut microbiome. Fares and Howlett showed that a 50% mortality rate occurs around 24 months in mice, which corresponds roughly to age 85 in humans . Our study did not reveal a strong indication that either FI or age was a better predictor of microbiome shifts due to the sample size of the study and the lack of observed variation in FI with respect to age in this particular subset of mice.
Taxonomic differences between young, middle, and old mouse groups
To determine if there are differences in the microbiomes of mice of different ages, we compared the taxonomic composition of the microbiome samples using several different methods. Initially, we extracted 16S rRNA ribosomal genes from the metagenomic sequences and assigned these to 97% operational taxonomic units (OTUs) from GreenGenes using the Quantitative Insights Into Microbial Ecology (QIIME) closed-reference OTU-picking protocol. Principal coordinate analysis of the samples showed significant separation by age groups using both weighted UniFrac (ANOSIM p value = 0.004) and unweighted UniFrac (ANOSIM p = 0.001) (Figure 2). Principal coordinate 1 (PC1) (percent variation explained: 15.2%) of the unweighted UniFrac separates all three sets of samples, while PC2 (10.8%) separates the old group from the other samples even further. In the weighted UniFrac plot, separation of the samples by age only appears after introducing PC2 (16.5%) and not with PC1 (54%), which appears to be driven by substantial changes in taxonomic relative abundance in two samples (Y7-Aug15 and 1E-May23). These two samples both have reduced levels of Bacteroidales relative to Clostridiales and vary greatly from all other samples including those taken from the same mouse within very close time spans (Y7-Aug15: 1 day and 1E-May23: 32 days). In spite of these differences, the clustering of samples by age group even for these two outliers by unweighted UniFrac suggests a strong relationship in qualitative terms. The clearer delineation of samples by age in the unweighted UniFrac analysis suggests that the presence or absence of particular OTUs is more important than the abundance of these OTUs in separating the samples by age. To further test if separation of samples was observed regardless of method, a PCA plot was created that used only the raw abundances of each OTU and did not use a beta diversity measure to relate the OTUs to each other. These results were consistent with our initial PCoA analysis as illustrated by the clear separation of the samples by age group when PC1 (40.0%) and PC2 (23.4%) were plotted (Additional file 3).
The second approach for determining taxonomic composition within the metagenomic samples used 40 conserved protein-coding genes that have previously been described  and were assigned using the PhyloSift package . Metagenomic reads matching to each of the conserved protein-coding genes were then inserted into their corresponding reference phylogenetic gene tree using Pplacer  and then compared using Pplacer’s edge principal component analysis (edge PCA) . The edge PCA method transforms the placed reads such that each edge in the tree becomes a variable of interest, which is then used for PCA. This has the added advantage of relating each of the principal components directly to the phylogenetic edges that are contributing to that signal. Similar to the weighted UniFrac PCoA, the first principal component (87.4%) from the edge PCA did not provide any separation of the samples, and after visualizing the drivers of PC1 on the phylogenetic tree provided by Pplacer, it seems to be driven by a trade-off between Firmicutes/Clostridia and Bacteroidetes . However, separation of samples by age did appear when using PC2 (6%) and PC3 (4%) (Additional file 4A). Separation between the old and middle samples occurred primarily by PC2 of the edge PCA plot and was driven by the middle group having more Rikenellaceae, Lactobacillaceae, Mycoplasmataceae, and Erysipelotrichaceae, and less Lachnospiraceae, Clostridiaceae, Ruminococcaceae, Prevotellaceae, and Porphyromonadaceae (Additional file 4B). PC3 primarily separated the young samples from the middle and old samples and was driven by the young having more Lactobacillaceae, Prevotellaceae, and Porphyromonadaceae and the middle and old groups having more Rikenellaceae, Lachnospiraceae, Ruminococcaceae, and Clostridiaceae (Additional file 4C). In all cases, the samples of young and old mice were more similar to each other than either was to the middle group. This observation agrees with studies in humans that also showed that middle-aged people had more distinct taxonomic  and functional  compositions then other adult age groups.
Some samples were taken from the same mouse over a period of time, raising the question of whether intra-mouse correlations artificially increased the apparent similarity within age groups. To address this question, we examined the average beta diversity between all samples taken from the same mouse in a given age group, versus all pairs of samples taken from different mice. We found no significant difference in means (Welch t test p < 0.05) between the two types of samples at all ages, for both weighted and unweighted UniFrac (Additional file 5). The average distances being compared never differed by more than 0.029 (middle mouse group; unweighted UniFrac), and in the young group, the samples from the same mice had greater beta diversity than those taken from different mice. This lack of difference suggests that there are substantial changes in the microbiome composition over short periods of time but that there is commonality in the microbiome as mice age under controlled laboratory conditions.
Taxonomic differences were compared across the different age groups to determine if particular taxa are associated with the gut microbiomes of aging and frail mice. The Rikenellaceae family, which contains the Alistipes genus and has previously been linked to microbiomes from elderly people , was the most significantly overrepresented family within the middle and old groups in comparison to young mice when using 16S data (Kruskal-Wallis H test, Benjamini-Hochberg FDR multiple test correction p = 0.007) (Figure 3; Additional files 6 and 7). To ensure that this taxonomic link between frailer mice and people was not an artifact of using the 16S rRNA gene as a marker, an independent analysis using protein markers (PhyloSift) and phylogenetic placement (Pplacer) was conducted. The latter method has the advantage of identifying the phylogenetic context of sequenced reads but does not provide statistical significance testing. However, the results from this method did provide additional support that the Alistipes were more abundant within the old and middle mouse groups (Figure 4). The PhyloSift approach indicated that the old mice have lower abundances of organisms from the Lachnospiraceae family, a group often associated with the production of beneficial SCFAs , even though the 16S-based method did not show a significant difference (p = 0.7) between the age groups. This contrasting result suggests that the type of phylogenetic marker used for taxonomic assessment can give conflicting results for some taxon families. Lastly, the PhyloSift approach indicated that some members of the Bacteroidaceae family were abundant in the old and middle groups, while fewer were abundant in the young group, which did agree with the 16S results (old: 9.2% ± 7.9%; middle: 10.3% ± 1.3%; young: 3.9% ± 1.4%, p = 0.06).
Previous research in humans showed a negative relationship between frailty and the abundance of Eubacteriaceae, Faecalibacterium, and Lactobacillus . We found virtually no Eubacteriaceae (eight sequences total) within any of our samples and no difference in proportions of Faecalibacterium between the age groups (p = 0.828; Additional file 6). These differing results could be due to unique differences between aging in mice and in humans or due to confounding signals from diet or habitat change in the elderly study. Although a recent finding did identify in aging mice the same decrease in Lactobacillus as in the elderly , no significant decrease in Lactobacillus was observed (old: 2.4% ± 2.4%; middle: 7.0% ± 4.6%; young: 4.9% ± 5.3%, p = 0.44; Additional file 6). Zhang et al. also reported a large phylum shift for all their aged mice from Firmicutes to Bacteroidetes, but we did not observe any change in the Firmicutes/Bacteroidetes ratio among our different age groups . Additionally, we noted that the genus Akkermansia was nearly absent from the old group (0.003% ± 0.004%), compared to the middle (0.39% ± 0.27%) and young mice (0.50% ± 0.8%), but the difference in means was again not significant at an alpha threshold of 0.05 (p = 0.058; Additional file 6). Akkermansia has been linked to healthy microbiomes in some studies ,, although members of this genus may also exacerbate certain types of infection via mucin degradation .
Functional differences among microbiomes
Annotation of metagenomic sequences using MG-RAST led to the identification of 99 SEED level 3 categories that had significant mean abundance differences across the young, middle, and old categories using STAMP  (Kruskal-Wallis p < 0.05 with Benjamini-Hochberg FDR correction). Detailed results are shown in Additional file 8; here, we focus on functions of particular interest due to their possible roles in aging and frailty.
Carbohydrate and lactate utilization
Most monosaccharide utilization categories were overrepresented in the old group compared to the young and middle groups, with the complex polysaccharide-xyloglucan utilization category underrepresented in the old group (0.112% ± 0.014%, p = 0.034) compared to the middle group (0.122% ± 0.016%) (Figure 5A). Monosaccharide utilization categories with statistically significant changes among age groups were L-fucose (old: 0.268% ± 0.015%; middle: 0.243% ± 0.015%; young: 0.172% ± 0.022%, p = 0.034), L-rhamnose (old: 0.347% ± 0.059%; middle: 0.146% ± 0.014%; young: 0.322% ± 0.064%, p = 0.040), D-galacturonate and D-glucuronate (old: 0.579% ± 0.032%; middle: 0.431% ± 0.052%; young: 0.492% ± 0.025%, p = 0.036), and xylose (old: 0.333% ± 0.040%; middle: 0.182% ± 0.040%; young: 0.317% ± 0.055%, p = 0.042). This observation agrees with earlier reports of decline in the saccharolytic potential of microbiome with age and further suggests a limited supply of simple sugars for gut epithelium function in the old age group.
Additionally, the D-galacturonate and D-glucuronate utilization category includes glucuronidases that catalyze the hydrolysis of D-glucuronic acid as a part of complex carbohydrate metabolism . Elevated levels of human β-glucuronidase activity in tumors compared to normal cells led to the development of several selective cancer chemotherapeutics . However, many microbial species contain β-glucuronidase homologs, which have been shown to reactivate the cancer drug irinotecan (e.g., CPT-11) leading to severe diarrhea and thus limiting the efficacy of the drug ,. Fructooligosaccharides (FOS) and raffinose are carbohydrates that are commonly used as prebiotics to promote growth of Lactobacilli and Bifidobacteria in the GI tract of older individuals ,. The underrepresentation of FOS and raffinose utilization in old mice (0.142% ± 0.037%, p = 0.031) compared to mice in the young (0.200% ± 0.044%) and middle (0.314% ± 0.048%) groups could be another contributing factor to perturbing beneficial microbe populations and an increased risk of opportunistic infection.
The concentration of lactate in the colon is typically maintained by lactate-utilizing bacteria and metabolism of D-lactate by host lactases . The reduced colonic ability to utilize lactate in the old group (0.018% ± 0.011%, p = 0.041) compared to the middle group (0.061% ± 0.014%) and the negative effect of aging on host lactase activity  can adversely affect host health with accumulation of lactate in the colon. High fecal lactate is associated with ulcerative colitis and other inflammatory bowel diseases ,.
Biosynthesis of vitamins and creatine degradation
Cobalamin (B12) biosynthesis was significantly underrepresented within the old group (0.037% ± 0.041%, p = 0.041) compared to the middle (0.129% ± 0.026%) and young (0.054% ± 0.018%) groups (Figure 5B). In a previous study, 55% of people aged 43 and under had a microbiome encoding a gene, cbiN (cobalamin biosynthesis protein; pfam ID PF02553), involved in B12 biosynthesis, while only 11% those over the age of 43 had a copy of it . Deficiency in B12 within elderly people has been well documented in the past ,, and the decreased potential of the microbiome to synthesize vitamin B12 has been suggested as a contributing factor . Conversely, creatine and creatinine degradation functionality was overrepresented in the gut microbiomes of the old group (0.0076% ± 0.002%, p = 0.030) compared to the middle (0.0031% ± 0.001%) and young (0.0018% ± 0.001%) groups. Creatine supplementation has been shown to promote muscle strength and hypertrophy in elderly people . Therefore, an increase in degradation of creatine by the microbiome may have negative consequences for the host if creatine abundance was limited and thus indirectly affected the FI by leading to reduced lean mass.
Another vitamin, biotin (vitamin B7), is synthesized in considerable amounts by microbiota of the large intestine (sometimes greater than the amount taken in the diet), and a significant portion of this is absorbed by human and animal colonocytes . Biotin synthesis was significantly underrepresented within the young (0.011% ± 0.022%, p = 0.040) and old (0.058% ± 0.103%) groups compared to the middle group (0.158% ± 0.044%). Younger individuals require less biotin compared to those of middle and old age, so low contributions from gut microbiota in the young group may not lead to biotin deficiencies under a normal diet. In the old group, however, low biotin can lead to biotin deficiency which in turn increases colon cancer risk .
DNA repair and persister cells
Bacterial SOS genes (including uvrABC and uvrD) play an essential function in bacterial nucleotide excision repair and decreasing the rate of mutagenesis . Although no observations have been made until now on the status of DNA repair in any microbiome studies, the significant underrepresentation of uvrD (0.0006% ± 0.0008%, p = 0.040) and uvrABC (0.57% ± 0.10%, p = 0.039) within the old group compared to the young group (uvrD: 0.061% ± 0.029%; uvrABC: 0.675% ± 0.066%) suggests potentially higher rates of mutagenesis and the possibility of acquiring antibiotic resistance within the gut microbiota of the old group. The rise in functional groups associated with persister cells from young (undetectable, p = 0.034) and middle (0.005% ± 0.002%) groups to the old group (0.007% ± 0.004%) points to the possibility of multidrug-resistant bacterial populations which may contribute to inflammation in old age .
Amino acid degradation
Rampelli et al.  suggested an association of aromatic amino acid degradation (except histidine) with aging. In the present study, there was no change in tryptophan and histidine metabolism among age groups. However, we found significant underrepresentation of other aromatic amino acid degradation within the old group (0.024% ± 0.003%, p = 0.049) compared to young (0.028% ± 0.004%) and middle (0.039% ± 0.006%) groups, although the effect size was very small.
Taxonomic assignments of identified functions
Using the RITA software package , we performed taxonomic assignment of sequence reads to determine if the shift in abundance was accompanied by a shift of a specific set of microbes. Many key functions showed a shift from genera Akkermansia and Parabacteroides to Bacteroides and Firmicutes. For example, 49.5% of cobalamin synthesis genes in the young and middle mouse groups were assigned to the genus Bacteroides, versus 71.1% in the old mice. Conversely, the contribution of Parabacteroides to this function dropped from 24.8% to 7.0%. Akkermansia was the predicted source of 74.4% of UvrD and related functional genes in the young and middle mouse groups, versus 13.3% in the old mouse cohort, which is instead dominated by Bacteroides and to a lesser extent Oscillibacter and Clostridium. Although the assignment of functions to taxonomic groups can be error-prone, the difference in these assignments between samples from mice of different ages indicates that differences exist, even if some assignments may be suspect.
Although previous microbiome studies in humans have shown changes with age, the results were confounded by other possible impacts on the microbiome such as diet, living conditions, and medications that also change with age. Here, we have used a mouse model that allows us to control for these factors, and we have observed taxonomic and functional changes in the microbiome that correlate with age and frailty. Consistent with an earlier study of frail humans, the genus Alistipes was found to be overrepresented in old mice, suggesting there may be some parallel shifts that occur in aging human and mouse populations. Further, we identify several additional bacterial taxa and functions that may be associated with the aging process. Our results suggest that the aging microbiome could have an effect on the availability of vitamins (B12 and B7) and creatine, DNA repair, and carbohydrate metabolism as well as the potential to interfere with some drug treatments.
Further work with larger mouse cohorts, particularly older mice with a range of frailty scores, will be needed to separate the effects of frailty and aging on the microbiome. Larger longitudinal studies would also clarify the relationship between specific clinical attributes and changes to the microbiome during aging. For example, previous work showed a steady increase in mouse weight up to a peak at 15–20 months, followed by a steady decline as the mice age . The microbiomes of obese individuals tend to have different taxonomic and metabolic properties ; there may be corresponding weight-associated microbiome shifts as mice age. Another important factor within this study is the possible role of cage effects on the microbiome. Hildebrand et al. showed that up to 30% of the variation in the microbiome can arise from cage effects, possibly due to coprophagy , and the authors recommended that groups of interest be housed in the same cage, or individually. However, neither of these approaches is appropriate since solitary cages for the entire life of a mouse can induce high levels of stress  which may in turn influence the microbiome, while a single cage for all mice is not practical for large-scale animal studies. Explicit tests of housing strategies are needed to examine the effects of this important factor. Although larger studies will be necessary to reveal the fine details of aging, frailty, and the microbiome, our results show that the mouse microbiome changes throughout various life stages and suggests that similar changes in humans may have a significant effect on health.
Subject mice and frailty assessment
Female C57BL/6 J mice purchased from Charles River (St. Constant, Quebec) at 3–4 weeks of age were housed in micro-isolator cages on a 12-h light/dark cycle in the Carlton Animal Care Facility at Dalhousie University. All mice in this study had free access to the same food (Prolab RMH 3000, LabDiet, St. Louis, MO) and water in their home cages. Three ages of mice were used in the experiments: a young adult group (174 ± 15 days old), a middle adult group (589 ± 18 days old), and old adult group (857 ± 16 days old). Study animals from the same age group were all housed in the same cage, with the exception of one old mouse which was housed together with another old mouse from which no data were obtained. Experiments were performed following the guidelines outlined by the Canadian Council on Animal Care Guide to the Care and Use of Experimental Animals (CCAC, Ottawa, ON: Vol. 1, 2nd edition, 1993; Vol. 2, 1984). All animal protocols were approved by the Dalhousie University Committee on Laboratory Animals.
Clinical examinations were performed between 10 am and 1 pm every day. Animals were first observed in their home cage and then taken to an assessment room. Mice were weighed, and body surface temperature was measured with an infrared temperature probe directed at the abdomen (Infrascan; La Crosse Technology). A brief clinical exam was then used to create a FI score for each mouse as described above.
Fecal sample collection and sequencing
Each mouse was placed in an empty, sterilized 43-cm-long, 21-cm-wide cage for 1.5 to 2 h to collect individual fecal samples in a manner similar to Lomansey et al. . A duration of 1.5 to 2 h of isolation may induce some stress in the mice, although much longer times seem to be required to induce major changes in behavior consistent with stress . Samples were collected between 10 am and 1 pm. To minimize sample contamination, mice had no access to food and water during this time and sterile forceps were used for sample handling. On average, mice produced between 100 and 200 mg of sample during each collection period. Immediately after collection, samples were frozen at −80°C.
Genomic DNA (gDNA) was isolated from fecal samples using the PowerSoil DNA isolation kit (MO BIO Laboratories, Inc.) according to instructions from the manufacturer. Once the isolation procedure was completed, the concentration of gDNA was determined for each sample by spectrophotometry and the sample was then frozen at −80°C. Samples were then shipped to Génome Québec (McGill University and Génome Québec Innovation Centre, Montreal, Québec, Canada) where they underwent quality control assessment to confirm sample integrity. Metagenomic sequencing of the microbial DNA extracted from the 21 stool samples was conducted on an Illumina HiSeq generating 100-bp paired-end reads. This produced a total of 85.61 Gbps (mean 2.04 ± 0.81 Gbps per sample) and 428.06 million reads (mean 9.96 ± 0.37 million reads per sample).
Raw FASTQ files for each paired end of the 21 samples were processed with MG-RAST  and made publicly available (project id: 3907). MG-RAST default pipeline options were used, including dereplication, screening against M. musculus NCBI v37 to remove host DNA sequences, and dynamic read trimming for quality control of sequence reads.
Taxonomic and functional analysis
In this study, three different methods were used to characterize the taxonomic diversity of samples: (i) using only the 16S rRNA gene (16S), (ii) using 40 different conserved protein genes , and (iii) assigning taxonomic classifications to metagenomic reads. Each of these methods is described in more detail below.
Files containing ribosomal 16S rRNA gene (16S) fragments as annotated by MG-RAST were downloaded for each paired end (e.g., ‘425.search.rna.fna’) and combined into their respective sample for further processing with QIIME . The putative 16S reads were assigned to existing OTUs using a closed picking protocol that considers only OTUs already present in a reference database. We used the QIIME-compatible UCLUST version 1.2.22q  to map the reads to version 13_5 reference package of GreenGenes  using a 97% identity threshold. Every sample was rarefied to 15,000 OTUs to remove possible bias from the variation in sequencing depth for different samples. To visualize and compare the differences between samples, principal coordinate analysis (PCoA) plots were generated with weighted and unweighted UniFrac beta diversity metrics . In addition, a PCA plot based only on the relative OTU abundances in each sample was generated using the R function ‘prcomp’. Taxonomic summary bar charts were created for 16S data by collapsing to different taxonomic ranks (level 5) using QIIME’s ‘summarize_taxa.py’ script, taking only the top ten most abundant families across all samples, renormalizing by sum samples, and then plotting using QIIME’s ‘plot_taxa_summary.py’ Python script.
PhyloSift was used as an alternative method to 16S for determining taxonomic composition . PhyloSift uses 40 conserved protein-coding genes that have previously been used for phylogenomic analysis  and places metagenomic read fragments onto a reference tree of sequenced genomes using Pplacer . Metagenomic reads passing MG-RAST quality control were downloaded (‘299.screen.passed.fna’), and paired ends were used as input to PhyloSift (git branch development version: 9a03023cb) with default options. The Pplacer ‘.jplace’ files representing the sequences placed on a concatenated 40-gene tree output by PhyloSift for each sample were used to create an ‘edge PCA’  and visualized in R. The individual sample ‘.jplace’ files were then merged into different age groupings using the ‘guppy merge’ command. Comparison of age groupings was conducted using the pairwise ‘guppy kr_heat’ program which allows visualization of the differences in organism abundance between two groupings on a phylogenetic tree.
Functional SEED  annotations including abundance information and actual annotated sequences were downloaded and formatted from the MG-RAST API using in-house scripts that are available at Github (https://github.com/mlangill/get_mgrast_data). Sequences with functions that statistically varied over the age groupings were filtered; a subset of these that yielded reliable sequence counts were subjected to taxonomic analysis. To determine and compare the taxonomic characterization of these particular functions, metagenomic reads annotated with these functions were analyzed using RITA . Briefly, RITA uses a reference database of sequenced and draft genomes to list the annotated proteins for each genome and construct compositional models of each genome. RITA uses a combination of homology search and compositional analysis to assign sequences to different classification groups. We used RITA v1.0.1 with a database of 2,987 finished and draft genomes, with USEARCH version 4.1.93  used for the homology search and FCP version 1.0.3  for compositional matching. Classes of RITA results considered for summary purposes were as follows: (i) cases where the top compositional and homology matches agreed with one another at the genus level and (ii) cases where the expectation value of the best-matching genome was at least 10 orders of magnitude better than the best-matching genome from another genus. Summaries were constructed for a small set of taxa that showed variation between a group comprising young and middle mice and old mice: Parabacteroides, Bacteroides, and Lactobacillus; all remaining genera in phyla Bacteroidetes and Firmicutes; and phylum Verrucomicrobia.
The statistical significance of samples between age groupings was calculated using ANOSIM as implemented in the R package vegan version 2.0-10 , with 999 permutations and the weighted and unweighted UniFrac distance matrices. All multiple-group comparisons were done using the Kruskal-Wallis H test with Benjamini-Hochberg FDR multiple test correction as calculated within STAMP version 2.0.7 . Where appropriate, reported p values are those corrected for multiple testing. Abundance values reported in the results and in the additional files are reported as mean relative abundances for the age groups along with the standard deviation.
Availability of supporting data
The data set supporting the results of this article is available in the MG-RAST repository, as project 3907, http://metagenomics.anl.gov/linkin.cgi?project=3907.
Principal component analysis
Principal coordinate analysis
Short-chain fatty acid
Cho I, Blaser MJ: The human microbiome: at the interface of health and disease. Nat Rev Genet. 2012, 13: 260-270.
Marchesi JR, Dutilh BE, Hall N, Peters WHM, Roelofs R, Boleij A, Tjalsma H: Towards the human colorectal cancer microbiome. PLoS One. 2011, 6: e20447-10.1371/journal.pone.0020447.
Fahlén A, Engstrand L, Baker BS, Powles A, Fry L: Comparison of bacterial microbiota in skin biopsies from normal and psoriatic skin. Arch Dermatol Res. 2012, 304: 15-22. 10.1007/s00403-011-1189-x.
Giongo A, Gano KA, Crabb DB, Mukherjee N, Novelo LL, Casella G, Drew JC, Ilonen J, Knip M, Hyöty H, Veijola R, Simell T, Simell O, Neu J, Wasserfall CH, Schatz D, Atkinson MA, Triplett EW: Toward defining the autoimmune microbiome for type 1 diabetes. ISME J. 2011, 5: 82-91. 10.1038/ismej.2010.92.
Major G, Spiller R: Irritable bowel syndrome, inflammatory bowel disease and the microbiome. Curr Opin Endocrinol Diabetes Obes. 2014, 21: 15-21. 10.1097/MED.0000000000000032.
Greenblum S, Turnbaugh PJ, Borenstein E: Metagenomic systems biology of the human gut microbiome reveals topological shifts associated with obesity and inflammatory bowel disease. Proc Natl Acad Sci U S A. 2012, 109: 594-599. 10.1073/pnas.1116053109.
Chang JY, Antonopoulos DA, Kalra A, Tonelli A, Khalife WT, Schmidt TM, Young VB: Decreased diversity of the fecal microbiome in recurrent Clostridium difficile-associated diarrhea. J Infect Dis. 2008, 197: 435-438. 10.1086/525047.
Petrof EO, Gloor GB, Vanner SJ, Weese SJ, Carter D, Daigneault MC, Brown EM, Schroeter K, Allen-Vercoe E: Stool substitute transplant therapy for the eradication of Clostridium difficile infection: ‘RePOOPulating’ the gut. Microbiome. 2013, 1: 3-10.1186/2049-2618-1-3.
Louis P, Young P, Holtrop G, Flint HJ: Diversity of human colonic butyrate-producing bacteria revealed by analysis of the butyryl-CoA:acetate CoA-transferase gene. Environ Microbiol. 2010, 12: 304-314. 10.1111/j.1462-2920.2009.02066.x.
Gill SR, Pop M, Deboy RT, Eckburg PB, Turnbaugh PJ, Samuel BS, Gordon JI, Relman DA, Fraser-Liggett CM, Nelson KE: Metagenomic analysis of the human distal gut microbiome. Science. 2006, 312: 1355-1359. 10.1126/science.1124234.
Ichinohe T, Pang IK, Kumamoto Y, Peaper DR, Ho JH, Murray TS, Iwasaki A: Microbiota regulates immune defense against respiratory tract influenza A virus infection. Proc Natl Acad Sci U S A. 2011, 108: 5354-5359. 10.1073/pnas.1019378108.
Wu H-J, Ivanov II, Darce J, Hattori K, Shima T, Umesaki Y, Littman DR, Benoist C, Mathis D: Gut-residing segmented filamentous bacteria drive autoimmune arthritis via T helper 17 cells. Immunity. 2010, 32: 815-827. 10.1016/j.immuni.2010.06.001.
Walls T, Power D, Tagg J: Bacteriocin-like inhibitory substance (BLIS) production by the normal flora of the nasopharynx: potential to protect against otitis media?. J Med Microbiol. 2003, 52 (Pt 9): 829-833. 10.1099/jmm.0.05259-0.
Liévin V, Peiffer I, Hudault S, Rochat F, Brassart D, Neeser JR, Servin AL: Bifidobacterium strains from resident infant human gastrointestinal microflora exert antimicrobial activity. Gut. 2000, 47: 646-652. 10.1136/gut.47.5.646.
Yatsunenko T, Rey FE, Manary MJ, Trehan I, Dominguez-Bello MG, Contreras M, Magris M, Hidalgo G, Baldassano RN, Anokhin AP, Heath AC, Warner B, Reeder J, Kuczynski J, Caporaso JG, Lozupone CA, Lauber C, Clemente JC, Knights D, Knight R, Gordon JI: Human gut microbiome viewed across age and geography. Nature. 2012, 486: 222-227.
Lozupone CA, Stombaugh J, Gonzalez A, Ackermann G, Wendel D, Vázquez-Baeza Y, Jansson JK, Gordon JI, Knight R: Meta-analyses of studies of the human microbiota. Genome Res. 2013, 23: 1704-1714. 10.1101/gr.151803.112.
Koenig JE, Spor A, Scalfone N, Fricker AD, Stombaugh J, Knight R, Angenent LT, Ley RE: Succession of microbial consortia in the developing infant gut microbiome. Proc Natl Acad Sci U S A. 2011, 108 (Suppl): 4578-4585. 10.1073/pnas.1000081107.
Kurokawa K, Itoh T, Kuwahara T, Oshima K, Toh H, Toyoda A, Takami H, Morita H, Sharma VK, Srivastava TP, Taylor TD, Noguchi H, Mori H, Ogura Y, Ehrlich DS, Itoh K, Takagi T, Sakaki Y, Hayashi T, Hattori M: Comparative metagenomics revealed commonly enriched gene sets in human gut microbiomes. DNA Res. 2007, 14: 169-181. 10.1093/dnares/dsm018.
Mariat D, Firmesse O, Levenez F, Guimarăes V, Sokol H, Doré J, Corthier G, Furet J-P: The Firmicutes/Bacteroidetes ratio of the human microbiota changes with age. BMC Microbiol. 2009, 9: 123-10.1186/1471-2180-9-123.
Biagi E, Nylund L, Candela M, Ostan R, Bucci L, Pini E, Nikkïla J, Monti D, Satokari R, Franceschi C, Brigidi P, De Vos W: Through ageing, and beyond: gut microbiota and inflammatory status in seniors and centenarians. PLoS One. 2010, 5: e10667-10.1371/journal.pone.0010667.
Van Tongeren SP, Slaets JPJ, Harmsen HJM, Welling GW: Fecal microbiota composition and frailty. Appl Environ Microbiol. 2005, 71: 6438-6442. 10.1128/AEM.71.10.6438-6442.2005.
Claesson MJ, Jeffery IB, Conde S, Power SE, O’Connor EM, Cusack S, Harris HMB, Coakley M, Lakshminarayanan B, O’Sullivan O, Fitzgerald GF, Deane J, O’Connor M, Harnedy N, O’Connor K, O’Mahony D, van Sinderen D, Wallace M, Brennan L, Stanton C, Marchesi JR, Fitzgerald AP, Shanahan F, Hill C, Ross RP, O’Toole PW: Gut microbiota composition correlates with diet and health in the elderly. Nature. 2012, 488: 178-184. 10.1038/nature11319.
Rampelli S, Candela M, Turroni S, Biagi E, Collino S, Franceschi C, O’Toole PW, Brigidi P: Functional metagenomic profiling of intestinal microbiome in extreme ageing. Aging (Albany NY). 2013, 5: 902-912.
Lutz CT, Quinn LS: Sarcopenia, obesity, and natural killer cell immune senescence in aging: altered cytokine levels as a common mechanism. Aging (Albany NY). 2012, 4: 535-546.
Fritz JV, Desai MS, Shah P, Schneider JG, Wilmes P: From meta-omics to causality: experimental models for human microbiome research. Microbiome. 2013, 1: 14-10.1186/2049-2618-1-14.
Kostic AD, Howitt MR, Garrett WS: Exploring host-microbiota interactions in animal models and humans. Genes Dev. 2013, 27: 701-718. 10.1101/gad.212522.112.
Parks RJ, Fares E, Macdonald JK, Ernst MC, Sinal CJ, Rockwood K, Howlett SE: A procedure for creating a frailty index based on deficit accumulation in aging mice. J Gerontol A Biol Sci Med Sci. 2012, 67: 217-227. 10.1093/gerona/glr193.
Rockwood K, Mitnitski A: Frailty in relation to the accumulation of deficits. J Gerontol A Biol Sci Med Sci. 2007, 62: 722-727. 10.1093/gerona/62.7.722.
Tachon S, Zhou J, Keenan M, Martin R, Marco ML: The intestinal microbiota in aged mice is modulated by dietary resistant starch and correlated with improvements in host responses. FEMS Microbiol Ecol. 2013, 83: 299-309. 10.1111/j.1574-6941.2012.01475.x.
Fares E, Howlett SE: Effect of age on cardiac excitation-contraction coupling. Clin Exp Pharmacol Physiol. 2010, 37: 1-7. 10.1111/j.1440-1681.2009.05276.x.
Wu M, Eisen JA: A simple, fast, and accurate method of phylogenomic inference. Genome Biol. 2008, 9: R151-10.1186/gb-2008-9-10-r151.
Darling AE, Jospin G, Lowe E, Matsen FA, Bik HM, Eisen JA: PhyloSift: phylogenetic analysis of genomes and metagenomes. PeerJ. 2014, 2: e243-10.7717/peerj.243.
Matsen FA, Kodner RB, Armbrust EV: pplacer: linear time maximum-likelihood and Bayesian phylogenetic placement of sequences onto a fixed reference tree. BMC Bioinformatics. 2010, 11: 538-10.1186/1471-2105-11-538.
Matsen FA, Evans SN: Edge principal components and squash clustering: using the special structure of phylogenetic placement data for sample comparison. PLoS One. 2013, 8: e56859-10.1371/journal.pone.0056859.
Huttenhower C, Gevers D, Huttenhower C, Gevers D, Knight R, Abubucker S, Badger JH, Chinwalla AT, Creasy HH, Earl AM, FitzGerald MG, Fulton RS, Giglio MG, Hallsworth-Pepin K, Lobos EA, Madupu R, Magrini V, Martin JC, Mitreva M, Muzny DM, Sodergren EJ, Versalovic J, Wollam AM, Worley KC, Wortman JR, Young SK, Zeng Q, Aagaard KM, Abolude OO, Allen-Vercoe E, et al: Structure, function and diversity of the healthy human microbiome. Nature. 2012, 486: 207-214. 10.1038/nature11234.
Lan Y, Kriete A, Rosen GL: Selecting age-related functional characteristics in the human gut microbiome. Microbiome. 2013, 1: 2-10.1186/2049-2618-1-2.
Meehan CJ, Beiko RG: A phylogenomic view of ecological specialization in the Lachnospiraceae, a family of digestive tract-associated bacteria. Genome Biol Evol. 2014, 6: 703-713. 10.1093/gbe/evu050.
Zhang C, Li S, Yang L, Huang P, Li W, Wang S, Zhao G, Zhang M, Pang X, Yan Z, Liu Y, Zhao L: Structural modulation of gut microbiota in life-long calorie-restricted mice. Nat Commun. 2013, 4: 2163-
Shin N-R, Lee J-C, Lee H-Y, Kim M-S, Whon TW, Lee M-S, Bae J-W: An increase in the Akkermansia spp. population induced by metformin treatment improves glucose homeostasis in diet-induced obese mice. Gut. 2014, 63: 727-735. 10.1136/gutjnl-2012-303839.
Kang C-S, Ban M, Choi E-J, Moon H-G, Jeon J-S, Kim D-K, Park S-K, Jeon SG, Roh T-Y, Myung S-J, Gho YS, Kim JG, Kim Y-K: Extracellular vesicles derived from gut microbiota, especially Akkermansia muciniphila, protect the progression of dextran sulfate sodium-induced colitis. PLoS One. 2013, 8: e76520-10.1371/journal.pone.0076520.
Ganesh BP, Klopfleisch R, Loh G, Blaut M: Commensal Akkermansia muciniphila exacerbates gut inflammation in Salmonella Typhimurium-infected gnotobiotic mice. PLoS One. 2013, 8: e74963-10.1371/journal.pone.0074963.
Parks DH, Beiko RG: Identifying biologically relevant differences between metagenomic communities. Bioinformatics. 2010, 26: 715-721. 10.1093/bioinformatics/btq041.
Vasella A, Davies GJ, Böhm M: Glycosidase mechanisms. Curr Opin Chem Biol. 2002, 6: 619-629. 10.1016/S1367-5931(02)00380-0.
Tranoy-Opalinski I, Legigan T, Barat R, Clarhaut J, Thomas M, Renoux B, Papot S: β-Glucuronidase-responsive prodrugs for selective cancer chemotherapy: an update. Eur J Med Chem. 2014, 74: 302-313. 10.1016/j.ejmech.2013.12.045.
Wallace BD, Wang H, Lane KT, Scott JE, Orans J, Koo JS, Venkatesh M, Jobin C, Yeh L-A, Mani S, Redinbo MR: Alleviating cancer drug toxicity by inhibiting a bacterial enzyme. Science. 2010, 330: 831-835. 10.1126/science.1191175.
Wilson ID, Nicholson JK: The role of gut microbiota in drug response. Curr Pharm Des. 2009, 15: 1519-1523. 10.2174/138161209788168173.
Dinoto A, Marques TM, Sakamoto K, Fukiya S, Watanabe J, Ito S, Yokota A: Population dynamics of Bifidobacterium species in human feces during raffinose administration monitored by fluorescence in situ hybridization-flow cytometry. Appl Environ Microbiol. 2006, 72: 7739-7747. 10.1128/AEM.01777-06.
Ewaschuk JB, Naylor JM, Zello GA: D-lactate in human and ruminant metabolism. J Nutr. 2005, 135: 1619-1625.
Lee MF, Russell RM, Montgomery RK, Krasinski SD: Total intestinal lactase and sucrase activities are reduced in aged rats. J Nutr. 1997, 127: 1382-1387.
Vernia P, Caprilli R, Latella G, Barbetti F, Magliocca FM, Cittadini M: Fecal lactate and ulcerative colitis. Gastroenterology. 1988, 95: 1564-1568.
Hove H, Nordgaard-Andersen I, Mortensen PB: Faecal DL-lactate concentration in 100 gastrointestinal patients. Scand J Gastroenterol. 1994, 29: 255-259. 10.3109/00365529409090473.
Markle HV: Cobalamin. Crit Rev Clin Lab Sci. 1996, 33: 247-356. 10.3109/10408369609081009.
Baik HW, Russell RM: Vitamin B12 deficiency in the elderly. Annu Rev Nutr. 1999, 19: 357-377. 10.1146/annurev.nutr.19.1.357.
Magrone T, Jirillo E: The interaction between gut microbiota and age-related changes in immune function and inflammation. Immun Ageing. 2013, 10: 31-10.1186/1742-4933-10-31.
Dalbo VJ, Roberts MD, Lockwood CM, Tucker PS, Kreider RB, Kerksick CM: The effects of age on skeletal muscle and the phosphocreatine energy system: can creatine supplementation help older adults. Dyn Med. 2009, 8: 6-10.1186/1476-5918-8-6.
Said HM, Ortiz A, McCloud E, Dyer D, Moyer MP, Rubin S: Biotin uptake by human colonic epithelial NCM460 cells: a carrier-mediated process shared with pantothenic acid. Am J Physiol. 1998, 275 (5 Pt 1): C1365-C1371.
O’Keefe SJD, Ou J, Aufreiter S, O’Connor D, Sharma S, Sepulveda J, Fukuwatari T, Shibata K, Mawhinney T: Products of the colonic microbiota mediate the effects of diet on colon cancer risk. J Nutr. 2009, 139: 2044-2048. 10.3945/jn.109.104380.
Reardon JT, Sancar A: Nucleotide excision repair. Prog Nucleic Acid Res Mol Biol. 2005, 79: 183-235. 10.1016/S0079-6603(04)79004-2.
MacDonald NJ, Parks DH, Beiko RG: Rapid identification of high-confidence taxonomic assignments for metagenomic data. Nucleic Acids Res. 2012, 40: e111-10.1093/nar/gks335.
Fahlström A, Yu Q, Ulfhake B: Behavioral changes in aging female C57BL/6 mice. Neurobiol Aging. 2011, 32: 1868-1880. 10.1016/j.neurobiolaging.2009.11.003.
Turnbaugh PJ, Ley RE, Mahowald MA, Magrini V, Mardis ER, Gordon JI: An obesity-associated gut microbiome with increased capacity for energy harvest. Nature. 2006, 444: 1027-1031. 10.1038/nature05414.
Hildebrand F, Nguyen TLA, Brinkman B, Yunta RG, Cauwe B, Vandenabeele P, Liston A, Raes J: Inflammation-associated enterotypes, host genotype, cage and inter-individual effects drive gut microbiota variation in common laboratory mice. Genome Biol. 2013, 14: R4-10.1186/gb-2013-14-1-r4.
Bartolomucci A, Palanza P, Sacerdote P, Ceresini G, Chirieleison A, Panerai AE, Parmigiani S: Individual housing induces altered immuno-endocrine responses to psychological stress in male mice. Psychoneuroendocrinology. 2003, 28: 540-558. 10.1016/S0306-4530(02)00039-2.
Lomasney KW, Cryan JF, Hyland NP: Converging effects of a Bifidobacterium and Lactobacillus probiotic strain on mouse intestinal physiology. Am J Physiol Gastrointest Liver Physiol. 2014, 307: G241-G247. 10.1152/ajpgi.00401.2013.
Martin AL, Brown RE: The lonely mouse: verification of a separation-induced model of depression in female mice. Behav Brain Res. 2010, 207: 196-207. 10.1016/j.bbr.2009.10.006.
Meyer F, Paarmann D, D’Souza M, Olson R, Glass EM, Kubal M, Paczian T, Rodriguez A, Stevens R, Wilke A, Wilkening J, Edwards RA: The metagenomics RAST server—a public resource for the automatic phylogenetic and functional analysis of metagenomes. BMC Bioinformatics. 2008, 9: 386-10.1186/1471-2105-9-386.
Caporaso JG, Kuczynski J, Stombaugh J, Bittinger K, Bushman FD, Costello EK, Fierer N, Peña AG, Goodrich JK, Gordon JI, Huttley GA, Kelley ST, Knights D, Koenig JE, Ley RE, Lozupone CA, McDonald D, Muegge BD, Pirrung M, Reeder J, Sevinsky JR, Turnbaugh PJ, Walters WA, Widmann J, Yatsunenko T, Zaneveld J, Knight R: QIIME allows analysis of high-throughput community sequencing data. Nat Methods. 2010, 7: 335-336. 10.1038/nmeth.f.303.
Edgar RC: Search and clustering orders of magnitude faster than BLAST. Bioinformatics. 2010, 26: 2460-2461. 10.1093/bioinformatics/btq461.
McDonald D, Price MN, Goodrich J, Nawrocki EP, DeSantis TZ, Probst A, Andersen GL, Knight R, Hugenholtz P: An improved Greengenes taxonomy with explicit ranks for ecological and evolutionary analyses of bacteria and archaea. ISME J. 2012, 6: 610-618. 10.1038/ismej.2011.139.
Lozupone C, Knight R: UniFrac: a new phylogenetic method for comparing microbial communities. Appl Environ Microbiol. 2005, 71: 8228-8235. 10.1128/AEM.71.12.8228-8235.2005.
Overbeek R, Olson R, Pusch GD, Olsen GJ, Davis JJ, Disz T, Edwards RA, Gerdes S, Parrello B, Shukla M, Vonstein V, Wattam AR, Xia F, Stevens R: The SEED and the Rapid Annotation of microbial genomes using Subsystems Technology (RAST). Nucleic Acids Res. 2014, 42 (Database issue): D206-D214. 10.1093/nar/gkt1226.
Parks DH, MacDonald NJ, Beiko RG: Classifying short genomic fragments from novel lineages using composition and homology. BMC Bioinformatics. 2011, 12: 328-10.1186/1471-2105-12-328.
Dixon P: VEGAN, a package of R functions for community ecology. J Veg Sci. 2003, 14: 927-930. 10.1111/j.1654-1103.2003.tb02228.x.
Whitehead JC, Hildebrand BA, Sun M, Rockwood MR, Rose RA, Rockwood K, Howlett SE: A clinical frailty index in aging mice: comparisons with frailty index data in humans. J Gerontol A Biol Sci Med Sci. 2014, 69: 621-632. 10.1093/gerona/glt136.
This work was supported by the Canadian Institutes of Health Research (CIHR) through an Emerging Team Grant to RGB, CIHR Operating Grants to SEH (MOP 126018) and RAR (MOP 93718), and a CIHR Fellowship to MGIL. Infrastructure was supported by the Canada Foundation for Innovation through a grant to RGB. RGB also acknowledges the support of the Canada Research Chairs program. The authors are grateful to Jocelyne Whitehead and Dr. Giban Ray for excellent technical assistance with the mouse frailty index and sample collection procedures.
The authors declare that they have no competing interests.
MGIL, RAR, SEH, and RGB designed the experiment; RAR and SEH performed the mouse frailty analysis, sample collection, and gDNA isolation; all authors analyzed the data; and all authors wrote the manuscript. All authors read and approved the final manuscript.