Metabolites

Mass Spectrometry Metabolomics Identification Connection Algorithm (MSMICA) is an algorithm for automated metabolite identification in untargeted liquid chromatography-high-resolution mass spectrometry (LC-HRMS) analyses. Limitations in metabolite identification can occur due to the availability and cost of standards and prevent recognition of metabolic factors impacting human health and disease. MSMICA performs mass-to-charge-ratio matching with chemical structures and clusters of LC-HRMS features for adduct and isotope forms. A local optimization is then used to integrate retention time prediction, metabolite precursor-product and transporter correlations, and biospecimen-specific abundance information for metabolite identification. Applying MSMICA to various internal and external mammalian datasets, validation results showed a 96.2 +- 5.1% correct rate of metabolite identification. When multiple LC-HRMS datasets were used, MSMICA enabled greater metabolite identifications, expanded metabolic pathway coverage, and data harmonization. Thus, MSMICA applies multiple pieces of evidence to substantially improve metabolite identification coverage and accuracy for known metabolites.

Isotopic Metabolic Flux Analysis (I-MFA) is a standard approach for estimating intracellular metabolic fluxes. I-MFA infers fluxes by comparing simulated and measured metabolite isotopologue distributions (MIDs) of metabolites from isotope labeling experiments. MIDs represent fractional abundances that strictly sum to one for any given metabolite, thus they are inherently compositional data. However, state-of-the-art estimation approaches rely on calculating standard Euclidean distances between MIDs in a non-compositional paradigm, introducing a systemic bias. To resolve this, our study proposes compositional I-MFA. We demonstrate how to construct a meaningful orthonormal basis for MIDs via ordered sequential binary partitioning, which can be used to perform isometric log-ratio (ILR) transformation. As a minimal change to existing I-MFA workflows, we suggest estimating fluxes by minimizing Euclidean distances between ILR-transformed MIDs. We validated this framework against traditional methods using both a toy model and a biologically realistic model, evaluating point estimates, sensitivity across varied true fluxes, and confidence intervals. In the two examples, compositional I-MFA consistently outperformed traditional approaches, reducing mean squared error of flux point estimates by an average of 42.6% and substantially narrowing confidence intervals. We conclude that compositional data analysis significantly improves I-MFA and can be implemented as a simple drop-in replacement for current pipelines. Graphical Abstract O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=156 SRC="FIGDIR/small/742769v1_ufig1.gif" ALT="Figure 1"> View larger version (34K): org.highwire.dtl.DTLVardef@1a53aa4org.highwire.dtl.DTLVardef@ad225aorg.highwire.dtl.DTLVardef@aa430eorg.highwire.dtl.DTLVardef@1880ca_HPS_FORMAT_FIGEXP M_FIG C_FIG HighlightsO_LINew compositional data approach improves metabolic flux estimation. C_LIO_LIThis data transformation requires minimal changes to existing workflows. C_LIO_LIThe new method reduced MSE of flux estimates by 42.6% in two examples tested. C_LIO_LIThe confidence intervals of the estimated fluxes were substantially narrowed. C_LIO_LIEstimation accuracy remained robust across a wide range of metabolic fluxes. C_LI

Metabolomics studies are increasingly being applied with hundreds to thousands, even tens of thousands of samples that demand rapid, automated data processing while maintaining analytical sensitivity or quantitative accuracy. A major computational bottleneck is feature finding, which is the transformation of LC-MS and LC-MS/MS data into a set of analyte signals aligned and quantified across samples. Feature finding can be computationally intensive and often requires manual iterative parameter optimization. To accelerate this process, we present the Everything Bagel (EB) feature finder, an ultra-fast automated feature finding tool that integrates feature detection, retention-time alignment, and gap filling designed for run-time and memory efficiency. We benchmarked EB against two automated feature finding methods on eight benchmarking datasets. Specifically, we evaluated these three feature finding methods by measuring spike-in standard detection coverage, dilution series quantification accuracy, and yeast 12C/13C credentialed features. In this evaluation, the EB feature finder achieved performance comparable to, and often exceeding, existing methods while requiring up to 150-fold lower CPU hours and up to 113-fold lower wall time. We further demonstrated the bioanalytical validity of EB by reanalyzing published datasets used for biomarker discovery and reproduced biologically significant features that matched the published findings using manually tuned feature finding settings. Taken along with the speed improvements, we anticipate EB will enhance the ability to automatically analyze datasets with thousands to tens of thousands of samples for the community.

