Atlas of Circadian Metabolism Reveals System-wide Coordination and Communication between Clocks.

Metabolic diseases are often characterized by circadian misalignment in different tissues, yet how altered coordination and communication among tissue clocks relate to specific pathogenic mechanisms remains largely unknown. Applying an integrated systems biology approach, we performed 24-hr metabolomics profiling of eight mouse tissues simultaneously. We present a temporal and spatial atlas of circadian metabolism in the context of systemic energy balance and under chronic nutrient stress (high-fat diet [HFD]). Comparative analysis reveals how the repertoires of tissue metabolism are linked and gated to specific temporal windows and how this highly specialized communication and coherence among tissue clocks is rewired by nutrient challenge. Overall, we illustrate how dynamic metabolic relationships can be reconstructed across time and space and how integration of circadian metabolomics data from multiple tissues can improve our understanding of health and disease.

Chemical Hybridization of Glucagon and Thyroid Hormone Optimizes Therapeutic Impact for Metabolic Disease.

Glucagon and thyroid hormone (T3) exhibit therapeutic potential for metabolic disease but also exhibit undesired effects. We achieved synergistic effects of these two hormones and mitigation of their adverse effects by engineering chemical conjugates enabling delivery of both activities within one precisely targeted molecule. Coordinated glucagon and T3 actions synergize to correct hyperlipidemia, steatohepatitis, atherosclerosis, glucose intolerance, and obesity in metabolically compromised mice. We demonstrate that each hormonal constituent mutually enriches cellular processes in hepatocytes and adipocytes via enhanced hepatic cholesterol metabolism and white fat browning. Synchronized signaling driven by glucagon and T3 reciprocally minimizes the inherent harmful effects of each hormone. Liver-directed T3 action offsets the diabetogenic liability of glucagon, and glucagon-mediated delivery spares the cardiovascular system from adverse T3 action. Our findings support the therapeutic utility of integrating these hormones into a single molecular entity that offers unique potential for treatment of obesity, type 2 diabetes, and cardiovascular disease.

Microglial phagolysosome dysfunction and altered neural communication amplify phenotypic severity in Prader-Willi Syndrome with larger deletion.

Prader-Willi Syndrome (PWS) is a rare neurodevelopmental disorder of genetic etiology, characterized by paternal deletion of genes located at chromosome 15 in 70% of cases. Two distinct genetic subtypes of PWS deletions are characterized, where type I (PWS T1) carries four extra haploinsufficient genes compared to type II (PWS T2). PWS T1 individuals display more pronounced physiological and cognitive abnormalities than PWS T2, yet the exact neuropathological mechanisms behind these differences remain unclear. Our study employed postmortem hypothalamic tissues from PWS T1 and T2 individuals, conducting transcriptomic analyses and cell-specific protein profiling in white matter, neurons, and glial cells to unravel the cellular and molecular basis of phenotypic severity in PWS sub-genotypes. In PWS T1, key pathways for cell structure, integrity, and neuronal communication are notably diminished, while glymphatic system activity is heightened compared to PWS T2. The microglial defect in PWS T1 appears to stem from gene haploinsufficiency, as global and myeloid-specific Cyfip1 haploinsufficiency in murine models demonstrated. Our findings emphasize microglial phagolysosome dysfunction and altered neural communication as crucial contributors to the severity of PWS T1's phenotype.

Skeletal muscle and intermuscular adipose tissue gene expression profiling identifies new biomarkers with prognostic significance for insulin resistance progression and intervention response.

