Structure and functionality of a multimeric human COQ7:COQ9 complex.

Coenzyme Q (CoQ) is a redox-active lipid essential for core metabolic pathways and antioxidant defense. CoQ is synthesized upon the mitochondrial inner membrane by an ill-defined "complex Q" metabolon. Here, we present structure-function analyses of a lipid-, substrate-, and NADH-bound complex comprising two complex Q subunits: the hydroxylase COQ7 and the lipid-binding protein COQ9. We reveal that COQ7 adopts a ferritin-like fold with a hydrophobic channel whose substrate-binding capacity is enhanced by COQ9. Using molecular dynamics, we further show that two COQ7:COQ9 heterodimers form a curved tetramer that deforms the membrane, potentially opening a pathway for the CoQ intermediates to translocate from the bilayer to the proteins' lipid-binding sites. Two such tetramers assemble into a soluble octamer with a pseudo-bilayer of lipids captured within. Together, these observations indicate that COQ7 and COQ9 cooperate to access hydrophobic precursors within the membrane and coordinate subsequent synthesis steps toward producing CoQ.

Defining mitochondrial protein functions through deep multiomic profiling.

Mitochondria are epicentres of eukaryotic metabolism and bioenergetics. Pioneering efforts in recent decades have established the core protein componentry of these organelles1 and have linked their dysfunction to more than 150 distinct disorders2,3. Still, hundreds of mitochondrial proteins lack clear functions4, and the underlying genetic basis for approximately 40% of mitochondrial disorders remains unresolved5. Here, to establish a more complete functional compendium of human mitochondrial proteins, we profiled more than 200 CRISPR-mediated HAP1 cell knockout lines using mass spectrometry-based multiomics analyses. This effort generated approximately 8.3 million distinct biomolecule measurements, providing a deep survey of the cellular responses to mitochondrial perturbations and laying a foundation for mechanistic investigations into protein function. Guided by these data, we discovered that PIGY upstream open reading frame (PYURF) is an S-adenosylmethionine-dependent methyltransferase chaperone that supports both complex I assembly and coenzyme Q biosynthesis and is disrupted in a previously unresolved multisystemic mitochondrial disorder. We further linked the putative zinc transporter SLC30A9 to mitochondrial ribosomes and OxPhos integrity and established RAB5IF as the second gene harbouring pathogenic variants that cause cerebrofaciothoracic dysplasia. Our data, which can be explored through the interactive online MITOMICS.app resource, suggest biological roles for many other orphan mitochondrial proteins that still lack robust functional characterization and define a rich cell signature of mitochondrial dysfunction that can support the genetic diagnosis of mitochondrial diseases.

PKA regulatory subunit Bcy1 couples growth, lipid metabolism, and fermentation during anaerobic xylose growth in Saccharomyces cerevisiae.

Organisms have evolved elaborate physiological pathways that regulate growth, proliferation, metabolism, and stress response. These pathways must be properly coordinated to elicit the appropriate response to an ever-changing environment. While individual pathways have been well studied in a variety of model systems, there remains much to uncover about how pathways are integrated to produce systemic changes in a cell, especially in dynamic conditions. We previously showed that deletion of Protein Kinase A (PKA) regulatory subunit BCY1 can decouple growth and metabolism in Saccharomyces cerevisiae engineered for anaerobic xylose fermentation, allowing for robust fermentation in the absence of division. This provides an opportunity to understand how PKA signaling normally coordinates these processes. Here, we integrated transcriptomic, lipidomic, and phospho-proteomic responses upon a glucose to xylose shift across a series of strains with different genetic mutations promoting either coupled or decoupled xylose-dependent growth and metabolism. Together, results suggested that defects in lipid homeostasis limit growth in the bcy1Δ strain despite robust metabolism. To further understand this mechanism, we performed adaptive laboratory evolutions to re-evolve coupled growth and metabolism in the bcy1Δ parental strain. The evolved strain harbored mutations in PKA subunit TPK1 and lipid regulator OPI1, among other genes, and evolved changes in lipid profiles and gene expression. Deletion of the evolved opi1 gene partially reverted the strain's phenotype to the bcy1Δ parent, with reduced growth and robust xylose fermentation. We suggest several models for how cells coordinate growth, metabolism, and other responses in budding yeast and how restructuring these processes enables anaerobic xylose utilization.

