Circulating cell-free DNA in health and disease - the relationship to health behaviours, ageing phenotypes and metabolomics.

Circulating cell-free DNA (cf-DNA) has emerged as a promising biomarker of ageing, tissue damage and cellular stress. However, less is known about health behaviours, ageing phenotypes and metabolic processes that lead to elevated cf-DNA levels. We sought to analyse the relationship of circulating cf-DNA level to age, sex, smoking, physical activity, vegetable consumption, ageing phenotypes (physical functioning, the number of diseases, frailty) and an extensive panel of biomarkers including blood and urine metabolites and inflammatory markers in three human cohorts (N = 5385; 17-82 years). The relationships were assessed using correlation statistics, and linear and penalised regressions (the Lasso), also stratified by sex.cf-DNA levels were significantly higher in men than in women, and especially in middle-aged men and women who smoke, and in older more frail individuals. Correlation statistics of biomarker data showed that cf-DNA level was higher with elevated inflammation (C-reactive protein, interleukin-6), and higher levels of homocysteine, and proportion of red blood cells and lower levels of ascorbic acid. Inflammation (C-reactive protein, glycoprotein acetylation), amino acids (isoleucine, leucine, tyrosine), and ketogenesis (3-hydroxybutyrate) were included in the cf-DNA level-related biomarker profiles in at least two of the cohorts.In conclusion, circulating cf-DNA level is different by sex, and related to health behaviour, health decline and metabolic processes common in health and disease. These results can inform future studies where epidemiological and biological pathways of cf-DNA are to be analysed in details, and for studies evaluating cf-DNA as a potential clinical marker.

Whole-genome sequencing analysis of semi-supercentenarians.

Extreme longevity is the paradigm of healthy aging as individuals who reached the extreme decades of human life avoided or largely postponed all major age-related diseases. In this study, we sequenced at high coverage (90X) the whole genome of 81 semi-supercentenarians and supercentenarians [105+/110+] (mean age: 106.6 ± 1.6) and of 36 healthy unrelated geographically matched controls (mean age 68.0 ± 5.9) recruited in Italy. The results showed that 105+/110+ are characterized by a peculiar genetic background associated with efficient DNA repair mechanisms, as evidenced by both germline data (common and rare variants) and somatic mutations patterns (lower mutation load if compared to younger healthy controls). Results were replicated in a second independent cohort of 333 Italian centenarians and 358 geographically matched controls. The genetics of 105+/110+ identified DNA repair and clonal haematopoiesis as crucial players for healthy aging and for the protection from cardiovascular events.

Genomic history of the Italian population recapitulates key evolutionary dynamics of both Continental and Southern Europeans.

BACKGROUND: The cline of human genetic diversity observable across Europe is recapitulated at a micro-geographic scale by variation within the Italian population. Besides resulting from extensive gene flow, this might be ascribable also to local adaptations to diverse ecological contexts evolved by people who anciently spread along the Italian Peninsula. Dissecting the evolutionary history of the ancestors of present-day Italians may thus improve the understanding of demographic and biological processes that contributed to shape the gene pool of European populations. However, previous SNP array-based studies failed to investigate the full spectrum of Italian variation, generally neglecting low-frequency genetic variants and examining a limited set of small effect size alleles, which may represent important determinants of population structure and complex adaptive traits. To overcome these issues, we analyzed 38 high-coverage whole-genome sequences representative of population clusters at the opposite ends of the cline of Italian variation, along with a large panel of modern and ancient Euro-Mediterranean genomes. RESULTS: We provided evidence for the early divergence of Italian groups dating back to the Late Glacial and for Neolithic and distinct Bronze Age migrations having further differentiated their gene pools. We inferred adaptive evolution at insulin-related loci in people from Italian regions with a temperate climate, while possible adaptations to pathogens and ultraviolet radiation were observed in Mediterranean Italians. Some of these adaptive events may also have secondarily modulated population disease or longevity predisposition. CONCLUSIONS: We disentangled the contribution of multiple migratory and adaptive events in shaping the heterogeneous Italian genomic background, which exemplify population dynamics and gene-environment interactions that played significant roles also in the formation of the Continental and Southern European genomic landscapes.

Insulin Resistance during normal child growth and development is associated with a distinct blood metabolic phenotype (Earlybird 72).

