Bioinformatics

Proteomics provides important information - that may not be inferable from indirect sources such as RNA or DNA - on key players in biological systems or disease states. However, it suffers from coverage and consistency problems. The advent of network-based analysis methods can help in overcoming these problems but requires careful application and interpretation. This review considers briefly current trends in proteomics technologies and understanding the causes of critical issues that need to be addressed - i.e., incomplete data coverage and inter-sample inconsistency. On the coverage issue, we argue that holistic analysis based on biological networks provides a suitable background on which more robust models and interpretations can be built upon; and we introduce some recently developed approaches. On consistency, group-based approaches based on identified clusters, as well as on properly integrated pathway databases, are particularly useful. Despite that protein interactions and pathway networks are still largely incomplete, given proper quality checks, applications and reasonably sized data sets, they yield valuable insights that greatly complement data generated from quantitative proteomics.

Traditional proteomics analysis is plagued by the use of arbitrary thresholds resulting in large loss of information. We propose here a novel method in proteomics that utilizes all detected proteins. We demonstrate its efficacy in a proteomics screen of 5 and 7 liver cancer patients in the moderate and late stage, respectively. Utilizing biological complexes as a cluster vector, and augmenting it with submodules obtained from partitioning an integrated and cleaned protein− protein interaction network, we calculate a Proteomics Signature Profile (PSP) for each patient based on the hit rates of their reported proteins, in the absence of fold change thresholds, against the cluster vector. Using this, we demonstrated that moderate-and late-stage patients segregate with high confidence. We also discovered a moderate-stage patient who displayed a proteomics profile similar to other poor-stage patients. We identified significant clusters using a modified version of the SNet approach. Comparing our results against the Proteomics Expansion Pipeline (PEP) on which the same patient data was analyzed, we found good correlation. Building on this finding, we report significantly more clusters (176 clusters here compared to 70 in PEP), demonstrating the sensitivity of this approach. Gene Ontology (GO) terms analysis also reveals that the significant clusters are functionally congruent with the liver cancer phenotype. PSP is a powerful and sensitive method for analyzing proteomics profiles even when sample sizes are small. It does not rely on the ratio scores but, rather, whether a protein is detected or not. Although consistency of individual proteins between patients is low, we found the reported proteins tend to hit clusters in a meaningful and informative manner. By extracting this information in the form of a Proteomics Signature Profile, we confirm that this information is conserved and can be used for (1) clustering of patient samples, (2) identification of significant clusters based on real biological complexes, and (3) overcoming consistency and coverage issues prevalent in proteomics data sets.

Objective To investigate the effect of serum amyloid A1 (SAA1) on global gene expression in macrophages derived from THP-1 monocytes. Materials and methods Global genetic expression in THP-1-derived macrophages was determined using Illumina HT-12 microarray chips and the results were validated by real-time PCR. Cytokine levels in cellular supernatant were quantified by ELISA. Results In total, 55 genes were upregulated with fold difference greater than two when THP-1-derived macrophages were incubated with SAA1 for 8 h. SAA1 is a strong cytokine inducer with significant upregulation of chemokines CCL1, CCL3, and CCL4 and this was confirmed by both real-time PCR and ELISA quantification. SAA1 also promotes the upregulation of genes involved in phagocytosis, anti-apoptosis, and tissue remodeling. Conclusions SAA1 appears to play an important role during the immune response and in chronic inflammatory diseases through the stimulation of genes involved in cytokine production, phagocytosis, and anti-apoptosis.

microRibonucleic acid (miRNAs) are small regulatory molecules that act by mRNA degradation or via translational repression. Although many miRNAs are ubiquitously expressed, a small subset have differential expression patterns that may give rise to tissue-specific complexes. Motivation: This work studies gene targeting patterns amongst miRNAs with differential expression profiles, and links this to control and regulation of protein complexes. Results: We find that, when a pair of miRNAs are not expressed in the same tissues, there is a higher tendency for them to target the direct partners of the same hub proteins. At the same time, they also avoid targeting the same set of hub-spokes. Moreover, the complexes corresponding to these hub-spokes tend to be specific and nonoverlapping. This suggests that the effect of miRNAs on the formation of complexes is specific.

