Published papers

The gut microbiota confers host protection against pathogen colonization early after infection. Several mechanisms underlying the protection have been described, but the contributions of nutrient competition versus direct inhibition are controversial. Using an ex vivo model of Salmonella growth in the mouse cecum with its indigenous microbes, we find that nutrient limitation and typical inhibitory factors alone cannot prevent pathogen growth. However, the addition of certain amino acids markedly reverses the microbiota's ability to suppress pathogen growth. Enhanced Salmonella colonization after antibiotic treatment is ablated by exclusion of dietary protein, which requires the presence of the microbiota. Thus, dietary protein and amino acids are important regulators of colonization resistance.

The Tsuruoka Metabolomics Cohort Study (TMCS) is an ongoing population-based cohort study being conducted in the rural area of Yamagata Prefecture, Japan. This study aimed to enhance the precision prevention of multi-factorial, complex diseases, including non-communicable and aging-associated diseases, by improving risk stratification and prediction measures. At baseline, 11,002 participants aged 35-74 years were recruited in Tsuruoka City, Yamagata Prefecture, Japan, between 2012 and 2015, with an ongoing follow-up survey. Participants underwent various measurements, examinations, tests, and questionnaires on their health, lifestyle, and social factors. This study uses an integrative approach with deep molecular profiling to identify potential biomarkers linked to phenotypes that underpin disease pathophysiology and provide better mechanistic insights into social health determinants. The TMCS incorporates multi-omics data, including genetic and metabolomic analyses of 10,933 participants, and comprehensive data collection ranging from physical, psychological, behavioral, and social to biological data. The metabolome is used as a phenotypic probe because it is sensitive to changes in physiological and external conditions. The TMCS focuses on collecting outcomes for cardiovascular disease, cancer incidence and mortality, disability and functional decline due to aging and disease sequelae, and the variation in health status within the body represented by omics analysis that lies between exposure and disease. It contains several sub-studies on aging, heated tobacco products, and women's health. This study is notable for its robust design, high participation rate (89%), and long-term repeated surveys. Moreover, it contributes to precision prevention in Japan and East Asia as a well-established multi-omics platform.

Plasma amino acids (AAs) have emerged as promising biomarkers for metabolic disorders, yet their causality remains unclear. We aimed to investigate the genetic determinants of AA levels in a cohort of 10,333 individuals and their causal effects on cardiometabolic traits using Mendelian randomization (MR). Plasma levels of 20 AAs were quantified using capillary electrophoresis mass spectrometry. Genome-wide association studies were conducted using BOLT-LMM and heritability estimation via LDSC analysis. Causal effects of AAs on 11 cardiometabolic traits were examined using two-sample MR analyses. We identified 85 locus-metabolite associations across 43 genes for 18 AAs, including 44 novel loci linked to metabolic genes. Heritability for AAs was estimated at 16%. MR analysis demonstrated cystine to positively associate with systolic blood pressure (SBP) (β = 0.056, SE = 0.010), while serine indicated protective effects on SBP (β = - 0.040, SE = 0.011), diastolic BP (β = - 0.044, SE = 0.010), and coronary artery disease (odds ratio 0.888, SE = 0.028). We identified potentially novel genetic loci associated with AA levels and demonstrated robust causal associations between several AAs and cardiometabolic traits. These findings reinforce the importance of AAs as potential biomarkers and therapeutic targets in cardiometabolic health.