BACKGROUND: While insulin resistance (IR) is associated with specific metabolite signatures in adults, there have been few truly longitudinal studies in healthy children, either to confirm which abnormalities are present, or to determine whether they precede or result from IR. Therefore, we investigated the association of serum metabolites with IR in childhood in the Earlybird cohort. METHODS: The Earlybird cohort is a well-characterized cohort of healthy children with annual measurements from age 5 to 16 years. For the first time, longitudinal association analyses between individual serum metabolites and homeostatic model assessment (HOMA) of insulin resistance (HOMA-IR) have been performed taking into account the effects of age, growth, puberty, adiposity, and physical activity. RESULTS: IR was higher in girls than in boys and was associated with increasing body mass index (BMI). In longitudinal analysis IR was associated with reduced concentrations of branched-chain amino acids (BCAA), 2-ketobutyrate, citrate and 3-hydroxybutyrate, and higher concentrations of lactate and alanine. These findings demonstrate the widespread biochemical consequences of IR for intermediary metabolism, ketogenesis, and pyruvate oxidation during normal child growth and development. CONCLUSIONS: Longitudinal analysis can differentiate metabolite signatures that precede or follow the development of greater levels of IR. In healthy normal weight children, higher levels of IR are associated with reduced levels of BCAA, ketogenesis, and fuel oxidation. In contrast, elevated lactate concentrations preceded the rise in IR. These changes reveal the metabolite signature of insulin action during normal growth, and they contrast with previous findings in obese children and adults that represent the consequences of IR and obesity.

Antioxidants linked with physical, cognitive and psychological frailty: Analysis of candidate biomarkers and markers derived from the MARK-AGE study.

Frailty among elderly people leads to an increased risk for negative health outcomes. To prevent frailty, we need a better understanding of the underlying mechanisms and early detection of individuals at risk. Both may be served by identifying candidate (bio)markers, i.e. biomarkers and markers, for the physical, cognitive, and psychological frailty domains. We used univariate (Rank-ANOVA) and multivariate (elastic net) approaches on the RASIG study population (age range: 35-74 years, n = 2220) of the MARK-AGE study to study up to 331 (bio)markers between individuals with and without frailty for each domain. Biomarkers and markers identified by both approaches were studied further regarding their association with frailty using logistic regression. Univariately, we found lower levels of antioxidants, including ß-cryptoxanthin and zeaxanthin, in those who were physically, cognitively or psychologically frail. Additionally, self-reported health was worse in these three frail groups. Multivariately, we observed lower levels of ß-cryptoxanthin and zeaxanthin in the cognitively frail. Levels of these carotenoids were inversely associated with the risk of being cognitively frail after adjusting for confounders. Antioxidants and self-reported health are potential (bio)markers to detect persons at risk of becoming frail. The biomarkers identified may indicate the involvement of inflammation in frailty, especially for physical and cognitive frailty.

Evaluation of Lymphocyte Response to the Induced Oxidative Stress in a Cohort of Ageing Subjects, including Semisupercentenarians and Their Offspring.

The production of reactive oxygen species (ROS) may promote immunosenescence if not counterbalanced by the antioxidant systems. Cell membranes, proteins, and nucleic acids become the target of ROS and progressively lose their structure and functions. This process could lead to an impairment of the immune response. However, little is known about the capability of the immune cells of elderly individuals to dynamically counteract the oxidative stress. Here, the response of the main lymphocyte subsets to the induced oxidative stress in semisupercentenarians (CENT), their offspring (OFF), elderly controls (CTRL), and young individuals (YO) was analyzed using flow cytometry. The results showed that the ratio of the ROS levels between the induced and noninduced (I/NI) oxidative stress conditions was higher in CTRL and OFF than in CENT and YO, in almost all T, B, and NK subsets. Moreover, the ratio of reduced glutathione levels between I/NI conditions was higher in OFF and CENT compared to the other groups in almost all the subsets. Finally, we observed significant correlations between the response to the induced oxidative stress and the degree of methylation in specific genes on the oxidative stress pathway. Globally, these data suggest that the capability to buffer dynamic changes in the oxidative environment could be a hallmark of longevity in humans.

Validation of the Brazilian Healthy Eating Index-Revised Using Biomarkers in Children and Adolescents.