AIMS/HYPOTHESIS: Although insulin resistance often leads to type 2 diabetes mellitus, its early stages are often unrecognised, thus reducing the probability of successful prevention and intervention. Moreover, treatment efficacy is affected by the genetics of the individual. We used gene expression profiles from a cross-sectional study to identify potential candidate genes for the prediction of diabetes risk and intervention response. METHODS: Using a multivariate regression model, we linked gene expression profiles of human skeletal muscle and intermuscular adipose tissue (IMAT) to fasting glucose levels and glucose infusion rate. Based on the expression patterns of the top predictive genes, we characterised and compared individual gene expression with clinical classifications using k-nearest neighbour clustering. The predictive potential of the candidate genes identified was validated using muscle gene expression data from a longitudinal intervention study. RESULTS: We found that genes with a strong association with clinical measures clustered into three distinct expression patterns. Their predictive values for insulin resistance varied substantially between skeletal muscle and IMAT. Moreover, we discovered that individual gene expression-based classifications may differ from classifications based predominantly on clinical variables, indicating that participant stratification may be imprecise if only clinical variables are used for classification. Of the 15 top candidate genes, ST3GAL2, AASS, ARF1 and the transcription factor SIN3A are novel candidates for predicting a refined diabetes risk and intervention response. CONCLUSION/INTERPRETATION: Our results confirm that disease progression and successful intervention depend on individual gene expression states. We anticipate that our findings may lead to a better understanding and prediction of individual diabetes risk and may help to develop individualised intervention strategies.

Diet triggers specific responses of hypothalamic astrocytes in time and region dependent manner.

Hypothalamic astrocytes are particularly affected by energy-dense food consumption. How the anatomical location of these glial cells and their spatial molecular distribution in the arcuate nucleus of the hypothalamus (ARC) determine the cellular response to a high caloric diet remains unclear. In this study, we investigated their distinctive molecular responses following exposure to a high-fat high-sugar (HFHS) diet, specifically in the ARC. Using RNA sequencing and proteomics, we showed that astrocytes have a distinct transcriptomic and proteomic profile dependent on their anatomical location, with a major proteomic reprogramming in hypothalamic astrocytes. By ARC single-cell sequencing, we observed that a HFHS diet dictates time- and cell- specific transcriptomic responses, revealing that astrocytes have the most distinct regulatory pattern compared to other cell types. Lastly, we topographically and molecularly characterized astrocytes expressing glial fibrillary acidic protein and/or aldehyde dehydrogenase 1 family member L1 in the ARC, of which the abundance was significantly increased, as well as the alteration in their spatial and molecular profiles, with a HFHS diet. Together, our results provide a detailed multi-omics view on the spatial and temporal changes of astrocytes particularly in the ARC during different time points of adaptation to a high calorie diet.

Transcriptional programming of lipid and amino acid metabolism by the skeletal muscle circadian clock.

Circadian clocks are fundamental physiological regulators of energy homeostasis, but direct transcriptional targets of the muscle clock machinery are unknown. To understand how the muscle clock directs rhythmic metabolism, we determined genome-wide binding of the master clock regulators brain and muscle ARNT-like protein 1 (BMAL1) and REV-ERBα in murine muscles. Integrating occupancy with 24-hr gene expression and metabolomics after muscle-specific loss of BMAL1 and REV-ERBα, here we unravel novel molecular mechanisms connecting muscle clock function to daily cycles of lipid and protein metabolism. Validating BMAL1 and REV-ERBα targets using luciferase assays and in vivo rescue, we demonstrate how a major role of the muscle clock is to promote diurnal cycles of neutral lipid storage while coordinately inhibiting lipid and protein catabolism prior to awakening. This occurs by BMAL1-dependent activation of Dgat2 and REV-ERBα-dependent repression of major targets involved in lipid metabolism and protein turnover (MuRF-1, Atrogin-1). Accordingly, muscle-specific loss of BMAL1 is associated with metabolic inefficiency, impaired muscle triglyceride biosynthesis, and accumulation of bioactive lipids and amino acids. Taken together, our data provide a comprehensive overview of how genomic binding of BMAL1 and REV-ERBα is related to temporal changes in gene expression and metabolite fluctuations.

Intermuscular adipose tissue directly modulates skeletal muscle insulin sensitivity in humans.