High dietary fructose consumption is a major contributor to the development of obesity and related cardiometabolic diseases, highlighting the need for accurate assessment of fructose metabolism in humans. Stable isotope tracer approaches, such as 13C6-fructose, require highly sensitive and specific analytical methods to quantify both concentrations and isotopic enrichments. In this study, we developed and validated a robust gas chromatography-triple quadrupole mass spectrometry (GC-QQQ)-based method for the simultaneous measurement of unlabeled and 13C6-fructose in human plasma. The method employs oximation and per-acetate derivatization, and demonstrates high specificity and accuracy. Intra- and inter-assay precision were below 10%, with no detectable carry-over, and a lower limit of quantification (LLOQ) of 0.1 nmol/mL for concentration and 0.02 molar percent excess (MPE%) for enrichment and no interference from glucose. We further compared data acquisition using multiple reaction monitoring (MRM) and selected ion monitoring (SIM). MRM showed superior performance at the low concentrations and enrichment levels characteristic of clinical plasma samples, resulting in improved sensitivity and lower LLOQs compared to SIM. Overall, this validated method provides a sensitive and reliable approach for fructose tracer studies in humans. Its application will facilitate robust investigations into fructose metabolism, and its role in metabolic dysregulation and obesity-related disease.

External biological knowledge provides valuable information about relationships among metabolites, yet this information is usually not incorporated directly into statistical estimation procedures. Most existing approaches estimate metabolite effects independently, ignoring known biochemical structure such as shared subclasses and pathway membership. We propose a Bayesian hierarchical framework that improves metabolite effect estimates by incorporating external biological information describing relationships among metabolites. The proposed method improves metabolite-specific estimates by allowing related metabolites to borrow information from one another while preserving metabolite-level inference. We evaluate the methodology using simulation studies across a range of sample sizes and heterogeneity regimes together with three metabolomics applications involving distinct biological annotation structures. Across both simulated and real datasets, incorporating external biological information consistently improves metabolite effect estimation. Gains are most pronounced when sample sizes are small and metabolite classes are informative, i.e. more homogenous within classes.

Background: Low-dose computed tomography (LDCT) screening reduces lung cancer mortality but is limited by low uptake and associated with high rates of false-positives and indeterminate-nodules. Breath volatile organic compound (VOC) analysis is a non-invasive candidate biomarker approach that could complement LDCT, but prior work has relied on laboratory-based high-resolution mass spectrometry (HRMS), limiting point-of-care deployment. Methods: In this pilot study, breath samples were collected from 40 patients with treatment-naive, pathologically confirmed non-small cell lung cancer (NSCLC) and 25 lung-cancer-screening-eligible healthy controls. Paired samples were analyzed via a compact point-of-care GC-MS platform (CLARION) and a laboratory HRMS reference. Diagnostic classification models were built independently for each platform using elastic net logistic regression with leave-one-out cross-validation, and performance was evaluated by area under the receiver operating characteristic curve (AUC). Results: CLARION identified 103 VOCs across breath specimens, compared to over 900 identified by HRMS. Despite this difference in panel size, CLARION achieved diagnostic performance nearly identical to HRMS for distinguishing NSCLC cases from controls (AUC 0.864 vs. 0.863). Compared to controls, performance statistics were similar for early-stage NSCLC (AUC 0.854 vs. 0.841) and adenocarcinoma (AUC 0.770 vs. 0.787). VOCs of interest include p-cymene, phenol, propylbenzene, tetradecane, {beta}-ocimene, 2,3-dihydro-indole, and 1-methylthio-(Z)-1-propene. Conclusion: A compact, point-of-care breath GC-MS platform achieved diagnostic performance for NSCLC detection comparable to a laboratory HRMS reference despite a substantially smaller detected VOC panel. These findings support continued development of point-of-care breath VOC testing as a non-invasive, field-deployable complement to LDCT-based lung cancer screening.