LipiDex 2 Integrates MSn Tree-Based Fragmentation Methods and Quality Control Modules to Improve Discovery Lipidomics.

As lipidomics experiments increase in scale and complexity, data processing tools must support workflows for new liquid chromatography-mass spectrometry (LC-MS) methods while simultaneously supporting quality controls to maximize the confidence in lipid identifications. LipiDex 2 improves lipidomics data processing algorithms from LipiDex 1 and introduces new tools for spectral matching and peak annotation functions, with improvements in speed and user-friendliness. In silico spectral library generation now supports tandem mass spectral (MSn) tree-based fragmentation methods, and the LipiDex 2 workflow fully integrates the fragmentation logic into the data processing steps to enable lipid identification at the appropriate level of structural resolution. Finally, LipiDex 2 features new modules for automated quality control checks that also allow users to visualize data quality in a data dashboard user interface.

The molecular and metabolic program by which white adipocytes adapt to cool physiologic temperatures.

Although visceral adipocytes located within the body's central core are maintained at approximately 37°C, adipocytes within bone marrow, subcutaneous, and dermal depots are found primarily within the peripheral shell and generally exist at cooler temperatures. Responses of brown and beige/brite adipocytes to cold stress are well studied; however, comparatively little is known about mechanisms by which white adipocytes adapt to temperatures below 37°C. Here, we report that adaptation of cultured adipocytes to 31°C, the temperature at which distal marrow adipose tissues and subcutaneous adipose tissues often reside, increases anabolic and catabolic lipid metabolism, and elevates oxygen consumption. Cool adipocytes rely less on glucose and more on pyruvate, glutamine, and, especially, fatty acids as energy sources. Exposure of cultured adipocytes and gluteal white adipose tissue (WAT) to cool temperatures activates a shared program of gene expression. Cool temperatures induce stearoyl-CoA desaturase-1 (SCD1) expression and monounsaturated lipid levels in cultured adipocytes and distal bone marrow adipose tissues (BMATs), and SCD1 activity is required for acquisition of maximal oxygen consumption at 31°C.

Mass Spectrometry-Based Multi-omics Integration with a Single Set of C. elegans Samples.

Mass spectrometry-based large-scale multi-omics research has proven to be powerful in answering biological questions; nonetheless, it faces many challenges from sample preparation to downstream data integration. To efficiently extract biomolecules of different physicochemical properties, preparation of various sample type needs specific tailoring, especially of difficult ones, such as Caenorhabditis elegans. In this study, we sought to develop a multi-omics sample preparation method starting with a single set ofC. elegans samples to save time, minimize variability, expand biomolecule coverage, and promote multi-omics integration. We investigated tissue disruption methods to effectively release biomolecules and optimized extraction strategies to achieve broader and more reproducible biomolecule coverage in proteomics, lipidomics, and metabolomics workflows. In our assessment, we also considered speediness and usability of the approaches. The developed method was validated through a study of 16C. elegans samples designed to shine light on mitochondrial unfolded protein response (UPRmt), induced by three unique stressors─knocking down electron transfer chain element cco-1, mitochondrial ribosome protein S5 mrps-5, and antibiotic treatment Doxycycline. Our findings suggested that the method achieved great coverage of proteome, lipidome, and metabolome with high reproducibility and validated that all stressors triggered UPRmt in C. elegans, although generating unique molecular signatures. Innate immune response was activated, and triglycerides were decreased under all three stressor conditions. Additionally, Doxycycline treatment elicited more distinct proteomic, lipidomic, and metabolomic response than the other two treatments. This method has been successfully used to process Saccharomyces cerevisiae (data not shown) and can likely be applied to other organisms for multi-omics research.

Rapid Multi-Omics Sample Preparation for Mass Spectrometry.

Multi-omics analysis is a powerful and increasingly utilized approach to gain insight into complex biological systems. One major hindrance with multi-omics, however, is the lengthy and wasteful sample preparation process. Preparing samples for mass spectrometry (MS)-based multi-omics involves extraction of metabolites and lipids with organic solvents, precipitation of proteins, and overnight digestion of proteins. These existing workflows are disparate and laborious. Here, we present a simple, efficient, and unified approach to prepare lipids, metabolites, and proteins for MS analysis. Our approach, termed the Bead-enabled Accelerated Monophasic Multi-omics (BAMM) method, combines an n-butanol-based monophasic extraction with unmodified magnetic beads and accelerated protein digestion. We demonstrate that the BAMM method affords comparable depth, quantitative reproducibility, and recovery of biomolecules as state-of-the-art multi-omics methods (e.g., Matyash extraction and overnight protein digestion). However, the BAMM method only requires about 3 h to perform, which saves 11 steps and 19 h on average compared to published multi-omics methods. Furthermore, we validate the BAMM method for multiple sample types and formats (biofluid, culture plate, and pellet) and show that in all cases, it produces high biomolecular coverage and data quality.