Intermuscular adipose tissue (IMAT) is negatively related to insulin sensitivity, but a causal role of IMAT in the development of insulin resistance is unknown. IMAT was sampled in humans to test for the ability to induce insulin resistance in vitro and characterize gene expression to uncover how IMAT may promote skeletal muscle insulin resistance. Human primary muscle cells were incubated with conditioned media from IMAT, visceral (VAT), or subcutaneous adipose tissue (SAT) to evaluate changes in insulin sensitivity. RNAseq analysis was performed on IMAT with gene expression compared with skeletal muscle and SAT, and relationships to insulin sensitivity were determined in men and women spanning a wide range of insulin sensitivity measured by hyperinsulinemic-euglycemic clamp. Conditioned media from IMAT and VAT decreased insulin sensitivity similarly compared with SAT. Multidimensional scaling analysis revealed distinct gene expression patterns in IMAT compared with SAT and muscle. Pathway analysis revealed that IMAT expression of genes in insulin signaling, oxidative phosphorylation, and peroxisomal metabolism related positively to donor insulin sensitivity, whereas expression of macrophage markers, inflammatory cytokines, and secreted extracellular matrix proteins were negatively related to insulin sensitivity. Perilipin 5 gene expression suggested greater IMAT lipolysis in insulin-resistant individuals. Combined, these data show that factors secreted from IMAT modulate muscle insulin sensitivity, possibly via secretion of inflammatory cytokines and extracellular matrix proteins, and by increasing local FFA concentration in humans. These data suggest IMAT may be an important regulator of skeletal muscle insulin sensitivity and could be a novel therapeutic target for skeletal muscle insulin resistance.

Epstein-Barr virus microRNAs reduce immune surveillance by virus-specific CD8+ T cells.

Infection with Epstein-Barr virus (EBV) affects most humans worldwide and persists life-long in the presence of robust virus-specific T-cell responses. In both immunocompromised and some immunocompetent people, EBV causes several cancers and lymphoproliferative diseases. EBV transforms B cells in vitro and encodes at least 44 microRNAs (miRNAs), most of which are expressed in EBV-transformed B cells, but their functions are largely unknown. Recently, we showed that EBV miRNAs inhibit CD4+ T-cell responses to infected B cells by targeting IL-12, MHC class II, and lysosomal proteases. Here we investigated whether EBV miRNAs also counteract surveillance by CD8+ T cells. We have found that EBV miRNAs strongly inhibit recognition and killing of infected B cells by EBV-specific CD8+ T cells through multiple mechanisms. EBV miRNAs directly target the peptide transporter subunit TAP2 and reduce levels of the TAP1 subunit, MHC class I molecules, and EBNA1, a protein expressed in most forms of EBV latency and a target of EBV-specific CD8+ T cells. Moreover, miRNA-mediated down-regulation of the cytokine IL-12 decreases the recognition of infected cells by EBV-specific CD8+ T cells. Thus, EBV miRNAs use multiple, distinct pathways, allowing the virus to evade surveillance not only by CD4+ but also by antiviral CD8+ T cells.

miR-335 promotes mesendodermal lineage segregation and shapes a transcription factor gradient in the endoderm.

Transcription factors (TFs) pattern developing tissues and determine cell fates; however, how spatio-temporal TF gradients are generated is ill defined. Here we show that miR-335 fine-tunes TF gradients in the endoderm and promotes mesendodermal lineage segregation. Initially, we identified miR-335 as a regulated intronic miRNA in differentiating embryonic stem cells (ESCs). miR-335 is encoded in the mesoderm-specific transcript (Mest) and targets the 3'-UTRs of the endoderm-determining TFs Foxa2 and Sox17. Mest and miR-335 are co-expressed and highly accumulate in the mesoderm, but are transiently expressed in endoderm progenitors. Overexpression of miR-335 does not affect initial mesendoderm induction, but blocks Foxa2- and Sox17-mediated endoderm differentiation in ESCs and ESC-derived embryos. Conversely, inhibition of miR-335 activity leads to increased Foxa2 and Sox17 protein accumulation and endoderm formation. Mathematical modeling predicts that transient miR-335 expression in endoderm progenitors shapes a TF gradient in the endoderm, which we confirm by functional studies in vivo. Taken together, our results suggest that miR-335 targets endoderm TFs for spatio-temporal gradient formation in the endoderm and to stabilize lineage decisions during mesendoderm formation.

