Search

Abstract Interactions between various molecular species in biological phenomena give rise to numerous networks. The investigation of these networks, including their statistical and biochemical interactions, supports a deeper understanding of biological phenomena. The clustering of nodes associated with molecular species and enrichment analysis is frequently applied to examine the biological significance of such network structures. However, these methods focus on delineating the function of a node. As such, in-depth investigations of the edges, which are the connections between the nodes, are rarely explored. In the current study, we aimed to investigate the functions of the edges rather than the nodes. To accomplish this, for each network, we categorized the edges and defined the edge type based on their biological annotations. Subsequently, we used the edge type to compare the network structures of the metabolome and transcriptome in the livers of healthy (wild-type) and obese (ob/ob) mice following oral glucose administration (OGTT). The findings demonstrate that the edge type can facilitate the characterization of the state of a network structure, thereby reducing the information available through datasets containing the OGTT response in the metabolome and transcriptome.

Abstract Oral glucose ingestion induces systemic changes of many blood metabolites related not only to glucose, but also other metabolites such as amino acids and lipids through many blood hormones. However, the detailed temporal changes in the concentrations of comprehensive metabolites and hormones over a long time by oral glucose ingestion are uncharacterized. We measured 83 metabolites and 7 hormones in 20 healthy human subjects in response to glucose ingestion. We characterized temporal patterns of blood molecules by four features: (i) the decomposability into “amplitude” and “rate” components, (ii) the similarity of temporal patterns among individuals, (iii) the relation of molecules over time among individuals, and (iv) the similarity of temporal patterns among molecules. Glucose and glucose metabolism-related hormones indicated a rapid increase, and citrulline and lipids, which indicated a rapid decrease, returned to fasting levels faster than amino acids. Compared to glucose metabolism-related molecules and lipids, amino acids showed similar temporal patterns among individuals. The four features of temporal patterns of blood molecules by oral glucose ingestion characterize the differences among individuals and among molecules.

Summary: An effective combination of multi-omic datasets can enhance our understanding of complex biological phenomena. To build a context-dependent network with multiple omic layers, i.e., a trans-omic network, we perform phosphoproteomics, transcriptomics, proteomics, and metabolomics of murine liver for 4 h after insulin administration and integrate the resulting time series. Structural characteristics and dynamic nature of the network are analyzed to elucidate the impact of insulin. Early and prominent changes in protein phosphorylation and persistent and asynchronous changes in mRNA and protein levels through non-transcriptional mechanisms indicate enhanced crosstalk between phosphorylation-mediated signaling and protein expression regulation. Metabolic response shows different temporal regulation with transient increases at early time points across categories and enhanced response in the amino acid and nucleotide categories at later time points as a result of process convergence. This extensive and dynamic view of the trans-omic network elucidates prominent regulatory mechanisms that drive insulin responses through intricate interlayer coordination.

Summary: 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.

Summary: 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. : Systems Biology; Omics; Metabolomics; Transcriptomics Subject Areas: Systems Biology, Omics, Metabolomics, Transcriptomics

Recent advances in measurement technology have enabled us to measure various omic layers, such as genome, transcriptome, proteome, and metabolome layers. The demand for data analysis to determine the network structure of the interaction between molecular species is increasing. The Gaussian graphical model is one method of estimating the network structure. However, biological omics data sets tend to include missing values, which is conventionally handled by preprocessing. We propose a novel method by which to estimate the network structure together with missing values by combining a sparse graphical model and matrix factorization. The proposed method was validated by artificial data sets and was applied to a signal transduction data set as a test run.

Summary: The thymic function to produce self-protective and self-tolerant T cells is chiefly mediated by cortical thymic epithelial cells (cTECs) and medullary TECs (mTECs). Recent studies including single-cell transcriptomic analyses have highlighted a rich diversity in functional mTEC subpopulations. Because of their limited cellularity, however, the biochemical characterization of TECs, including the proteomic profiling of cTECs and mTECs, has remained unestablished. Utilizing genetically modified mice that carry enlarged but functional thymuses, here we show a combination of proteomic and transcriptomic profiles for cTECs and mTECs, which identified signature molecules that characterize a developmental and functional contrast between cTECs and mTECs. Our results reveal a highly specific impact of the thymoproteasome on proteasome subunit composition in cTECs and provide an integrated trans-omics platform for further exploration of thymus biology. : Ohigashi et al. show that the use of cyclin D1-transgenic mice allows quantitative proteomic analysis of cortical and medullary thymic epithelial cells (TECs). Results provide a trans-omics platform for further exploration of TEC biology and reveal the specific impact of the thymoproteasome on proteasome subunit composition in cortical TECs. Keywords: cortical thymic epithelial cells, medullary thymic epithelial cells, proteome, trans-omics, thymoproteasome, Psmb11, cyclin D1

