Trans-omic analysis reveals opposite metabolic dysregulation between feeding and fasting in liver associated with obesity.

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.

DNA hypomethylation characterizes genes encoding tissue-dominant functional proteins in liver and skeletal muscle.

Each tissue has a dominant set of functional proteins required to mediate tissue-specific functions. Epigenetic modifications, transcription, and translational efficiency control tissue-dominant protein production. However, the coordination of these regulatory mechanisms to achieve such tissue-specific protein production remains unclear. Here, we analyzed the DNA methylome, transcriptome, and proteome in mouse liver and skeletal muscle. We found that DNA hypomethylation at promoter regions is globally associated with liver-dominant or skeletal muscle-dominant functional protein production within each tissue, as well as with genes encoding proteins involved in ubiquitous functions in both tissues. Thus, genes encoding liver-dominant proteins, such as those involved in glycolysis or gluconeogenesis, the urea cycle, complement and coagulation systems, enzymes of tryptophan metabolism, and cytochrome P450-related metabolism, were hypomethylated in the liver, whereas those encoding-skeletal muscle-dominant proteins, such as those involved in sarcomere organization, were hypomethylated in the skeletal muscle. Thus, DNA hypomethylation characterizes genes encoding tissue-dominant functional proteins.

A stepwise and digital pattern of RSK phosphorylation determines the outcome of thymic selection.

Developing CD4+CD8+ double-positive (DP) thymocytes with randomly generated T cell receptors (TCRs) undergo positive (maturation) or negative (apoptosis) selection on the basis of the strength of TCR stimulation. Selection fate is determined by engagement of TCR ligands with a subtle difference in affinity, but the molecular details of TCR signaling leading to the different selection outcomes have remained unclear. We performed phosphoproteome analysis of DP thymocytes and found that p90 ribosomal protein kinase (RSK) phosphorylation at Thr562 was induced specifically by high-affinity peptide ligands. Such phosphorylation of RSK triggered its translocation to the nucleus, where it phosphorylated the nuclear receptor Nur77 and thereby promoted its mitochondrial translocation for apoptosis induction. Inhibition of RSK activity protected DP thymocytes from antigen-induced cell death. We propose that RSK phosphorylation constitutes a mechanism by which DP thymocytes generate a stepwise and binary signal in response to exposure to TCR ligands with a graded affinity.

RPL3L-containing ribosomes determine translation elongation dynamics required for cardiac function.

Although several ribosomal protein paralogs are expressed in a tissue-specific manner, how these proteins affect translation and why they are required only in certain tissues have remained unclear. Here we show that RPL3L, a paralog of RPL3 specifically expressed in heart and skeletal muscle, influences translation elongation dynamics. Deficiency of RPL3L-containing ribosomes in RPL3L knockout male mice resulted in impaired cardiac contractility. Ribosome occupancy at mRNA codons was found to be altered in the RPL3L-deficient heart, and the changes were negatively correlated with those observed in myoblasts overexpressing RPL3L. RPL3L-containing ribosomes were less prone to collisions compared with RPL3-containing canonical ribosomes. Although the loss of RPL3L-containing ribosomes altered translation elongation dynamics for the entire transcriptome, its effects were most pronounced for transcripts related to cardiac muscle contraction and dilated cardiomyopathy, with the abundance of the encoded proteins being correspondingly decreased. Our results provide further insight into the mechanisms and physiological relevance of tissue-specific translational regulation.

The ASC-1 complex promotes translation initiation by scanning ribosomes.

Translation initiates when the eIF4F complex binds the 5' mRNA cap, followed by 5' untranslated region scanning for the start codon by scanning ribosomes. Here, we demonstrate that the ASC-1 complex (ASCC), which was previously shown to promote the dissociation of colliding 80S ribosomes, associates with scanning ribosomes to regulate translation initiation. Selective translation complex profiling (TCP-seq) analysis revealed that ASCC3, a helicase domain-containing subunit of ASCC, localizes predominantly to the 5' untranslated region of mRNAs. Ribo-seq, TCP-seq, and luciferase reporter analyses showed that ASCC3 knockdown impairs 43S preinitiation complex loading and scanning dynamics, thereby reducing translation efficiency. Whereas eIF4A, an RNA helicase in the eIF4F complex, is important for global translation, ASCC was found to regulate the scanning process for a specific subset of transcripts. Our results have thus revealed that ASCC is required not only for dissociation of colliding 80S ribosomes but also for efficient translation initiation by scanning ribosomes at a subset of transcripts.