Time-matched analysis of DNA adduct formation and early gene expression as predictive tool for renal carcinogenesis in methylazoxymethanol acetate treated Eker rats.

Genotoxic carcinogens pose great hazard to human health. Uncertainty of current risk assessment strategies and long latency periods between first carcinogen exposure and diagnosis of tumors have raised interest in predictive biomarkers. Initial DNA adduct formation is a necessary step for genotoxin induced carcinogenesis. However, as DNA adducts not always translate into tumorigenesis, their predictive value is limited. Here we hypothesize that the combined analysis of pro-mutagenic DNA adducts along with time-matched gene expression changes could serve as a superior prediction tool for genotoxic carcinogenesis. Eker rats, heterozygous for the tuberous sclerosis (Tsc2) tumor suppressor gene and thus highly susceptible towards genotoxic renal carcinogens, were continuously treated with the DNA alkylating carcinogen methylazoxymethanol acetate (MAMAc). Two weeks of MAMAc treatment resulted in a time-dependent increase of O6-methylguanine and N7-methylguanine adducts in the kidney cortex, which was however not reflected by significant expression changes of cyto-protective genes involved in DNA repair, cell cycle arrest or apoptosis. Instead, we found a transcriptional regulation of genes involved in the tumor-related MAPK, FoxO and TGF-beta pathways. Continuous MAMAc treatment for up to 6 months resulted in a mild but significant increase of cancerous lesions. In summary, the combined analysis of DNA adducts and early gene expression changes could serve as a suitable predictive tool for genotoxicant-induced carcinogenesis.

K1 and K15 of Kaposi's Sarcoma-Associated Herpesvirus Are Partial Functional Homologues of Latent Membrane Protein 2A of Epstein-Barr Virus.

UNLABELLED: The human herpesviruses Epstein-Barr virus (EBV) and Kaposi's sarcoma-associated herpesvirus (KSHV) are associated with Hodgkin's lymphoma (HL) and Primary effusion lymphomas (PEL), respectively, which are B cell malignancies that originate from germinal center B cells. PEL cells but also a quarter of EBV-positive HL tumor cells do not express the genuine B cell receptor (BCR), a situation incompatible with survival of normal B cells. EBV encodes LMP2A, one of EBV's viral latent membrane proteins, which likely replaces the BCR's survival signaling in HL. Whether KSHV encodes a viral BCR mimic that contributes to oncogenesis is not known because an experimental model of KSHV-mediated B cell transformation is lacking. We addressed this uncertainty with mutant EBVs encoding the KSHV genes K1 or K15 in lieu of LMP2A and infected primary BCR-negative (BCR(-)) human B cells with them. We confirmed that the survival of BCR(-) B cells and their proliferation depended on an active LMP2A signal. Like LMP2A, the expression of K1 and K15 led to the survival of BCR(-) B cells prone to apoptosis, supported their proliferation, and regulated a similar set of cellular target genes. K1 and K15 encoded proteins appear to have noncomplementing, redundant functions in this model, but our findings suggest that both KSHV proteins can replace LMP2A's key activities contributing to the survival, activation and proliferation of BCR(-) PEL cells in vivo. IMPORTANCE: Several herpesviruses encode oncogenes that are receptor-like proteins. Often, they are constitutively active providing important functions to the latently infected cells. LMP2A of Epstein-Barr virus (EBV) is such a receptor that mimics an activated B cell receptor, BCR. K1 and K15, related receptors of Kaposi's sarcoma-associated herpesvirus (KSHV) expressed in virus-associated tumors, have less obvious functions. We found in infection experiments that both viral receptors of KSHV can replace LMP2A and deliver functions similar to the endogenous BCR. K1, K15, and LMP2A also control the expression of a related set of cellular genes in primary human B cells, the target cells of EBV and KSHV. The observed phenotypes, as well as the known characteristics of these genes, argue for their contributions to cellular survival, B cell activation, and proliferation. Our findings provide one possible explanation for the tumorigenicity of KSHV, which poses a severe problem in immunocompromised patients.