Proteomics Signature Profiling (PSP) is a novel hit-rate based method that proved useful in resolving consistency and coverage issues in proteomics. As a follow-up study, several points need to be addressed: 1/ PSP's generalisability to pathways, 2/ understand the biological interplay between significant complexes and pathway subnets co-located on the same pathways on our liver cancer dataset, 3/ understand PSP's false positive rate and 4/ demonstration that PSP works on other suitable proteomics datasets as well as expanding PSP's analytical resolution via the use of specialised ontologies.

Despite its prominence for characterization of complex mixtures, LC−MS/MS frequently fails to identify many proteins. Networkbased analysis methods, based on protein−protein interaction networks (PPINs), biological pathways, and protein complexes, are useful for recovering non-detected proteins, thereby enhancing analytical resolution. However, network-based analysis methods do come in varied flavors for which the respective efficacies are largely unknown. We compare the recovery performance and functional insights from three distinct instances of PPIN-based approaches, viz., Proteomics Expansion Pipeline (PEP), Functional Class Scoring (FCS), and Maxlink, in a test scenario of valproic acid (VPA)-treated mice. We find that the most comprehensive functional insights, as well as best non-detected protein recovery performance, are derived from FCS utilizing real biological complexes. This outstrips other network-based methods such as Maxlink or Proteomics Expansion Pipeline (PEP). From FCS, we identified known biological complexes involved in epigenetic modifications, neuronal system development, and cytoskeletal rearrangements. This is congruent with the observed phenotype where adult mice showed an increase in dendritic branching to allow the rewiring of visual cortical circuitry and an improvement in their visual acuity when tested behaviorally. In addition, PEP also identified a novel complex, comprising YWHAB, NR1, NR2B, ACTB, and TJP1, which is functionally related to the observed phenotype. Although our results suggest different network analysis methods can produce different results, on the whole, the findings are mutually supportive. More critically, the non-overlapping information each provides can provide greater holistic understanding of complex phenotypes.

Troglitazone, a first-generation thiazolidinedione of antihyperglycaemic properties, was withdrawn from the market due to unacceptable idiosyncratic hepatotoxicity. Despite intensive research, the underlying mechanism of troglitazone-induced liver toxicity remains unknown. Here we report the use of the Sod2 +/− mouse model of silent mitochondrial oxidative-stress-based and quantitative mass spectrometry-based proteomics to track the mitochondrial proteome changes induced by physiologically relevant troglitazone doses. By quantitative untargeted proteomics, we first globally profiled the Sod2 +/− hepatic mitochondria proteome and found perturbations including GSH metabolism that enhanced the toxicity of the normally nontoxic troglitazone. Short-and long-term troglitazone administration in Sod2 +/− mouse led to a mitochondrial proteome shift from an early compensatory response to an eventual phase of intolerable oxidative stress, due to decreased mitochondrial glutathione (mGSH) import protein, decreased dicarboxylate ion carrier (DIC), and the specific activation of ASK1-JNK and FOXO3a with prolonged troglitazone exposure. Furthermore, mapping of the detected proteins onto mouse specific proteincentered networks revealed lipid-associated proteins as contributors to overt mitochondrial and liver injury when under prolonged exposure to the lipid-normalizing troglitazone. By integrative toxicoproteomics, we demonstrated a powerful systems approach in identifying the collapse of specific fragile nodes and activation of crucial proteome reconfiguration regulators when targeted by an exogenous toxicant.

The proper combination of proteomics with bioinformatics is highly synergistic, capable of propelling proteomics into a truly high-throughput platform. However, appropriate experimental design and analytical considerations are needed to maximize analytical outcome. This review highlights key issues and caveats in converting raw data to protein identifications, and subsequently biological insight. It offers some insights on how the establishment of highly robust proteomics pipelines can be used for studying novel areas such as computational epigenetics, high performance functional studies and new discovery paradigms for drug targets and biomarkers Reviews INFORMATICS