There is currently a global epidemic of obesity and obesity-related diseases such as type 2 diabetes due to decreased physical activity, excessive food intake, and/or genetic predisposition. The Hippo-YAP1 pathway has attracted attention as a potential therapeutic target because YAP1/TAZ activation in murine immature adipocytes in vitro suppresses their differentiation by inhibiting PPARγ activity. However, the role of YAP1 activation in mature adipocytes in vivo remains unclear. MOB1, whose expression is increased in obesity, is the hub of the Hippo core molecule complex and negatively regulates YAP1/TAZ activation. Therefore, we generated aMob1DKO mutant mice, which feature deficiency of Mob1a/b specifically in mature adipocytes. Compared to controls, aMob1DKO mice subjected to a high-fat diet showed beneficial changes consistent with resistance to diet-induced obesity. The mutants exhibited increases in basal lipolysis, "beiging," and energy expenditure, as well as suppression of ROS production and inflammation in white adipose tissue. Insulin sensitivity and glucose tolerance were improved, and ectopic fat accumulation was reduced. Most of these changes were dependent on the YAP1 activation observed in mature white adipose tissue of aMob1DKO mice. FGF21, which improves lipid metabolism, was upregulated directly via YAP1 activation, and many of the phenotypes seen in aMob1DKO mice were also dependent on FGF21. Thus, the aMob1DKO mouse is an interesting model for the study of the metabolic effects of diet-induced obesity and protection against diabetes. Our work suggests that a YAP1-FGF21 axis exists in adipocytes that may be a potential therapeutic target for obesity.

Adaptation to starvation is a multimolecular and temporally ordered process. We sought to elucidate how the healthy liver regulates various molecules in a temporally ordered manner during starvation and how obesity disrupts this process. We used multiomic data collected from the plasma and livers of wild-type and leptin-deficient obese (ob/ob) mice at multiple time points during starvation to construct a starvation-responsive metabolic network that included responsive molecules and their regulatory relationships. Analysis of the network structure showed that in wild-type mice, the key molecules for energy homeostasis, ATP and AMP, acted as hub molecules to regulate various metabolic reactions in the network. Although neither ATP nor AMP was responsive to starvation in ob/ob mice, the structural properties of the network were maintained. In wild-type mice, the molecules in the network were temporally ordered through metabolic processes coordinated by hub molecules, including ATP and AMP, and were positively or negatively coregulated. By contrast, both temporal order and coregulation were disrupted in ob/ob mice. These results suggest that the metabolic network that responds to starvation was structurally robust but temporally disrupted by the obesity-associated loss of responsiveness of the hub molecules. In addition, we propose how obesity alters the response to intermittent fasting.

Although mutations in human patatin-like phospholipase PNPLA6 are associated with hereditary retinal degenerative diseases, its mechanistic action in the retina is poorly understood. Here, we uncover the molecular mechanism by which PNPLA6 dysfunction disturbs retinal homeostasis and visual function. PNPLA6, by acting as a phospholipase B, regulates choline mobilization from phosphatidylcholine and subsequent choline turnover for phosphatidylcholine regeneration in retinal pigment epithelial cells. PNPLA6-driven choline is supplied from retinal pigment epithelial cells to adjacent photoreceptor cells to support their survival. Inhibition of this pathway results in abnormal morphology, proliferation, metabolism, and functions of retinal pigment epithelial and photoreceptor cells, and mice with retina-specific PNPLA6 deletion exhibit retinitis pigmentosa-like retinal degeneration. Notably, these abnormalities are entirely rescued by choline supplementation. Thus, PNPLA6 plays an essential role in retinal homeostasis by controlling choline availability for phospholipid recycling and provide a framework for the development of an ophthalmic drug target for retinal degeneration.

Recent advances in cancer therapy have underscored the importance of targeting specific metabolic pathways. In this study, we propose a precision nutrition approach aimed at lysosomal function in glioblastoma multiforme (GBM). Using patient-derived GBM cells, we identify lysosomal activity as a unique metabolic biomarker of tumorigenesis, controlling the efficacy of temozolomide (TMZ), a standard GBM therapy. Employing combined analyses of clinical patient samples and xenograft models, we further elucidate the pivotal role of Transcription Factor Binding To IGHM Enhancer 3 (TFE3), a master regulator of lysosomal biogenesis, in modulating malignant properties, particularly TMZ tolerance, by regulating peroxisome proliferator-activated receptor-gamma coactivator 1-alpha (PGC1α)-mediated mitochondrial activity. Notably, we find that lysine protects GBM cells from lysosomal stress by counteracting arginine's effects on nitric oxide production. The lysine restriction mimetic, homoarginine administration, significantly enhances the efficacy of anticancer therapies through lysosomal dysfunction. This study underscores the critical role of lysosomal function modulated by amino acid metabolism in GBM pathogenesis and treatment.

