Gene and protein expression and metabolic flux analysis reveals metabolic scaling in liver ex vivo and in vivo
Abstract
Metabolic scaling, the inverse correlation of metabolic rates to body mass, has been appreciated for more than 80 years. Studies of metabolic scaling have largely been restricted to mathematical modeling of caloric intake and oxygen consumption, and mostly rely on computational modeling. The possibility that other metabolic processes scale with body size has not been comprehensively studied. To address this gap in knowledge, we employed a systems approach including transcriptomics, proteomics, and measurement of in vitro and in vivo metabolic fluxes. Gene expression in livers of five species spanning a 30,000-fold range in mass revealed differential expression according to body mass of genes related to cytosolic and mitochondrial metabolic processes, and to detoxication of oxidative damage. To determine whether flux through key metabolic pathways is ordered inversely to body size, we applied stable isotope tracer methodology to study multiple cellular compartments, tissues, and species. Comparing C57BL/6 J mice with Sprague-Dawley rats, we demonstrate that while ordering of metabolic fluxes is not observed in in vitro cell-autonomous settings, it is present in liver slices and in vivo. Together, these data reveal that metabolic scaling extends beyond oxygen consumption to other aspects of metabolism, and is regulated at the level of gene and protein expression, enzyme activity, and substrate supply.
Editor's evaluation
Key metabolic processes have been shown to scale inversely with the body mass of different animals. This study provides direct evidence for metabolic scaling of key metabolic fluxes in the livers of mice and rats, as well as species-specific differences in the transcription and expression of enzymes involved in energy metabolism that could contribute to metabolic scaling. The finding suggests that metabolic scaling likely reflects multiple levels of regulation and have broad implications for studying animal metabolism and physiology.
https://doi.org/10.7554/eLife.78335.sa0Introduction
In 1932, Max Kleiber published a seminal study (Kleiber, 1932), integrating prior reports demonstrating a phenomenon that came to be termed ‘Kleiber’s law,’ or the principle of metabolic scaling. Metabolic scaling refers to the phenomenon that the metabolic processes in many animals, if not all, scale inversely to three-quarters of their body mass (West et al., 1997). In simpler terms, there is a reduction in metabolic rate as body size increases. For example, an elephant is 25 million times larger than a fruit fly, yet its energy expenditure is only 20 thousand times higher; thus, from the fruit fly to elephant, the metabolic rate per gram of body weight scales down 1250 times. While there is experimental evidence for metabolic scaling from bacteria to large mammals, data have been generated almost exclusively from observations of caloric intake and oxygen consumption, with gene and protein expression, and substrate fluxes almost entirely unexplored.
The concept of hierarchical regulation, whereby gene expression initiates the cascade that allows for the flux of metabolic pathways (Rossell et al., 2005; Suarez and Moyes, 2012), provides a systems framework to begin to understand scaling. Beginning at the transcriptional level, we studied liver gene expression across five species: mice (Mus musculus), rats (Rattus norvegicus), monkeys (Macaca mulatta), humans (Homo sapiens), and cattle (bos taurus), species with a 30,000-fold range of average body weight in adults (from 30 g in mice, to 900 kg in cattle). Numerous metabolic genes related to glycolysis, gluconeogenesis, fatty acid metabolism, oxygen consumption, electron transport, and redox function, and detoxification of oxidative damage, were expressed at levels inverse to body size. Further analysis of liver proteomics revealed that approximately half of the genes in the liver that were expressed inversely proportionally to body size at the transcriptional level, were also expressed at levels inversely proportional to body size at the level of protein expression. To determine if gene and protein expression would correlate with enzyme activity and metabolic flux, we performed a comprehensive assessment of liver metabolism in vivo and in vitro using modified Positional Isotopomer NMR Tracer Analysis (PINTA) (Perry et al., 2017b) and stable isotope-derived turnover (Perry et al., 2015) methods. Our analysis shows that rats exhibit lower metabolic rates when compared to mice in and ex vivo; however, no significant differences were observed when we isolated hepatocytes and cultured them in vitro under identical conditions. Taken together, this study demonstrates the variation of metabolic fluxes according to body size, extending prior studies of metabolic scaling, and provides unique insight into the regulation of metabolic flux across species.
Results
Genes within the liver that are expressed inversely proportional to body weight are predominantly metabolic genes
We examined gene expression in livers from mice (Mus musculus), rats (Rattus norvegicus), monkeys (Macaca mulatta), humans (Homo sapiens), and cattle (Bos taurus). Using recent advances in high throughput mRNA sequencing and bioinformatics tools that allow for intra-species data preprocessing (Bray et al., 2016; Conesa et al., 2016; Ritchie et al., 2015), we searched for a set of genes in the liver, the metabolic hub of mammals, whose expression correlates inversely with body mass. After normalizing for differences in transcript length and abundance across species, we filtered out genes that followed the pattern of mouse >rat > monkey >human > cow. The genes that met these criteria were predominantly related to metabolic pathways, including pyruvate metabolism, amino acid metabolism, and glucose metabolism (Figure 1A). Genes from this list were further restricted to genes involved in amino acid, carbohydrate, energy, lipid, vitamin, and TCA cycle metabolism, and demonstrated a range of degrees of inverse correlation with body mass, with only TCA cycle genes clustering together (Figure 1B).
Genes encoding enzymes involved in hepatic metabolism are expressed inversely proportionally to body mass and involve metabolite detoxification, intertissue metabolism, substrate metabolism, electron transport, and NAD metabolism
In order to further understand the functional aspects of the metabolic genes that are expressed inversely proportionally to body size, the gene list from Figure 1B was categorized into several functional categories, converging on optimizing energy provision, oxidative metabolism, and damage control from oxidative stress and ammonia (Figure 2). Furthermore, to understand whether or not certain genes that are expressed inversely proportionally to body size involved anabolic or catabolic processes, they were further classified by their properties to be energy suppliers or consumers. Eleven of sixteen critical metabolic enzymes that scaled required molecular oxygen, NAD+/NADH, or ATP/ADP for function, possibly indicating exquisite regulation of energy-consuming processes at the individual gene level. Genes involved in the detoxication of lipid peroxidation-derived aldehydes (ALDH7A1), hydrogen peroxide (GPX1), and ammonia (GLUL) suggest scaling of damage control mechanisms that are associated with increased oxidative metabolism across species (Figure 2A). The inverse correlation between body size and expression of genes that are associated with interorgan crosstalk is consistent with scaling in vivo which would not be expected in plated cells. For example, the differentially expressed genes include GPT1, which is involved in recycling skeletal muscle-derived alanine back to liver-derived glucose (Felig and Wahren, 1971; Petersen et al., 2019), and the adiponectin receptor (ADIPOR1), which binds an adipose tissue-derived hormone that regulates gluconeogenesis and fatty acid oxidation (Lin and Accili, 2011; Li et al., 2020 Figure 2B). Genes involved in fatty acid metabolism included the rate-limiting steps of the synthesis of CoA (PANK1), of de novo fatty acid synthesis (ACACA), and of fatty acid elongation (ELOVL2), in addition to the oxidation of diacylglycerols (LIPE) (Figure 2C). NAD and ATP-dependent genes involved in glycolysis (PDK1), fructose/glucose metabolism (SORD1), and DLST of the TCA cycle also correlated inversely with body size (Figure 2D–E). Differentially regulated genes also couple oxygen consumption to NAD provision (MDH1, TDO2), and are involved with the function of the electron transport chain (subunits of complex I, NDUFA9, and complex IV, COX8A, which catalyzes oxygen accepting the final electrons of the electron transport chain) (Figure 2F–G).
To examine the possibility that the inverse correlation between body mass and gene expression observed in the transcriptomics analysis could be a consequence of global alterations in mRNA (for example, as a consequence of alterations in RNA turnover rates), we performed targeted quantitative polymerase chain reaction (qPCR), measuring in liver tissue abundance of mRNA encoding several enzymes that were found to scale in the five-species transcriptomics analysis, relative to the common housekeeping gene β-actin (Actb). We found that all three enzymes (Glul, Lipe, and Dlst) scaled relative to Actb (Figure 2—figure supplement 1A–C), whereas structural genes (collagenase 3 [Mmp3] and Larp1) did not (Figure 2—figure supplement 1D–E), indicating that the differences in metabolic gene expression observed across species is likely not a result of global changes in RNA levels.
In addition to transcriptomics, we assessed proteomics data to evaluate the protein levels corresponding to the genes that were found to be expressed inversely proportionally to body size at the level of mRNA expression. Our proteomics data were limited to mouse, rat, and human, as all the open-source proteomic databases that we identified lacked data from monkey or cow. An important limitation for finding such data is that even with careful post-processing, we cannot combine data from different studies, because differences in methods of tissue preparation may influence results. Therefore, we were limited to a single experiment that had generated proteomics data for mouse, rat, and human using the same experimental procedures. The dataset contained protein expression corresponding to eight of the twenty genes identified to scale in our transcriptomics data analysis. Of these, three (GLUL, GPX1, and MDH1) were found to follow a reverse correlation with body size (Figure 3A–C). Interestingly, one of these proteins (GLUL) was also found to be expressed inversely proportionally to body size in the left ventricle of the heart (Figure 3D). Additionally, we measured liver transaminase concentrations and observed that both alanine aminotransferase (ALT) and aspartate aminotransferase (AST) exhibited lower concentrations in humans as compared to rats and rats as compared to mice (Figure 3E–F), consistent with scaling at the level of protein expression as well as mRNA expression. Finally, we utilized established enzymatic assays to measure the activity of peroxidase and pyruvate carboxylase in the livers of mice and rats. 30–40% lower activity of each enzyme per mg tissue was observed in rats as compared to mice (Figure 3—figure supplement 1A–B), suggesting scaling at the level of metabolic enzyme activity.
Metabolic rates of mouse vs. rat hepatocytes in vitro are not significantly different
Considering prior data reporting higher oxygen consumption per unit body mass in smaller as compared to larger animals (Gilman et al., 2013; Brody, 1945; Urbina and Glover, 2013), we first asked whether these differences were cell-intrinsic, or whether in vivo or hepatocyte-extrinsic signals are required. We incubated plated hepatocytes in [3-13C] lactate and first validated that the data met the assumptions of PINTA, including reaching steady-state in [13C] lactate and glucose enrichment, and producing glucose at a linear rate throughout the 6 hr incubation (Figure 4—figure supplement 1A–C). Consistent with the possibility that hepatocyte-extrinsic signals are primarily responsible for metabolic scaling, when we used PINTA to assess cytosolic and mitochondrial fluxes, we observed no significant differences between species in any of the fluxes measured in plated hepatocytes: glucose production, VPC, VCS, the contribution of glucose or fatty acids to the tricarboxylic acid (TCA) cycle, or lipolysis (Figure 4A–H, Figure 4—figure supplement 1D–F). Similarly, a mitochondrial stress test in plated hepatocytes revealed no difference in any parameter: neither basal mitochondrial and non-mitochondrial respiration, ATP production, maximal (uncoupled) respiration, spare respiratory capacity, nor proton leak differed between plated hepatocytes from mice and rats (Figure 4I). Previous studies have demonstrated scaling in vitro in cell suspensions only when analyzed immediately after hepatocyte isolation (Porter and Brand, 1995), and have suggested that the phenomenon of scaling gradually disappears around 24 hr post removal (Brown et al., 2007), similar to the conditions in which we performed these studies. Most prior in vitro studies have also demonstrated an absence of scaling, in contrast to in vivo (Glazier, 2015), and we extend these results to gluconeogenic and lipolytic fluxes in hepatocytes, glucose production in liver slices, and multimodal flux analysis in vivo.
Glucose production per gram tissue is higher ex vivo in liver slices from mice than in rats
Next, considering that hepatocytes comprise approximately 70–80% of liver mass and that their culture in vitro does not replicate in vivo conditions (Krebs, 1950), we asked whether glucose production would be different between mice and rats in slices of liver. Indeed, we found that liver glucose production per gram liver mass was threefold greater in mouse liver slices as compared to rats (Figure 5A–B), suggesting that hepatocyte-extrinsic signals (for example, from other liver cell types) are involved in liver metabolic scaling.
Metabolic rates in multiple tissue types are higher in vivo in mice relative to rats
We utilized multimodal stable isotope metabolic flux analysis to compare rats and mice with respect to a panel of metabolic fluxes (Figure 6A). First, we validated tracer assumptions in vivo, including the metabolic and isotopic steady state in plasma and negligible liver glycogen concentrations, although in the recently fed state, hepatic glycogenolysis was higher in mice than that in rats (Figure 6—figure supplement 1A–G). Using PINTA, we found that both endogenous glucose production and gluconeogenesis from pyruvate (VPC) per gram liver were more than twofold higher in mice than rats (Figure 6B–C), although the fractional contribution of pyruvate to gluconeogenesis did not differ between mice and rats (Figure 6—figure supplement 1H). Mitochondrial oxidation scaled similarly, increasing threefold in mice as compared to rats studied under the same conditions, due to increases in both glucose oxidation (pyruvate dehydrogenase flux, VPDH) and fatty acid oxidation (Figure 6D–F), associated with an increase in the ratio of pyruvate carboxylase anaplerosis to citrate synthase flux (VPC/VCS) without any difference in the fraction of VCS flux fueled by glucose through PDH (Figure 6—figure supplement 1I–J). While we did not have the capacity to measure liver fluxes in larger mammals in the current study, endogenous glucose production, VPC, and VCS previously measured using PINTA were 50–60% lower in overnight fasted humans than in rats (Petersen et al., 2019), assuming a liver size of 1500 g in humans. These differences in metabolic fluxes according to body size applied not only to rodent liver metabolism but also to adipose tissue metabolism: whole-body fatty acid turnover, reflecting lipolysis, was 2.5-fold higher in mice than in rats (Figure 6G). No sex differences were observed in any of the measured fluxes (Figure 6—figure supplement 1K–P). Taken together, these data emphasize the inadequacy of common in vitro methods as a readout of in vivo metabolism: whereas in vivo mitochondrial oxidation (TCA cycle flux) was threefold higher in mice than in rats, in vitro measurements of oxygen consumption throughout a mitochondrial stress test, TCA cycle flux, and glucose production were not different between the species (Figure 4).
Classical clustering defined species-specific clusters based on in vivo metabolic fluxes but not in vitro fluxes
A clustering dendrogram was applied to our in vitro flux data and showed no distinct clustering between species (Figure 7A–B). However, the in vivo metabolic flux data led to distinct clustering of rats and mice (Figure 7C), providing a classical clustering-based objective analysis of in vitro versus in vivo metabolic flux.
Discussion
Oxygen consumption has been shown to scale inversely with body mass in species ranging in mass across 20 orders of magnitude, from 10–14 to 106 grams (Ernest et al., 2003; Gillooly et al., 2001; Kleiber, 1932; Makarieva et al., 2005a; Makarieva et al., 2008; Savage et al., 2004; West et al., 2002). This phenomenon has been most studied in mammals, but is highly conserved, having also been shown to occur in prokaryotes (Fenchel and Finlay, 1983; Makarieva et al., 2005a; Makarieva et al., 2005a; Makarieva et al., 2008; Moses et al., 2008), plants (Makarieva et al., 2005b; Mori et al., 2010; Reich et al., 2006), insects (Chown et al., 2007; Maino and Kearney, 2014; Makarieva et al., 2005a), fish (Clarke and Johnston, 1999; Gjoni et al., 2020; Rubalcaba et al., 2020), and birds (Glazier, 2008; Hudson et al., 2013; Makarieva et al., 2005b). However, a major limitation of prior studies in this field has been that observations have been largely limited to oxygen consumption and caloric intake, leaving other metabolic processes unexplored. This study sought to address this issue by examining the generalizability of the inverse relationship between body mass and metabolic rates, using both experimental measurements and previously assembled databases that have not previously been employed in this context.
It is important to note that the metabolic processes which we observed to be higher in mice as compared to rats did not necessarily adhere quantitatively to the classic metabolic scaling relationship, with metabolic rates proportional to three-quarters of body mass. This speaks to the idea that the scaling relationship is multidimensional: it is entirely conceivable that whole-body oxygen consumption could be proportional to three-quarters of body mass, while other metabolic processes may exhibit a different scaling relationship. Further studies across species beyond rodents will be required to address this question.
The possibility that gene expression, as reflected by mRNA abundance, may also scale with body mass has not been previously addressed. We observed that the expression of key genes in glycolysis, gluconeogenesis, fatty acid metabolism, NAD synthesis and transport, mitochondrial oxygen consumption, and protection from oxidative damage scale with body mass. More compelling, however, is the observation that those genes for which an inverse relationship of expression with body mass is observed, are not randomly distributed across the genome. Rather, the collection of genes whose expression is inversely correlated with body mass is enriched for genes related to metabolic processes, and whose corresponding proteins’ enzymatic action are constrained by the supply of substrate, NAD, ATP, or oxygen.
The notion that body mass is a variable related to the level of expression of certain genes has not previously been considered as an aspect of metabolic scaling. However, it should be noted that metabolic scaling cannot fully be explained at the transcriptional level, because many rate-limiting enzymes in the metabolic processes measured in vivo did not scale at the transcriptional level, and only approximately half of genes that scaled at the level of mRNA scaled at the level of protein. Thus, it is likely that both transcriptional and other mechanisms – such as enzyme activity – are responsible for variations in metabolic flux per unit mass, inversely proportionally to body size. Additionally, the currently available data do not allow us to assess whether the expression of certain isoforms of key metabolic enzymes scales differentially across species.
It is also informative to contrast the lack of an inverse relationship between body size and metabolic fluxes per tissue weight (oxygen consumption, mitochondrial oxidation, lipolysis, and glucose production in hepatocytes) measured in the in vitro setting, to our findings in vivo (where all fluxes proceeded faster in mice than in rats). This emphasizes the need to employ tracer methods in vivo to generate a comprehensive picture of differences in metabolic fluxes between species.
Our findings emphasize that measurement of oxygen consumption in vitro may fail to detect any influence of scaling processes present in vivo. Glucose production was threefold higher in mouse liver slices relative to rat liver slices, but did not significantly differ between plated hepatocytes from mice and rats. Future studies using metabolic flux analysis may have the further capacity to generate new insights as to the scope and mechanism of metabolic scaling (Wiechert, 2001). For example, our data do not allow us to ascertain whether differences in oxygen consumption orchestrate metabolic alterations, as has been suggested in the setting of cancer cells (Nakazawa et al., 2016), or metabolic alterations require changes in oxygen consumption. Additionally, there are limitations to the fact that metabolic flux studies were performed only in the two arguably most related species included in the transcriptomics analysis (with monkey and human perhaps being similarly related). Our laboratory does not have the capacity to perform flux analysis in larger or smaller species, but we fully recognize that widening our scope beyond the 10-fold range in body size between mice and rats may yield different results. For example, it is possible that the mathematical relationship between body size and metabolic flux may reach a threshold of body size, after which this relationship may change. However, experimental and computational data suggest that the metabolic scaling relationship with regard to oxygen consumption holds at least from birds to cattle, and that a similar – but not identical – scaling relationship exists in plants (Niklas and Kutschera, 2015), so we posit that it is likely that if such a threshold exists, it occurs at a body size orders of magnitude larger than the rodents included in our study.
Taken together, the findings of this study show that the phenomenon of metabolic scaling extends to processes beyond oxygen consumption and caloric intake, and is bookended by scaling not only of gene and protein expression, but also of enzyme activity and cytosolic and mitochondrial fluxes in vivo. It may be that if one were to alter the activity of an enzyme in a larger species using a mutation to increase it to the level of activity of the same enzyme in a smaller species, that the phenomenon of metabolic scaling might disappear. Future work will be required to examine this possibility. The mechanisms underlying the phenomenon of metabolic scaling remain obscure, and deserve investigation as a fundamental question in biology (Hatton et al., 2019; Kolokotrones et al., 2010; White and Seymour, 2003). Our data support the notion that cell-extrinsic regulatory factors are involved.
Ideas and speculation
These data have implications for metabolic regulation during hibernation, in which metabolic rates are reproducibly and markedly but reversibly suppressed (Jansen et al., 2021; Jansen et al., 2019; Tøien et al., 2011). Metabolic scaling phenomena may also have medical relevance. Recent reports provide data consistent with ‘Peto’s paradox’ by showing that cancer is not more prevalent in larger, long-lived organisms than smaller ones (Makarieva et al., 2005a; Vincze et al., 2022), despite the fact that more cells are at risk for transformation, and over longer periods of time. Metabolic scaling provides a potential explanation: the slower rates of metabolic processes (including oxygen consumption) and reduced oxidative damage in larger animals imply lower rates of proliferation and carcinogenic DNA damage, so that while more cells are at risk in large animals, carcinogenesis proceeds more rapidly in smaller ones.
Metabolic scaling phenomena that lead to different rates of metabolic fluxes according to organism size deserve consideration in the context of murine models of human diseases. These differences may contribute to the faster rate of aging and shorter life expectancy of rodents as compared to humans. A separate example concerns cancer models, where species-specific different rates of in vivo metabolic fluxes might be relevant to carcinogenesis or efficacy of treatments.
Materials and methods
All animal studies were approved by the Yale University Institutional Animal Care and Use Committee (protocol 2019–20290). Source data are found in Supplementary file 1.
Rodents
Request a detailed protocolHealthy male (or, when designated, female), wild-type C57BL/6 J mice (catalog number 000664) and Sprague-Dawley rats (strain code 400) were obtained at 8 weeks of age from Jackson Laboratories and Charles River Laboratories, respectively, and given ad lib access to regular chow and water. Rodents underwent surgery under isoflurane anesthesia to place catheters in the jugular vein (mice) and in both the jugular vein and carotid artery, with the tip of the arterial catheter advanced into the right atrium of the heart (rats). After a week of recovery and confirmation that the animals had regained their pre-surgical body weight and following a 24 hr fast, rodents underwent the in vivo tracer studies described below. Rodents used for hepatocyte studies were fed ad lib until isoflurane euthanasia and liver isolation as described in the in vitro studies below.
Human participants
Request a detailed protocolPlasma samples were obtained, deidentified, from healthy, lean participants (all male, in order to maintain consistency with the rodent studies; age 26±3, body mass index 24.2 ± 0.8; ethnicity data not collected) in a Yale University IRB-approved (protocol 0108012609, NCT03738852, beginning 11/7/2018) clinical study. Eligibility included age >18 years and body mass index >18.0, without evidence of renal, thyroid, neurologic, or psychiatric dysfunction. Pre-defined outcomes were focused on the counterregulatory response to hypoglycemia, but baseline samples could be used in a deidentified manner for other analyses. All volunteers provided informed consent to participate in the study, and for the data to be published. Data were collected at Yale University. Blood was obtained from a venous catheter in the left forearm, centrifuged, and stored at –80 °C to await further analysis. No harms were observed.
Gene expression analysis
Request a detailed protocolAll liver raw transcriptomics data were obtained from Array Express (Brazma et al., 2003) (https://doi.org/10.1093/nar/gkg091), and were all preprocessed using the same methods. Cattle, monkey, and rat were obtained from E-MTAB-4550, mouse from E-MTAB-5166, and human from E-MTAB-6814. Two replicates from each species were used. Raw counts from each species were normalized to counts per million (CPM) and were then TMM-normalized to account for differences in sequencing depth as well as transcript length across species and scanners using the R package edgeR (Robinson et al., 2010). All transcript homologs were converted to human gene names using ENSEMBL in the biomaRt R package (Durinck et al., 2009). For the analyses in Figure 5, genes were filtered to those that followed the allometric scaling pattern mouse > rat > monkey > human > cattle. This gene list was put into EnrichR and KEGG pathway enrichment analysis was performed (Kuleshov et al., 2016). Genes in this list that followed the scaling pattern were then filtered based on whether they met one of six Reactome metabolic superpathways related to the metabolism of amino acids, carbohydrates, energy, lipids, tricarboxylic acid cycle, and vitamin cofactors (Jassal et al., 2020). A clustering heatmap was performed on these genes using the Seaborn Python package (Waskom, 2021). Genes were plotted again upon this map using the Seaborn Python package. To assess for differences in gene expression in genes displayed on the metabolic map, a one-way ANOVA with a Bonferroni correction for multiple comparisons (thus we used an adjusted p-value threshold for this test of 0.01, corresponding to an alpha level of 0.05, for significance) using the Python Scipy package (Virtanen et al., 2020).
Targeted gene expression analysis was performed by qPCR in our laboratory. mRNA was prepared using TRIzol (Life Technologies), purified using RNeasy mini-columns (Qiagen), and used to prepare cDNA. Quantitative real-time PCR was performed using SYBR Green (Applied Biosystems). Primers were obtained from IDT with the sequences shown in Supplementary file 1.
Protein expression analysis
Request a detailed protocolAll liver raw mass spectrometry proteomics data were obtained from the Proteomics Identifications Database (PRIDE). Mouse liver S9 fraction (PXD000733), rat liver S9 fraction (PXD000717), and human liver S9 fraction (PXD00734) raw data from a previously published study (Golizeh et al., 2015) were processed using MaxQuant software (version 2.1.4.0, Max-Planck Institute of Biochemistry, Munich, Germany), and proteins were identified with built-in Andromeda search engine based on Uniprot canonical protein collections for each species. The false discovery rate cut-offs were set to 1% on peptide, protein, and site decoy levels (default), thereby only allowing the high-quality identifications to pass. Raw intensities of all species showed similar distributions; therefore, data were normalized across species by median centering, and we used the EggNOG database to match orthologous proteins between species. Perseus was used to perform data analysis. We looked exclusively at the twenty genes that were identified with the pattern mouse > rat > human in our transcriptomics data analysis. We used iBAQ intensity (intensity Based Absolute Quantification) as a proxy for protein expression because iBAQ intensities are the raw intensities divided by the number of theoretical peptides. Therefore, iBAQ values are proportional to the molar quantities of the proteins. Additionally, the iBAQ algorithm can roughly estimate the relative abundance of proteins within a sample (Krey et al., 2018).
Because the primary goal of Golizeh et al., 2015 , was to evaluate protein/peptide-level ion exchange fractionation and proteome coverage in liver microsomes versus S9 fractions, the protein data did not include protein abundance. Consequently, the dataset did not have replicates as two livers from each species were digested differently (one with trypsin and the second with pepsin). Since the different digestion methods could alter both protein identification and abundance, we evaluated only the data from the liver S9 fraction digested with trypsin. Since we had no replicates, statistical analysis could not be carried out.
In vitro tracer analysis
Request a detailed protocolPrimary hepatocytes were isolated by the Yale Liver Center’s Cell Isolation Core, plated in a six-well plate (4.0 × 105 cells per well), and allowed to recover for 6 hr at 37 °C in substrate replete media (DMEM high glucose containing 10% FBS, 2% penicillin–streptomycin, 100 nM dexamethasone, 10 mM HEPES, and 1 nM insulin). The attached cells were then washed once in PBS and incubated overnight in low-glucose culture media (DMEM low glucose containing 10% FBS, 2% penicillin–streptomycin and 10 mM HEPES) for glucose production assays, or serum-free low-glucose culture medium (DMEM low glucose supplemented with 0.5% fatty acid free BSA, 2% penicillin–streptomycin, and 10 mM HEPES) for lipolysis assays. Following the overnight incubation, for glucose production assays, cells were washed twice with PBS and the media was replaced with 2 ml of substrate-replete glucose production medium (DMEM no-glucose base media containing 0.5% fatty acid-free BSA, 20 mM sodium lactate (50% 3-13C), 2 mM sodium pyruvate, 2 mM GlutaMAX, MEM non-essential amino acids, and 10 mM HEPES). After a s6 hr incubation at 37 °C, the media was collected. Glucose concentrations were measured as described in the Biochemical Analysis section, and normalized to total protein measured using the bicinchoninic acid (BCA) assay. For lipolysis assays, cells were washed twice with PBS and the media was replaced with a fresh 2 ml volume of low-glucose culture medium. Cells were incubated at 37 °C for 6 hr, after which the media was collected and non-esterified fatty acid concentrations measured as described in the Biochemical Analysis section. The lipolysis rate was calculated after normalizing to total protein concentrations determined using the BCA assay. Flux ratios were measured using Equations 5-13 in Supplementary file 2, and back-calculated from net glucose production determined by measuring the glucose concentration in the media using the Sekisui Glucose Assay and assuming a linear rate of net glucose production during the 6 hr incubation.
Glucose production in liver slices
Request a detailed protocolMice and rats were fasted overnight (16 hr) and sacrificed. Rodent livers were extracted and washed in Krebs-Henseleit buffer (KHB) containing 550 mM sodium chloride, 23 mM potassium chloride, 6.3 mM calcium dichloride, 10 mM magnesium sulfate, and 6.9 mM sodium phosphate monobasic. Rodent livers were cored into smaller bits using the Alabama R&D tissue coring press. Cored livers were sliced to the thickness of 230 microns, the lowest setting on the Alabama R&D tissue slicer. Liver slices were transferred to 24-well plates containing gluconeogenesis media (DMEM without glucose but with 20 mM lactate, 2 mM pyruvate, 10 mM HEPES, and 44 mM sodium bicarbonate) (400 ul for mouse livers and 500 ul for rat livers). The 24-well plates were placed in a tissue culture incubator (5% CO2) and shaken at 80 RPM for 6 hr. At the end of the 6 hr incubation, liver slices and media were collected. Liver slice weights were measured on a scale. Glucose concentrations from liver slice media were measured using the Sekisui Glucose Assay. Glucose concentrations were normalized to liver weights measured in milligrams and microliters of gluconeogenesis media.
In vitro oxygen consumption analysis
Request a detailed protocolPrimary hepatocytes were isolated from ad lib fed mice and rats by the Yale Liver Center’s Cell Isolation core and plated recovery media as described previously Camporez et al., 2013 in 24-well XF24 V7 cell culture plates coated with type I collagen. After 6–8 hr of recovery at 37 °C in 5% CO2, cells were washed twice with PBS and the media was replaced with low-glucose culture media (DMEM base medium containing 5 mM glucose, 2 mM glutamine, and non-essential amino acids, pH 7.4), in which cells were cultured overnight. The next morning, as we have previously described (Perry et al., 2020), cells were washed twice with PBS and the media was replaced with 500 μL XF24 assay medium (DMEM base medium containing 5.5 mM glucose, 1 mM pyruvate, and 2 mM glutamine, pH 7.4) and equilibrated at 37 °C for 60 min. The Seahorse XFe 24 Analyzer was used to perform a mitochondrial stress test: after three baseline measurements of O2 consumption (10 min apart), oligomycin (an inhibitor of ATP synthase) was injected, and three subsequent measures of O2 consumption was performed using a 4 min mix/2 min wait/4 min measure protocol. Next, the uncoupler 2-[2-[4-(trifluoromethoxy)phenyl]hydrazinylidene]-propanedinitrile (FCCP) was injected to dissipate the proton gradient, with three O2 consumption measurements taken as described above. Finally, rotenone (0.5 μM) and antimycin (10 μM) were injected to inhibit Complexes I and III, respectively. Oxygen consumption was normalized to total protein measured using the Pierce BCA Protein Assay.
In vivo tracer analysis
Request a detailed protocolMice and rats received a 3 X (5 min) primed-continuous infusion of [3-13C] sodium lactate (4.5 mg/kg body weight/min), [1, 2, 3, 4, 5, 6, 6-2H7] glucose (1.0 mg/kg/min), and [U-13C16] potassium palmitate (0.8 mg/kg/min) for 120 min. Tracers were infused into the jugular vein (mice) or the right atrium (rats) to ensure systemic delivery. After 100, 110, and 120 min, blood was collected from the tail vein (mice) or from the jugular venous catheter (rats) and centrifuged to obtain plasma, then animals were sacrificed with IV pentobarbital and their livers freeze-clamped in liquid nitrogen.
Biochemical analysis
Request a detailed protocolIn plated hepatocyte and liver slice studies, glucose concentrations in the cell media were determined using the Sekisui Glucose Assay. Non-esterified fatty acid concentrations in the hepatocyte media were measured using the Sekisui Non-Esterified Fatty Acid kit.
Glucose production in hepatocytes was determined by measuring the glucose concentration in media after 1, 2, 3, 4, 5, and 6 hr of incubation. Glucose enrichment in hepatocytes, plasma, and livers was measured by gas chromatography/mass spectrometry (GC/MS). Samples were deprotonized using 1:1 barium hydroxide:zinc sulfate (100 μl for cells and 300 μl for ~100 mg liver samples, which were subsequently homogenized using a TissueLyser) and derivatized with 50 μl 1:1 acetic anhydride:pyridine. After 20 min heating to 65 °C, 50 μl methanol was added, and glucose enrichment ([13C1], [13C2] and, in plasma, [2H7] (all using the Chemical Ionization mode), [4, 5, 6-13C1] and [4, 5, 6-13C2] (both using the electron ionization mode)) were measured by GC/MS.
[13C] alanine enrichment was measured by GC/MS (Perry et al., 2016). We have previously shown that flux through pyruvate kinase and malic enzyme – which would add [13C] label to carbon 2 of alanine – is minimal under fasting conditions in normal rats Perry et al., 2016; therefore, the measured [13C] alanine enrichment can be attributed entirely to [3-13C] lactate enrichment. [13C] malate enrichment was also measured by GC/MS (Perry et al., 2017b). Total malate enrichment and C1C2C3 malate enrichment were measured and the C2 + C3 malate enrichment was determined according to Supplementary file 2, equations 1 and 2, which relies on the assumption that the C4 enrichment is approximately equal to the C1 enrichment of malate (Perry et al., 2017b). Total and C4C5 glutamate enrichment was measured by LC-MS/MS (Perry et al., 2018). With 13C lactate infusion, no label enters glutamate carbon 5, so all label in the C4C5 fragment was assumed to be labeled in glutamate C4.
Plasma transaminase concentrations were measured by ELISA (Abcam; ALT: ab282882 [mouse]/ab234579 [rat], AST: ab263882 [mouse]/ab263883 [rat]). Liver glycogen concentrations were measured in 0 and 4 hr fasted mice and rats using the phenol-sulfuric acid method reported by Schaubroeck et al., 2022. PC activity (Vatner et al., 2013) and peroxidase activity (Sigma ‘Enzymatic Assay of Peroxidase (EC 1.11.1.7)’ protocol) were measured enzymatically in the liver.
In vivo flux analysis
Request a detailed protocolEndogenous palmitate and glucose production were determined using equation 3 in Supplementary file 2; Perry et al., 2017a comparing plasma enrichment to the infused tracer enrichment, measured using gas chromatography/mass spectrometry (GC/MS) as described in the Biochemical Analysis section. In the fasted/substrate-depleted and, therefore, glycogen-depleted state (Perry et al., 2018), endogenous glucose production can be attributed entirely to gluconeogenesis (equation 4). Based on equations previously described (Perry et al., 2017b), and after verifying minimal renal and hepatic bicarbonate enrichment by GC/MS (Perry et al., 2020), we measured the whole-body ratio of phosphoenolpyruvate carboxykinase (PEPCK) flux (i.e. gluconeogenesis from pyruvate) to total gluconeogenesis by mass isotopomer distribution analysis (equations 5 and 6), correcting for any [13C2] glucose synthesized from [13C2] trioses (equation 7). At steady state, VPEPCK is equal to the sum of pyruvate kinase flux and pyruvate carboxylase flux (VPK +VPC); based on our previously published data, under fasting conditions pyruvate kinase flux is minimal, less than 10% of pyruvate carboxylase flux (Perry et al., 2016). Therefore, we can assume that the rate of gluconeogenesis from pyruvate (i.e. VPEPCK) is approximately equal to VPC.
Next, we measured the ratio of pyruvate carboxylase anaplerosis to citrate synthase flux (VPC/VCS) using the enrichment of liver alanine, malate, and glutamate (equation 8). This equation, in which pyruvate cycling is again assumed minimal, is derived in detail in our recent publication (Perry et al., 2017b). We then measured the fractional contribution of glycolytic carbons to the TCA cycle (i.e. pyruvate dehydrogenase flux relative to citrate synthase flux, VPDH/VCS) (equation 9) using the model and assumptions we (Befroy et al., 2014; Perry et al., 2018; Perry et al., 2016; Song et al., 2020) and others Alves et al., 2011; Petersen et al., 2016; Petersen et al., 2015 have described. Finally, absolute turnover rates (equations 10-13) were determined, utilizing the ratios measured with equations 5, 8, and 9, and the absolute gluconeogenesis rate measured using equations 3 and 4.
Statistical considerations, power calculations, and statistical analysis
Request a detailed protocolThe sample size (n=4–6 in vitro replicates, 2–4 slices, or 4–6 animals in vivo) was calculated to supply 80% power at α=0.05 to detect the expected twofold difference with 50% variance. Power calculations were performed using the ClinCalc online calculator. Groups in the in vivo and in vitro studies were compared using the two-tailed unpaired Student’s t-test after confirming that the data met the assumptions of the test. The in vivo studies compared 4–6 biological replicates (unique animals), the liver slice studies compared 2–4 biological replicates, and the in vitro studies compared two biological replicates, with three technical replicates (separate wells) from each. No samples were excluded from analysis in the in vitro or ex vivo studies. Randomization and blinding were not possible during the in vivo studies because of the readily apparent differences between mice and rats, but all analyses were performed by investigators who were blinded as to species. No adverse events occurred. ARRIVE reporting guidelines were followed.
Data availability
Source data generated for this manuscript are included in Source Data files 1-5. Source data for the gene expression analyses are in the original paper (DOI is provided in the manuscript: https://doi.org/10.1093/nar/gkg091). The source of the proteomics data is also provided in the manuscript: All liver raw mass spectrometry proteomics data were obtained from the Proteomics Identifications Database (PRIDE): mouse liver S9 fraction (PXD000733), rat liver S9 fraction (PXD000717) and human liver S9 fraction (PXD00734) from a previously published study (Golizeh et al., 2015).
References
-
Near-optimal probabilistic RNA-Seq QuantificationNature Biotechnology 34:525–527.https://doi.org/10.1038/nbt0816-888d
-
Arrayexpress--a public repository for Microarray gene expression data at the EBINucleic Acids Research 31:68–71.https://doi.org/10.1093/nar/gkg091
-
Metabolic rate does not scale with body mass in cultured mammalian cellsAmerican Journal of Physiology. Regulatory, Integrative and Comparative Physiology 292:R2115–R2121.https://doi.org/10.1152/ajpregu.00568.2006
-
Scaling of metabolic rate with body mass and temperature in Teleost fishJournal of Animal Ecology 68:893–905.https://doi.org/10.1046/j.1365-2656.1999.00337.x
-
BookInterrelationship Between Amino Acid and Carbohydrate Metabolism During Exercise: The Glucose Alanine-Cycle InIn: Pernow B, Saltin B, editors. Muscle Metabolism During Exercise: Proceedings of a Karolinska Institutet Symposium Held in Stockholm, Sweden, September 6–9, 1970 Honorary Guest: E Hohwü Christensen, Advances in Experimental Medicine and Biology. Boston, MA: Springer US. pp. 205–214.https://doi.org/10.1007/978-1-4613-4609-8_19
-
Respiration rates in Heterotrophic, free-living protozoaMicrobial Ecology 9:99–122.https://doi.org/10.1007/BF02015125
-
Oxygen consumption in relation to body size, wave exposure, and Cirral beat behavior in the Barnacle Balanus GlandulaJournal of Crustacean Biology 33:317–322.https://doi.org/10.1163/1937240X-00002147
-
Effects of metabolic level on the body size Scaling of metabolic rate in birds and mammalsProceedings. Biological Sciences 275:1405–1410.https://doi.org/10.1098/rspb.2008.0118
-
Linking scaling laws across eukaryotesPNAS 116:21616–21622.https://doi.org/10.1073/pnas.1900492116
-
The relationship between body mass and field metabolic rate among individual birds and mammalsThe Journal of Animal Ecology 82:1009–1020.https://doi.org/10.1111/1365-2656.12086
-
Can offsetting the energetic cost of Hibernation restore an active season phenotype in grizzly bears (Ursus Arctos Horribilis)The Journal of Experimental Biology 224:jeb242560.https://doi.org/10.1242/jeb.242560
-
The Reactome pathway KnowledgebaseNucleic Acids Research 48:D498–D503.https://doi.org/10.1093/nar/gkz1031
-
Body size and tissue respirationBiochimica et Biophysica Acta 4:249–269.https://doi.org/10.1016/0006-3002(50)90032-1
-
Enrichr: A comprehensive Gene set enrichment analysis web server 2016 updateNucleic Acids Research 44:W90–W97.https://doi.org/10.1093/nar/gkw377
-
Ontogenetic and Interspecific metabolic Scaling in insectsThe American Naturalist 184:695–701.https://doi.org/10.1086/678401
-
Energetics of the smallest: Do bacteria breathe at the same rate as Whales?Proceedings. Biological Sciences 272:2219–2224.https://doi.org/10.1098/rspb.2005.3225
-
Revisiting a model of Ontogenetic growth: Estimating model parameters from theory and dataThe American Naturalist 171:632–645.https://doi.org/10.1086/587073
-
Oxygen availability and metabolic adaptationsNature Reviews. Cancer 16:663–673.https://doi.org/10.1038/nrc.2016.84
-
Kleiber’s law: How the fire of life ignited debate, fueled theory, and neglected plants as model organismsPlant Signaling & Behavior 10:e1036216.https://doi.org/10.1080/15592324.2015.1036216
-
Propionate increases hepatic pyruvate Cycling and Anaplerosis and alters mitochondrial metabolismThe Journal of Biological Chemistry 291:12161–12170.https://doi.org/10.1074/jbc.M116.720631
-
Mechanism for Leptin’s acute insulin-independent effect to reverse diabetic KetoacidosisThe Journal of Clinical Investigation 127:657–669.https://doi.org/10.1172/JCI88477
-
Regulation of hepatic mitochondrial oxidation by glucose-alanine Cycling during starvation in humansThe Journal of Clinical Investigation 129:4671–4675.https://doi.org/10.1172/JCI129913
-
Causes of differences in respiration rate of hepatocytes from mammals of different body massAmerican Journal of Physiology-Regulatory, Integrative and Comparative Physiology 269:R1213–R1224.https://doi.org/10.1152/ajpregu.1995.269.5.R1213
-
The predominance of quarter-power Scaling in biologyFunctional Ecology 18:257–282.https://doi.org/10.1111/j.0269-8463.2004.00856.x
-
An optimized method for tissue glycogen quantificationPhysiological Reports 10:e15195.https://doi.org/10.14814/phy2.15195
-
Metabolism in the age of ‘OmesThe Journal of Experimental Biology 215:2351–2357.https://doi.org/10.1242/jeb.059725
-
Relationship between fish size and metabolic rate in the Oxyconforming Inanga Galaxias Maculatus reveals size-dependent strategies to withstand hypoxiaPhysiological and Biochemical Zoology 86:740–749.https://doi.org/10.1086/673727
-
Thyroid hormone receptor-Β agonists prevent hepatic steatosis in fat-Fed rats but impair insulin sensitivity via discrete pathwaysAmerican Journal of Physiology. Endocrinology and Metabolism 305:E89–E100.https://doi.org/10.1152/ajpendo.00573.2012
-
Seaborn: Statistical data visualizationJournal of Open Source Software 6:3021.https://doi.org/10.21105/joss.03021
Article and author information
Author details
Funding
National Institutes of Health (R00CA215315)
- Rachel J Perry
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
The authors thank Traci LaMoia for her assistance in performing the liver slice experiments, Kathy Harry of the Yale Liver Center for isolating hepatocytes, and members of the Perry and Pollak labs for helpful discussions. This study was funded by grants from the U.S. Public Health Service (K99/R00 R00CA215315 [to R.J.P.], T32-GM0007324 [supporting N.D.A.], T32GM136651 [supporting B.P.L.], and P30DK034989 [supporting the Yale Liver Center]).
Ethics
Clinical trial registration NCT03738852.
Plasma samples were obtained, deidentified, from healthy, lean participants (all male, in order to maintain consistency with the rodent studies; age 26±3, body mass index 24.2±0.8; ethnicity data not collected) in a Yale University IRB-approved (protocol 0108012609) clinical study. All volunteers provided informed consent to participate in the study, and for the data to be published.
All animal studies were approved by the Yale University Institutional Animal Care and Use Committee (protocol 2019 and later 2022-20290).
Copyright
© 2023, Akingbesote et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 1,803
- views
-
- 120
- downloads
-
- 1
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Biochemistry and Chemical Biology
Ubiquitination typically involves covalent linking of ubiquitin (Ub) to a lysine residue on a protein substrate. Recently, new facets of this process have emerged, including Ub modification of non-proteinaceous substrates like ADP-ribose by the DELTEX E3 ligase family. Here, we show that the DELTEX family member DTX3L expands this non-proteinaceous substrate repertoire to include single-stranded DNA and RNA. Although the N-terminal region of DTX3L contains single-stranded nucleic acid binding domains and motifs, the minimal catalytically competent fragment comprises the C-terminal RING and DTC domains (RD). DTX3L-RD catalyses ubiquitination of the 3’-end of single-stranded DNA and RNA, as well as double-stranded DNA with a 3’ overhang of two or more nucleotides. This modification is reversibly cleaved by deubiquitinases. NMR and biochemical analyses reveal that the DTC domain binds single-stranded DNA and facilitates the catalysis of Ub transfer from RING-bound E2-conjugated Ub. Our study unveils the direct ubiquitination of nucleic acids by DTX3L, laying the groundwork for understanding its functional implications.
-
- Biochemistry and Chemical Biology
Escherichia coli ClpB and Saccharomyces cerevisiae Hsp104 are AAA+ motor proteins essential for proteome maintenance and thermal tolerance. ClpB and Hsp104 have been proposed to extract a polypeptide from an aggregate and processively translocate the chain through the axial channel of its hexameric ring structure. However, the mechanism of translocation and if this reaction is processive remains disputed. We reported that Hsp104 and ClpB are non-processive on unfolded model substrates. Others have reported that ClpB is able to processively translocate a mechanically unfolded polypeptide chain at rates over 240 amino acids (aa) per second. Here, we report the development of a single turnover stopped-flow fluorescence strategy that reports on processive protein unfolding catalyzed by ClpB. We show that when translocation catalyzed by ClpB is challenged by stably folded protein structure, the motor enzymatically unfolds the substrate at a rate of ~0.9 aa s−1 with a kinetic step-size of ~60 amino acids at sub-saturating [ATP]. We reconcile the apparent controversy by defining enzyme catalyzed protein unfolding and translocation as two distinct reactions with different mechanisms of action. We propose a model where slow unfolding followed by fast translocation represents an important mechanistic feature that allows the motor to rapidly translocate up to the next folded region or rapidly dissociate if no additional fold is encountered.