Improved Structural Characterization of Glycerophospholipids and Sphingomyelins with Real-Time Library Searching.

In mass spectrometry-based lipidomics, complex lipid mixtures undergo chromatographic separation, are ionized, and are detected using tandem MS (MSn) to simultaneously quantify and structurally characterize eluting species. The reported structural granularity of these identified lipids is strongly reliant on the analytical techniques leveraged in a study. For example, lipid identifications from traditional collisionally activated data-dependent acquisition experiments are often reported at either species level or molecular species level. Structural resolution of reported lipid identifications is routinely enhanced by integrating both positive and negative mode analyses, requiring two separate runs or polarity switching during a single analysis. MS3+ can further elucidate lipid structure, but the lengthened MS duty cycle can negatively impact analysis depth. Recently, functionality has been introduced on several Orbitrap Tribrid mass spectrometry platforms to identify eluting molecular species on-the-fly. These real-time identifications can be leveraged to trigger downstream MSn to improve structural characterization with lessened impacts on analysis depth. Here, we describe a novel lipidomics real-time library search (RTLS) approach, which utilizes the lipid class of real-time identifications to trigger class-targeted MSn and to improve the structural characterization of phosphotidylcholines, phosphotidylethanolamines, phosphotidylinositols, phosphotidylglycerols, phosphotidylserine, and sphingomyelins in the positive ion mode. Our class-based RTLS method demonstrates improved selectivity compared to the current methodology of triggering MSn in the presence of characteristic ions or neutral losses.

Proteomics, Lipidomics, Metabolomics, and 16S DNA Sequencing of Dental Plaque From Patients With Diabetes and Periodontal Disease.

Oral microbiome influences human health, specifically prediabetes and type 2 diabetes (Pre-DM/DM) and periodontal diseases (PDs), through complex microbial interactions. To explore these relations, we performed 16S rDNA sequencing, metabolomics, lipidomics, and proteomics analyses on supragingival dental plaque collected from individuals with Pre-DM/DM (n = 39), Pre-DM/DM and PD (n = 37), PD alone (n = 11), or neither (n = 10). We identified on average 2790 operational taxonomic units and 2025 microbial and host proteins per sample and quantified 110 metabolites and 415 lipids. Plaque samples from Pre-DM/DM patients contained higher abundance of Fusobacterium and Tannerella than plaques from metabolically healthy patients. Phosphatidylcholines, plasmenyl phosphatidylcholines, ceramides containing non-OH fatty acids, and host proteins related to actin filament rearrangement were elevated in plaques from PD versus non-PD samples. Cross-omic correlation analysis enabled the detection of a strong association between Lautropia and monomethyl phosphatidylethanolamine (PE-NMe), which is striking because synthesis of PE-NMe is uncommon in oral bacteria. Lipidomics analysis of in vitro cultures of Lautropia mirabilis confirmed the synthesis of PE-NMe by the bacteria. This comprehensive analysis revealed a novel microbial metabolic pathway and significant associations of host-derived proteins with PD.

Maximizing MS/MS Acquisition for Lipidomics Using Capillary Separation and Orbitrap Tribrid Mass Spectrometer.

Liquid chromatography-mass spectrometry (LC-MS) is a typical strategy for lipidomics analysis. Although capillary LC-MS is a common analytical technique for proteomics analysis, its application to lipidomics has been limited. In this study, we aim at improving lipid identifications achieved in a single LC-MS analysis by a 3-fold approach: capillary LC and nanoelectrospray for enhanced ionization, ion trap for higher sensitivity tandem MS, and parallelization of mass analyzers for increased speed of acquisition on an Orbitrap hybrid system. By applying the methods to a complex lipid mixture of human plasma, we identified and performed relative quantification on over 1500 lipids within a 60 min capillary LC-MS analysis.

Multi-Omic Single-Shot Technology for Integrated Proteome and Lipidome Analysis.