Background/aim: 18F-fluorodeoxyglucose positron emission tomography/computed tomography (FDG-PET/CT) is reportedly associated with the malignant potential of cancer. This study aimed to evaluate the association between FDG accumulation and tumor metabolism in pancreatic ductal adenocarcinoma (PDAC). Patients and methods: A prognostic analysis of data from 131 patients with PDAC who underwent FDG-PET/CT before curative-intent pancreatic surgery was performed. Capillary electrophoresis-mass spectrometry (CE-MS) was used to analyze the metabolome of tumor and non-neoplastic pancreas from 80 patients. These patients were divided into two groups: low SUVmax group (SUVmax <6.09) and high SUVmax group (SUVmax ≥6.09). Results: Carbohydrate antigen 19-9 (CA19-9), maximum standardized uptake value (SUVmax) of PET, N stage, and postoperative chemotherapy were identified as significant prognostic factors by univariate analysis. SUVmax emerged as an independent prognostic factor for overall survival [hazard ratio (HR)=1.88, p<0.05] and disease-free survival (HR=2.01, p<0.05) in multivariate analysis. Metabolic analyses confirmed that 43 metabolites significantly differed depending on the accumulation of SUV in tumors. Metabolites involved in the removal of reactive oxygen species (e.g., hypotaurine, glutathione, Met), treatment resistance (UDP-N-acetylglucosamine), and proliferation (e.g., choline, leucine, isoleucine) were increased in the high SUVmax group. Conclusion: FDG accumulation is an important independent prognostic factor reflecting tumor activity associated with metabolic changes in cancer cells.

Metabolic abnormalities play a pivotal role in various pathological conditions, necessitating the quantification of specific metabolites for diagnosis. While mass spectrometry remains the primary method for metabolite measurement, its limited throughput underscores the need for biosensors capable of rapid detection. Previously, we reported that pillar[6]arene with 12 carboxylate groups (P6AC) forms host-guest complexes with 1-methylnicotinamide (1-MNA), which is produced in vivo by nicotinamide N-methyltransferase (NNMT). P6AC acts as a biosensor by measuring the fluorescence quenching caused by photoinduced electron transfer upon 1-MNA binding. However, the low sensitivity of P6AC makes it impractical for detecting 1-MNA in unpurified biological samples. In this study, we found that P6A with 12 sulfonate groups (P6AS) is a specific and potent supramolecular host for 1-MNA interactions even in biological samples. The 1-MNA binding affinity of P6AS in water was found to be (5.68 ± 1.02) × 106 M-1, which is approximately 700-fold higher than that of P6AC. Moreover, the 1-MNA detection limit of P6AS was determined to be 2.84 × 10-7 M, which is substantially lower than that of P6AC. Direct addition of P6AS to culture medium was sufficient to quantify 1-MNA produced by cancer cells. Furthermore, this sensor was able to specifically detect 1-MNA even in unpurified human urine. P6AS therefore enables rapid and high-throughput quantification of 1-MNA, and further improvement of our strategy will contribute to the establishment of high-throughput screening of NNMT inhibitors, diagnosis of liver diseases, and imaging of human cancer cells in vivo.

Hepatic glucose metabolism serves dual purposes: maintaining glucose homeostasis and converting glucose into energy sources; however, the underlying mechanisms are unclear. We quantitatively measured liver metabolites, gene expression, and phosphorylated insulin signaling molecules in mice orally administered varying doses of glucose, and constructed a transomic network. Rapid phosphorylation of insulin signaling molecules in response to glucose intake was observed, in contrast to the more gradual changes in gene expression. Glycolytic and gluconeogenic metabolites and expression of genes involved in glucose metabolism including glucose-6-phosphate, G6pc, and Pck1, demonstrated high glucose dose sensitivity. Whereas, glucokinase expression and glycogen accumulation showed low glucose dose sensitivity. During the early phase after glucose intake, metabolic flux was geared towards glucose homeostasis regardless of the glucose dose but shifted towards energy conversion during the late phase at higher glucose doses. Our research provides a comprehensive view of time- and dose-dependent selective glucose metabolism.