The integration of networks with genomics (network genomics) is a familiar field. Conventional network analysis takes advantage of the larger coverage and relative stability of gene expression measurements. Network proteomics on the other hand has to develop further on two critical factors: (1) expanded data coverage and consistency, and (2) suitable reference network libraries, and data mining from them. Concerning (1) we discuss several contemporary themes that can improve data quality, which in turn will boost the outcome of downstream network analysis. For (2), we focus on network analysis developments, specifically, the need for context-specific networks and essential considerations for localized network analysis.

The recruitment and migration of macrophages and neutrophils is an important process during the early stages of the innate immune system in response to acute injury. Transgenic pu.1:EGFP zebrafish permit the acquisition of leukocyte migration trajectories during inflammation. Currently, these high-quality live-imaging data are mainly analysed using general statistics, for example, cell velocity. Here, we present a spatio-temporal analysis of the cell dynamics using transition matrices, which provide information of the type of cell migration. We find evidence that leukocytes exhibit types of migratory behaviour, which differ from previously described random walk processes. Dimethyl sulfoxide treatment decreased the level of persistence at early time points after wounding and ablated temporal dependencies observed in untreated embryos. We then use pharmacological inhibition of p38 and c-Jun N-terminal kinase mitogen-activated protein kinases to determine their effects on in vivo leukocyte migration patterns and discuss how they modify the characteristics of the cell migration process. In particular, we find that their respective inhibition leads to decreased and increased levels of persistent motion in leukocytes following wounding. This example shows the high level of information content, which can be gained from live-imaging data if appropriate statistical tools are used.

When analysing gene expression time series data an often overlooked but crucial aspect of the model is that the regulatory network structure may change over time. Whilst some approaches have addressed this problem previously in the literature, many are not well suited to the sequential nature of the data. Here we present a method that allows us to infer regulatory network structures that may vary between time points, utilising a set of hidden states that describe the network structure at a given time point. To model the distribution of the hidden states we have applied the Hierarchical Dirichlet Process Hidden Markov Model, a nonparametric extension of the traditional Hidden Markov Model, that does not require us to fix the number of hidden states in advance. We apply our method to existing microarray expression data as well as demonstrating is efficacy on simulated test data.

The 40-fold increase in childhood megakaryocyte-erythroid and B-cell leukemia in Down syndrome implicates trisomy 21 (T21) in perturbing fetal hematopoiesis. Here, we show that compared with primary disomic controls, primary T21 fetal liver (FL) hematopoietic stem cells (HSC) and megakaryocyte-erythroid progenitors are markedly increased, whereas granulocyte-macrophage progenitors are reduced. Commensurately, HSC and megakaryocyte-erythroid progenitors show higher clonogenicity, with increased megakaryocyte, megakaryocyte-erythroid, and replatable blast colonies. Biased megakaryocyte-erythroid-primed gene expression was detected as early as the HSC compartment. In lymphopoiesis, T21 FL lymphoid-primed multipotential progenitors and early lymphoid progenitor numbers are maintained, but there was a 10-fold reduction in committed PreproB-lymphoid progenitors and the functional B-cell potential of HSC and early lymphoid progenitor is severely impaired, in tandem with reduced early lymphoid gene expression. The same pattern was seen in all T21 FL samples and no samples had GATA1 mutations. Therefore, T21 itself causes multiple distinct defects in FL myelo-and lymphopoiesis. transient myeloproliferative disorder | aneuploidy | human fetus C onstitutional trisomy 21 (T21) causes Down syndrome (DS), the most common syndrome-associated chromosomal anomaly in humans (1). As well as with neurodevelopmental, cardiac, and gut anomalies (2), there is a striking increase in childhood acute leukemia in DS, even though the risk of solid tumors is much lower than with the general population (3). Intriguingly, this susceptibility to hematopoietic tumors manifests as an increased risk both of acute megakaryocyte (MK)-erythroid leukemia (known as ML-DS) by 150-fold and of acute B-lymphoblastic leukemia (B-ALL) by 33-fold (3, 4).