In situ digestion of alcohol-fixed cells for quantitative proteomics.

Currently, the bottom-up approach, in which proteins are digested by enzymes such as trypsin prior to mass spectrometry, is the mainstream approach in mass spectrometer-based proteomics. In this approach, the enzymatic digestion process strongly affects the reproducibility of protein identification and quantification. Here, we quantitatively evaluated the enzymatic digestion of proteins under various conditions by quantitative proteomics using data-independent acquisition and found that proteins precipitated with acetone after solubilization with SDS were fully digestible without re-solubilization. This result implies that organic solvent treatment makes cells amenable to trypsin digestion. Direct trypsin digestion of methanol-fixed cells achieved the same digestion efficiency and quantitative reproducibility as the conventional method. Furthermore, this method was found to be equally applicable to mouse liver samples. The establishment of this method indicates that the sample preparation process in bottom-up proteomics can be simplified while maintaining high digestion efficiency and is expected to become a general method for sample preparation in bottom-up proteomics in the future.

In vivo transomic analyses of glucose-responsive metabolism in skeletal muscle reveal core differences between the healthy and obese states.

Metabolic regulation in skeletal muscle is essential for blood glucose homeostasis. Obesity causes insulin resistance in skeletal muscle, leading to hyperglycemia and type 2 diabetes. In this study, we performed multiomic analysis of the skeletal muscle of wild-type (WT) and leptin-deficient obese (ob/ob) mice, and constructed regulatory transomic networks for metabolism after oral glucose administration. Our network revealed that metabolic regulation by glucose-responsive metabolites had a major effect on WT mice, especially carbohydrate metabolic pathways. By contrast, in ob/ob mice, much of the metabolic regulation by glucose-responsive metabolites was lost and metabolic regulation by glucose-responsive genes was largely increased, especially in carbohydrate and lipid metabolic pathways. We present some characteristic metabolic regulatory pathways found in central carbon, branched amino acids, and ketone body metabolism. Our transomic analysis will provide insights into how skeletal muscle responds to changes in blood glucose and how it fails to respond in obesity.

Multi-omics-based label-free metabolic flux inference reveals obesity-associated dysregulatory mechanisms in liver glucose metabolism.

Glucose homeostasis is maintained by modulation of metabolic flux. Enzymes and metabolites regulate the involved metabolic pathways. Dysregulation of glucose homeostasis is a pathological event in obesity. Analyzing metabolic pathways and the mechanisms contributing to obesity-associated dysregulation in vivo is challenging. Here, we introduce OMELET: Omics-Based Metabolic Flux Estimation without Labeling for Extended Trans-omic Analysis. OMELET uses metabolomic, proteomic, and transcriptomic data to identify relative changes in metabolic flux, and to calculate contributions of metabolites, enzymes, and transcripts to the changes in metabolic flux. By evaluating the livers of fasting ob/ob mice, we found that increased metabolic flux through gluconeogenesis resulted primarily from increased transcripts, whereas that through the pyruvate cycle resulted from both increased transcripts and changes in substrates of metabolic enzymes. With OMELET, we identified mechanisms underlying the obesity-associated dysregulation of metabolic flux in the liver.

Meta-analysis of transcriptional regulatory networks for lipid metabolism in neural cells from schizophrenia patients based on an open-source intelligence approach.

There have been a number of reports about the transcriptional regulatory networks in schizophrenia. However, most of these studies were based on a specific transcription factor or a single dataset, an approach that is inadequate to understand the diverse etiology and underlying common characteristics of schizophrenia. Here we reconstructed and compared the transcriptional regulatory network for lipid metabolism enzymes using 15 public transcriptome datasets of neural cells from schizophrenia patients. Since many of the well-known schizophrenia-related SNPs are in enhancers, we reconstructed a network including enhancer-dependent regulation and found that 53.3 % of the total number of edges (7,577 pairs) involved regulation via enhancers. By examining multiple datasets, we found common and unique transcriptional modes of regulation. Furthermore, enrichment analysis of SNPs that were connected with genes in the transcriptional regulatory networks by eQTL suggested an association with hematological cell counts and some other traits/diseases, whose relationship to schizophrenia was either not or insufficiently reported in previous studies. Based on these results, we suggest that in future studies on schizophrenia, information on genotype, comorbidities and hematological cell counts should be included, along with the transcriptome, for a more detailed genetic stratification and mechanistic exploration of schizophrenia.