Epithelial-to-mesenchymal transition (EMT) contributes to the poor prognosis of patients with cancer by promoting distant metastasis and anti-cancer drug resistance. Several distinct metabolic alterations have been identified as key EMT phenotypes. In the present study, we further characterize the role of transforming growth factor-β (TGF-β)-induced EMT in non-small-cell lung cancer. Our study revealed that TGF-β plays a role in EMT functions by upregulation of cytidine 5'-triphosphate synthetase 1 (CTPS), a vital enzyme for CTP biosynthesis in the pyrimidine metabolic pathway. Both knockdown and enzymatic inhibition of CTPS reduced TGF-β-induced changes in EMT marker expression, chemoresistance and migration in vitro. Moreover, CTPS knockdown counteracted the TGF-β-mediated downregulation of UDP-glucuronate, glutarate, creatine, taurine and nicotinamide. These findings indicate that CTPS plays a multifaceted role in EMT metabolism, which is crucial for the malignant transformation of cancer through EMT, and underline its potential as a promising therapeutic target for preventing drug resistance and metastasis in non-small-cell lung cancer.

Gut bacteria are linked to neurodegenerative diseases but the risk factors beyond microbiota composition are limited. Here we used a pre-clinical model of multiple sclerosis (MS), experimental autoimmune encephalomyelitis (EAE), to identify microbial risk factors. Mice with different genotypes and complex microbiotas or six combinations of a synthetic human microbiota were analysed, resulting in varying probabilities of severe neuroinflammation. However, the presence or relative abundances of suspected microbial risk factors failed to predict disease severity. Akkermansia muciniphila, often associated with MS, exhibited variable associations with EAE severity depending on the background microbiota. Significant inter-individual disease course variations were observed among mice harbouring the same microbiota. Evaluation of microbial functional characteristics and host immune responses demonstrated that the immunoglobulin A coating index of certain bacteria before disease onset is a robust individualized predictor of disease development. Our study highlights the need to consider microbial community networks and host-specific bidirectional interactions when aiming to predict severity of neuroinflammation.

Accumulating evidence demonstrates clear correlation between the gut microbiota and sporadic colorectal cancer (CRC). Despite this, there is limited understanding of the association between the gut microbiota and CRC in Lynch Syndrome (LS), a hereditary type of CRC. Here, we analyzed fecal shotgun metagenomic and targeted metabolomic of 71 Japanese LS subjects. A previously published Japanese sporadic CRC cohort, which includes non-LS controls, was utilized as a non-LS cohort (n = 437). LS subjects exhibited reduced microbial diversity and low-Faecalibacterium enterotypes compared to non-LS. Patients with LS-CRC had higher levels of Fusobacterium nucleatum and fap2. Differential fecal metabolites and functional genes suggest heightened degradation of lysine and arginine in LS-CRC. A comparison between LS and non-LS subjects prior to adenoma formation revealed distinct fecal metabolites of LS subjects. These findings suggest that the gut microbiota plays a more responsive role in CRC tumorigenesis in patients with LS than those without LS.

Dysregulation of liver metabolism associated with obesity during feeding and fasting leads to the breakdown of metabolic homeostasis. However, the underlying mechanism remains unknown. Here, we measured multi-omics data in the liver of wild-type and leptin-deficient obese (ob/ob) mice at ad libitum feeding and constructed a differential regulatory trans-omic network of metabolic reactions. We compared the trans-omic network at feeding with that at 16 h fasting constructed in our previous study. Intermediate metabolites in glycolytic and nucleotide metabolism decreased in ob/ob mice at feeding but increased at fasting. Allosteric regulation reversely shifted between feeding and fasting, generally showing activation at feeding while inhibition at fasting in ob/ob mice. Transcriptional regulation was similar between feeding and fasting, generally showing inhibiting transcription factor regulations and activating enzyme protein regulations in ob/ob mice. The opposite metabolic dysregulation between feeding and fasting characterizes breakdown of metabolic homeostasis associated with obesity.