Acetyl-CoA and other acyl-CoA thioesters are critical intermediates in the metabolic reactions that cells rely on to produce energy and carry out biosynthesis. Therefore, the analysis of acyl-CoA provides valuable information about the metabolic activity of cells, especially when combined with stable isotope tracing. Acyl-CoA species are routinely monitored by reversed-phase liquid chromatography coupled to tandem mass spectrometry (RPLC-MS/MS). However, drastic differences in the hydrophobicity of short-chain versus long-chain acyl-CoA species have been challenging to accommodate with a single set of RPLC conditions. Here, we describe a convenient method based on hydrophilic interaction liquid chromatography (HILIC-MS/MS) for the concurrent detection of both short-chain and long-chain acyl-CoA and their corresponding acyl-carnitine species. Using this strategy, we tracked the metabolism of isotope-labeled fatty acids in multiple cell lines, which revealed differences in their propensities for fatty acid oxidation and the extent to which isotope incorporation into acyl-CoA mirrored that of acyl-carnitine. We also applied the HILIC-MS/MS workflow to the analysis of NADH and ATP, making it a useful technique for gauging cellular bioenergetics as reflected by the acetyl-CoA/CoA, NADH/NAD+, and ATP/ADP ratios. Altogether, this HILIC-MS/MS platform enables a streamlined analysis of acyl-CoA species and other key intermediates in cell metabolism.

Toxicological omics studies require comprehensive metadata to support reproducibility, interoperability, and regulatory reuse. However, metadata requirements differ across public repositories, reporting frameworks, and laboratory workflows, resulting in inconsistent annotation and limited data integration. To address this challenge, we developed MeSyTo (Metadata for Systems Toxicology), an ontology-driven framework for harmonizing metadata across toxicological omics. Metadata concepts from public repositories, the OECD Omics Reporting Framework (OORF), community standards, and institutional workflows were semantically aligned and implemented as the MeSyTo Metadata Model (MMM). The MMM serves as the basis for the automatic generation of SHACL validation shapes and framework-specific metadata profiles, while curated value sets are represented as SKOS controlled vocabularies to support metadata collection and validation. The current implementation comprises 105 ontology classes and 527 data properties and supports transcriptomics, proteomics, and metabolomics. A prototype web application demonstrates ontology-driven metadata collection with integrated semantic validation and ontology-based term resolution. The ontology, validation shapes, controlled vocabularies, generation scripts, and software are publicly available as open-source resources. MeSyTo provides a reusable semantic foundation for harmonized, machine-actionable metadata and facilitates repository submission, regulatory reporting, and interoperable data exchange across toxicological omics studies.

BackgroundEmerging evidence suggests that PFAS can cross blood-brain barrier and lead to neurotoxicity. Recent evidence also suggest that PFAS can bioaccumulate in gut microbiota resident in the gut. However, how gut microbes influence PFAS-driven reorganization of metabolic homeostasis in the brain is poorly understood. MethodsTo address this gap, we investigated how gut microbiota influences brain metabolomic and lipidomic responses to PFAS exposure. Specific pathogen-free (SPF) and germ-free (GF) mice were fed an obesogenic diet for 8 weeks to promote metabolic disturbance. After 1 week of acclimation, half received control water and half received water containing a PFAS mixture (PFHxS, GenX, PFOA, PFOS, and FTOH mixture). Plasma and brain samples (cortex, subcortex, cerebellum, olfactory bulb, and brainstem) were collected after 8 weeks. Untargeted analyses were performed for lipidomic, metabolomic and PFAS using high resolution liquid chromatography tandem mass spectrometry (LC-MS/MS). Data was processed using MassCube with open-sources libraries. ResultsPFHxS, GenX, PFOA, PFOS, PFDA, and PFDS were detected in plasma. PFHxS, PFOA, PFOS, and PFDS were detected across all five brain regions, with PFOS as the predominant brain-enriched species. Pathway analysis identified nicotinate and nicotinamide metabolism as the most consistently PFAS-altered pathway in both SPF and GF mice. PFAS exposure induced region-specific metabolic remodeling, with gut microbiota differentially modulating responses in the cortex, cerebellum, and brainstem, whereas the olfactory bulb showed a largely microbiota-independent response. In addition to local effects within individual brain regions, plasma-brain analysis suggested systemic metabolic responses across tissues, with association strength varying by brain region and microbiome status. Gut microbiota also shaped PFAS-induced lipid dysregulation in the brain, and methylnicotinamide and delta-valerobetaine were among the most responsive metabolites. ConclusionThis study is the first to demonstrate that resident microbiota impact PFAS-associated metabolic remodeling across the gut-plasma-brain axis. HighlightsO_LIPFAS-induced metabolic remodeling in the brain is modified by gut microbiota. C_LIO_LIPFAS exposure alters nicotinate and nicotinamide metabolism throughout the brain. C_LIO_LIPFAS-induced brain metabolic responses are region specific and microbiota dependent. C_LIO_LIPlasma-brain analysis suggests potential systemic metabolic disruption by PFAS. C_LI Graphical abstract O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=104 SRC="FIGDIR/small/743341v1_ufig1.gif" ALT="Figure 1"> View larger version (38K): org.highwire.dtl.DTLVardef@15301deorg.highwire.dtl.DTLVardef@9fac0aorg.highwire.dtl.DTLVardef@d7f0f4org.highwire.dtl.DTLVardef@10c29c2_HPS_FORMAT_FIGEXP M_FIG C_FIG