Cells decode information of signaling activation at a scale of tens of minutes by downstream gene expression with a scale of hours to days, leading to cell fate decisions such as cell differentiation. However, no system identification method with such different time scales exists. Here we used compressed sensing technology and developed a system identification method using data of different time scales by recovering signals of missing time points. We measured phosphorylation of ERK and CREB, immediate early gene expression products, and mRNAs of decoder genes for neurite elongation in PC12 cell differentiation and performed system identification, revealing the input-output relationships between signaling and gene expression with sensitivity such as graded or switch-like response and with time delay and gain, representing signal transfer efficiency. We predicted and validated the identified system using pharmacological perturbation. Thus, we provide a versatile method for system identification using data with different time scales.

Cellular homeostasis is regulated by signals through multiple molecular networks that include protein phosphorylation and metabolites. However, where and when the signal flows through a network and regulates homeostasis has not been explored. We have developed a reconstruction method for the signal flow based on time-course phosphoproteome and metabolome data, using multiple databases, and have applied it to acute action of insulin, an important hormone for metabolic homeostasis. An insulin signal flows through a network, through signaling pathways that involve 13 protein kinases, 26 phosphorylated metabolic enzymes, and 35 allosteric effectors, resulting in quantitative changes in 44 metabolites. Analysis of the network reveals that insulin induces phosphorylation and activation of liver-type phosphofructokinase 1, thereby controlling a key reaction in glycolysis. We thus provide a versatile method of reconstruction of signal flow through the network using phosphoproteome and metabolome data.

Signaling networks are made up of limited numbers of molecules and yet can code information that controls different cellular states through temporal patterns and a combination of signaling molecules. In this study, we used a data-driven modeling approach, the Laguerre filter with partial least square regression, to describe how temporal and combinatorial patterns of signaling molecules are decoded by their downstream targets. The Laguerre filter is a time series model used to represent a nonlinear system based on Volterra series expansion. Furthermore, with this approach, each component of the Volterra series expansion is expanded by Laguerre basis functions. We combined two approaches, application of a Laguerre filter and partial least squares (PLS) regression, and applied the combined approach to analysis of a signal transduction network. We applied the Laguerre filter with PLS regression to identify input and output (IO) relationships between MAP kinases and the products of immediate early genes (IEGs). We found that Laguerre filter with PLS regression performs better than Laguerre filter with ordinary regression for the reproduction of a time series of IEGs. Analysis of the nonlinear characteristics extracted using the Laguerre filter revealed a priming effect of ERK and CREB on c-FOS induction. Specifically, we found that the effects of a first pulse of ERK enhance the subsequent effects on c-FOS induction of treatment with a second pulse of ERK, a finding consistent with prior molecular biological knowledge. The variable importance of projections and output loadings in PLS regression predicted the upstream dependency of each IEG. Thus, a Laguerre filter with partial least square regression approach appears to be a powerful method to find the processing mechanism of temporal patterns and combination of signaling molecules by their downstream gene expression.

Homeostatic control of blood glucose is regulated by a complex feedback loop between glucose and insulin, of which failure leads to diabetes mellitus. However, physiological and pathological nature of the feedback loop is not fully understood. We made a mathematical model of the feedback loop between glucose and insulin using time course of blood glucose and insulin during consecutive hyperglycemic and hyperinsulinemic-euglycemic clamps in 113 subjects with variety of glucose tolerance including normal glucose tolerance (NGT), impaired glucose tolerance (IGT) and type 2 diabetes mellitus (T2DM). We analyzed the correlation of the parameters in the model with the progression of glucose intolerance and the conserved relationship between parameters. The model parameters of insulin sensitivity and insulin secretion significantly declined from NGT to IGT, and from IGT to T2DM, respectively, consistent with previous clinical observations. Importantly, insulin clearance, an insulin degradation rate, significantly declined from NGT, IGT to T2DM along the progression of glucose intolerance in the mathematical model. Insulin clearance was positively correlated with a product of insulin sensitivity and secretion assessed by the clamp analysis or determined with the mathematical model. Insulin clearance was correlated negatively with postprandial glucose at 2h after oral glucose tolerance test. We also inferred a square-law between the rate constant of insulin clearance and a product of rate constants of insulin sensitivity and secretion in the model, which is also conserved among NGT, IGT and T2DM subjects. Insulin clearance shows a conserved relationship with the capacity of glucose disposal among the NGT, IGT and T2DM subjects. The decrease of insulin clearance predicts the progression of glucose intolerance.