The Brazilian Healthy Eating Index-Revised (BHEI-R) can be used to determine overall dietary patterns. We assessed the BHEI-R scores in children and adolescents, aged from 9 to 13 years old, and associated its component scores with biomarkers of health and dietary exposure. Three 24-h recalls were used to generate BHEI-R. Biomarkers were analyzed in plasma and red blood cells. Correlation tests, agreement, and covariance analyses were used to associate BHEI-R components with biomarkers. Data from 167 subjects were used. The strongest correlations were between fruits, vegetables and legumes with omega-6 and omega-3 fatty acids, and ß-carotene intakes. Milk and dairy correlated with plasma retinol and pyridoxine. All components rich in vegetable and animal protein sources correlated with plasma creatine. Total BHEI-R scores were positively associated with intakes of omega-6, omega-3, fiber and vitamin C, and inversely associated with energy and saturated fat intakes of individuals. Plasma ß-carotene and riboflavin biomarkers were positively associated with total BHEI-R. An inadequate food consumption pattern was captured by both biomarkers of health and dietary exposure. BHEI-R was validated for the above dietary components and can be associated with metabolomics and nutritional epidemiological data in future pediatric studies.

Circadian and Feeding Rhythms Orchestrate the Diurnal Liver Acetylome.

Lysine acetylation is involved in various biological processes and is considered a key reversible post-translational modification in the regulation of gene expression, enzyme activity, and subcellular localization. This post-translational modification is therefore highly relevant in the context of circadian biology, but its characterization on the proteome-wide scale and its circadian clock dependence are still poorly described. Here, we provide a comprehensive and rhythmic acetylome map of the mouse liver. Rhythmic acetylated proteins showed subcellular localization-specific phases that correlated with the related metabolites in the regulated pathways. Mitochondrial proteins were over-represented among the rhythmically acetylated proteins and were highly correlated with SIRT3-dependent deacetylation. SIRT3 activity being nicotinamide adenine dinucleotide (NAD)+ level-dependent, we show that NAD+ is orchestrated by both feeding rhythms and the circadian clock through the NAD+ salvage pathway but also via the nicotinamide riboside pathway. Hence, the diurnal acetylome relies on a functional circadian clock and affects important diurnal metabolic pathways in the mouse liver.

One-carbon metabolism, cognitive impairment and CSF measures of Alzheimer pathology: homocysteine and beyond.

BACKGROUND: Hyperhomocysteinemia is a risk factor for cognitive decline and dementia, including Alzheimer disease (AD). Homocysteine (Hcy) is a sulfur-containing amino acid and metabolite of the methionine pathway. The interrelated methionine, purine, and thymidylate cycles constitute the one-carbon metabolism that plays a critical role in the synthesis of DNA, neurotransmitters, phospholipids, and myelin. In this study, we tested the hypothesis that one-carbon metabolites beyond Hcy are relevant to cognitive function and cerebrospinal fluid (CSF) measures of AD pathology in older adults. METHODS: Cross-sectional analysis was performed on matched CSF and plasma collected from 120 older community-dwelling adults with (n = 72) or without (n = 48) cognitive impairment. Liquid chromatography-mass spectrometry was performed to quantify one-carbon metabolites and their cofactors. Least absolute shrinkage and selection operator (LASSO) regression was initially applied to clinical and biomarker measures that generate the highest diagnostic accuracy of a priori-defined cognitive impairment (Clinical Dementia Rating-based) and AD pathology (i.e., CSF tau phosphorylated at threonine 181 [p-tau181]/ß-Amyloid 1-42 peptide chain [Aß1-42] >0.0779) to establish a reference benchmark. Two other LASSO-determined models were generated that included the one-carbon metabolites in CSF and then plasma. Correlations of CSF and plasma one-carbon metabolites with CSF amyloid and tau were explored. LASSO-determined models were stratified by apolipoprotein E (APOE) ε4 carrier status. RESULTS: The diagnostic accuracy of cognitive impairment for the reference model was 80.8% and included age, years of education, Aß1-42, tau, and p-tau181. A model including CSF cystathionine, methionine, S-adenosyl-L-homocysteine (SAH), S-adenosylmethionine (SAM), serine, cysteine, and 5-methyltetrahydrofolate (5-MTHF) improved the diagnostic accuracy to 87.4%. A second model derived from plasma included cystathionine, glycine, methionine, SAH, SAM, serine, cysteine, and Hcy and reached a diagnostic accuracy of 87.5%. CSF SAH and 5-MTHF were associated with CSF tau and p-tau181. Plasma one-carbon metabolites were able to diagnose subjects with a positive CSF profile of AD pathology in APOE ε4 carriers. CONCLUSIONS: We observed significant improvements in the prediction of cognitive impairment by adding one-carbon metabolites. This is partially explained by associations with CSF tau and p-tau181, suggesting a role for one-carbon metabolism in the aggregation of tau and neuronal injury. These metabolites may be particularly critical in APOE ε4 carriers.