A ubiquitin-like protein encoded by the "noncoding" RNA TINCR promotes keratinocyte proliferation and wound healing.

Although long noncoding RNAs (lncRNAs) are transcripts that do not encode proteins by definition, some lncRNAs actually contain small open reading frames that are translated. TINCR (terminal differentiation-induced ncRNA) has been recognized as a lncRNA that contributes to keratinocyte differentiation. However, we here show that TINCR encodes a ubiquitin-like protein that is well conserved among species and whose expression was confirmed by the generation of mice harboring a FLAG epitope tag sequence in the endogenous open reading frame as well as by targeted proteomics. Forced expression of this protein promoted cell cycle progression in normal human epidermal keratinocytes, and mice lacking this protein manifested a delay in skin wound healing associated with attenuated cell cycle progression in keratinocytes. We termed this protein TINCR-encoded ubiquitin-like protein (TUBL), and our results reveal a role for TINCR in the regulation of keratinocyte proliferation and skin regeneration that is dependent on TUBL.

Trans-omic analysis reveals obesity-associated dysregulation of inter-organ metabolic cycles between the liver and skeletal muscle.

Systemic metabolic homeostasis is regulated by inter-organ metabolic cycles involving multiple organs. Obesity impairs inter-organ metabolic cycles, resulting in metabolic diseases. The systemic landscape of dysregulated inter-organ metabolic cycles in obesity has yet to be explored. Here, we measured the transcriptome, proteome, and metabolome in the liver and skeletal muscle and the metabolome in blood of fasted wild-type and leptin-deficient obese (ob/ob) mice, identifying components with differential abundance and differential regulation in ob/ob mice. By constructing and evaluating the trans-omic network controlling the differences in metabolic reactions between fasted wild-type and ob/ob mice, we provided potential mechanisms of the obesity-associated dysfunctions of metabolic cycles between liver and skeletal muscle involving glucose-alanine, glucose-lactate, and ketone bodies. Our study revealed obesity-associated systemic pathological mechanisms of dysfunction of inter-organ metabolic cycles.

Transomics analysis reveals allosteric and gene regulation axes for altered hepatic glucose-responsive metabolism in obesity.

Impaired glucose tolerance associated with obesity causes postprandial hyperglycemia and can lead to type 2 diabetes. To study the differences in liver metabolism in healthy and obese states, we constructed and analyzed transomics glucose-responsive metabolic networks with layers for metabolites, expression data for metabolic enzyme genes, transcription factors, and insulin signaling proteins from the livers of healthy and obese mice. We integrated multiomics time course data from wild-type and leptin-deficient obese (ob/ob) mice after orally administered glucose. In wild-type mice, metabolic reactions were rapidly regulated within 10 min of oral glucose administration by glucose-responsive metabolites, which functioned as allosteric regulators and substrates of metabolic enzymes, and by Akt-induced changes in the expression of glucose-responsive genes encoding metabolic enzymes. In ob/ob mice, the majority of rapid regulation by glucose-responsive metabolites was absent. Instead, glucose administration produced slow changes in the expression of carbohydrate, lipid, and amino acid metabolic enzyme-encoding genes to alter metabolic reactions on a time scale of hours. Few regulatory events occurred in both healthy and obese mice. Thus, our transomics network analysis revealed that regulation of glucose-responsive liver metabolism is mediated through different mechanisms in healthy and obese states. Rapid changes in allosteric regulators and substrates and in gene expression dominate the healthy state, whereas slow changes in gene expression dominate the obese state.

Trans-omic Analysis Reveals ROS-Dependent Pentose Phosphate Pathway Activation after High-Frequency Electrical Stimulation in C2C12 Myotubes.