Background: The application of metabolomics-based profiles in environmental epidemiological studies is a promising approach to refine the process of health risk assessment. We aimed to identify potential metabolomics-based profiles in urine and plasma for the detection of relatively low-level cadmium (Cd) exposure in large population-based studies. Method: We analyzed 123 urinary metabolites and 94 plasma metabolites detected in fasting urine and plasma samples collected from 1,412 men and 2,022 women involved in the Tsuruoka Metabolomics Cohort Study. Regression analysis was performed for urinary N-acetyl-beta-D-glucosaminidase (NAG), plasma, and urinary metabolites as dependent variables, and urinary Cd (U-Cd, quartile) as an independent variable. The multivariable regression model included age, gender, systolic blood pressure, smoking, rice intake, BMI, glycated hemoglobin, low-density lipoprotein cholesterol, alcohol consumption, physical activity, educational history, dietary energy intake, urinary Na/K ratio, and uric acid. Pathway-network analysis was carried out to visualize the metabolite networks linked to Cd exposure. Result: Urinary NAG was positively associated with U-Cd, but not at lower concentrations (Q2). Among urinary metabolites in the total population, 45 metabolites showed associations with U-Cd in the unadjusted and adjusted models after adjusting for the multiplicity of comparison with FDR. There were 12 urinary metabolites which showed consistent associations between Cd exposure from Q2 to Q4. Among plasma metabolites, six cations and one anion were positively associated with U-Cd, whereas alanine, creatinine, and isoleucine were negatively associated with U-Cd. Our results were robust by statistical adjustment of various confounders. Pathway-network analysis revealed metabolites and upstream regulator changes associated with mitochondria (ACACB, UCP2, and metabolites related to the TCA cycle). Conclusion: These results suggested that U-Cd was associated with metabolites related to upstream mitochondrial dysfunction in a dose-dependent manner. Our data will help develop environmental Cd exposure profiles for human populations.

Small-cell lung cancer (SCLC) has a poor prognosis, emphasizing the necessity for developing new therapies. The de novo synthesis pathway of purine nucleotides, which is involved in the malignant growth of SCLC, has emerged as a novel therapeutic target. Purine nucleotides are supplied by two pathways: de novo and salvage. However, the role of the salvage pathway in SCLC and the differences in utilization and crosstalk between the two pathways remain largely unclear. Here, we found that deletion of the HPRT1 gene, which codes for the rate-limiting enzyme of the purine salvage pathway, significantly suppressed tumor growth in vivo in several SCLC cells. We also demonstrated that HPRT1 expression confers resistance to lemetrexol (LMX), an inhibitor of the purine de novo pathway. Interestingly, HPRT1-knockout had less effect on SCLC SBC-5 cells, which are more sensitive to LMX than other SCLC cell lines, suggesting that a preference for either the purine de novo or salvage pathway occurs in SCLC. Furthermore, metabolome analysis of HPRT1-knockout cells revealed increased intermediates in the pentose phosphate pathway and elevated metabolic flux in the purine de novo pathway, indicating compensated metabolism between the de novo and salvage pathways in purine nucleotide biosynthesis. These results suggest that HPRT1 has therapeutic implications in SCLC and provide fundamental insights into the regulation of purine nucleotide biosynthesis.

Nicotinamide phosphoribosyltransferase (NAMPT) plays a major role in NAD biosynthesis in many cancers and is an attractive potential cancer target. However, factors dictating therapeutic efficacy of NAMPT inhibitors (NAMPTi) are unclear. We report that neuroendocrine phenotypes predict lung and prostate carcinoma vulnerability to NAMPTi, and that NAMPTi therapy against those cancers is enhanced by dietary modification. Neuroendocrine differentiation of tumor cells is associated with down-regulation of genes relevant to quinolinate phosphoribosyltransferase-dependent de novo NAD synthesis, promoting NAMPTi susceptibility in vitro. We also report that circulating nicotinic acid riboside (NAR), a non-canonical niacin absent in culture media, antagonizes NAMPTi efficacy as it fuels NAMPT-independent but nicotinamide riboside kinase 1-dependent NAD synthesis in tumors. In mouse transplantation models, depleting blood NAR by nutritional or genetic manipulations is synthetic lethal to tumors when combined with NAMPTi. Our findings provide a rationale for simultaneous targeting of NAR metabolism and NAMPT therapeutically in neuroendocrine carcinoma.

