Journal of Nutrition and Metabolism

Journal of Nutrition and Metabolism / 2017 / Article

Research Article | Open Access

Volume 2017 |Article ID 4535710 |

Rondo P. Middleton, Sebastien Lacroix, Marie-Pier Scott-Boyer, Nikola Dordevic, Adam D. Kennedy, Amanda R. Slusky, Jerome Carayol, Christina Petzinger-Germain, Alison Beloshapka, Jim Kaput, "Metabolic Differences between Dogs of Different Body Sizes", Journal of Nutrition and Metabolism, vol. 2017, Article ID 4535710, 11 pages, 2017.

Metabolic Differences between Dogs of Different Body Sizes

Academic Editor: Phillip B. Hylemon
Received13 Jul 2017
Accepted25 Sep 2017
Published26 Oct 2017


Introduction. The domesticated dog, Canis lupus familiaris, has been selectively bred to produce extreme diversity in phenotype and genotype. Dogs have an immense diversity in weight and height. Specific differences in metabolism have not been characterized in small dogs as compared to larger dogs. Objectives. This study aims to identify metabolic, clinical, and microbiota differences between small and larger dogs. Methods. Gas chromatography/mass spectrometry, liquid chromatography/tandem mass spectrometry, clinical chemistry analysis, dual-energy X-ray absorptiometry, and 16S pyrosequencing were used to characterize blood metabolic, clinical, and fecal microbiome systems, respectively. Eighty-three canines from seven different breeds, fed the same kibble diet for 5 weeks, were used in the study. Results. 449 metabolites, 16 clinical parameters, and 6 bacteria (at the genus level) were significantly different between small and larger dogs. Hierarchical clustering of the metabolites yielded 8 modules associated with small dog size. Conclusion. Small dogs had a lower antioxidant status and differences in circulating amino acids. Some of the amino acid differences could be attributed to differences in microflora. Additionally, analysis of small dog metabolites and clinical parameters reflected a network which strongly associates with kidney function.

1. Introduction

Significant genetic and metabolic variation occurs within the Canis lupus species that spans from the wolf through all domesticated canines. Artificial selection for phenotypic traits generated profound genetic differences within Canis lupus familiaris, the domesticated dog that is widespread in all human cultures. Dogs can vary in size between roughly 2 and 90 Kg. Segmenting dogs by size reveals some inherent characteristics. Small dogs are more likely to suffer from integumentary, cardiovascular, and dental diseases and have a higher incidence of endocrine-related deaths compared to larger dogs [1, 2]. Small dogs have a lower basal metabolic rate and a higher mass-specific metabolic rate than larger dogs [3]. Interestingly, small dogs have a longer lifespan than larger dogs [4], which differs from what is observed in other mammalian species.

We conducted a diet- and environment-controlled study in canines to understand metabolic, clinical, and microbiota differences between small and larger dogs. Multiple differences were found, specifically in blood concentrations of antioxidants and amino acids, as well as in microbiota composition. To better understand these differences, we identified modules of highly cooccurring metabolites and further analyzed correlations between metabolic and clinical data.

2. Materials and Methods

2.1. Cohort Study Design

Eighty-three (83) canines from seven different breeds were all fed the same dry extruded kibble diet for 5 weeks. Small dogs (34 individuals; Beagle, Small Fox Terrier, and Miniature Schnauzer) had a mean weight of 9.3 kg (range of 6.1–15.6 kg). Larger dogs (49 individuals; Labrador Retriever, English Setter, Siberian Husky, and Rottweiler) had a mean weight of 31.5 kg (range of 18.4–54.4 kg). Small dogs had a mean age of 6.8 years (range of 2.4–13.3 years) while larger dogs had a mean age of 6.2 years (range of 2.2–9.8 years). All dogs except, Rottweilers, were housed in the same location. Samples were handled and processed in the same manner to avoid technical variability or bias. Plasma-EDTA and serum samples were taken after overnight fasting during the fifth week of feeding. Fecal samples were collected during the fifth week of feeding. Plasma, serum, and fecal samples not immediately analyzed were frozen at −80°C.

2.2. Canine Plasma Metabolites

Metabolite profiling was performed as described previously [5]. In summary, metabolites from each blood sample were extracted and analyzed by GC/MS (Thermo Fischer, DSQ mass spectrometer) and LC/MS (Thermo Fischer, LTQ mass spectrometers). We carried out chromatographic separation followed by full-scan MS to record and quantify all detectable ions in the samples [6]. All metabolites with known chemical structure were identified by matching the ions’ chromatographic retention index and MS fragmentation signatures with reference library entries created from authentic standard metabolites [7]. The reported masses within the supplemental table represent the respective masses of the ion features utilized to identify each molecule. These ion features could be the mass of the parent molecule, the mass of an ion feature/fragment of the molecule, or an adduct (e.g., sodium adduct) of the parental molecule or an adduct of an ion fragment. Additional library entries were added for ions that were not covered by the standards based on their unique ion signatures (chromatographic and mass spectral) so these ions could then be routinely detected and quantified. For quality control and run-day performance analysis, labeled internal standards were spiked into all samples at different stages of the data acquisition process. The median relative standard deviations were 7% and 12% for the internal standards and endogenous biochemicals, 589 detected in total, respectively.

We identified 589 metabolites in plasma prior to removal of some for low percentage presence across all breeds. Of these, 401 biochemicals matched a named structure in the reference library (named). The remaining 188 biochemicals represented distinct chemical entities that represent a single molecule of discrete molecular formula and structure but that do not currently match a named chemical in the reference library (unnamed).

2.3. Canine Clinical Parameters

Most clinical parameters were analyzed in all dogs. Some anthropometric analyses using dual-energy X-ray absorptiometry (DEXA) and thyroid hormones triiodothyronine (T3) and thyroxine (T4) were evaluated only in 34 small dogs and 35 larger dogs (not assessed in Rottweilers). Serum creatinine, creatine kinase, potassium, total bilirubin, total antioxidant status (TAS), aspartate transaminase (AST), gamma-glutamyltransferase (GGT), T3, and T4 were measured using the Cobas c311 or e411 clinical chemistry analyzer, according to manufacturer’s directions. Protein digestibility was determined based on amount of protein in food consumed, amount in feces, and corrected for microbial nitrogen [8]. DEXA was performed according to manufacturer’s directions.