External quality assessment (EQA) of multianalyte assays is commonly interpreted analyte by analyte, although many panels contain known relations among measured features that may reveal joint quality patterns. We propose PathEQA, a feature-graph-guided random forest framework in which a user-supplied graph can represent biochemical pathways, molecular interactions, shared measurement processes, or other domain relations. The same graph is allowed to influence feature representation, node-level candidate generation, and split selection, with an optional local grouped decision. We evaluated the framework in graph-aligned and graph-misspecified simulations and used a six-analyte catecholamine-related liquid chromatography-tandem mass spectrometry EQA data set as an illustrative case study (929 records from 58 laboratories and 117 complete multianalyte panels). In graph-aligned simulations, the grouped variant reduced test root mean squared error by 7.4-9.4% relative to ordinary random forest across training sizes of 60-240, whereas graph misspecification could worsen prediction. In the catecholamine case study, full PathEQA was comparable with ordinary random forest in laboratory-grouped cross-validation (RMSE 0.570 versus 0.569) and modestly better in the final-round temporal holdout (0.307 versus 0.318); a simpler static network-sampling baseline performed best. Dopamine-norepinephrine was the strongest pair, whereas dopamine-norepinephrine-epinephrine best estimated multianalyte failure burden. These results support a general conclusion: feature-graph guidance can improve small-sample multivariate quality assessment when the supplied structure is outcome-relevant, but graph relevance must be tested rather than assumed. Catecholamines serve here as a worked example rather than a restriction of the framework.

BackgroundChronic obstructive pulmonary disease, primarily caused by exposure to cigarette smoke, is a heterogeneous lung condition characterized by complex metabolic alterations. The metabolic changes associated with smoking status have not been thoroughly investigated. Our study aims to explore the metabolite profile of COPD patients categorised by their smoking habits, including smokers, ex-smokers, and non-smokers. MethodsIn this study, the plasma metabolome of smoking stratified COPD patients were assessed using gas chromatography coupled to mass spectrometry. We applied multivariate and univariate statistical analysis to identify the differentially abundant metabolites. ResultsWe identified 23 altered metabolites in the smokers and 36 in the ex-smokers COPD subgroups. Interestingly, in comparison to the control group, no significant alteration was observed in the plasma of non-smoker COPD patients. Additionally, pathway enrichment analysis revealed top dysregulated metabolic pathways, including biosynthesis of unsaturated fatty acids, galactose metabolism, phenylalanine, tyrosine, and tryptophan biosynthesis, and glycosylphosphatidylinositol (GPI)-anchor biosynthesis. The receiver operating characteristic curve screened five metabolites, such as tetradecanoic acid, 2,4-di-tert-butylphenol, chloroxylenol, tetradecanal, and 1-dodecene, with the highest diagnostic performance (AUC > 0.8). ConclusionThis study reveals distinct plasma metabolic signatures across COPD subgroups categorized by cigarette smoking history.

Background: Evidence suggests chronic distress influences ovarian cancer (OC) etiology and metabolomic profiles. Here, we evaluated the association of a metabolite-based distress score (MDS) and OC risk. Methods: We included two matched case-control studies nested within the Nurses' Health Studies (N=584) and the Prostate, Lung, Colorectal, and Ovarian Cancer Screening Trial (N=348). Metabolites were measured 3-27 years before diagnosis using liquid-chromatography tandem mass spectrometry. We examined the association of quintiles of MDS and 19 constituent metabolites with OC risk using unconditional logistic regression and stratified by tumor histotype, menopausal status, and age at diagnosis. Results: We observed women in the highest versus lowest quintile of MDS had an increased OC risk (OR=1.62,95%CI=1.03-2.54,ptrend=0.07), and type 2 tumors (OR=1.71,95%CI=1.03-2.83,ptrend=0.11). Associations were suggestively stronger for premenopausal and <69-year-old women, and driven by pseudouridine, and N2,N2-dimethylguanosine. Conclusion: Our findings suggest chronic distress-associated metabolic dysregulation may represent a novel OC risk factor, especially among younger women.

