Ontogeny Related Changes in the Pediatric Liver Metabolome

Introduction

Most drugs prescribed for children have not been studied in the relevant pediatric patient population to determine the appropriate dosing regimen, with ~20% of drugs approved by the FDA being labeled for use in children as of 2016 (from FDA website https://www.fda.gov/drugs/resourcesforyou/consumers/ucm143565.htm). Hence, the majority of drug dosing regimens for use in pediatric populations are based on dosing guidelines developed for adults and modified for use in children based on body weight and size of the child. However, children differ from adults in other ways beyond just size or weight, including body composition and organ development (1, 2). The disposition of many drugs is dependent on hepatic factors, such as blood flow and activity of drug-metabolizing enzymes and transporters.

In particular, the activity of many cytochrome P450 (CYP) isoforms, such as CYP3A4, CYP1A2, and CYP2D6, has been shown to be lower in newborns compared to adults, with each having a distinct developmental trajectory (3, 4). Similar findings have been reported for glucuronosyl transferases (5) and the expression of liver transporters (6). In addition to the critical need to understand the ontogeny of drug disposition to aid in age-appropriate dose selection, there needs to be recognition that a child is a dynamically changing biological system. For example, age-dependent changes in organic acid profiles imply that mitochondrial function may change during growth and development, especially relevant in the context of valproate hepatotoxicity (7). Less studied are the developmental trajectories of hepatic pathways that serve as targets of drug action, such as cholesterol biosynthesis (statins) and glucose homeostasis (metformin).

Technological advances in measuring metabolites and the rapid commercialization of novel instrumentation have sped up the adoption of metabolomics in all aspects of basic, population and clinical biomedical research (8, 9). Metabolites are the substrates, cofactors, and products needed for biological pathways and essential for cellular functions. In addition, many endogenous compounds are also substrates for “drug” metabolizing enzymes and transporters. Therefore, the purpose of this study was to characterize the biochemical changes occurring in liver between birth and 18 years of age and gain initial insight into ages/developmental stages that may be associated with altered drug response or increased susceptibility for age-related drug toxicity not apparent from adult data.

Materials and Methods

Liver Samples

Postmortem pediatric human liver tissue samples were obtained through the Brain and Tissue Bank for Developmental Disorders at the University of Maryland (Baltimore, MD), the Liver Tissue Cell Distribution System (LTCDS; University of Pittsburgh and University of Minnesota), and XenoTech LLC (Lenexa, KS). The use of these tissues was classified as non-human subject research by the Children's Mercy Hospital Pediatric Institutional Review Board. A replication set of post-mortem liver tissue samples from autopsies of fetuses (from therapeutic abortions or stillbirths) and infants was provided by the Erasmus Medical Center Tissue Bank, Sophia Children's Hospital, Rotterdam, the Netherlands. Tissue was procured at the time of autopsy within 24 h after death, snap-frozen in liquid nitrogen and stored at −80 °C for later research use. The Erasmus Medical Center Research Ethics Board waived the need for formal ethics approval according to the Dutch Law on Medical Research in Humans. Tissue was collected when parental written informed consent for both autopsy and the explicit use of the tissue for research was present. Samples were selected based on the absence of a clinical diagnosis or medications affecting the liver (CMH and Erasmus Medical Center), and tissue that was histologically normal (Erasmus Medical Center). Samples were stratified into four age groups: <1 year of age (age group 1), 1 to <6 years (age group 2), 6 to <12 years (age group 3), and 12–18 years of age (age group 4). In total 98 liver samples were available for metabolomic analysis. Characteristics of the study group are presented in Table 1.

TABLE 1

Table 1. Summary of pediatric liver samples included in the study.

Metabolomic Analysis

Untargeted metabolomic profiling was conducted by Metabolon Inc. (Durham, NC). Samples in two sets of experiments, Experiment 1 and Experiment 2, as described below and depicted in Figure 1A.

FIGURE 1

Figure 1. (A) Description of the two procedures used to obtain metabolite concentrations and how the samples were allocated in the three batches. (B) Flowchart of the entire analysis batch including how the quality control was conducted for the raw data, and summarizes the metabolite validation.

Experiment 1

The first experiment was completed using the first set of samples (N = 48) (referred to as “batch 1”).

Metabolite extraction and detection as previously described (10). Briefly, liver sample preparation was conducted using a proprietary series of organic and aqueous extractions to remove the protein fraction and optimize recovery of small molecules through the automated MicroLab STAR^® system (Hamilton Company, UT, USA), centrifuged, and the resulting supernatants were analyzed by analyzed by ultra-performance liquid chromatography mass spectrometry (UPLC-MS/MS) in a positive and negative ion mode (UPLC: Waters, Milford, MA; mass spectrometer: Thermo-Finnigan LTQ, Thermo Fisher Scientific, Waltham, MA, scan range, 80–1,000 m/z) and by GC-MS (Thermo-Finnigan Trace DSQ fast-scanning single-quadrupole mass spectrometer, scan range 50–750 m/z). The final experiment 1 metabolomic dataset comprised a total of 751 biochemicals, 478 compounds of known identity (named biochemicals) and 273 compounds of unknown structural identity. As initial statistical analysis revealed an age-dependent effect that could not be distinguished from a tissue source-related effect, a replication set of group 1 samples was obtained through collaboration with the Erasmus Medical Center/Sophia Children's Hospital.

Experiment 2

Given that the metabolomic platform changed between the first analysis and the sample set containing the replication samples, the second experiment examined the entire set of 98 samples. The same 48 samples previously processed in Experiment 1 and designated as “batch 1” above were re-analyzed on the new platform, with the results designated “batch 2.” The replication samples from the Erasmus Medical Center/Sophia Children's Hospital and additional samples from CMH (N = 50) are designated as “batch 3.” Following the sample extraction, the resulting extract was analyzed using a Waters ACQUITY UPLC and a Thermo Scientific Q-Exactive high resolution/accurate mass spectrometer interfaced with a heated electrospray ionization (HESI-II) source and Orbitrap mass analyzer operated at 35,000 mass resolution (11). Four methods were utilized: two separate reverse phase (RP)/UPLC-MS/MS methods with positive ion mode electrospray ionization (ESI), RP/UPLC-MS/MS with negative ion mode ESI, and HILIC/UPLC-MS/MS with negative ion mode ESI. The MS analysis alternated between MS and data-dependent MSⁿ scans using dynamic exclusion. The scan range varied slightly between methods but covered 70–1,000 m/z. The final experiment 2 metabolomic dataset comprised a total of 971 biochemicals, 779 compounds of known identity (named biochemicals) and 192 compounds of unknown structural identity.

Data Processing

Raw data was extracted, peak-identified and QC processed using proprietary hardware and software. Compounds were identified by comparison to library entries of purified standards or recurrent unknown entities (12). Peaks were quantified using area-under-the-curve. As the metabolomic assays span multiple days, a data normalization step was performed to correct variation resulting from instrument inter-day tuning differences. Essentially, each compound was corrected in run-day blocks by registering the medians to equal one (1.00) and normalizing each data point proportionately.

Statistical Analyses

Additional quality control measures were taken and are outlined in Figure 1B, such as imputation of missing data, and normalization. First, metabolites with more than 75% missing data were removed from subsequent analysis. After removing metabolites from analysis data set, data were log₂ transformed and scaled followed by imputation of missing values with the minimum detected level for the specific metabolite divided by $\sqrt{2}$. Latent effects are unmeasurable and unobservable factors can bias results, some examples of latent effects are batch effects, day-to-day variations in instrument performance. Latent factors were identified using principal component analysis (PCA), where the top 2 principal components (PCs) were adjusted for by using a linear model for each batch. Removing the top 2 PCs strikes a balance between removing latent technical factors, while the biologically relevant factors are still present. On the other hand it is possible that some of the estimated latent factors are not technical artifacts but rather represent true biology presented in the data.

It is unreasonable to assume that the abundance of each metabolite is normally distributed, hence the rank-based non-parametric Kruskal-Wallis (KW) test was performed, for each batch and metabolite individually, to determine if there were any differences in metabolite levels among the four age groups (AG1—children <1 year; AG2—children with age between 1 and 6 years; AG3—children with age between 6 and 12 years; AG4—children with age between 12 and 18 years). Linear associations between age group and metabolite levels were investigated using Spearman correlation test. To adjust for multiple testing a Bonferroni correction was applied. Statistical analysis was restricted to metabolites in common between the three batches of samples to enable assessment of both technical and biological validation. Biological validation was achieved by determining which metabolites had a significant p-value from KW test in both batches 2 and 3. This set of metabolites was then compared to the set of metabolites with a significant p-value from KW test (p < 0.05) in Batch 1 for technical validation. Pairwise differences were assessed for technically validated metabolite using Mann-Whitney tests.

The implication of an association between the concentrations of a single metabolite and age can be difficult to relate to important biological processes, while grouping metabolites into biologically meaningful sets can help understand an entire metabolite set. To determine association between age and a metabolite set we used the Gamma Method, a variation of Fisher's method (13). In the Gamma method, p-values from the metabolite-age association analysis were combined into a test statistic defined as ${\displaystyle {\sum }_{i=1}^k}h(p_i)$, where k represents the number of metabolites in a pathway, p_i represents the p-value for the association between the ith metabolite and age group, and h() represent the inverse gamma cumulative distribution function [i.e., ${G_{\omega ,1}}^{-1}(1-p_i)$]. This transformation gives more weight with p-values below a soft truncation threshold value determined by the shape parameter ω (13, 14). For this analysis we set the shape parameter to 0.0382, which gives more weight to p < 0.1. Due to the correlation between metabolites, the sum of transformed p-values does not necessarily follows a known distribution; therefore empirical p-values were computing using permutation methods, with p-values estimated from 100,000 permutations. Definition and mapping of metabolites to the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways and drug metabolite sets was completed using pathway mapping data from MetaboAnalyst (15). Metabolite pathways with at least 5 metabolites and coverage of at least 50% were included in the analysis, excluding pathways related to disease processes, resulting in 2 KEGG and 30 drug metabolite pathways being analyzed. Only metabolites that were present all three batches were considered for pathway analysis. Pathways were considered enriched if the empirical p-value derived using the Gamma method is <0.01 in all three batches.

Results

Metabolites with more than 75% missing data (across all samples) were removed from the analyses; 88 out of 751 metabolites were removed from consideration from batch 1, while 66 and 57 out of 971 metabolites were not considered in batches 2 and 3, respectively. After the removal of metabolites with >75% missing values, there were 322 metabolites in common between the 3 batches for comparative analysis (884 in common between batch 2 and batch 3). First, the effect due to age was removed from the data and then PCA was conducted to estimate latent features. Then, the first 2 principal component removed for each batch individually, with the resulting data presented in Figure 2. As Figure 2 illustrates, there is a large difference in global metabolites levels between the infant (AG1) and the early childhood/late childhood/adolescent age groups (AG2–AG4). The 2 principal components that were removed consisted of 29.83, 26.34, and 27.33% of the total variance in each batch, respectively. There were 16, 21, and 11 metabolites that were present in only AG1 in batch 1, 2, and 3, respectively. However, there was only one metabolite, estriol 3-sulfate, that was only present in AG1 for all batches.

FIGURE 2

Figure 2. Principal Component Analysis of each metabolite for each batch. The metabolomic liver profiles for children <1 year old clearly separate from the profiles for the children older than 1 year of age in PCA plot.

Assessment of the association of each of the 884 metabolites with age group in Batch 2 and Batch 3 utilizing the KW test detected 70 and 197 metabolites in batch 2 and batch 3, respectively (Bonferroni adjusted p < 0.05/884 = 5.7e-5) (Supplemental Tables 1, 2). Spearman correlation tests to determine linear trends between metabolites and the four age groups detected 95 and 286 metabolites in batch 2 and 3, respectively (p < 0.05/884). Under the null hypothesis that a metabolite is not associated with age and independent between tests (i.e., metabolites), we would have expected to have 0 metabolites detected at the 5.7e-5 significance level. Thus, there appears to be a significant departure from the null hypothesis. Forty-one of these metabolites were detected in both batch 2 and batch 3 and had similar direction of effect as measured by Spearman correlation (Rho), were thus considered biologically validated.

For technical validation, 25 of these 41 biologically-validated metabolites were also associated with age in batch 1 (p < 0.05), of which 24 are known metabolites (Supplemental Table 3). The 24 known metabolites are presented in Table 2, with a heatmap of the 25 metabolite levels presented in Figure 3 and Supplemental Figure 1. Of the 25 metabolites, 20 were found to have increasing levels over childhood development with the lowest levels present in the infant age group (AG1) and little difference in metabolite levels between the three older age groups. Additionally, pairwise comparison between the four age groups revealed that in many cases the distribution of metabolites levels varied significantly between infants and all other age groups (p < 0.001; Table 2). Of the 24 known metabolites found to be associated with age, 3 metabolites were in the sub-pathways “purine metabolism (hypo)xanthine/Inosine containing,” “pyrimidine metabolism, uracil containing,” and “long fatty acids.”

TABLE 2

Table 2. Twenty-four known metabolites out of the 25 biologically and technically replicated metabolites.

FIGURE 3

Figure 3. Heatmaps for the 25 biologically and technically validated metabolites for each of the three batches.

Analysis of the pathways/metabolite sets (removing pathways related to disease etiology) with more than 5 metabolites and at least 50% coverage of the pathway resulted in a total of 32 significant metabolite sets (p < 0.001 in all 3 batches) (Supplemental Table 4). Based on the use of 100,000 permutations, the smallest observable empirical p < 1/100,000 = 1.0 × 10⁻⁵. The most significant pathway in all three batches was the “Betaine Metabolism” pathway (Table 3). For this pathway, 12 of the 21 metabolites in this pathway were present in our study, resulting in a p < 1.0 × 10⁻⁵ all three batches. The 12 metabolites included in our analysis are: 5-methyltetrahydrofolate (5MeTHF), adenosine, betaine, choline, dimethylglycine, flavin adenine dinucleotide (FAD), homocysteine, methionine, nicotinamide adenine dinucleotide (NAD+), phosphate, S-adenosylhomocysteine (SAH), S-adenosylmethionine (SAM). The majority of these metabolites were individually significant associated with age. To put these results in reference to the biological pathway, we have plotted them with the pathway in Figure 4A and have denoted the metabolites measured in our study and the results for association with age. The other KEGG pathway associated with age was alpha-linolenic Acid and linoleic acid metabolism (Table 3 and Figure 4B), for which we observed 5 metabolites with p < 0.05. These 5 metabolites (arachidonate (20:4n6), docosapentaenoate (n3 DPA; 22:5n3), linoleate (18:2n6), linolenate [alpha or gamma; (18:3n3 or 6)]) were all observed to be increasing in abundance age. For 30 of the 32 metabolism and/or drug related pathways (Supplemental Table 4), the same 10 metabolites (adenosine monophosphate, L-alanine, L-methionine, L-leucine, L-histidine, L-proline, L-asparagine, L-valine, L-threonine, L-soleucine) were observed in the pathways resulting in the same pathway level results for batch 1, batch 2 and batch 3 of p = 0.00011, 0.00017, and <1.0 × 10⁻⁵, respectively (Table 3 and Supplemental Figure 2).

TABLE 3

Table 3. Summary of results for the measured metabolites in the each of the ontogeny related pathways.

FIGURE 4

Figure 4. (A) Betaine pathway and (B) alpha linolenic acid and Linoleic acid metabolism. The results for metabolites measured in this study are indicated with red indicating the metabolite was increasing with age and blue indicating the metabolite was decreasing in abundance with age. Metabolites in the pathway not measured are indicated in black. The shape corresponds to the significance of the Spearman correlation of metabolite's concentration and age, where a square indicates that the p < 0.05 and a circle shows that the p > 0.05.

Discussion

This study set out to investigate age-specific patterns of cellular metabolites in liver between birth and 18 years of age that may be indicative of altered drug response or susceptibility to drug toxicity unique to pediatric patient populations. Using an untargeted mass-spectroscopy based metabolomics analysis; we assessed hundreds of metabolites in a set of pediatric liver samples, with both technical and biological validation built into the study design (Figure 1A). A major challenge for investigations like the one we have conducted is access to a large number of tissue samples of sufficiently high quality to obtain interpretable results; differences in retrieval and preservation methods across various publicly funded and commercial providers has considerable potential to not only contribute to variability in sample quality (16), but may also confound data interpretation. For example, in the United States, availability of tissue samples from infants <1 year of age is relatively infrequent, and those that are available for research tend to be procured from the NIH-funded University of Maryland Brain and Tissue Bank for Developmental Disorders. This period of time also tends to be a period of rapid growth and can be accompanied by developmental patterns of gene expression as we noted in a previous study involving RNA-Seq analysis of a similar set of samples (17). Thus, it is difficult to differentiate true developmental differences in RNA or metabolite expression from confounding effects related to tissue source. Given that initial statistical analysis of the metabolomic data corroborated the source effect in the RNA-Seq data (17), the current investigation included samples <1 year of age from a second, independent source.

This experimental approach detected 24 metabolites that changed in a set of liver samples spanning birth to late adolescence (18 years of age; Table 2). Of the 24 metabolites, 21 were found to have increasing levels over childhood development with the lowest levels present in the infant age group (children <1 year of age) and little difference in metabolite levels between the three older age groups (children >1 year of age) (Figure 3 and Supplemental Figure 1). Key sub-pathways were observed in the list of significant known metabolites, including “purine metabolism (hypo)xanthine/Inosine containing,” “pyrimidine metabolism, uracil containing,” and “long fatty acids”—all three replicated metabolites in long fatty acids sub-pathway increasing in abundance throughout childhood development (Table 2). Lastly, we observed consistent age-related differences in KEGG pathways “Betaine Metabolism” (pathway p-value of < 1.0 × 10⁻⁵ for all 3 batches) and “Alpha-Linolenic Acid and Linoleic Acid Metabolism” (pathway p-value of 0.002, 0.007, and < 1.0 × 10⁻⁵ for batch 1, batch 2 and batch 3, respectively). A limitation to our analysis is that metabolites may not be considered simply because their concentration is natural concentration is near the lower limit of detection. There were 353 metabolites that were measured in all three batches and of these 31 metabolites were removed because there was 75% missing in at least one batch (Supplemental Figure 3). Interestingly, X – 11,795, glycerol 2–phosphate, gamma–glutamylhistidine, formiminoglutamate, and disulfide, benzoate, may have been removed due to the different technologies being used in Batch 1, and Batches 2 and 3. On the other hand, many of the other metabolites have a consistently high percentage of missing observations across all batches (Supplemental Figure 4).

Individual analytes and those assigned to specific pathways may be derived from both endogenous as well as exogenous (dietary) sources. For example, the betaine metabolism pathway contained 12 measured metabolites, with all but two metabolites found to be related to ontogeny (p < 0.10). Eight out of the 10 ontogeny related metabolites increased in abundance increased with age (Table 3 and Figure 4A). In humans, betaine itself is found in food and can also be formed endogenously from choline. Betaine is vital in transmethylation and provides control of hepatocellular hydration and provides protection of the liver from various forms of stress, including osmotic stress (18). Similarly, within the “Alpha-Linolenic Acid and Linoleic Acid Metabolism” pathway, 4 metabolites consistently were observed to increase with increasing age across the batches with p < 0.10: docosapentaenoic acid, docosapentaenoic acid (22n-6), linoleic acid, gamma-linolenic acid and alpha-linolenic acid, whereas arachidonic acid was present at higher concentrations in children <1 year of age relative to the older age groups (Table 3 and Figure 4B). Alpha-linolenic acid (ALA) is an essential omega-3 fatty acid found in many nuts and vegetable oils, and thus it is difficult to differentiate analyte changes due to biological maturation with external factors that also change with increasing age.

One obvious difference between infants <1 year of age and older children/adolescents is diet, with breast milk and formula, either cow's milk- or soy-based formulas, representing the primary source of nutrition until solid oral foods are introduced later in the first year of life. However, it is also now well-recognized that the intestinal microbiota of newborns is different from older children and adults, being established after birth and influenced not only by mode-of delivery, vaginal vs. cesarean, but also to a considerable extent by the mode of feeding. For example, several studies now report higher levels of Lactobacilli and bifidobacteria in stool from breast-fed infants compared to stool from formula-fed infants that is dominated by a more diverse variety of species, Bacteroides, Clostridia, Staphylococci, enterobacteria, Enterococci, and Atopobium (19–21). Introduction of solid foods also results in increased diversity of the intestinal microbiota over time. Thus, it is not unreasonable to expect that concentrations of dietary constituents and the products of metabolism of those constituents by a changing intestinal microbiota might result in a unique infant metabolome.

Developmental differences in the gut microbiome between infants <1 year of age and older children, or more specifically establishment of the gut microbiome after birth, would be expected to result in changes in bile acid composition. For example, primary bile acids derived from endogenous liver metabolism would be expected to be higher and the concentrations of gut microbiome-derived secondary bile acids to be lower in the younger group, with a shift toward accumulation of predominantly secondary bile acids in the older age groups following maturation of the intestinal microbiome. Higher concentrations of primary bile acids in premature infants and term newborns compared to older ages is a well-known phenomenon (22, 23), and our observation of lower concentrations of primary and secondary bile acids taurodeoxycholate and glycolithocholate sulfate in the <1 year old group is consistent with colonization and maturation of the gut microbiome over the first year of life.

Studies in which germ-free mice are transferred to a conventional environment simulate the change in environment experienced during birth when newborns enter their extrauterine environment. Claus et al. (24) applied metabolomic approaches to characterize the metabolic adaptation to bacterial colonization of the gut by transferring germ-free mice to conventional environment in which bedding used by conventional mice was provided to expose the germ-free animals to the same microbial ecosystem. They observed rapid increases in weight over the first few days after colonization as well as increases in gluconeogenesis that were followed by increases in hepatic triglyceride synthesis and alterations in bile acid metabolites. It is interesting to note that five analytes in our dataset that were differentially present between the <1 year old and older age groups, heme, taurodexycholate, 10-heptadecenoate (17:1n7), 10-non-adecenoate (19:1n9), and margarate (17:0), were also observed to change by Claus et al., supporting the hypothesis that some of the observed “developmental” changes in metabolome are a consequence of bacterial colonization of the gut after birth.

In conclusion, this limited dataset illustrates the complexity underlying observed “developmental” changes in the hepatic metabolome as assessed by the snapshot of developmental changes provided by metabolomic analysis of liver tissue. The first year of life is characterized by a velocity of change in height and weight that exceeds that observed in older children, and thus, some changes in the metabolome reflect the processes of growth and development. However, it is apparent that other factors, such as change in diet and gut microbiome with increasing age, also contribute to the overall picture. The interplay between diet, gut microbiome and hepatic physiology can be expected to underlie developmental changes in expression of drug metabolizing enzymes and transporters involved in the absorption, distribution, metabolism and excretion of drugs, nutrients, toxicants, and other foreign compounds. For example, bile salts are essential for absorption of lipid soluble vitamins and poorly water-soluble medications. Furthermore, expansion of the bile acid pool is also accompanied by changes in expression of hepatocellular and biliary uptake and efflux transporters (25) and likely contributes to the developmental trajectories of transporters that is of considerable interest for modeling and simulation of drug disposition in pediatrics (26). Studies with germ-free mice introduced into a conventional environment also resulted in increased expression of drug metabolizing enzymes, such as CYP2C29 and CYP3A11 that are similar to the developmental trajectories of important drug biotransformation pathways in humans, such as CYP2C9 (27) and CYP3A4 (28). Thus, improved understanding of the factors contributing to the developmental processes governing the role of the liver and its interaction with other systems will improve our understanding of not only drug disposition in a vulnerable, understudied patient population, but will also aid in identification of factors that influence drug disposition and response in individual children.

Data Availability Statement

The normalized and imputed data used in the analysis along with the analysis results can be found at http://explorerpedpgx.moffitt.org using Google Chrome. The data after QC is also available for download on the Metabolomic workbench, https://www.metabolomicsworkbench.org/. These data can be found with the following study numbers are ST001402, ST001403, and ST001404 corresponding to Batch 1, 2, 3 respectively.

Ethics Statement

The studies involving human participants were reviewed and approved by the use of these tissues was classified as non-human subject research by the Children's Mercy Hospital Pediatric Institutional Review Board. A replication set of post-mortem liver tissue samples from autopsies of fetuses (from therapeutic abortions or stillbirths) and infants was provided by the Erasmus Medical Center Tissue Bank, Sophia Children's Hospital, Rotterdam, the Netherlands. The Erasmus Medical Center Research Ethics Board waived the need for formal ethics approval according to the Dutch Law on Medical Research in Humans. Tissue was collected when parental written informed consent for both autopsy and the explicit use of the tissue for research was present. Written informed consent to participate in this study was provided by the participants' legal guardian/next of kin.

Author Contributions

SW and JL conceived and planned the experiments. CW, QL, CB, and BF completed the statistical and bioinformatics analysis. JL and BF supervised the project. RG, SW, JL, and CB completed experiments. CW, JL, and BF wrote the manuscript. All authors reviewed and edited the manuscript.

Funding

This research was provided by NIH/NICHD U54 HD090258, Children's Mercy Hospital, and Moffitt Cancer Center.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

Liver tissue was obtained from the University of Maryland Brain and Tissue Bank for Developmental Disorders (Baltimore, MD; funded by NIH Contract HHSN275200900011C, Reference No. N01-HD-9- 0011) and the Liver Tissue Cell Distribution System (funded by NIH Contract No. N01-DK-7-0004/HHSN267200700004C). We thank the Department of Pathology at Erasmus Medical Center for sample collection for the Erasmus Medical Center Tissue Biobank.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fped.2020.00549/full#supplementary-material

Supplemental Figure 1. Side by side boxplots of metabolite concentration and age for the 25 biologically and technically validated metabolites for each of the three batches.

Supplemental Figure 2. Drug action pathways. The results for metabolites measured in this study are indicated with red indicating the metabolite was increasing with age and blue indicating the metabolite was decreasing in abundance with age. Metabolites in the pathway not measured are indicated in black. The shape corresponds to the significance of the Spearman correlation of metabolite's concentration and age, where a square indicates that the p < 0.05 and a circle shows that the p > 0.05. For the drug action pathway, the one individual metabolite that is unique to a drug action pathway is denoted by a purple trapezoid.

Supplemental Figure 3. Venn diagram of the number of metabolites removed from analysis, by batch, because 75% of the observations were missing in at least one batch (31 metabolites).

Supplemental Figure 4. Heatmap of the percent of observations missing for metabolites that were eliminated from analysis due to more than 75% of observations were missing for at most 2 batches (20 metabolites).

Supplemental Table 1. Results for all 905 metabolites in batch 2. N infant = 13, N early childhood = 14, N late childhood = 10, N adolescent = 11.

Supplemental Table 2. Results for all 914 metabolites in batch 3. N infant = 21, N early childhood = 7, N late childhood = 11, N adolescent = 11.

Supplemental Table 3. Results for all 663 metabolites in batch 1. N infant = 13, N early childhood = 14, N late childhood = 10, N adolescent = 11.

Supplemental Table 4. Results for all 32 KEGG and drug pathways with at least 5 metabolites and pathway coverage of at least 50%.

References

1. Kearns GL, Abdel-Rahman SM, Alander SW, Blowey DL, Leeder JS, Kauffman RE. Developmental pharmacology–drug disposition, action, and therapy in infants and children. N Engl J Med. (2003) 349:1157–67. doi: 10.1056/NEJMra035092

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Rodieux F, Wilbaux M, Van Den Anker JN, Pfister M. Effect of kidney function on drug kinetics and dosing in neonates, infants, and children. Clin Pharmacokinet. (2015) 54:1183–204. doi: 10.1007/s40262-015-0298-7

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Hines RN. The ontogeny of drug metabolism enzymes and implications for adverse drug events. Pharmacol Ther. (2008) 118:250–67. doi: 10.1016/j.pharmthera.2008.02.005

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Hines RN. Developmental expression of drug metabolizing enzymes: impact on disposition in neonates and young children. Int J Pharm. (2013) 452:3–7. doi: 10.1016/j.ijpharm.2012.05.079

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Bhatt DK, Mehrotra A, Gaedigk A, Chapa R, Basit A, Zhang H, et al. Age- and genotype-dependent variability in the protein abundance and activity of six major uridine diphosphate-glucuronosyltransferases in human liver. Clin Pharmacol Ther. (2018) 105:131–41. doi: 10.1002/cpt.1109

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Prasad B, Gaedigk A, Vrana M, Gaedigk R, Leeder JS, Salphati L, et al. Ontogeny of hepatic drug transporters as quantified by LC-MS/MS proteomics. Clin Pharmacol Ther. (2016) 100:362–70. doi: 10.1002/cpt.409

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Price KE, Pearce RE, Garg UC, Heese BA, Smith LD, Sullivan JE, et al. Effects of valproic acid on organic acid metabolism in children: a metabolic profiling study. Clin Pharmacol Ther. (2011) 89:867–74. doi: 10.1038/clpt.2011.47

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Kell DB, Oliver SG. The metabolome 18 years on: a concept comes of age. Metabolomics. (2016) 12:148. doi: 10.1007/s11306-016-1108-4

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Azad RK, Shulaev V. Metabolomics technology and bioinformatics for precision medicine. Brief Bioinform. (2019) 20:1957–71. doi: 10.1093/bib/bbx170

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Evans AM, Dehaven CD, Barrett T, Mitchell M, Milgram E. 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. Anal Chem. (2009) 81:6656–67. doi: 10.1021/ac901536h

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Evans AM, Bridgewater BR, Liu Q, Mitchell MW, Robinson RJ, Dai H, et al. High resolution mass spectrometry improves data quantity and quality as compared to unit mass resolution mass spectrometry in high- throughput profiling metabolomics. Metabolomics Open Access. (2014) 4:1–7. doi: 10.4172/2153-0769.1000132

CrossRef Full Text | Google Scholar

12. Dehaven CD, Evans AM, Dai H, Lawton KA. Organization of GC/MS and LC/MS metabolomics data into chemical libraries. J Cheminform. (2010) 2:9. doi: 10.1186/1758-2946-2-9

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Zaykin DV, Zhivotovsky LA, Czika W, Shao S, Wolfinger RD. Combining p-values in large-scale genomics experiments. Pharm Stat. (2007) 6:217–26. doi: 10.1002/pst.304

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Biernacka JM, Jenkins GD, Wang L, Moyer AM, Fridley BL. Use of the gamma method for self-contained gene-set analysis of SNP data. Eur J Hum Genet. (2012) 20:565–71. doi: 10.1038/ejhg.2011.236

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Xia J, Sinelnikov IV, Han B, Wishart DS. MetaboAnalyst 3.0–making metabolomics more meaningful. Nucleic Acids Res. (2015) 43:W251–7. doi: 10.1093/nar/gkv380

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Leeder JS, Meibohm B. Challenges and opportunities for increasing the knowledge base related to drug biotransformation and pharmacokinetics during growth and development. Drug Metab Dispos. (2016) 44:916–23. doi: 10.1124/dmd.116.071159

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Meier R, Bi C, Gaedigk R, Heruth DP, Ye SQ, Leeder JS, et al. Ontogeny-related pharmacogene changes in the pediatric liver transcriptome. Pharmacogenet Genomics. (2018) 28:86–94. doi: 10.1097/FPC.0000000000000326

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Kempson SA, Vovor-Dassu K, Day C. Betaine transport in kidney and liver: use of betaine in liver injury. Cell Physiol Biochem. (2013) 32:32–40. doi: 10.1159/000356622

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Hill CJ, Lynch DB, Murphy K, Ulaszewska M, Jeffery IB, O'shea CA, et al. Evolution of gut microbiota composition from birth to 24 weeks in the INFANTMET cohort. Microbiome. (2017) 5:4. doi: 10.1186/s40168-016-0213-y

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Stewart CJ, Ajami NJ, O'brien JL, Hutchinson DS, Smith DP, Wong MC, et al. Temporal development of the gut microbiome in early childhood from the TEDDY study. Nature. (2018) 562:583–8. doi: 10.1038/s41586-018-0617-x

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Zhuang L, Chen H, Zhang S, Zhuang J, Li Q, Feng Z. Intestinal microbiota in early life and its implications on childhood health. Genomics Proteomics Bioinformatics. (2019) 17:13–25. doi: 10.1016/j.gpb.2018.10.002

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Suchy FJ, Balistreri WF, Heubi JE, Searcy JE, Levin RS. Physiologic cholestasis: elevation of the primary serum bile acid concentrations in normal infants. Gastroenterology. (1981) 80:1037–41. doi: 10.1016/0016-5085(81)90078-0

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Abuhelwa AY, Foster DJR, Upton RN. A quantitative review and meta-models of the variability and factors affecting oral drug absorption-part I: gastrointestinal pH. AAPS J. (2016) 18:1309–21. doi: 10.1208/s12248-016-9952-8

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Claus SP, Ellero SL, Berger B, Krause L, Bruttin A, Molina J, et al. Colonization-induced host-gut microbial metabolic interaction. mBio. (2011) 2:e00271-10. doi: 10.1128/mBio.00271-10

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Suchy FJ, Bucuvalas JC, Novak DA. Determinants of bile formation during development: ontogeny of hepatic bile acid metabolism and transport. Semin Liver Dis. (1987) 7:77–84. doi: 10.1055/s-2008-1040567

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Elmorsi Y, Barber J, Rostami-Hodjegan A. Ontogeny of hepatic drug transporters and relevance to drugs used in pediatrics. Drug Metab Dispos. (2016) 44:992–8. doi: 10.1124/dmd.115.067801

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Koukouritaki SB, Manro JR, Marsh SA, Stevens JC, Rettie AE, Mccarver DG, et al. Developmental expression of human hepatic CYP2C9 and CYP2C19. J Pharmacol Exp Ther. (2004) 308:965–74. doi: 10.1124/jpet.103.060137

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Stevens JC, Hines RN, Gu C, Koukouritaki SB, Manro JR, Tandler PJ, et al. Developmental expression of the major human hepatic CYP3A enzymes. J Pharmacol Exp Ther. (2003) 307:573–82. doi: 10.1124/jpet.103.054841

PubMed Abstract | CrossRef Full Text | Google Scholar