High-throughput method for the quantitation of metabolites and co-factors from homocysteine-methionine cycle for nutritional status assessment.

AIM: There is increasing interest in the profiling and quantitation of methionine pathway metabolites for health management research. Currently, several analytical approaches are required to cover metabolites and co-factors. RESULTS: We report the development and the validation of a method for the simultaneous detection and quantitation of 13 metabolites in red blood cells. The method, validated in a cohort of healthy human volunteers, shows a high level of accuracy and reproducibility. CONCLUSION: This high-throughput protocol provides a robust coverage of central metabolites and co-factors in one single analysis and in a high-throughput fashion. In large-scale clinical settings, the use of such an approach will significantly advance the field of nutritional research in health and disease.

Systems medicine of inflammaging.

Systems Medicine (SM) can be defined as an extension of Systems Biology (SB) to Clinical-Epidemiological disciplines through a shifting paradigm, starting from a cellular, toward a patient centered framework. According to this vision, the three pillars of SM are Biomedical hypotheses, experimental data, mainly achieved by Omics technologies and tailored computational, statistical and modeling tools. The three SM pillars are highly interconnected, and their balancing is crucial. Despite the great technological progresses producing huge amount of data (Big Data) and impressive computational facilities, the Bio-Medical hypotheses are still of primary importance. A paradigmatic example of unifying Bio-Medical theory is the concept of Inflammaging. This complex phenotype is involved in a large number of pathologies and patho-physiological processes such as aging, age-related diseases and cancer, all sharing a common inflammatory pathogenesis. This Biomedical hypothesis can be mapped into an ecological perspective capable to describe by quantitative and predictive models some experimentally observed features, such as microenvironment, niche partitioning and phenotype propagation. In this article we show how this idea can be supported by computational methods useful to successfully integrate, analyze and model large data sets, combining cross-sectional and longitudinal information on clinical, environmental and omics data of healthy subjects and patients to provide new multidimensional biomarkers capable of distinguishing between different pathological conditions, e.g. healthy versus unhealthy state, physiological versus pathological aging.

Decreased epigenetic age of PBMCs from Italian semi-supercentenarians and their offspring.

Given the dramatic increase in ageing populations, it is of great importance to understand the genetic and molecular determinants of healthy ageing and longevity. Semi-supercentenarians (subjects who reached an age of 105-109 years) arguably represent the gold standard of successful human ageing because they managed to avoid or postpone the onset of major age-related diseases. Relatively few studies have looked at epigenetic determinants of extreme longevity in humans. Here we test whether families with extreme longevity are epigenetically distinct from controls according to an epigenetic biomarker of ageing which is known as "epigenetic clock". We analyze the DNA methylation levels of peripheral blood mononuclear cells (PBMCs) from Italian families constituted of 82 semi-supercentenarians (mean age: 105.6 ± 1.6 years), 63 semi-supercentenarians' offspring (mean age: 71.8 ± 7.8 years), and 47 age-matched controls (mean age: 69.8 ± 7.2 years). We demonstrate that the offspring of semi-supercentenarians have a lower epigenetic age than age-matched controls (age difference=5.1 years, p=0.00043) and that centenarians are younger (8.6 years) than expected based on their chronological age. By contrast, no significant difference could be observed for estimated blood cell counts (such as naïve or exhausted cytotoxic T cells or helper T cells). Future studies will be needed to replicate these findings in different populations and to extend them to other tissues. Overall, our results suggest that epigenetic processes might play a role in extreme longevity and healthy human ageing.

The Use of Non-Variant Sites to Improve the Clinical Assessment of Whole-Genome Sequence Data.