Mass spectrometry (MS) serves as the centerpiece technology for proteome, lipidome, and metabolome analysis. To gain a better understanding of the multifaceted networks of myriad regulatory layers in complex organisms, integration of different multiomic layers is increasingly performed, including joint extraction methods of diverse biomolecular classes and comprehensive data analyses of different omics. Despite the versatility of MS systems, fractured methodology drives nearly all MS laboratories to specialize in analysis of a single ome at the exclusion of the others. Although liquid chromatography-mass spectrometry (LC-MS) analysis is similar for different biomolecular classes, the integration on the instrument level is lagging behind. The recent advancements in high flow proteomics enable us to take a first step towards integration of protein and lipid analysis. Here, we describe a technology to achieve broad and deep coverage of multiple molecular classes simultaneously through multi-omic single-shot technology (MOST), requiring only one column, one LC-MS instrument, and a simplified workflow. MOST achieved great robustness and reproducibility. Its application to a Saccharomyces cerevisiae study consisting of 20 conditions revealed 2842 protein groups and 325 lipids and potential molecular relationships.

Unique inflammatory profile is associated with higher SARS-CoV-2 acute respiratory distress syndrome (ARDS) mortality.

The COVID19 pandemic has caused more than a million of deaths worldwide, primarily due to complications from COVID19-associated acute respiratory distress syndrome (ARDS). Controversy surrounds the circulating cytokine/chemokine profile of COVID19-associated ARDS, with some groups suggesting that it is similar to patients without COVID19 ARDS and others observing substantial differences. Moreover, although a hyperinflammatory phenotype associates with higher mortality in non-COVID19 ARDS, there is little information on the inflammatory landscape's association with mortality in patients with COVID19 ARDS. Even though the circulating leukocytes' transcriptomic signature has been associated with distinct phenotypes and outcomes in critical illness including ARDS, it is unclear whether the mortality-associated inflammatory mediators from patients with COVID19 are transcriptionally regulated in the leukocyte compartment. Here, we conducted a prospective cohort study of 41 mechanically ventilated patients with COVID19 infection using highly calibrated methods to define the levels of plasma cytokines/chemokines and their gene expressions in circulating leukocytes. Plasma IL1RA and IL8 were found positively associated with mortality, whereas RANTES and EGF negatively associated with that outcome. However, the leukocyte gene expression of these proteins had no statistically significant correlation with mortality. These data suggest a unique inflammatory signature associated with severe COVID19.

Orosomucoid-like 3 Supports Rhinovirus Replication in Human Epithelial Cells.

Polymorphism at the 17q21 gene locus and wheezing responses to rhinovirus (RV) early in childhood conspire to increase the risk of developing asthma. However, the mechanisms mediating this gene-environment interaction remain unclear. In this study, we investigated the impact of one of the 17q21-encoded genes, ORMDL3 (orosomucoid-like 3), on RV replication in human epithelial cells. ORMDL3 knockdown inhibited RV-A16 replication in HeLa, BEAS-2B, A549, and NCI-H358 epithelial cell lines and primary nasal and bronchial epithelial cells. Inhibition varied by RV species, as both minor and major group RV-A subtypes RV-B52 and RV-C2 were inhibited but not RV-C15 or RV-C41. ORMDL3 siRNA did not affect expression of the major group RV-A receptor ICAM-1 or initial internalization of RV-A16. The two major outcomes of ORMDL3 activity, SPT (serine palmitoyl-CoA transferase) inhibition and endoplasmic reticulum (ER) stress induction, were further examined: silencing ORMDL3 decreased RV-induced ER stress and IFN-ß mRNA expression. However, pharmacologic induction of ER stress and concomitant increased IFN-ß inhibited RV-A16 replication. Conversely, blockade of ER stress with tauroursodeoxycholic acid augmented replication, pointing to an alternative mechanism for the effect of ORMDL3 knockdown on RV replication. In comparison, the SPT inhibitor myriocin increased RV-A16 but not RV-C15 replication and negated the inhibitory effect of ORMDL3 knockdown. Furthermore, lipidomics analysis revealed opposing regulation of specific sphingolipid species (downstream of SPT) by myriocin and ORMDL3 siRNA, correlating with the effect of these treatments on RV replication. Together, these data revealed a requirement for ORMDL3 in supporting RV replication in epithelial cells via SPT inhibition.