Gastroesophageal reflux disease (GERD) is a prevalent chronic disorder where dietary modifications, particularly reducing spicy foods, are a primary management strategy. Salsa, a widely consumed condiment whose spiciness comes from capsaicin, lacks standardized heat labelling, potentially leading to inconsistent capsaicin exposure for consumers. To address this, our study aimed to develop and apply an LC-MS workflow for accurate capsaicin quantification in commercially available salsas. This approach seeks to provide objective "reflux-conscious" spice classification, supporting evidence-based dietary recommendations for individuals with GERD. In the eight commercial brands we examined, we found that products labelled "mild" had significantly lower capsaicin levels as compared to "medium" or "hot", but that there was an almost 15-fold range of capsaicin within this group. Surprisingly, there was no statistical difference in capsaicin content between those groups labelled "medium" or "hot" facilitating unambiguous assignment to either category, revealing that product labelling alone is insufficient to guide consumers seeking to control capsaicin exposure in their food. The results in this study enable improved brand-specific recommendations for GERD symptom management.

Colorectal cancer cells increase glycolysis to help meet the metabolic demands required for cell growth. Many factors, both endogenous and exogenous, likely drive cellular metabolism and enhance glycolytic flux in colorectal cells. Interleukin-1 beta (IL-1{beta}) is a pro-inflammatory cytokine that is elevated in colorectal cancer. In this study, we investigated the effect of IL-1{beta} toward driving the cancer cell to increase glycolysis, while also suppressing the oxidation of the fiber-derived nutrient butyrate. The results presented here demonstrate that IL-1{beta} stimulated glycolysis and inhibited maximal mitochondrial respiration. IL-1{beta} also increased the phosphorylation of AKT and hypoxia-inducible factor 1 alpha (HIF1) levels. Utilizing colorectal cancer cells with AKT1/2 or HIF1 knocked out showed the requirement of these proteins in mediating the increase in glycolysis following IL-1{beta} treatment. Importantly, AKT1/2 was identified as upstream of HIF1, as IL-1{beta} still increased phosphorylation of AKT even in the absence of HIF1. However, loss of AKT1/2 completely abolished the ability of IL-1{beta} to increase HIF1 protein levels. Tumor necrosis factor alpha (TNF), another cytokine found to be elevated in colorectal cancer, also increased glycolysis in an AKT and HIF1-dependent manner. Our data point to a common pathway through AKT activation and HIF1 upregulation, by which pro-inflammatory cytokines increase glycolysis in colorectal cancer cells to help promote cancer progression.

Hepatocellular Carcinoma (HCC) arising from Metabolic Dysfunction-Associated Steatotic Liver Disease (MASLD) is an increasing public health burden with high mortality, highlighting the need for improved early detection strategies. Current surveillance tools, including Alpha-fetoprotein (AFP) and ultrasound, lack sufficient sensitivity for early-stage HCC detection. We analyzed serum samples from 131 patients, including 58 with cirrhosis and 73 with MASLD-related HCC (42 early-stage, 31 late-stage), using an nLC-stepped HCD-PRM-MS/MS workflow for targeted N-glycome profiling of glycopeptides derived from haptoglobin and vitronectin. Combining targeted glycopeptides with AFP significantly improved HCC detection compared with AFP alone. The optimal panel for all HCC versus cirrhosis (AFP + VTNC_169_A2G2F0S1 + VTNC_242_A3G3F2S2) achieved an AUC of 0.859 and 76.7% sensitivity at 90% specificity. For early-stage HCC, AFP + HP_184_A3G3F1S3 + VTNC_169_A2G2F0S1 yielded an AUC of 0.890 with 66.7% sensitivity at 1% specificity. A SHAP-selected Gaussian Naive Bayes model based on seven molecular/glycopeptide features, without demographic variables, further improved performance, achieving ROC-AUC values of 0.9985 in training and 1.0000 in independent testing cohorts, with accuracies of 98.1% and 100.0%, respectively.