2.4. Gut Microbiota Composition

The sequences for 16S amplicon PCR forward and reverse primers for the variable regions V4 to V6 (V456) were 5′AGGCCAGCAGCCGCGGTAA and 5′GCCRRCACGAGCTGACGAC, respectively. The pyrosequencing was performed using Roche 454 GS-FLX Pyrosequencer. 489,290 sequences were generated for 98 samples (83 samples were used in this analysis). Data quality control and sequence trimming were performed using QIIME’s python script “split libraries” [9] with default settings except the following parameters: (i) no barcode mismatches were allowed, (ii) maximum sequence length was set to 520 bp, and (iii) a sliding window of 50 nucleotides was used with average quality score ≥ 25. Pyrosequencing error was removed using flowgram clustering. The chimeric sequences were detected and removed using UCHIME [10]. A total of 265,401 high quality sequences were obtained with the average of 3,198 sequences per sample for the 83 samples used in this analysis. The cleaned sequences were then clustered into operational taxonomic units (OTU) using a closed reference-based OTU picking method with similarity threshold of 97% [10, 11], where the reference data file was obtained from the greengenes website (, August 2013 release) [11]. A consensus taxonomic assignment for each OTU was performed using the ribosomal database project naïve Bayesian classifier [12] at a minimum confidence interval of 80%. Welch’s -test was used to calculate values.

2.5. Imputation of Missing Value and Outlier Detection

Parameters (metabolites or clinical) missing more than 20% values across all breeds were discarded while missing values were imputed in cases where less than 20% values were absent. Imputed metabolites were randomly determined from a uniform distribution between 0 and the lowest measured value. Imputed clinical parameters were replaced by the breed-specific median value of corresponding clinical parameter.

Outlier detection was performed using robust principal component analysis (PCA) as described previously [13] and implemented in the rrcov package for R [14]. PCA was performed separately for metabolites and clinical variables. One small dog (Small Fox Terrier) was identified as an outlier and removed from subsequent metabolomics analyses. After preprocessing, 449 metabolites (of which 131 were unnamed chemicals) were included in the final metabolomics dataset. No clinical outliers were detected.

2.6. Evaluation of Variable Importance

The ability of metabolomics and clinical variables to distinguish between body sizes was evaluated using random forest (RF) with the rfPermute package for R [15]. This method not only is an appropriate approach for variable selection [16] but also provides a measure of variable importance. RF model is initially created to calculate variable importance and the response variable is then permuted 1000 times to estimate values for RF importance.

Mann–Whitney tests were also used to compare differences in distribution of measured variables between body size groups in addition to RF. Biomarker studies show that two types of methods may provide complementary results [17]. Metabolite ratio between small and larger dogs was determined by dividing the metabolite value by its median value across all samples, then dividing the small dog mean by the larger dog mean.

2.7. Detection of Metabolite Modules and Correlation with Clinical Parameters

Weighted gene coexpression network analysis (WGCNA) implemented in the WGCNA R package [18] was used for determining coexpression patterns between genes. We applied such methodology to metabolomics data to identify groups or modules of highly cooccurring metabolites. Modules were identified using unsupervised clustering Dynamic Branch Cut method with an optimal value of 4 for soft threshold and a minimal number of 10 metabolites per module. Note that modules were randomly assigned colors, while unassigned metabolites were grouped in the grey module. A variable representative of all metabolites within a module—equivalent to the first eigenvalue of the PCA—was obtained (eigengene) and Spearman rank pairwise correlations between each eigengene and clinical parameters were calculated with unadjusted and Bonferroni Hochberg [19] adjusted value significance thresholds of 0.05 and 0.1, respectively.

To assess the ability of the modules to distinguish between body sizes, a hierarchical tree was constructed with the metabolites from larger dogs and metabolites were colored according to the small dog’s module memberships. -scores were then calculated to evaluate how well a module identified with one population of reference (in this case small dogs) was preserved when the reference was modified (i.e., larger dogs) [20]. Greater -scores show higher preservation and thus lower ability to distinguish between body sizes. All analyses were conducted with R software version 3.0.1 [21].

Associations with kidney function were performed using Pathway Studio Mammalian, ChemEffect and DiseaseFX databases (Elsevier). Kidney function was represented by the following cell process and clinical parameter terms: kidney function, kidney elimination, kidney excretion, renal reabsorption, kidney tubule function, kidney filtration, kidney blood flow, kidney vascularization, renal vasodilation, renal acidification, renin-angiotensin system, kidney development, renal water reabsorption, renal tubular secretion, renal clearance, and glomerular filtration rate.

3. Results

3.1. Plasma Circulating Metabolites and Clinical Measures Are Different between Dogs of Different Body Size