Studies examining long-term longitudinal metabolomic data and their reliability in large-scale populations are limited. Therefore, we aimed to evaluate the reliability of repeated measurements of plasma metabolites in a prospective cohort setting and to explore intra-individual concentration changes at three time points over a 6-year period. The study participants included 2999 individuals (1317 men and 1682 women) from the Tsuruoka Metabolomics Cohort Study, who participated in all three surveys-at baseline, 3 years, and 6 years. In each survey, 94 plasma metabolites were quantified for each individual and quality control (QC) sample. The coefficients of variation of QC, intraclass correlation coefficients, and change rates of QC were calculated for each metabolite, and their reliability was classified into three categories: excellent, fair to good, and poor. Seventy-six percent (71/94) of metabolites were classified as fair to good or better. Of the 39 metabolites grouped as excellent, 29 (74%) in men and 26 (67%) in women showed significant intra-individual changes over 6 years. Overall, our study demonstrated a high degree of reliability for repeated metabolome measurements. Many highly reliable metabolites showed significant changes over the 6-year period, suggesting that repeated longitudinal metabolome measurements are useful for epidemiological studies.

Objective: A significant number of patients experience early recurrence after surgical resection for pancreatic ductal adenocarcinoma (PDAC), negating the benefit of surgery. The present study conducted clinicopathologic and metabolomic analyses to explore the factors associated with the early recurrence of PDAC. Materials and methods: Patients who underwent pancreatectomy for PDAC at Kagawa University Hospital between 2011 and 2020 were enrolled. Tissue samples of PDAC and nonneoplastic pancreas were collected and frozen immediately after resection. Charged metabolites were quantified by capillary electrophoresis-mass spectrometry. Patients who relapsed within 1 year were defined as the early recurrence group. Results: Frozen tumor tissue and nonneoplastic pancreas were collected from 79 patients. The clinicopathologic analysis identified 11 predictive factors, including preoperative carbohydrate antigen 19-9 levels. The metabolomic analysis revealed that only hypotaurine was a significant risk factor for early recurrence. A multivariate analysis, including clinical and metabolic factors, showed that carbohydrate antigen 19-9 and hypotaurine were independent risk factors for early recurrence (P = 0.045 and P = 0.049, respectively). The recurrence-free survival rate 1 year after surgery with both risk factors was only 25%. Conclusions: Our results suggested that tumor hypotaurine is a potential metabolite associated with early recurrence. Carbohydrate antigen 19-9 and hypotaurine showed a vital utility for predicting early recurrence.

Malignant tumors are characterized by a hypoxic microenvironment, and metabolic reprogramming is necessary to ensure energy production and oxidative stress resistance. Although the microenvironmental properties of tumors vary under acute and chronic hypoxia, studies on chronic hypoxia-induced metabolic changes are limited. In the present study, we performed a comprehensive metabolic analysis in a chronic hypoxia model using colorectal cancer (CRC) organoids, and identified an amino acid supply system through the γ-glutamyl cycle, a glutathione recycling pathway. We analyzed the metabolic changes caused by hypoxia over time and observed that chronic hypoxia resulted in an increase in 5-oxoproline and a decrease in oxidized glutathione (GSSG) compared to acute hypoxia. These findings suggest that chronic hypoxia induces metabolic changes in the γ-glutamyl cycle. Moreover, inhibition of the γ-glutamyl cycle via γ-glutamyl cyclotransferase (GGCT) and γ-glutamyl transferase 1 (GGT1) knockdown significantly reversed chronic hypoxia-induced upregulation of 5-oxoproline and several amino acids. Notably, GGT1 knockdown downregulated the intracellular levels of γ-glutamyl amino acids. Conclusively, these results indicate that the γ-glutamyl cycle serves as an amino acid supply system in CRC under chronic hypoxia, which provides fresh insight into cancer metabolism under chronic hypoxia.