Genetic testing, which is now a routine part of clinical practice and disease management protocols, is often based on the assessment of small panels of variants or genes. On the other hand, continuous improvements in the speed and per-base costs of sequencing have now made whole exome sequencing (WES) and whole genome sequencing (WGS) viable strategies for targeted or complete genetic analysis, respectively. Standard WGS/WES data analytical workflows generally rely on calling of sequence variants respect to the reference genome sequence. However, the reference genome sequence contains a large number of sites represented by rare alleles, by known pathogenic alleles and by alleles strongly associated to disease by GWAS. It's thus critical, for clinical applications of WGS and WES, to interpret whether non-variant sites are homozygous for the reference allele or if the corresponding genotype cannot be reliably called. Here we show that an alternative analytical approach based on the analysis of both variant and non-variant sites from WGS data allows to genotype more than 92% of sites corresponding to known SNPs compared to 6% genotyped by standard variant analysis. These include homozygous reference sites of clinical interest, thus leading to a broad and comprehensive characterization of variation necessary to an accurate evaluation of disease risk. Altogether, our findings indicate that characterization of both variant and non-variant clinically informative sites in the genome is necessary to allow an accurate clinical assessment of a personal genome. Finally, we propose a highly efficient extended VCF (eVCF) file format which allows to store genotype calls for sites of clinical interest while remaining compatible with current variant interpretation software.

MARK-AGE biomarkers of ageing.

Many candidate biomarkers of human ageing have been proposed in the scientific literature but in all cases their variability in cross-sectional studies is considerable, and therefore no single measurement has proven to serve a useful marker to determine, on its own, biological age. A plausible reason for this is the intrinsic multi-causal and multi-system nature of the ageing process. The recently completed MARK-AGE study was a large-scale integrated project supported by the European Commission. The major aim of this project was to conduct a population study comprising about 3200 subjects in order to identify a set of biomarkers of ageing which, as a combination of parameters with appropriate weighting, would measure biological age better than any marker in isolation.

Blood plasma lipidomic signature of epicardial fat in healthy obese women.

OBJECTIVES: A lipidomic approach was employed in a clinically well-defined cohort of healthy obese women to explore blood lipidome phenotype ascribed to body fat deposition, with emphasis on epicardial adipose tissue (EAT). METHODS: The present investigation delivered a lipidomics signature of epicardial adiposity under healthy clinical conditions using a cohort of 40 obese females (age: 25-45 years, BMI: 28-40 kg/m(2) ) not showing any metabolic disease traits. Lipidomics analysis of blood plasma was employed in combination with in vivo quantitation of mediastinal fat depots by computerized tomography. RESULTS: All cardiac fat depots correlated to indicators of hepatic dysfunctions (ALAT and ASAT), which describe physiological connections between hepatic and cardiac steatosis. Plasma lipidomics encompassed overall levels of lipid classes, fatty acid profiles, and individual lipid species. EAT and visceral fat associated with diacylglycerols (DAG), triglycerides, and distinct phospholipid and sphingolipid species. A pattern of DAG and phosphoglycerols was specific to EAT. CONCLUSIONS: Human blood plasma lipidomics appears to be a promising clinical and potentially diagnostic readout for patient stratification and monitoring. Association of blood lipidomics signature to regio-specific mediastinal and visceral adiposity under healthy clinical conditions may help provide more biological insights into obese patient stratification for cardiovascular disease risks.

The three genetics (nuclear DNA, mitochondrial DNA, and gut microbiome) of longevity in humans considered as metaorganisms.

Usually the genetics of human longevity is restricted to the nuclear genome (nDNA). However it is well known that the nDNA interacts with a physically and functionally separated genome, the mitochondrial DNA (mtDNA) that, even if limited in length and number of genes encoded, plays a major role in the ageing process. The complex interplay between nDNA/mtDNA and the environment is most likely involved in phenomena such as ageing and longevity. To this scenario we have to add another level of complexity represented by the microbiota, that is, the whole set of bacteria present in the different part of our body with their whole set of genes. In particular, several studies investigated the role of gut microbiota (GM) modifications in ageing and longevity and an age-related GM signature was found. In this view, human being must be considered as "metaorganism" and a more holistic approach is necessary to grasp the complex dynamics of the interaction between the environment and nDNA-mtDNA-GM of the host during ageing. In this review, the relationship between the three genetics and human longevity is addressed to point out that a comprehensive view will allow the researchers to properly address the complex interactions that occur during human lifespan.

Genome-wide association study of metabolic traits reveals novel gene-metabolite-disease links.