1. Since plants interact with their environment through complex combinations of phytochemicals rather than metabolites in isolation, quantifying plant chemodiversity is of increasing interest. Although functionally relevant, structural disparity is mostly neglected in measures of chemodiversity because its integration relies on known compound identity, limiting broad application. 2. We established an approach deriving compound dissimilarity from UV-Vis spectra in HPLC-DAD datasets, evaluated how it relates to structure- or biosynthesis-based approaches, and examined whether it provides meaningful contributions to chemodiversity measures. For this, we applied it in an experiment full-factorially testing the effects of drought and herbivory on Solanum dulcamara leaf chemodiversity. 3. UV-Vis spectral dissimilarity aligned well with fMCS-based structural dissimilarity and reflected structural relationships within a set of standards. Incorporating it into chemodiversity analysis improved separation of the effects of drought and herbivory. Especially in the combination of both stresses, their distinct effects on different leaf metabolites were only fully reflected when accounting for compound disparity. 4. Hence, UV-Vis spectral dissimilarity captures chemically meaningful compound relatedness without requiring compound identity. Using it as a proxy for compound disparity adds a biologically relevant dimension to measures of chemodiversity with broader implications when assessing functional consequences of phytochemical diversity.

Metabolomics uses tandem mass spectrometry (MS/MS) data to gain structural insights of small molecules that play biological roles, generating datasets whose size and complexity demand systematic organisation. Molecular networking addresses this by representing MS/MS spectra as nodes and their pairwise similarity as edges, but its output is critically sensitive to user-defined parameters: similarity score cut-off, maximum component size, maximum links and minimum matching peaks. These parameters are routinely left at default values, which can either collapse interpretable molecular families into entangled "hairballs" or fragment them into disconnected singletons. In the absence of ground truth, no standardised framework exists to evaluate molecular networks or to assess whether their connections are robust to run-to-run variability present in metabolomic experiments. Here, we introduce arteMIS (Accelerated Ranking and Tuning using Multi-metric Interpretability across Scores), a framework for systematic parameter optimisation that uses Latin Hypercube Sampling to efficiently cover the four-dimensional parameter space and ranks candidate networks through a user-tuneable composite Z-score, combining topology- and chemistry-based metrics. This framework supports three complementary modes: global, seed, and target-class, adapting optimisation to fully unannotated datasets, curated subset of reference features or class-focused discovery, respectively. Benchmarking across four spectral libraries (~600 to ~13,000 spectra) and four scoring methods (Cosine, Modified Cosine, Spec2Vec, MS2DeepScore), we provide practical guidance for parameter selection as a function of scoring method and dataset size and show that optimal settings do not transfer between them. Top-ranked arteMIS configurations outperformed GNPS defaults in chemistry and topology metrics and produced networks with higher edge-stability under subsampling. Applied to actinobacteria and fungal samples, arteMIS rescued structurally meaningful families that remained fragmented under default settings. We conclude that arteMIS reframes molecular network construction from a default-driven step into a task-customisable optimisation.

Brook trout (Salvelinus fontinalis) exhibit strain-level variation in growth performance, environmental tolerance, and survival, yet the biochemical mechanisms underlying these differences remain poorly understood. We developed and applied a ratiometric biochemical framework integrating the pentose-phosphate pathway (PPP) and glutathione metabolism to characterize strain-specific hepatic metabolic organization in brook trout. Five strains reared under standardized conditions differed significantly in hepatic soluble protein density, glutathione pool size, total NADP(H) concentration, and activities of glucose-6-phosphate dehydrogenase (G6PDH), glutathione reductase (GR), and transketolase (TKT). These differences were not uniformly coordinated across pathways, demonstrating that metabolic phenotype cannot be inferred from individual biomarkers alone. Derived ratios describing oxidative-to-non-oxidative PPP capacity (G6PDH/TKT) and glutathione buffering relative to recycling capacity ((GSH+GSSG)/GR) resolved distinct patterns of metabolic allocation among strains. Despite shared ancestry, the Temiscamie (TEM) strain and its domestic x TEM hybrid (TXD) exhibited markedly divergent metabolic phenotypes, demonstrating that closely related strains can differ substantially in hepatic metabolic organization. Together, these findings identify relative allocation among interconnected metabolic pathways as an axis of physiologic diversity and establish a ratiometric approach for comparing metabolic organization across populations and species. Graphical abstractHepatic metabolic phenotypes of brook trout strains were characterized by integrating pentose phosphate pathway enzyme capacities, glutathione metabolism, NADP(H) availability, and soluble protein into a ratiometric framework. Ratios distinguish investment in oxidative versus non-oxidative PPP capacity (G6PDH/TKT), antioxidant buffering versus glutathione recycling capacity (total glutathione/GR), and hepatic protein density (soluble protein/liver mass), revealing distinct metabolic organization among strains. O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=88 SRC="FIGDIR/small/743818v1_ufig1.gif" ALT="Figure 1"> View larger version (25K): org.highwire.dtl.DTLVardef@1694676org.highwire.dtl.DTLVardef@90f2d4org.highwire.dtl.DTLVardef@365327org.highwire.dtl.DTLVardef@8d56ca_HPS_FORMAT_FIGEXP M_FIG C_FIG HighlightsO_LIA ratiometric framework was developed to characterize hepatic metabolic organization in brook trout C_LIO_LIGlutathione buffering and recycling capacity distinguish alternative redox phenotypes C_LIO_LIInvestment in oxidative and non-oxidative PPP capacity varies independently among strains C_LIO_LIG6PDH/TKT and total glutathione (GSH+GSSG)/GR reveal distinct metabolic phenotypes C_LIO_LIRatiometric indices provide a framework for interpreting redox metabolism and carbon allocation C_LI