High-content phenotypic analysis of a C. elegans recombinant inbred population identifies genetic and molecular regulators of lifespan.

Lifespan is influenced by complex interactions between genetic and environmental factors. Studying those factors in model organisms of a single genetic background limits their translational value for humans. Here, we mapped lifespan determinants in 85 genetically diverse C. elegans recombinant intercross advanced inbred lines (RIAILs). We assessed molecular profiles - transcriptome, proteome, and lipidome - and life-history traits, including lifespan, development, growth dynamics, and reproduction. RIAILs exhibited large variations in lifespan, which positively correlated with developmental time. Among the top candidates obtained from multi-omics data integration and QTL mapping, we validated known and novel longevity modulators, including rict-1, gfm-1 and mltn-1. We translated their relevance to humans using UK Biobank data and showed that variants in RICTOR and GFM1 are associated with an elevated risk of age-related heart disease, dementia, diabetes, kidney, and liver diseases. We organized our dataset as a resource (https://lisp-lms.shinyapps.io/RIAILs/) that allows interactive explorations for new longevity targets.

Lysosomal storage disease proteo/lipidomic profiling using nMOST links ferritinophagy with mitochondrial iron deficiencies in cells lacking NPC2.

Lysosomal storage diseases (LSDs) comprised ~50 monogenic diseases characterized by the accumulation of cellular material in lysosomes and associated defects in lysosomal function, but systematic molecular phenotyping is lacking. Here, we develop a nanoflow-based multi-omic single-shot technology (nMOST) workflow allowing simultaneously quantify HeLa cell proteomes and lipidomes from more than two dozen LSD mutants, revealing diverse molecular phenotypes. Defects in delivery of ferritin and its autophagic receptor NCOA4 to lysosomes (ferritinophagy) were pronounced in NPC2-/- cells, which correlated with increased lyso-phosphatidylcholine species and multi-lamellar membrane structures visualized by cryo-electron-tomography. Ferritinophagy defects correlated with loss of mitochondrial cristae, MICOS-complex components, and electron transport chain complexes rich in iron-sulfur cluster proteins. Strikingly, mitochondrial defects were alleviated when iron was provided through the transferrin system. This resource reveals how defects in lysosomal function can impact mitochondrial homeostasis in trans and highlights nMOST as a discovery tool for illuminating molecular phenotypes across LSDs.

Scd1 and monounsaturated lipids are required for autophagy and survival of adipocytes.

OBJECTIVE: Exposure of adipocytes to 'cool' temperatures often found in the periphery of the body induces expression of Stearoyl-CoA Desaturase-1 (Scd1), an enzyme that converts saturated fatty acids to monounsaturated fatty acids. The goal of this study is to further investigate the roles of Scd in adipocytes. METHOD: In this study, we employed Scd1 knockout cells and mouse models, along with pharmacological Scd1 inhibition to dissect the enzyme's function in adipocyte physiology. RESULTS: Our study reveals that production of monounsaturated lipids by Scd1 is necessary for fusion of autophagosomes to lysosomes and that with a Scd1-deficiency, autophagosomes accumulate. In addition, Scd1-deficiency impairs lysosomal and autolysosomal acidification resulting in vacuole accumulation and eventual cell death. Blocking autophagosome formation or supplementation with monounsaturated fatty acids maintains vitality of Scd1-deficient adipocytes. CONCLUSION: This study demonstrates the indispensable role of Scd1 in adipocyte survival, with its inhibition in vivo triggering autophagy-dependent cell death and its depletion in vivo leading to the loss of bone marrow adipocytes.

Hem25p is a mitochondrial IPP transporter.

Coenzyme Q (CoQ, ubiquinone) is an essential cellular cofactor comprised of a redox-active quinone head group and a long hydrophobic polyisoprene tail. How mitochondria access cytosolic isoprenoids for CoQ biosynthesis is a longstanding mystery. Here, via a combination of genetic screening, metabolic tracing, and targeted uptake assays, we reveal that Hem25p-a mitochondrial glycine transporter required for heme biosynthesis-doubles as an isopentenyl pyrophosphate (IPP) transporter in Saccharomyces cerevisiae. Mitochondria lacking Hem25p fail to efficiently incorporate IPP into early CoQ precursors, leading to loss of CoQ and turnover of CoQ biosynthetic proteins. Expression of Hem25p in Escherichia coli enables robust IPP uptake demonstrating that Hem25p is sufficient for IPP transport. Collectively, our work reveals that Hem25p drives the bulk of mitochondrial isoprenoid transport for CoQ biosynthesis in yeast.

A transcriptional network governing ceramide homeostasis establishes a cytokine-dependent developmental process.

Transcriptional mechanisms controlling developmental processes establish and maintain proteomic networks, which can govern the levels of intracellular small molecules. Although dynamic changes in bioactive small molecules can link transcription factor and genome activity with cell state transitions, many mechanistic questions are unresolved. Using quantitative lipidomics and multiomics, we discover that the hematopoietic transcription factor GATA1 establishes ceramide homeostasis during erythroid differentiation by regulating genes encoding sphingolipid metabolic enzymes. Inhibiting a GATA1-induced sphingolipid biosynthetic enzyme, delta(4)-desaturase, or disrupting ceramide homeostasis with cell-permeable dihydroceramide or ceramide is detrimental to erythroid, but not myeloid, progenitor activity. Coupled with genetic editing-based rewiring of the regulatory circuitry, we demonstrate that ceramide homeostasis commissions vital stem cell factor and erythropoietin signaling by opposing an inhibitory protein phosphatase 2A-dependent, dual-component mechanism. Integrating bioactive lipids as essential components of GATA factor mechanisms to control cell state transitions has implications for diverse cell and tissue types.

Hem25p is required for mitochondrial IPP transport in fungi.

Coenzyme Q (CoQ, ubiquinone) is an essential cellular cofactor composed of a redox-active quinone head group and a long hydrophobic polyisoprene tail. How mitochondria access cytosolic isoprenoids for CoQ biosynthesis is a longstanding mystery. Here, via a combination of genetic screening, metabolic tracing and targeted uptake assays, we reveal that Hem25p-a mitochondrial glycine transporter required for haem biosynthesis-doubles as an isopentenyl pyrophosphate (IPP) transporter in Saccharomyces cerevisiae. Mitochondria lacking Hem25p failed to efficiently incorporate IPP into early CoQ precursors, leading to loss of CoQ and turnover of CoQ biosynthetic proteins. Expression of Hem25p in Escherichia coli enabled robust IPP uptake and incorporation into the CoQ biosynthetic pathway. HEM25 orthologues from diverse fungi, but not from metazoans, were able to rescue hem25∆ CoQ deficiency. Collectively, our work reveals that Hem25p drives the bulk of mitochondrial isoprenoid transport for CoQ biosynthesis in fungi.

Plasma metabolomics supports non-fasted sampling for metabolic profiling across a spectrum of glucose tolerance in the Nile rat model for type 2 diabetes.

Type 2 diabetes is a challenge in modern healthcare, and animal models are necessary to identify underlying mechanisms. The Nile rat (Arvicanthis niloticus) develops diet-induced diabetes rapidly on a conventional rodent chow diet without genetic or chemical manipulation. Unlike common laboratory models, the outbred Nile rat model is diurnal and has a wide range of overt diabetes onset and diabetes progression patterns in both sexes, better mimicking the heterogeneous diabetic phenotype in humans. While fasted blood glucose has historically been used to monitor diabetic progression, postprandial blood glucose is more sensitive to the initial stages of diabetes. However, there is a long-held assumption that ad libitum feeding in rodent models leads to increased variance, thus masking diabetes-related metabolic changes in the plasma. Here we compared repeatability within triplicates of non-fasted or fasted plasma samples and assessed metabolic changes relevant to glucose tolerance in fasted and non-fasted plasma of 8-10-week-old male Nile rats. We used liquid chromatography-mass spectrometry lipidomics and polar metabolomics to measure relative metabolite abundances in the plasma samples. We found that, compared to fasted metabolites, non-fasted plasma metabolites are not only more strongly associated with glucose tolerance on the basis of unsupervised clustering and elastic net regression model, but also have a lower replicate variance. Between the two sampling groups, we detected 66 non-fasted metabolites and 32 fasted metabolites that were associated with glucose tolerance using a combined approach with multivariable elastic net and individual metabolite linear models. Further, to test if metabolite replicate variance is affected by age and sex, we measured non-fasted replicate variance in a cohort of mature 30-week-old male and female Nile rats. Our results support using non-fasted plasma metabolomics to study glucose tolerance in Nile rats across the progression of diabetes.