c-Jun/c-Fos heterodimers regulate cellular genes via a newly identified class of methylated DNA sequence motifs.

CpG methylation in mammalian DNA is known to interfere with gene expression by inhibiting the binding of transactivators to their cognate sequence motifs or recruiting proteins involved in gene repression. An Epstein-Barr virus-encoded transcription factor, Zta, was the first example of a sequence-specific transcription factor that preferentially recognizes and selectively binds DNA sequence motifs with methylated CpG residues, reverses epigenetic silencing and activates gene transcription. The DNA binding domain of Zta is homologous to c-Fos, a member of the cellular AP-1 (activator protein 1) transcription factor family, which regulates cell proliferation and survival, apoptosis, transformation and oncogenesis. We have identified a novel AP-1 binding site termed meAP-1, which contains a CpG dinucleotide. If methylated, meAP-1 sites are preferentially bound by the AP-1 heterodimer c-Jun/c-Fos in vitro and in cellular chromatin in vivo. In activated human primary B cells, c-Jun/c-Fos locates to these methylated elements in promoter regions of transcriptionally activated genes. Reminiscent of the viral Zta protein, c-Jun/c-Fos is the first identified cellular member of the AP-1 family of transactivators that can induce expression of genes with methylated, hence repressed promoters, reversing epigenetic silencing.

Torpor at high ambient temperature in a neotropical didelphid, the grey short-tailed opossum (Monodelphis domestica).

The grey short-tailed opossum, Monodelphis domestica, has been an established research animal for more than five decades, but relatively, little is known about its thermophysiology. Here we studied core body temperature (T b) and metabolic rate (MR) of female adult M. domestica housed in the laboratory at an ambient temperature (T a) of 26 °C. In expanding previous reports, the average recorded core T b of M. domestica was 34.3 °C. The T b of an individual M. domestica can drop below 30 °C (minimal T b: 28.6 °C) accompanied by a reduction in MR of up to 52 % even while having ad libitum access to food. These findings demonstrate for the first time the presence of spontaneous torpor in M. domestica. Metabolic suppression at relatively high T a and T b furthermore broadens our perspective on the use of torpor as a metabolic strategy not just restricted to cold climates.

Two-stage evolution of mammalian adipose tissue thermogenesis.

Brown adipose tissue (BAT) is a heater organ that expresses thermogenic uncoupling protein 1 (UCP1) to maintain high body temperatures during cold stress. BAT thermogenesis is considered an overarching mammalian trait, but its evolutionary origin is unknown. We show that adipose tissue of marsupials, which diverged from eutherian mammals ~150 million years ago, expresses a nonthermogenic UCP1 variant governed by a partial transcriptomic BAT signature similar to that found in eutherian beige adipose tissue. We found that the reconstructed UCP1 sequence of the common eutherian ancestor displayed typical thermogenic activity, whereas therian ancestor UCP1 is nonthermogenic. Thus, mammalian adipose tissue thermogenesis may have evolved in two distinct stages, with a prethermogenic stage in the common therian ancestor linking UCP1 expression to adipose tissue and thermal stress. We propose that in a second stage, UCP1 acquired its thermogenic function specifically in eutherians, such that the onset of mammalian BAT thermogenesis occurred only after the divergence from marsupials.

Class 3 PI3K coactivates the circadian clock to promote rhythmic de novo purine synthesis.