Metabolic traits are molecular phenotypes that can drive clinical phenotypes and may predict disease progression. Here, we report results from a metabolome- and genome-wide association study on (1)H-NMR urine metabolic profiles. The study was conducted within an untargeted approach, employing a novel method for compound identification. From our discovery cohort of 835 Caucasian individuals who participated in the CoLaus study, we identified 139 suggestively significant (P<5×10(-8)) and independent associations between single nucleotide polymorphisms (SNP) and metabolome features. Fifty-six of these associations replicated in the TasteSensomics cohort, comprising 601 individuals from São Paulo of vastly diverse ethnic background. They correspond to eleven gene-metabolite associations, six of which had been previously identified in the urine metabolome and three in the serum metabolome. Our key novel findings are the associations of two SNPs with NMR spectral signatures pointing to fucose (rs492602, P = 6.9×10(-44)) and lysine (rs8101881, P = 1.2×10(-33)), respectively. Fine-mapping of the first locus pinpointed the FUT2 gene, which encodes a fucosyltransferase enzyme and has previously been associated with Crohn's disease. This implicates fucose as a potential prognostic disease marker, for which there is already published evidence from a mouse model. The second SNP lies within the SLC7A9 gene, rare mutations of which have been linked to severe kidney damage. The replication of previous associations and our new discoveries demonstrate the potential of untargeted metabolomics GWAS to robustly identify molecular disease markers.

Reprint of: Musculoskeletal system in the old age and the demand for healthy ageing biomarkers.

Population ageing has emerged as a major demographic trend worldwide due to improved health and longevity. This global ageing phenomenon will have a major impact on health-care systems worldwide due to increased morbidity and greater needs for hospitalization/institutionalization. As the ageing population increases worldwide, there is an increasing awareness not only of increased longevity but also of the importance of "healthy ageing" and "quality of life". Yet, the age related chronic inflammation is believed to be pathogenic with regards to its contribution to frailty and degenerative disorders. In particular, the frailty syndrome is increasingly being considered as a key risk indicator of adverse health outcomes. In addition, elderly may be also prone to be resistant to anabolic stimuli which is likely a key factor in the loss of skeletal muscle mass with ageing. Vital to understand these key biological processes is the development of biological markers, through system biology approaches, aiding at strategies for tailored therapeutic and personalized nutritional program. Overall aim is to prevent or attenuate decline of key physiological functions required to live an active, independent life. This review focus on core indicators of health and functions in older adults, where nutrition and tailored personalized programs could exhibit preventive roles, and where the aid of metabolomics technologies are increasingly displaying potential in revealing key molecular mechanisms/targets linked to specific ageing and/or healthy ageing processes.

Serum profiling of healthy aging identifies phospho- and sphingolipid species as markers of human longevity.

As centenarians well represent the model of healthy aging, there are many important implications in revealing the underlying molecular mechanisms behind such successful aging. By combining NMR metabonomics and shot-gun lipidomics in serum we analyzed metabolome and lipidome composition of a group of centenarians with respect to elderly individuals. Specifically, NMR metabonomics profiling of serum revealed that centenarians are characterized by a metabolic phenotype distinct from that of elderly subjects, in particular regarding amino acids and lipid species. Shot- gun lipidomics approach displays unique changes in lipids biosynthesis in centenarians, with 41 differently abundant lipid species with respect to elderly subjects. These findings reveal phospho/sphingolipids as putative markers and biological modulators of healthy aging, in humans. Considering the particular actions of these metabolites, these data are suggestive of a better counteractive antioxidant capacity and a well-developed membrane lipid remodelling process in the healthy aging phenotype.

Impact of breast-feeding and high- and low-protein formula on the metabolism and growth of infants from overweight and obese mothers.

BACKGROUND: The combination of maternal obesity in early pregnancy and high protein intake in infant formula feeding might predispose to obesity risk in later life. METHODS: This study assesses the impact of breast- or formula-feeding (differing in protein content by 1.65 or 2.7 g/100 kcal) on the metabolism of term infants from overweight and obese mothers. From birth to 3 mo of age, infants received exclusively either breast- or starter formula-feeding and until 6 mo, exclusively either a formula designed for this study or breast-feeding. From 6 to 12 mo, infants received complementary weaning food. Metabonomics was conducted on the infants' urine and stool samples collected at the age of 3, 6, and 12 mo. RESULTS: Infant formula-feeding resulted in higher protein-derived short-chain fatty acids and amino acids in stools. Urine metabonomics revealed a relationship between bacterial processing of dietary proteins and host protein metabolism stimulated with increasing protein content in the formula. Moreover, formula-fed infants were metabolically different from breast-fed infants, at the level of lipid and energy metabolism (carnitines, ketone bodies, and Krebs cycle). CONCLUSION: Noninvasive urine and stool metabolic monitoring of responses to early nutrition provides relevant readouts to assess nutritional requirements for infants' growth.