Purpose: To evaluate whether fluorescence lifetime imaging ophthalmoscopy (FLIO) combined with deep learning can detect metabolic signatures for classification of type 2 diabetes mellitus (T2DM). Design: Cross-sectional analysis of participants included AI-READI dataset (version 3) with FLIO imaging and and hemoglobin A1c (HbA1c) measurement. Subjects: 1,783 participants from the AI-READI dataset (version 3) with HbA1c measurements and FLIO imaging scans (6,912 total): 671 normoglycemic, 726 prediabetic, and 386 diabetic. Methods: Mean fluorescence lifetime maps were generated using a center-of-mass approach and used as inputs to AI models. We trained convolutional neural networks (CNNs), ResNet-18, and XGBoost under three-class (normal, prediabetic, diabetic) and two binary (normal vs. impaired; normal vs. diabetic) classification schemes, using nested 5-fold cross-validation with participant-level grouping. Main Outcome Measures: Macro-averaged accuracy, F1 score, area under the receiver operating characteristic curve (AUROC), sensitivity, specificity, and positive predictive value (PPV). Results: Group-averaged lifetime maps demonstrated consistent spatial differences across glycemic groups, with progressively longer lifetimes from normal to diabetic participants. The CNN achieved the best overall performance in the 3-class classification (accuracy 0.41 +/- 0.03, F1 score 0.39 +/- 0.02, AUROC 0.58 +/- 0.02), compared to the random classifier for 3-class classification (AUROC = 0.50; accuracy = F1 = 0.33). ResNet-18 and XGBoost showed similar performance (AUROC 0.53-0.58). Confusion matrices revealed substantial overlap between classes, with frequent misclassification toward the prediabetes group. Binary reformulation (normal vs. diabetic) improved performance substantially, with the CNN resulting in AUROC 0.63 +/- 0.02 and XGBoost 0.67 +/- 0.07. Conclusions: FLIO-derived lifetime maps capture metabolic signals associated with glycemic status but yield modest classification performance with current AI models. These findings highlight both the potential and the challenges of using FLIO for early metabolic screening and monitoring, informing future development of clinically applicable imaging biomarkers.

Long-term in vivo transcriptomic time courses are costly, limiting assessment of chronic molecular responses from short studies. We developed a pharmacokinetics-informed transcriptomic ordinary differential equation model (PKT-ODE) that links an oral pharmacokinetic profile and Hill drug-effect function to first-order turnover of co-expression modules. The model was fitted to rat liver responses to fenofibrate at three doses in Open TG-GATEs through day 8. At the held-out day-29 endpoint, PKT-ODE achieved Pearson r = 0.960 and mean squared error (MSE) = 0.148. In this dataset, these values achieved lower prediction error and higher correlation than four statistical baselines and validation-selected linear and multilayer-perceptron transition models. Literature-curated peroxisome proliferator-activated receptor target genes occurred only in modules with positive fitted drug effects. These results provide a proof of concept for pharmacokinetics-informed transcriptomic extrapolation; cross-compound, cross-organ and alternative-regimen performance remain to be tested.