Skeletal muscle adaptation is mediated by cooperative regulation of metabolism, signal transduction, and gene expression. However, the global regulatory mechanism remains unclear. To address this issue, we performed electrical pulse stimulation (EPS) in differentiated C2C12 myotubes at low and high frequency, carried out metabolome and transcriptome analyses, and investigated phosphorylation status of signaling molecules. EPS triggered extensive and specific changes in metabolites, signaling phosphorylation, and gene expression during and after EPS in a frequency-dependent manner. We constructed trans-omic network by integrating these data and found selective activation of the pentose phosphate pathway including metabolites, upstream signaling molecules, and gene expression of metabolic enzymes after high-frequency EPS. We experimentally validated that activation of these molecules after high-frequency EPS was dependent on reactive oxygen species (ROS). Thus, the trans-omic analysis revealed ROS-dependent activation in signal transduction, metabolome, and transcriptome after high-frequency EPS in C2C12 myotubes, shedding light on possible mechanisms of muscle adaptation.

Chemical acetylation of mitochondrial transcription factor A occurs on specific lysine residues and affects its ability to change global DNA topology.

Chemical acetylation is postulated to occur in mitochondria. Mitochondrial transcription factor A (TFAM or mtTFA), a mitochondrial transcription initiation factor as well as the major mitochondrial nucleoid protein coating the entire mitochondrial genome, is proposed to be acetylated in animals and cultured cells. This study investigated the properties of human TFAM, in conjunction with the mechanism and effects of TFAM acetylation in vitro. Using highly purified recombinant human TFAM and 3 kb circular DNA as a downsized mtDNA model, we studied how the global TFAM-DNA interaction is affected/regulated by the quantitative TFAM-DNA relationship and TFAM acetylation. Results showed that the TFAM-DNA ratio strictly affects the TFAM property to unwind circular DNA in the presence of topoisomerase I. Mass spectrometry analysis showed that in vitro chemical acetylation of TFAM with acetyl-coenzyme A occurs preferentially on specific lysine residues, including those reported to be acetylated in exogenously expressed TFAM in cultured human cells, indicating that chemical acetylation plays a crucial role in TFAM acetylation in mitochondria. Intriguingly, the modification significantly decreased TFAM's DNA-unwinding ability, while its DNA-binding ability was largely unaffected. Altogether, we propose TFAM is chemically acetylated in vivo, which could change mitochondrial DNA topology, leading to copy number and gene expression modulation.

Cell cycle-dependent localization of the proteasome to chromatin.

An integrative understanding of nuclear events including transcription in normal and cancer cells requires comprehensive and quantitative measurement of protein dynamics that underlie such events. However, the low abundance of most nuclear proteins hampers their detailed functional characterization. We have now comprehensively quantified the abundance of nuclear proteins with the use of proteomics approaches in both normal and transformed human diploid fibroblasts. We found that subunits of the 26S proteasome complex were markedly down-regulated in the nuclear fraction of the transformed cells compared with that of the wild-type cells. The intranuclear proteasome abundance appeared to be inversely related to the rate of cell cycle progression, with restraint of the cell cycle being associated with an increase in the amount of proteasome subunits in the nucleus, suggesting that the nuclear proteasome content is dependent on the cell cycle. Furthermore, chromatin enrichment for proteomics (ChEP) analysis revealed enrichment of the proteasome in the chromatin fraction of quiescent cells and its apparent dissociation from chromatin in transformed cells. Our results thus suggest that translocation of the nuclear proteasome to chromatin may play an important role in control of the cell cycle and oncogenesis through regulation of chromatin-associated transcription factors.

Reconstruction of global regulatory network from signaling to cellular functions using phosphoproteomic data.

Cellular signaling regulates various cellular functions via protein phosphorylation. Phosphoproteomic data potentially include information for a global regulatory network from signaling to cellular functions, but a procedure to reconstruct this network using such data has yet to be established. In this paper, we provide a procedure to reconstruct a global regulatory network from signaling to cellular functions from phosphoproteomic data by integrating prior knowledge of cellular functions and inference of the kinase-substrate relationships (KSRs). We used phosphoproteomic data from insulin-stimulated Fao hepatoma cells and identified protein phosphorylation regulated by insulin specifically over-represented in cellular functions in the KEGG database. We inferred kinases for protein phosphorylation by KSRs, and connected the kinases in the insulin signaling layer to the phosphorylated proteins in the cellular functions, revealing that the insulin signal is selectively transmitted via the Pi3k-Akt and Erk signaling pathways to cellular adhesions and RNA maturation, respectively. Thus, we provide a method to reconstruct global regulatory network from signaling to cellular functions based on phosphoproteomic data.