Cells use common signaling molecules for the selective control of downstream gene expression and cell-fate decisions. The relationship between signaling molecules and downstream gene expression and cellular phenotypes is a multiple-input and multiple-output (MIMO) system and is difficult to understand due to its complexity. For example, it has been reported that, in PC12 cells, different types of growth factors activate MAP kinases (MAPKs) including ERK, JNK, and p38, and CREB, for selective protein expression of immediate early genes (IEGs) such as c-FOS, c-JUN, EGR1, JUNB, and FOSB, leading to cell differentiation, proliferation and cell death; however, how multiple-inputs such as MAPKs and CREB regulate multiple-outputs such as expression of the IEGs and cellular phenotypes remains unclear. To address this issue, we employed a statistical method called partial least squares (PLS) regression, which involves a reduction of the dimensionality of the inputs and outputs into latent variables and a linear regression between these latent variables. We measured 1,200 data points for MAPKs and CREB as the inputs and 1,900 data points for IEGs and cellular phenotypes as the outputs, and we constructed the PLS model from these data. The PLS model highlighted the complexity of the MIMO system and growth factor-specific input-output relationships of cell-fate decisions in PC12 cells. Furthermore, to reduce the complexity, we applied a backward elimination method to the PLS regression, in which 60 input variables were reduced to 5 variables, including the phosphorylation of ERK at 10 min, CREB at 5 min and 60 min, AKT at 5 min and JNK at 30 min. The simple PLS model with only 5 input variables demonstrated a predictive ability comparable to that of the full PLS model. The 5 input variables effectively extracted the growth factor-specific simple relationships within the MIMO system in cell-fate decisions in PC12 cells.

A wide range of growth factors encode information into specific temporal patterns of MAP kinase (MAPK) and CREB phosphorylation, which are further decoded by expression of immediate early gene products (IEGs) to exert biological functions. However, the IEG decoding system remain unknown. We built a data-driven based on time courses of MAPK and CREB phosphorylation and IEG expression in response to various growth factors to identify how signal is processed. We found that IEG expression uses common decoding systems regardless of growth factors and expression of each IEG differs in upstream dependency, switch-like response, and linear temporal filters. Pulsatile ERK phosphorylation was selectively decoded by expression of EGR1 rather than c-FOS. Conjunctive NGF and PACAP stimulation was selectively decoded by synergistic JUNB expression through switch-like response to c-FOS. Thus, specific temporal patterns and combinations of MAPKs and CREB phosphorylation can be decoded by selective IEG expression via distinct temporal filters and switch-like responses. The data-driven modeling is versatile for analysis of signal processing and does not require detailed prior knowledge of pathways.

BACKGROUND: Modeling of cellular functions on the basis of experimental observation is increasingly common in the field of cellular signaling. However, such modeling requires a large amount of quantitative data of signaling events with high spatio-temporal resolution. A novel technique which allows us to obtain such data is needed for systems biology of cellular signaling. METHODOLOGY/PRINCIPAL FINDINGS: We developed a fully automatable assay technique, termed quantitative image cytometry (QIC), which integrates a quantitative immunostaining technique and a high precision image-processing algorithm for cell identification. With the aid of an automated sample preparation system, this device can quantify protein expression, phosphorylation and localization with subcellular resolution at one-minute intervals. The signaling activities quantified by the assay system showed good correlation with, as well as comparable reproducibility to, western blot analysis. Taking advantage of the high spatio-temporal resolution, we investigated the signaling dynamics of the ERK pathway in PC12 cells. CONCLUSIONS/SIGNIFICANCE: The QIC technique appears as a highly quantitative and versatile technique, which can be a convenient replacement for the most conventional techniques including western blot, flow cytometry and live cell imaging. Thus, the QIC technique can be a powerful tool for investigating the systems biology of cellular signaling.