Metabolic demands fluctuate rhythmically and rely on coordination between the circadian clock and nutrient-sensing signalling pathways, yet mechanisms of their interaction remain not fully understood. Surprisingly, we find that class 3 phosphatidylinositol-3-kinase (PI3K), known best for its essential role as a lipid kinase in endocytosis and lysosomal degradation by autophagy, has an overlooked nuclear function in gene transcription as a coactivator of the heterodimeric transcription factor and circadian driver Bmal1-Clock. Canonical pro-catabolic functions of class 3 PI3K in trafficking rely on the indispensable complex between the lipid kinase Vps34 and regulatory subunit Vps15. We demonstrate that although both subunits of class 3 PI3K interact with RNA polymerase II and co-localize with active transcription sites, exclusive loss of Vps15 in cells blunts the transcriptional activity of Bmal1-Clock. Thus, we establish non-redundancy between nuclear Vps34 and Vps15, reflected by the persistent nuclear pool of Vps15 in Vps34-depleted cells and the ability of Vps15 to coactivate Bmal1-Clock independently of its complex with Vps34. In physiology we find that Vps15 is required for metabolic rhythmicity in liver and, unexpectedly, it promotes pro-anabolic de novo purine nucleotide synthesis. We show that Vps15 activates the transcription of Ppat, a key enzyme for the production of inosine monophosphate, a central metabolic intermediate for purine synthesis. Finally, we demonstrate that in fasting, which represses clock transcriptional activity, Vps15 levels are decreased on the promoters of Bmal1 targets, Nr1d1 and Ppat. Our findings open avenues for establishing the complexity for nuclear class 3 PI3K signalling for temporal regulation of energy homeostasis.

Adipocyte-derived extracellular vesicles increase insulin secretion through transport of insulinotropic protein cargo.

Adipocyte-derived extracellular vesicles (AdEVs) are membranous nanoparticles that convey communication from adipose tissue to other organs. Here, to delineate their role as messengers with glucoregulatory nature, we paired fluorescence AdEV-tracing and SILAC-labeling with (phospho)proteomics, and revealed that AdEVs transfer functional insulinotropic protein cargo into pancreatic ß-cells. Upon transfer, AdEV proteins were subjects for phosphorylation, augmented insulinotropic GPCR/cAMP/PKA signaling by increasing total protein abundances and phosphosite dynamics, and ultimately enhanced 1st-phase glucose-stimulated insulin secretion (GSIS) in murine islets. Notably, insulinotropic effects were restricted to AdEVs isolated from obese and insulin resistant, but not lean mice, which was consistent with differential protein loads and AdEV luminal morphologies. Likewise, in vivo pre-treatment with AdEVs from obese but not lean mice amplified insulin secretion and glucose tolerance in mice. This data suggests that secreted AdEVs can inform pancreatic ß-cells about insulin resistance in adipose tissue in order to amplify GSIS in times of increased insulin demand.

High-calorie diets uncouple hypothalamic oxytocin neurons from a gut-to-brain satiation pathway via κ-opioid signaling.

Oxytocin-expressing paraventricular hypothalamic neurons (PVNOT neurons) integrate afferent signals from the gut, including cholecystokinin (CCK), to adjust whole-body energy homeostasis. However, the molecular underpinnings by which PVNOT neurons orchestrate gut-to-brain feeding control remain unclear. Here, we show that mice undergoing selective ablation of PVNOT neurons fail to reduce food intake in response to CCK and develop hyperphagic obesity on a chow diet. Notably, exposing wild-type mice to a high-fat/high-sugar (HFHS) diet recapitulates this insensitivity toward CCK, which is linked to diet-induced transcriptional and electrophysiological aberrations specifically in PVNOT neurons. Restoring OT pathways in diet-induced obese (DIO) mice via chemogenetics or polypharmacology sufficiently re-establishes CCK's anorexigenic effects. Last, by single-cell profiling, we identify a specialized PVNOT neuronal subpopulation with increased κ-opioid signaling under an HFHS diet, which restrains their CCK-evoked activation. In sum, we document a (patho)mechanism by which PVNOT signaling uncouples a gut-brain satiation pathway under obesogenic conditions.