449 metabolites were used for statistical and computational analysis (see Supplementary Table in the Supplementary Material available online at Of these, 131 represent unnamed metabolites. Levels of 66 metabolites, of which 14 were unnamed, were significantly different between small and larger dogs (adjusted value < 0.05; Table 1). 69 clinical parameters were measured in this study (Supplementary Table ) with 16 being significantly different between small and larger dogs (adjusted value < 0.05; Table 2). Many of these parameters were based on weight (overall weight, DEXA parameters, etc.), which were expected to be different since this study compared two groups of dogs based on weight.

Metabolite nameRatio (small/larger)ClassAdj. PRIMass

Creatine2.81Amino acid758132.1
Creatinine0.78Amino acid730114.1
5-Oxoproline0.78Amino acid744128.2
Hydroxyproline0.72Amino acid705132.1
Phenylalanine0.84Amino acid2056166.1
p-Cresol sulfate1.48Amino acid2896187.1
Phenol sulfate1.74Amino acid2150173.1
C-Glycosyltryptophan0.79Amino acid1912367.1
Glutamine0.89Amino acid0.002684147.2
Glycolate (hydroxyacetate)0.87Xenobiotics0.0021119177
Hydroquinone sulfate1.57Xenobiotics0.0051383189
4-Ethylphenyl sulfate1.50Xenobiotics0.0063570201.1
Arabonate0.69Cofactors and vitamins0.0061736292.1
Pantothenate (Vitamin B5)1.32Cofactors and vitamins0.0062218220.1
Pyroglutamine0.74Amino acid0.006764129.2
Citrulline1.21Amino acid0.010715176.1
Tyrosine0.86Amino acid0.0101516182.1
Gulono-1,4-lactone0.68Cofactors and vitamins0.0111862333.1
Methylpalmitate (15 or 2)1.23Lipid0.0115698269.4
Arginine1.15Amino acid0.015728173.2
Ophthalmate0.47Amino acid0.0231457290.1
Homocitrulline0.77Amino acid0.024832190.1
N-Formylmethionine0.89Amino acid0.0291541176.1
Bilirubin (E,E)0.50Cofactors and vitamins0.0314625585.2
Palmitoyl sphingomyelin0.84Lipid0.0322524311.3
Cysteine-glutathione disulfide0.81Amino acid0.038821427.1
4-Vinylphenol sulfate1.29Xenobiotics0.0403323199.1
Dihomolinolenate (20:3n3 or 3n6)1.19Lipid0.0435600305.4
Anthranilate1.36Amino acid0.0493213138.1
Lysine0.75Amino acid0.0491836.7317.2
Threonate0.74Cofactors and vitamins0.0491560.7292.1

Ratios of median transformed values, Benjamini-Hochberg adjusted values, metabolite biological class, retention index (RI), and mass (Da) are shown. X-##### metabolites represent distinct chemical entities of discrete molecular formula and structure but do not currently match a named biochemical in our reference library.

Clinical measureMean smallSEM smallMean largerSEM largerAdj. P

DEXA tissue (gm)9135.12401.8725392569.82
DEXA lean (gm)6884.76308.5617942.89491.91
DEXA bone mineral content (gm)379.5815.991009.0624.88
Weight (gm)9486.15421.5831768.721302.21
DEXA fat (gm)2250.42145.257449.06427.25
DEXA total mass (kg)9.510.4226.400.59
Bone density (gm/cm2)0.650.010.780.01
Serum creatinine (mg/dL)0.640.020.860.02
Serum creatine kinase (IU/L)263.8833.51154.0421.37
Serum potassium (mmol/L)4.310.054.650.05
Serum total bilirubin (mg/dL)
Serum aspartate transaminase (U/L)35.681.9927.911.070.02
Serum gamma-glutamyltransferase (g/dL)6.930.435.940.360.02
Serum total antioxidant status (mmol/L)1.490.031.600.020.03
Serum T3 (nmol/L)
Protein digestibility (percentage)88.000.3886.000.570.01

Mean, standard error of the mean (SEM), and adjusted P values are shown. DEXA, dual-energy X-ray absorptiometry; T3, total triiodothyronine.
3.2. Antioxidant Status Is Lower in Small Dogs

Total antioxidant status, bilirubin, glutathione metabolites, and urate (measures of overall antioxidant status) were lower in small dogs compared to larger dogs (Tables 1 and 2). Total antioxidant status is representative of antioxidants present in the blood that inhibit the oxidation reaction assay. Inhibition can be caused by many types of antioxidants and is thus not specific for an individual antioxidant. Bilirubin is associated with heme breakdown but is also a powerful antioxidant [22]. Glutathione metabolites, including multiple gamma-glutamyl amino acids, cysteine-glutathione, and 5-oxoproline, were lower in small dogs. Glutathione (GSH) is a powerful antioxidant present at high concentrations in cells. Urate, another powerful antioxidant, was also lower in small dogs. Urate, initially produced from xanthine by xanthine oxidase in the metabolism of purines, is converted to allantoin by the uricase enzyme in most species. In humans and great apes, however, the uricase enzyme activity is absent and urate concentrations in these species are elevated [23].

3.3. Circulating Levels of Amino Acids Are Different between Dogs of Different Body Size

Several amino acids and their metabolites differed between small and larger dogs (Table 1). Circulating levels of the essential amino acids phenylalanine, tyrosine, lysine, and the nonessential amino acids glutamine, hydroxyproline, and prolylhydroxyproline were lower in small dogs compared with larger dogs. Circulating levels of the essential amino acid arginine were higher in small dogs. Protein digestibility (Table 2) was unexpectedly higher in small dogs even though levels of many circulating amino acids were lower. The amino acid metabolites phenol sulfate and p-cresol-sulfate levels were higher in small dogs. These metabolites are formed in the liver by sulfation of bacterial-derived tyrosine metabolites [24].

3.4. Other Plasma Metabolites and Clinical Measures Are Different between Dogs of Different Body Size

Other metabolites and clinical measures involved in biological processes differed between small dogs and larger dogs (Tables 1 and 2). While most of these processes were represented by only one metabolite or clinical measure, they have physiological significance since they function across processes. Creatine levels were higher in small dogs as were two of its building blocks, citrulline and arginine (Table 1). Creatine can be endogenously synthesized by the kidneys, pancreas, and liver, or it can be ingested [25]. Most creatine is phosphorylated and stored as phosphocreatine in skeletal muscle, where it is used to replenish ATP concentrations. Creatine is recycled into creatinine, which is then excreted by the kidney. Creatinine levels were lower in small dogs.

Homocitrulline concentrations were lower in larger dogs. Homocitrulline is formed by the carbamylation of lysine residues and results in loss of the biological activity of lysine [26]. Triiodothyronine (T3) was higher in small dogs. T3 is a thyroid hormone that affects basal metabolic rate and increases metabolic consumption affecting both fatty acid and carbohydrate metabolism. Pantothenic acid (vitamin B5), an essential component of coenzyme A, was higher in small dogs. Bone mineral density was lower in small dogs consistent with the results of others [27].

3.5. The Fecal Microbiome Differed between Dogs of Different Body Size

The fecal microbiome was analyzed to determine if small dogs had differences in bacterial populations compared to larger dogs. Since bacterial populations can change in as little as one day in response to a macronutrient change [28] and all dogs were fed identical diets, fecal samples were taken after 5 weeks, based on the assumption that bacterial populations had normalized. 16S analysis identified 6 bacteria ( value < 0.05) at the genus level that differed between small dogs and larger dogs. Of these, Bacteroides, Faecalibacterium (higher in small dogs), and Collinsella and Lactobacillus (higher in larger dogs) had the highest mean proportion difference (mean relative frequency > 1%) (Figure 1).

3.6. Plasma Metabolites and Clinical Measures Associate with a Kidney Function

Visualization with PCA performed using only selected clinical or metabolite parameters simultaneously confirmed clear discrimination between small and larger dogs (Figures 2(a) and 2(b)). In order to understand relationships between metabolites and clinical parameters, metabolite small dog modules were first identified with WGCNA using the 449 metabolites (Figure 3). 8 modules were identified (Table 3). The blue module containing 64 named metabolites (Supplementary Table ) had the highest ability to distinguish between small dogs and larger dogs (lowest -score and thus lowest preservation when the reference group was changed). A correlation analysis was then performed between the clinical parameters and metabolomics modules eigengenes (Figure 4). Age (negative), creatinine, and bilirubin (both positive) among other clinical parameters showed significant correlations (Figure 4) with blue module metabolites. Metabolites and clinical parameters were analyzed for associations with kidney function. Greater than 30% of the metabolites (Supplementary Table ) and nearly all blue module associated clinical parameters (Supplementary Table ) have associations with kidney function.



-scores associated with each metabolite module. Greater -scores depict greater preservation when the reference is changed from small to larger dogs. Thus, modules with greater -scores have a lower ability to distinguish between groups.

4. Discussion

Although Monod’s dictum that anything true of E. coli is also true of the elephant describes well the commonality of regulatory and biochemical reactions across the tree of life, subtle and important differences in transcriptional and metabolic processes occur between species and within individuals of the same species. Analysis of plasma metabolites and the microbiome revealed differences based on body size between small and larger dogs. Hierarchical clustering of the final, preprocessed metabolites identified metabolites with similar concentrations (i.e., the blue module) that had the highest ability to distinguish between small and larger dogs. This module contained 64 named metabolites in different metabolic classes (amino acid, nucleotides, carbohydrates, xenobiotics, and cofactors; see Supplementary Table ).

Small dogs have lower plasma levels of multiple antioxidants, including total antioxidant status, urate, glutathione metabolites, and bilirubin as compared to larger dogs suggesting a higher mass-specific metabolic rate [3] in tissues. Lower levels of antioxidants may be an adaptation to a higher production of free radicals because of the higher basal metabolic rate in small dogs. Bilirubin exerts its most potent antioxidant effects against lipid oxidation, while water-soluble urate and glutathione have much more potent antioxidant effects on protein [29], which when considered together indicated a systemic difference between small and larger dogs.

Urate concentrations were higher in larger dogs compared with small dogs. Urate has both anti- and prooxidant characteristics. Urate utilizes glutathione to reduce brain free radicals, by increasing cysteine uptake via EAAT-1 transporters in neurons of the hippocampus [30]. Urate has specifically demonstrated protective effects against peroxide [31], 1-methyl-4-phenyl-pyridinium [32], the reaction products of peroxynitrite and CO2, , and NO2 [33], as well as 6-hydroxydopamine (6-OHDA) [34]. Appropriate urate levels have protective effects against oxidative diseases such as Alzheimer’s Disease [35], Parkinson’s disease [36], and Multiple Sclerosis [37].

However, urate can also be prooxidant. A 1 mg/dL increase in serum urate concentrations has been associated with a significant increased incidence of hypertension [38]. Higher urate concentrations have also been shown to predict worsening of renal disease in patients with renal failure [39] and affect glomerular filtration rate in healthy study participants [40]. Proper regulation of urate concentration is thus important to mediate its anti- or prooxidant properties.

Glutathione is present in all cells in high concentrations and acts as a buffer in redox reactions. GSH is the most prevalent antioxidant in the liver, is stored in its reduced state, and protects hepatocytes against oxidative damage. GSH is involved in maintaining thiol disulfide balance, peroxide detoxification, leukotriene biosynthesis, and amino acid transport and has been shown to play a role in multiple cellular processes such as transcription, proliferation, and apoptosis [41]. Cardiac disease is associated with decreased GSH levels in dogs [42]. Compared with healthy dogs, GSH concentrations in clinically ill dogs are significantly decreased in erythrocytes, and GSH depletion correlates with severity of illness and mortality [43]. We showed here that multiple metabolites of GSH metabolism were lower in the plasma of small dogs suggesting utilization for GSH production. However, the true response of the GSH system in small dogs is not known since oxidative states of the GSH were not measured in this study.

We have shown that bilirubin concentrations were lower in small dogs. Bilirubin represents a lipophilic cytoprotectant antioxidant which is complementary to water-soluble antioxidants, such as GSH [29]. The clinical implications of low bilirubin in small dogs are unclear, although lower concentrations suggest increased metabolic consumption and the possibility of ongoing oxidative damage in small dogs.

Small dogs had lower levels of multiple amino acids, including the essential amino acids phenylalanine, tyrosine, and lysine. Arginine was an exception since it was found at higher levels in small dogs. We previously determined that protein digestibility was higher in small dogs suggesting that metabolic differences produced the lower levels of amino acids in small dogs.

Phenylalanine is involved in the synthesis of tyrosine, which is then used as a precursor for the synthesis of adrenaline, noradrenaline, and dopamine [44]. Tyrosine can be fermented in the large intestine, resulting in the microbial formation of p-cresol and indoxyl (a phenolic compound), which after hepatic sulfation become p-cresyl sulfate and phenol sulfate [24, 45]. Concentrations of both of these compounds were higher in small dogs suggesting that phenylalanine and tyrosine were lower due to fermentation and sulfation, perhaps through increased microbial metabolism. Levels of certain microbial populations affect formation of p-cresyl sulfate indicative of differences in the gut microflora between small and larger dogs [45]. 16S analysis revealed higher levels of Bacteroides in small dogs. Bacteroides fragilis has been shown to produce phenols from tyrosine [45]. Previous studies have shown differences between fecal versus illeal protein digestibility which takes into account microbial contributions to amino acid metabolism [46, 47] consistent with the results presented here.

We observed higher levels of arginine and lower levels of lysine in small dogs. Lysine and arginine share the cationic amino acid transporter system (CAT) mechanism, and concentrations of one competitively inhibit transport of the other which has been observed in mammalian cells [48, 49]. The clinical implications of differences in concentrations of these amino acids cannot be interpreted with the current data, although arginine supplementation is known to be beneficial in wound healing and is essential for maintaining intestinal integrity [50]. Additionally, the lower levels of plasma glutamine in small dogs may be due to absorption into enterocytes, since enterocytes use glutamine as their primary fuel source [51].

Creatine concentrations in plasma were lower in larger dogs. In addition to dietary source, creatine can be synthesized endogenously and is stored as phosphocreatine in muscles as an ATP buffer. Creatine is excreted as creatinine in urine and often serves as a marker of kidney function. Dogs with increased body weight have increased urinary creatinine excretion, suggesting increased creatine turnover, consistent with the higher levels of creatinine in larger dogs noted in our study [52]. However, creatinine may also be associated with increased muscle mass or muscle volume. Sighthounds have more muscle volume than breeds of similar body weight with concomitant increased creatinine concentrations when compared with other breeds [53, 54]. Reduced concentrations of creatine might be associated with increased storage as phosphocreatine, increased usage in skeletal muscle, and consequent increased creatinine excretion.

Plasma homocitrulline was higher in larger dogs. Elevated homocitrulline concentrations have been correlated with the progression of renal disease [55, 56]. Homocitrulline concentration is inversely associated with estimated glomerular filtration rate [55, 57]. Kidney disease has a prevalence of 0.5–1.5% in the canine population [58]. In a retrospective study in Sweden an average of 15.8 cases of kidney disease were diagnosed per 10,000 dog years at risk (DYAR), with a mortality of 9.7 deaths per 10,000 DYAR [59]. Of the breeds exceeding the mean incidence of kidney disease, 23 are considered breeds of larger body size. However, the three breeds with the lowest incidence of mortality due to kidney disease were all small breeds. This is especially interesting because a different survey-based study found that smaller dogs are at greater risk for kidney disease: with every decrease of 10 kg in body weight, a 50% increased risk of kidney disease was found [60].

Many other metabolites in the blue module were associated with kidney function (Supplementary Table ). Other clinical parameters correlated with metabolites in this module also have associations with kidney function (Supplementary Table ). Based on these data, this metabolite module and its clinical correlates may contribute to a kidney metabolite network. Many of the differences noted here between small and larger dogs can be attributed to functional differences rather than disease since the animals used in this study did not have clinical kidney conditions. However, no predispositions are known between the small and large breed dog populations, but more research needs to be done here to determine the clinical significance.

5. Conclusions

We have described here metabolic differences between small dogs and dogs of larger body sizes. This included the difference in circulating metabolites, clinical parameters, and microbiota from dogs in a diet and environment-controlled study. Small dogs had a lower antioxidant status as measured by multiple metabolites and clinical parameters. Differences were also shown in circulating amino acids, some of which could be tied to variations in specific bacteria of the microbiota. Additionally, analysis of small dog metabolites and clinical parameters reflected a network which strongly associates with kidney function. This analysis represents a unique and initial view of metabolic differences between body sizes within a mammalian species. These differences reflect not only morphometric induced variability, but also metabolic-specific genetic differences associated with the creation of breeds due to artificial selection. More will need to be done in order to identify other unique metabolic characteristics of small dogs including studies representing more breeds belonging to the body size classes of Canis lupus familiaris.


AST:Aspartate transaminase
DEXA:Dual-energy X-ray absorptiometry
DYAR:Dog years at risk
GC/MS:Gas chromatography/mass spectrometry
LC/MS:Liquid chromatography/tandem mass spectrometry
OTU:Operational taxonomic units
PCA:Principal component analysis
RF:Random forest
TAS:Total antioxidant status
WGCNA:Weighted gene coexpression network analysis.

Ethical Approval

The Nestlé Purina Animal Care and Use Committee approved the study protocol.


Present address of Jim Kaput is Vydiant, Inc., Gold River, CA. R. Middleton, C. Petzinger-Germain, and A. Beloshapka are employees of Nestlé Purina Research. J. Carayol is an employee of Nestlé Institute of Health Sciences and J. Kaput was employed by Nestlé Institute of Health Sciences during the project work reflected in this manuscript.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this article.


The authors would like to thank Rachel Anderson, Patricia Turpin, Cathy Chase, Aaron Parker, and Carolyn Cupp for their help with conducting the animal trial and sample collection. They thank Qinghong Li, Jeff Labuda, and Ziad Ramadan for their help with the fecal microbiome analysis. They thank Wendell Kerr and Xuemei Si for their help with the clinical statistics analysis. This work was funded by Nestlé Purina.

Supplementary Materials

Supplemental Table 1 provides the information for all 449 metabolites measured. Supplemental Table 2 provides information for clinical parameters measures. Supplemental Table 3 provides the information for metabolites in the blue module. Supplemental Table 4 provides the association of the blue module metabolites with kidney function. Supplemental Table 5 provides the association of blue module associated-clinical parameters with kidney function.

  1. Supplementary Material


  1. C. L. Turner and D. Christian, “Health Risks in Small Dogs,” Supplement to Compendium on Continuing Education for the Practicing Veterinarian, vol. 24, p. 76, 2002. View at: Google Scholar
  2. J. M. Fleming, K. E. Creevy, and D. E. L. Promislow, “Mortality in north american dogs from 1984 to 2004: an investigation into age-, size-, and breed-related causes of death,” Journal of Veterinary Internal Medicine, vol. 25, no. 2, pp. 187–198, 2011. View at: Publisher Site | Google Scholar
  3. A. G. Jimenez, “Physiological underpinnings in life-history trade-offs in man’s most popular selection experiment: the dog,” Journal of Comparative Physiology B: Biochemical, Systemic, and Environmental Physiology, vol. 186, no. 7, pp. 813–827, 2016. View at: Publisher Site | Google Scholar
  4. V. J. Adams, K. M. Evans, J. Sampson, and J. L. Wood, “Methods and mortality results of a health survey of purebred dogs in the UK,” J Small Anim Pract, vol. 51, no. 10, pp. 512–524, 2010. View at: Publisher Site | Google Scholar
  5. Q. Li, L. M. Freeman, J. E. Rush et al., “Veterinary medicine and multi-omics research for future nutrition targets: metabolomics and transcriptomics of the common degenerative mitral valve disease in dogs,” OMICS: A Journal of Integrative Biology, vol. 19, no. 8, pp. 461–470, 2015. View at: Publisher Site | Google Scholar
  6. A. M. Evans, C. D. DeHaven, T. Barrett, M. Mitchell, and E. Milgram, “Integrated, nontargeted ultrahigh performance liquid chromatography/electrospray ionization tandem mass spectrometry platform for the identification and relative quantification of the small-molecule complement of biological systems,” Analytical Chemistry, vol. 81, no. 16, pp. 6656–6667, 2009. View at: Publisher Site | Google Scholar
  7. C. D. DeHaven, A. M. Evans, H. Dai, and K. A. Lawton, “Organization of GC/MS and LC/MS metabolomics data into chemical libraries,” J Cheminformatics, vol. 2, pp. 9–21, 2010. View at: Google Scholar
  8. L. K. Karr-Lilienthal, C. M. Grieshop, J. K. Spears et al., “Estimation of the proportion of bacterial nitrogen in canine feces using diaminopimelic acid as an internal bacterial marker,” Journal of Animal Science, vol. 82, no. 6, pp. 1707–1712, 2004. View at: Publisher Site | Google Scholar
  9. J. Reeder and R. Knight, “Rapidly denoising pyrosequencing amplicon reads by exploiting rank-abundance distributions,” Nature Methods, vol. 7, no. 9, pp. 668-669, 2010. View at: Publisher Site | Google Scholar
  10. R. C. Edgar, B. J. Haas, J. C. Clemente, C. Quince, and R. Knight, “UCHIME improves sensitivity and speed of chimera detection,” Bioinformatics, vol. 27, no. 16, pp. 2194–2200, 2011. View at: Publisher Site | Google Scholar
  11. T. Z. DeSantis, P. Hugenholtz, N. Larsen et al., “Greengenes, a chimera-checked 16S rRNA gene database and workbench compatible with ARB,” Applied and Environmental Microbiology, vol. 72, no. 7, pp. 5069–5072, 2006. View at: Publisher Site | Google Scholar
  12. Q. Wang, G. M. Garrity, J. M. Tiedje, and J. R. Cole, “Naive Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy,” Appl Environ Microbiol, vol. 73, pp. 5261–5267, 2007. View at: Google Scholar
  13. M. Hubert, P. J. Rousseeuw, and K. V. Branden, “ROBPCA: a new approach to robust principal component analysis,” Technometrics, vol. 47, no. 1, pp. 64–79, 2005. View at: Publisher Site | Google Scholar
  14. V. Todorov and P. Filzmoser, “An object-oriented framework for robust multivariate analysis,” Journal of Statistical Software , vol. 32, no. 3, pp. 1–47, 2009. View at: Google Scholar
  15. E. Archer, rfPermute: Estimate Permutation p-Values for Random Forest Importance Metrics, R package version 2.1.5, 2016,
  16. L. Breiman, “Machine Learning,” vol. 45, 2001. View at: Publisher Site | Google Scholar
  17. E. Saccenti, H. C. J. Hoefsloot, A. K. Smilde, J. A. Westerhuis, and M. M. W. B. Hendriks, “Reflections on univariate and multivariate analysis of metabolomics data,” Metabolomics, vol. 10, no. 3, pp. 361–374, 2014. View at: Publisher Site | Google Scholar
  18. P. Langfelder and S. Horvath, “WGCNA: an R package for weighted correlation network analysis,” BMC Bioinformatics, vol. 9, article 559, 2008. View at: Publisher Site | Google Scholar
  19. Y. Benjamini and Y. Hochberg, “Controlling the false discovery rate: a practical and powerful approach to multiple testing,” Journal of the Royal Statistical Society: Series B (Statistical Methodology), vol. 57, pp. 289–300, 1995. View at: Google Scholar
  20. P. Langfelder, R. Luo, M. C. Oldham, and S. Horvath, “Is My Network Module Preserved and Reproducible?” PLoS Comput Biol, vol. 7, Article ID 1001057, p. 10, 2011. View at: Google Scholar
  21. R Core Team, R: A language and environment for statistical computing, R Foundation for Statistical Computing, Vienna, Austria, 2016,
  22. D. E. Barañano, M. Rao, C. D. Ferris, and S. H. Snyder, “Biliverdin reductase: a major physiologic cytoprotectant,” Proceedings of the National Acadamy of Sciences of the United States of America, vol. 99, no. 25, pp. 16093–16098, 2002. View at: Publisher Site | Google Scholar
  23. M. Oda, Y. Satta, O. Takenaka, and N. Takahata, “Loss of urate oxidase activity in hominoids and its evolutionary implications,” Molecular Biology and Evolution, vol. 19, no. 5, pp. 640–653, 2002. View at: Publisher Site | Google Scholar
  24. W. R. Wikoff, A. T. Anfora, and J. Liu, “Metabolomics analysis reveals large effects of gut microflora on mammalian blood metabolites,” Proceedings of the National Acadamy of Sciences of the United States of America, vol. 106, no. 10, pp. 3698–3703, 2009. View at: Publisher Site | Google Scholar
  25. M. G. Bemben and H. S. Lamont, “Creatine supplementation and exercise performance: Recent findings,” Sports Medicine, vol. 35, no. 2, pp. 107–125, 2005. View at: Publisher Site | Google Scholar
  26. Z. Wang, S. J. Nicholls, E. R. Rodriguez et al., “Protein carbamylation links inflammation, smoking, uremia and atherogenesis,” Nature Medicine, vol. 13, no. 10, pp. 1176–1184, 2007. View at: Publisher Site | Google Scholar
  27. C. J. Cupp and G. Czarnecki-Maulden, “Effect of breed and size of dog on body composition,” J An Sci, vol. 82, supplement 1, p. 244, 2004. View at: Google Scholar
  28. L. A. David, C. F. Maurice, R. N. Carmody et al., “Diet rapidly and reproducibly alters the human gut microbiome,” Nature, vol. 505, no. 7484, pp. 559–563, 2014. View at: Publisher Site | Google Scholar
  29. T. W. Sedlak, M. Saleh, D. S. Higginson, B. D. Paul, K. R. Juluri, and S. H. Snyder, “Bilirubin and glutathione have complementary antioxidant and cytoprotective roles,” Proceedings of the National Acadamy of Sciences of the United States of America, vol. 106, no. 13, pp. 5171–5176, 2009. View at: Publisher Site | Google Scholar
  30. K. Aoyama, N. Matsumura, M. Watabe, F. Wang, K. Kikuchi-Utsumi, and T. Nakaki, “Caffeine and uric acid mediate glutathione synthesis for neuroprotection,” Neuroscience, vol. 181, pp. 206–215, 2011. View at: Publisher Site | Google Scholar
  31. B. N. Ames, R. Cathcart, E. Schwiers, and P. Hochstein, “Uric acid provides an antioxidant defense in humans against oxidant- and radical-caused aging and cancer: a hypothesis,” Proceedings of the National Acadamy of Sciences of the United States of America, vol. 78, no. 11, pp. 6858–6862, 1981. View at: Publisher Site | Google Scholar
  32. S. Cipriani, C. A. Desjardins, T. C. Burdett, Y. Xu, K. Xu, and M. A. Schwartzchild, “Urate and its transgenic depletion modulate neuronal vulnerability in a cellular model of Parkinsons disease,” PLoS One, vol. 7, Article ID e37331, 2012. View at: Google Scholar
  33. G. L. Squadrito, R. Cueto, A. E. Splenser et al., “Reaction of uric acid with peroxynitrite and implications for the mechanism of neuroprotection by uric acid,” Archives of Biochemistry and Biophysics, vol. 376, no. 2, pp. 333–337, 2000. View at: Publisher Site | Google Scholar
  34. L. Gong, Q.-L. Zhang, N. Zhang et al., “Neuroprotection by urate on 6-OHDA-lesioned rat model of Parkinson's disease: Linking to Akt/GSK3β signaling pathway,” Journal of Neurochemistry, vol. 123, no. 5, pp. 876–885, 2012. View at: Publisher Site | Google Scholar
  35. T.-S. Kim, C.-U. Pae, S.-J. Yoon et al., “Decreased plasma antioxidants in patients with Alzheimer's disease,” International Journal of Geriatric Psychiatry, vol. 21, no. 4, pp. 344–348, 2006. View at: Publisher Site | Google Scholar
  36. J. W. Davis, A. Grandinetti, C. I. Waslien, G. W. Ross, L. R. White, and D. M. Morens, “Observations on serum uric acid levels and the risk of idiopathic Parkinson's disease,” American Journal of Epidemiology, vol. 144, no. 5, pp. 480–484, 1996. View at: Publisher Site | Google Scholar
  37. D. C. Hooper, G. S. Scott, A. Zborek et al., “Uric acid, a peroxynitrite scavenger, inhibits CNS inflammation, blood- CNS barrier permeability changes, and tissue damage in a mouse model of multiple sclerosis,” The FASEB Journal, vol. 14, no. 5, pp. 691–698, 2000. View at: Google Scholar
  38. P. C. Grayson, S. Young Kim, M. Lavalley, and H. K. Choi, “Hyperuricemia and incident hypertension: a systematic review and meta-analysis,” Arthritis Care & Research, vol. 63, no. 1, pp. 102–110, 2011. View at: Publisher Site | Google Scholar
  39. G. Bellomo, S. Venanzi, C. Verdura, P. Saronio, A. Esposito, and M. Timio, “Association of uric acid with change in kidney function in healthy normotensive individuals,” American Journal of Kidney Diseases, vol. 56, no. 2, pp. 264–272, 2010. View at: Publisher Site | Google Scholar
  40. C.-Y. Hsu, C. Iribarren, C. E. McCulloch, J. Darbinian, and A. S. Go, “Risk factors for end-stage renal disease - 25-year follow-up,” JAMA Internal Medicine, vol. 169, no. 4, pp. 342–350, 2009. View at: Publisher Site | Google Scholar
  41. E. V. Kalinina, N. N. Chernov, and M. D. Novichkova, “Role of glutathione, glutathione transferase, and glutaredoxin in regulation of redox-dependent processes,” Biochemistry (Moscow), vol. 79, no. 13, pp. 1562–1583, 2014. View at: Publisher Site | Google Scholar
  42. L. M. Freeman, J. E. Rush, P. E. Milbury, and J. B. Blumberg, “Antioxidant status and biomarkers of oxidative stress in dogs with congestive heart failure,” Journal of Veterinary Internal Medicine, vol. 19, no. 4, pp. 537–541, 2005. View at: Publisher Site | Google Scholar
  43. K. R. Viviano, S. N. Lavergne, L. Goodman et al., “Glutathione, cysteine, and ascorbate concentrations in clinically ill dogs and cats,” Journal of Veterinary Internal Medicine, vol. 23, no. 2, pp. 250–257, 2009. View at: Publisher Site | Google Scholar
  44. P. F. Fitzpatrick, “The aromatic amino acid hydroxylases,” Advances in Enzymology and Related Areas of Molecular Biology, vol. 74, pp. 235–294, 2000. View at: Publisher Site | Google Scholar
  45. E. Bone, A. Tamm, and M. Hill, “The production of urinary phenols by gut bacteria and their possible role in the causation of large bowel cancer,” American Journal of Clinical Nutrition, vol. 29, no. 12, pp. 1448–1454, 1976. View at: Google Scholar
  46. M. Fuller, “Determination of protein and amino acid digestibility in foods including implications of gut microbial amino acid synthesis,” British Journal of Nutrition, vol. 108, no. 2, pp. S238–S246, 2012. View at: Publisher Site | Google Scholar
  47. P. Reeds, G. Schaafsma, D. Tomé, and V. Young, “Criteria and significance of dietary protein sources in humans. Summary of the workshop with recommendations,” J Nutr, vol. 130, pp. 1874S–1876S, 2000. View at: Google Scholar
  48. R. L. O’Kane, J. R. Vina, I. Simpson, R. Zaragoza, A. Mokashi, and R. A. Hawkins, “Cationic amino acid transport across the blood-brain barrier is mediated exclusively by system y+,” American Journal of Physiology Endocrinology and Metabolism Clinics of North America, vol. 291, pp. E412–E419, 2009. View at: Google Scholar
  49. B. Stieger, G. Stange, J. Biber, and H. Murer, “Transport of l-lysine by rat renal brush border membrane vesicles,” Pflügers Archiv - European Journal of Physiology, vol. 397, no. 2, pp. 106–113, 1983. View at: Publisher Site | Google Scholar
  50. K. Singh, L. A. Coburn, D. P. Barry, J.-L. Boucher, R. Chaturvedi, and K. T. Wilson, “L-arginine uptake by cationic amino acid transporter 2 is essential for colonic epithelial cell restitution,” American Journal of Physiology-Gastrointestinal and Liver Physiology, vol. 302, no. 9, pp. G1061–G1073, 2012. View at: Publisher Site | Google Scholar
  51. S. T. O'Dwyer, R. J. Smith, T. L. Hwang, and D. W. Wilmore, “Maintenance of Small Bowel Mucosa with Glutamine-Enriched Parenteral Nutrition,” Journal of Parenteral and Enteral Nutrition, vol. 13, no. 6, pp. 579–585, 1989. View at: Publisher Site | Google Scholar
  52. S. A. Center, E. Wilkinson, C. A. Smith, H. Erb, and R. M. Lewis, “24-Hour urine protein/creatinine ratio in dogs with protein-losing nephropathies.,” Journal of the American Veterinary Medical Association, vol. 187, no. 8, pp. 820–824, 1985. View at: Google Scholar
  53. W. E. Feeman III, C. G. Couto, and T. L. Gray, “Serum creatinine concentrations in retired racing greyhounds,” Veterinary Clinical Pathology, vol. 32, no. 1, pp. 40–42, 2003. View at: Publisher Site | Google Scholar
  54. M. Hillpö, “Some haematological and clinical-chemical parameters of sight hounds (Afghan hound, saluki and whippet),” Nordisk Veterinærmedicin, vol. 38, no. 3, pp. 148–155, 1986. View at: Google Scholar
  55. W. H. W. Tang, K. Shrestha, Z. Wang et al., “Protein carbamylation in chronic systolic heart failure: Relationship with renal impairment and adverse long-term outcomes,” Journal of Cardiac Failure, vol. 19, no. 4, pp. 219–224, 2013. View at: Publisher Site | Google Scholar
  56. L. M. Kraus, L. Gaber, C. R. Handorf, H. Marti, Krais. Jr, and A. P, “Carbamylation of glomerular and tubular proteins in patients with kidney failure: a potential mechanism of ongoing renal damage,” Swiss Medical Weekly, vol. 131, pp. 139–145, 2001. View at: Google Scholar
  57. A. S. Levey, J. P. Bosch, J. B. Lewis, T. Greene, N. Rogers, and D. Roth, “A more accurate method to estimate glomerular filtration rate from serum creatinine: a new prediction equation,” Annals of Internal Medicine, vol. 130, no. 6, pp. 461–470, 1999. View at: Publisher Site | Google Scholar
  58. D. J. Polzin, “Chronic kidney disease,” in Textbook of Veterinary Internal Medicine, S. J. Ettinger and E. C. Feldman, Eds., pp. 1990–2021, Elsevier Saunders, St. Louis, Missouri, MO, USA, 7th edition, 2010. View at: Google Scholar
  59. L. Pelander, I. Ljungvall, A. Egenvall, H. Syme, J. Elliott, and J. Häggström, “Incidence of and mortality from kidney disease in over 600,000 insured Swedish dogs,” Veterinary Record, vol. 176, no. 25, p. 656, 2015. View at: Publisher Site | Google Scholar
  60. P. C. Bartlett, J. W. Van Buren, A. D. Bartlett, and C. Zhou, “Case-control study of risk factors associated with feline and canine chronic kidney disease,” Veterinary Medicine International, vol. 2010, Article ID 957570, 9 pages, 2010. View at: Publisher Site | Google Scholar

Copyright © 2017 Rondo P. Middleton et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Related articles

No related content is available yet for this article.
 PDF Download Citation Citation
 Download other formatsMore
 Order printed copiesOrder

Related articles

No related content is available yet for this article.

Article of the Year Award: Outstanding research contributions of 2021, as selected by our Chief Editors. Read the winning articles.