BACKGROUND: Continuous NGF stimulation induces PC12 cell differentiation. However, why continuous NGF stimulation is required for differentiation is unclear. In this study, we investigated the underlying mechanisms of the timing-dependent requirement of NGF action for cell differentiation. METHODOLOGY/PRINCIPAL FINDINGS: To address the timing-dependency of the NGF action, we performed a discontinuous stimulation assay consisting of a first transient stimulation followed by an interval and then a second sustained stimulation and quantified the neurite extension level. Consequently, we observed a timing-dependent action of NGF on cell differentiation, and discontinuous NGF stimulation similarly induced differentiation. The first stimulation did not induce neurite extension, whereas the second stimulation induced fast neurite extension; therefore, the first stimulation is likely required as a prerequisite condition. These observations indicate that the action of NGF can be divided into two processes: an initial stimulation-driven latent process and a second stimulation-driven extension process. The latent process appears to require the activities of ERK and transcription, but not PI3K, whereas the extension-process requires the activities of ERK and PI3K, but not transcription. We also found that during the first stimulation, the activity of NGF can be replaced by PACAP, but not by insulin, EGF, bFGF or forskolin; during the second stimulation, however, the activity of NGF cannot be replaced by any of these stimulants. These findings allowed us to identify potential genes specifically involved in the latent process, rather than in other processes, using a microarray. CONCLUSIONS/SIGNIFICANCE: These results demonstrate that NGF induces the differentiation of PC12 cells via mechanically distinct processes: an ERK-driven and transcription-dependent latent process, and an ERK- and PI3K-driven and transcription-independent extension process.

Three-dimensional density measurement of unsteady flow field around a projectile (φ = 8mm) is carried out in the ballistic range at Institute of Fluid Science, Tohoku University. The purpose of this study is to obtain the density distribution including its wake region in non-axisymmetric unsteady flow. The projectile’s Mach number in the experiment is 1.35. Simultaneous multi-angle BOS measurement system using twelve digital cameras and pulsed LEDs is installed in the test chamber of the ballistic range. The Color-Grid Background Oriented Schlieren (CGBOS) technique is used in the measurement system to obtain the projection data of density gradient. Algebraic Reconstruction Technique (ART) method and Successive Over Relaxation (SOR) method are used to obtain the density distribution from projection data. To improve the accuracy of 3D reconstruction, we evaluate the view-angle of the measurement system in the present study. The result shows that the bow shock, expansion wave and recirculation zone can be confirmed by the proposed method. However, bow shock in each method are thicker than the schlieren image due to the blur of the CGBOS images. Even then, the 3D shape of density distribution reconstructed by proposed method can be reliable, and the relative error of obtained density between the theoretical value for normal shock relation is about 5% at the bow shock.

Robust transmission of information despite the presence of variation is a fundamental problem in cellular functions. However, the capability and characteristics of information transmission in signaling pathways remain poorly understood. We describe robustness and compensation of information transmission of signaling pathways at the cell population level. We calculated the mutual information transmitted through signaling pathways for the growth factor-mediated gene expression. Growth factors appeared to carry only information sufficient for a binary decision. Information transmission was generally more robust than average signal intensity despite pharmacological perturbations, and compensation of information transmission occurred. Information transmission to the biological output of neurite extension appeared robust. Cells may use information entropy as information so that messages can be robustly transmitted despite variation in molecular activities among individual cells. [ABSTRACT FROM AUTHOR]

A latent process involving signal transduction and gene expression is needed as a preparation step for cellular function. We previously found that nerve growth factor (NGF)-induced cell differentiation has a latent process, which is dependent on ERK activity and gene expression and required for subsequent neurite extension. A latent process can be considered as a preparation step that decodes extracellular stimulus information into cellular functions; however, molecular mechanisms of this process remain unknown. We identified Metrnl, Dclk1 and Serpinb1a as genes that are induced during the latent process (LP) with distinct temporal expression profiles and are required for subsequent neurite extension in PC12 cells. The LP genes showed distinct dependency on the duration of ERK activity, and they were also induced during the latent process of PACAP- and forskolin-induced cell differentiation. Regardless of neurotrophic factors, expression levels of the LP genes during the latent process (0-12 hours), but not phosphorylation levels of ERK, always correlated with subsequent neurite extension length (12-24 hours). Overexpression of all LP genes together, but not of each gene separately, enhanced NGF-induced neurite extension. The LP gene products showed distinct spatial localization. Thus, the LP genes appear to be the common decoders for neurite extension length regardless of neurotrophic factors, and they might function in distinct temporal and spatial manners during the latent process. Our findings provide molecular insight into the physiological meaning of the latent process as the preparation step for decoding information for future phenotypic change. [ABSTRACT FROM AUTHOR]