An ensemble approach for inferring semi-quantitative regulatory dynamics for the differentiation of mouse embryonic stem cells using prior knowledge.

The process of differentiation of embryonic stem cells (ESCs) is currently becoming the focus of many systems biologists not only due to mechanistic interest but also since it is expected to play an increasingly important role in regenerative medicine, in particular with the advert to induced pluripotent stem cells. These ESCs give rise to the formation of the three germ layers and therefore to the formation of all tissues and organs. Here, we present a computational method for inferring regulatory interactions between the genes involved in ESC differentiation based on time resolved microarray profiles. Fully quantitative methods are commonly unavailable on such large-scale data; on the other hand, purely qualitative methods may fail to capture some of the more detailed regulations. Our method combines the beneficial aspects of qualitative and quantitative (ODE-based) modeling approaches searching for quantitative interaction coefficients in a discrete and qualitative state space. We further optimize on an ensemble of networks to detect essential properties and compare networks with respect to robustness. Applied to a toy model our method is able to reconstruct the original network and outperforms an entire discrete boolean approach. In particular, we show that including prior knowledge leads to more accurate results. Applied to data from differentiating mouse ESCs reveals new regulatory interactions, in particular we confirm the activation of Foxh1 through Oct4, mediating Nodal signaling.

Correlation-guided Network Integration (CoNI), an R package for integrating numerical omics data that allows multiform graph representations to study molecular interaction networks.

Summary: Today's immense growth in complex biological data demands effective and flexible tools for integration, analysis and extraction of valuable insights. Here, we present CoNI, a practical R package for the unsupervised integration of numerical omics datasets. Our tool is based on partial correlations to identify putative confounding variables for a set of paired dependent variables. CoNI combines two omics datasets in an integrated, complex hypergraph-like network, represented as a weighted undirected graph, a bipartite graph, or a hypergraph structure. These network representations form a basis for multiple further analyses, such as identifying priority candidates of biological importance or comparing network structures dependent on different conditions. Availability and implementation: The R package CoNI is available on the Comprehensive R Archive Network (https://cran.r-project.org/web/packages/CoNI/) and GitLab (https://gitlab.com/computational-discovery-research/coni). It is distributed under the GNU General Public License (version 3). Supplementary information: Supplementary data are available at Bioinformatics Advances online.

Untargeted and Targeted Circadian Metabolomics Using Liquid Chromatography-Tandem Mass Spectrometry (LC-MS/MS) and Flow Injection-Electrospray Ionization-Tandem Mass Spectrometry (FIA-ESI-MS/MS).

A diverse array of 24-h oscillating hormones and metabolites direct and reflect circadian clock function. Circadian metabolomics uses advanced high-throughput analytical chemistry techniques to comprehensively profile these small molecules (<1.5 kDa) across 24 h in cells, media, body fluids, breath, tissues, and subcellular compartments. The goals of circadian metabolomics experiments are often multifaceted. These include identifying and tracking rhythmic metabolic inputs and outputs of central and peripheral circadian clocks, quantifying endogenous free-running period, monitoring relative phase alignment between clocks, and mapping pathophysiological consequences of clock disruption or misalignment. Depending on the particular experimental question, samples are collected under free-running or entrained conditions. Here we describe both untargeted and targeted liquid chromatography-tandem mass spectrometry (LC-MS/MS) and flow injection-electrospray ionization-tandem mass spectrometry (FIA-ESI-MS/MS) based assays we have used for circadian metabolomics studies. We discuss tissue homogenization, chemical derivatization, measurement, and tips for data processing, normalization, scaling, how to handle outliers, and imputation of missing values.