Sparse Gaussian graphical model with missing values

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

Timing-Dependent Actions of NGF Required for Cell Differentiation

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

Logical design of oral glucose ingestion pattern minimizing blood glucose in humans

Excessive increase in blood glucose level after eating increases the risk of macroangiopathy, and a method for not increasing the postprandial blood glucose level is desired. However, a logical design method of the dietary ingestion pattern controlling the postprandial blood glucose level has not yet been established. We constructed a mathematical model of blood glucose control by oral glucose ingestion in three healthy human subjects, and predicted that intermittent ingestion 30 min apart was the optimal glucose ingestion patterns that minimized the peak value of blood glucose level. We confirmed with subjects that this intermittent pattern consistently decreased the peak value of blood glucose level. We also predicted insulin minimization pattern, and found that the intermittent ingestion 30 min apart was optimal, which is similar to that of glucose minimization pattern. Taken together, these results suggest that the glucose minimization is achieved by suppressing the peak value of insulin concentration, rather than by enhancing insulin concentration. This approach could be applied to design optimal dietary ingestion patterns

Trans-omics Impact of Thymoproteasome in Cortical Thymic Epithelial Cells

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

Reconstruction of Insulin Signal Flow from Phosphoproteome and Metabolome Data

SummaryCellular 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

A Quantitative Image Cytometry Technique for Time Series or Population Analyses of Signaling Networks

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

Reconstruction of Insulin Signal Flow from Phosphoproteome and Metabolome Data

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.UTokyo Research掲載「細胞内のビッグデータから大規模ネットワークの再構築に成功」URI: http://www.u-tokyo.ac.jp/ja/utokyo-research/research-news/reconstruction-of-molecular-network-from-cellular-big-data/UTokyo Research "Reconstruction of molecular network from cellular big data" URI: http://www.u-tokyo.ac.jp/en/utokyo-research/research-news/reconstruction-of-molecular-network-from-cellular-big-data

Minimizing the cross validation error to mix kernel matrices of heterogeneous biological data,” Neural Process

Abstract. In biological data, it is often the case that objects are described in two or more representations. In order to perform classification based on such data, we have to combine them in a certain way. In the context of kernel machines, this task amounts to mix several kernel matrices into one. In this paper, we present two ways to mix kernel matrices, where the mixing weights are optimized to minimize the cross validation error. In bacteria classification and gene function prediction experiments, our methods significantly outperformed single kernel classifiers in most cases. Key words. bacteria classification, bioinformatics, kernel machines, mixing kernel matrices 1