We present an analysis of protein interaction network data via the comparison of models of network evolution to the observed data. We take a Bayesian approach and perform posterior density estimation using an approximate Bayesian computation with sequential Monte Carlo method. Our approach allows us to perform model selection over a selection of potential network growth models. The methodology we apply uses a distance defined in terms of graph spectra which captures the network data more naturally than previously used summary statistics such as the degree distribution. Furthermore, we include the effects of sampling into the analysis, to properly correct for the incompleteness of existing datasets, and have analysed the performance of our method under various degrees of sampling. We consider a number of models focusing not only on the biologically relevant class of duplication models, but also including models of scale-free network growth that have previously been claimed to describe such data. We find a preference for a duplication-divergence with linear preferential attachment model in the majority of the interaction datasets considered. We also illustrate how our method can be used to perform multi-model inference of network parameters to estimate properties of the full network from sampled data.

The mitogen-activated protein kinase (MAPK) family of proteins is involved in regulating cellular fates such as proliferation, differentiation and apoptosis. In particular, the dynamics of the Erk/Mek system, which has become the canonical example for MAPK signaling systems, have attracted considerable attention. Erk is encoded by two genes, Erk1 and Erk2, that until recently had been considered equivalent as they differ only subtly at the sequence level. However, these proteins exhibit radically different trafficking between cytoplasm and nucleus and this fact may have functional implications. Here we use spatially resolved data on Erk1/2 to develop and analyze spatio-temporal models of these cascades, and we discuss how sensitivity analysis can be used to discriminate between mechanisms. Our models elucidate some of the factors governing the interplay between signaling processes and the Erk1/2 localization in different cellular compartments, including competition between Erk1 and Erk2. Our approach is applicable to a wide range of signaling systems, such as activation cascades, where translocation of molecules occurs. Our study provides a first model of Erk1 and Erk2 activation and their nuclear shuttling dynamics, revealing a role in the regulation of the efficiency of nuclear signaling.

The scientific literature contains a wealth of information about biological systems. Manual curation lacks the scalability to extract this information due to the ever-increasing numbers of papers being published. The development and application of text mining technologies has been proposed as a way of dealing with this problem. However, the inter-species ambiguity of the genomic nomenclature makes mapping of gene mentions identified in text to their corresponding Entrez gene identifiers an extremely difficult task. We propose a novel method, which transforms a MEDLINE record into a mixture of adjacency matrices; by performing a random walkover the resulting graph, we can perform multi-class supervised classification allowing the assignment of taxonomy identifiers to individual gene mentions. The ability to achieve good performance at this task has a direct impact on the performance of normalizing gene mentions to Entrez gene identifiers. Such graph mixtures add flexibility and allow us to generate probabilistic classification schemes that naturally reflect the uncertainties inherent, even in literature-derived data. Results: Our method performs well in terms of both microand macro-averaged performance, achieving micro-F 1 of 0.76 and macro-F 1 of 0.36 on the publicly available DECA corpus. Re-curation of the DECA corpus was performed, with our method achieving 0.88 micro-F 1 and 0.51 macro-F 1 . Our method improves over standard classification techniques [such as support vector machines (SVMs)] in a number of ways: flexibility, interpretability and its resistance to the effects of class bias in the training data. Good performance is achieved without the need for computationally expensive parse tree generation or 'bag of words classification'.

Background: Human T lymphotropic virus Type 1 (HTLV-1) causes a chronic inflammatory disease of the central nervous system known as HTLV-1-associated myelopathy/tropical spastic paraparesis (HAM) which resembles chronic spinal forms of multiple sclerosis (MS). The pathogenesis of HAM remains uncertain. To aid in the differential diagnosis of HAM and to identify pathogenetic mechanisms, we analysed the plasma proteome in asymptomatic HTLV-1 carriers (ACs), patients with HAM, uninfected controls, and patients with MS. We used surface-enhanced laser desorption-ionization (SELDI) mass spectrometry to analyse the plasma proteome in 68 HTLV-1-infected individuals (in two non-overlapping sets, each comprising 17 patients with HAM and 17 ACs), 16 uninfected controls, and 11 patients with secondary progressive MS. Candidate biomarkers were identified by tandem Q-TOF mass spectrometry.