Maintaining homeostasis in higher organisms involves a complex interplay of multiple ubiquitous and organ-specific molecular mechanisms that can be characterized using functional genomics technologies such as transcriptomics, proteomics, and metabonomics and dissected out through genetic investigations in healthy and diseased individuals. We characterized the genomic, metabolic, and physiological divergence of several inbred rat strains-Brown Norway, Lewis, Wistar Kyoto, Fisher (F344)-frequently used as healthy controls in genetic studies of the cardiometabolic syndrome. Hierarchical clustering of H-1 NMR-based metabolic profiles (n = 20 for urine, n = 16 for plasma) identified metabolic phenotype (metabotype) divergence patterns similar to the phylogenetic variability based on single nucleotide polymorphisms. However, the observed urinary metabotype variation exceeded that explainable by genetic polymorphisms. To understand further this natural variation, we used an integrative, knowledge-based network biology metabolic pathway analysis approach, coined Metabolite-Set Enrichment Analysis (MSEA). MSEA reveals that homeostasis and physiological plasticity can be achieved despite widespread divergences in glucose, lipid, amino acid, and energy metabolism in the host, together with different gut microbiota contributions suggestive of strain-specific transgenomic interactions. This work illustrates the concept of natural metabolomic variation, leading to physiologically stable albeit diverse strategies within the range of normality, all of which are highly relevant to animal model physiology, genetical genomics, and patient stratification in personalized healthcare.