Trans-omic Analysis Reveals Selective Responses to Induced and Basal Insulin across Signaling, Transcriptional, and Metabolic Networks.

The concentrations of insulin selectively regulate multiple cellular functions. To understand how insulin concentrations are interpreted by cells, we constructed a trans-omic network of insulin action in FAO hepatoma cells using transcriptomic data, western blotting analysis of signaling proteins, and metabolomic data. By integrating sensitivity into the trans-omic network, we identified the selective trans-omic networks stimulated by high and low doses of insulin, denoted as induced and basal insulin signals, respectively. The induced insulin signal was selectively transmitted through the pathway involving Erk to an increase in the expression of immediate-early and upregulated genes, whereas the basal insulin signal was selectively transmitted through a pathway involving Akt and an increase of Foxo phosphorylation and a reduction of downregulated gene expression. We validated the selective trans-omic network in vivo by analysis of the insulin-clamped rat liver. This integrated analysis enabled molecular insight into how liver cells interpret physiological insulin signals to regulate cellular functions.

Selective control of up-regulated and down-regulated genes by temporal patterns and doses of insulin.

Secretion of insulin transiently increases after eating, resulting in a high circulating concentration. Fasting limits insulin secretion, resulting in a low concentration of insulin in the circulation. We analyzed transcriptional responses to different temporal patterns and doses of insulin in the hepatoma FAO cells and identified 13 up-regulated and 16 down-regulated insulin-responsive genes (IRGs). The up-regulated IRGs responded more rapidly than did the down-regulated IRGs to transient stepwise or pulsatile increases in insulin concentration, whereas the down-regulated IRGs were repressed at lower concentrations of insulin than those required to stimulate the up-regulated IRGs. Mathematical modeling of the insulin response as two stages-(i) insulin signaling to transcription and (ii)transcription and mRNA stability-indicated that the first stage was the more rapid stage for the down-regulated IRGs, whereas the second stage of transcription was the more rapid stage for the up-regulated IRGs. A subset of the IRGs that were up-regulated or down-regulated in the FAO cells was similarly regulated in the livers of rats injected with a single dose of insulin. Thus, not only can cells respond to insulin but they can also interpret the intensity and pattern of signal to produce distinct transcriptional responses. These results provide insight that may be useful in treating obesity and type 2 diabetes associated with aberrant insulin production or tissue responsiveness.

Phosphoproteomics analyses show subnetwork systems in T-cell receptor signaling.

A key issue in the study of signal transduction is how multiple signaling pathways are systematically integrated into the cell. We have now performed multiple phosphoproteomics analyses focused on the dynamics of the T-cell receptor (TCR) signaling network and its subsystem mediated by the Ca2+ signaling pathway. Integration of these phosphoproteomics data sets and extraction of components of the TCR signaling network dependent on Ca2+ signaling showed unexpected phosphorylation kinetics for candidate substrates of the Ca2+ -dependent phosphatase calcineurin (CN) during TCR stimulation. Detailed characterization of the TCR-induced phosphorylation of a novel CN substrate, Itpkb, showed that phosphorylation of this protein is regulated by both CN and the mitogen-activated protein kinase Erk in a competitive manner. Phosphorylation of additional CN substrates was also found to be regulated by Erk and CN in a similar manner. The combination of multiple phosphoproteomics approaches thus showed two major subsystems mediated by Erk and CN in the TCR signaling network, with these subsystems regulating the phosphorylation of a group of proteins in a competitive manner.

Trans-Omics: How To Reconstruct Biochemical Networks Across Multiple 'Omic' Layers.

We propose 'trans-omic' analysis for reconstructing global biochemical networks across multiple omic layers by use of both multi-omic measurements and computational data integration. We introduce technologies for connecting multi-omic data based on prior knowledge of biochemical interactions and characterize a biochemical trans-omic network by concepts of a static and dynamic nature. We introduce case studies of metabolism-centric trans-omic studies to show how to reconstruct a biochemical trans-omic network by connecting multi-omic data and how to analyze it in terms of the static and dynamic nature. We propose a trans-ome-wide association study (trans-OWAS) connecting phenotypes with trans-omic networks that reflect both genetic and environmental factors, which can characterize several complex lifestyle diseases as breakdowns in the trans-omic system.