Intuitive Visualization and Analysis of Multi-Omics Data and Application to Escherichia coli Carbon Metabolism

Combinations of ‘omics’ investigations (i.e, transcriptomic, proteomic, metabolomic and/or fluxomic) are increasingly applied to get comprehensive understanding of biological systems. Because the latter are organized as complex networks of molecular and functional interactions, the intuitive interpretation of multi-omics datasets is difficult. Here we describe a simple strategy to visualize and analyze multi-omics data. Graphical representations of complex biological networks can be generated using Cytoscape where all molecular and functional components could be explicitly represented using a set of dedicated symbols. This representation can be used i) to compile all biologically-relevant information regarding the network through web link association, and ii) to map the network components with multi-omics data. A Cytoscape plugin was developed to increase the possibilities of both multi-omic data representation and interpretation. This plugin allowed different adjustable colour scales to be applied to the various omics data and performed the automatic extraction and visualization of the most significant changes in the datasets. For illustration purpose, the approach was applied to the central carbon metabolism of Escherichia coli. The obtained network contained 774 components and 1232 interactions, highlighting the complexity of bacterial multi-level regulations. The structured representation of this network represents a valuable resource for systemic studies of E. coli, as illustrated from the application to multi-omics data. Some current issues in network representation are discussed on the basis of this work

A fragile metabolic network adapted for cooperation in the symbiotic bacterium Buchnera aphidicola

<p>Abstract</p> <p>Background</p> <p><it>In silico </it>analyses provide valuable insight into the biology of obligately intracellular pathogens and symbionts with small genomes. There is a particular opportunity to apply systems-level tools developed for the model bacterium <it>Escherichia coli </it>to study the evolution and function of symbiotic bacteria which are metabolically specialised to overproduce specific nutrients for their host and, remarkably, have a gene complement that is a subset of the <it>E. coli </it>genome.</p> <p>Results</p> <p>We have reconstructed and analysed the metabolic network of the γ-proteobacterium <it>Buchnera aphidicola </it>(symbiont of the pea aphid) as a model for using systems-level approaches to discover key traits of symbionts with small genomes. The metabolic network is extremely fragile with > 90% of the reactions essential for viability <it>in silico</it>; and it is structured so that the bacterium cannot grow without producing the essential amino acid, histidine, which is released to the insect host. Further, the amount of essential amino acid produced by the bacterium <it>in silico </it>can be controlled by host supply of carbon and nitrogen substrates.</p> <p>Conclusion</p> <p>This systems-level analysis predicts that the fragility of the bacterial metabolic network renders the symbiotic bacterium intolerant of drastic environmental fluctuations, whilst the coupling of histidine production to growth prevents the bacterium from exploiting host nutrients without reciprocating. These metabolic traits underpin the sustained nutritional contribution of <it>B. aphidicola </it>to the host and, together with the impact of host-derived substrates on the profile of nutrients released from the bacteria, point to a dominant role of the host in controlling the symbiosis.</p

Use of reconstituted metabolic networks to assist in metabolomic data visualization and mining

Metabolomics experiments seldom achieve their aim of comprehensively covering the entire metabolome. However, important information can be gleaned even from sparse datasets, which can be facilitated by placing the results within the context of known metabolic networks. Here we present a method that allows the automatic assignment of identified metabolites to positions within known metabolic networks, and, furthermore, allows automated extraction of sub-networks of biological significance. This latter feature is possible by use of a gap-filling algorithm. The utility of the algorithm in reconstructing and mining of metabolomics data is shown on two independent datasets generated with LC–MS LTQ-Orbitrap mass spectrometry. Biologically relevant metabolic sub-networks were extracted from both datasets. Moreover, a number of metabolites, whose presence eluded automatic selection within mass spectra, could be identified retrospectively by virtue of their inferred presence through gap filling

Industrial Systems Biology of Saccharomyces cerevisiae Enables Novel Succinic Acid Cell Factory.

Saccharomyces cerevisiae is the most well characterized eukaryote, the preferred microbial cell factory for the largest industrial biotechnology product (bioethanol), and a robust commerically compatible scaffold to be exploitted for diverse chemical production. Succinic acid is a highly sought after added-value chemical for which there is no native pre-disposition for production and accmulation in S. cerevisiae. The genome-scale metabolic network reconstruction of S. cerevisiae enabled in silico gene deletion predictions using an evolutionary programming method to couple biomass and succinate production. Glycine and serine, both essential amino acids required for biomass formation, are formed from both glycolytic and TCA cycle intermediates. Succinate formation results from the isocitrate lyase catalyzed conversion of isocitrate, and from the alpha-keto-glutarate dehydrogenase catalyzed conversion of alpha-keto-glutarate. Succinate is subsequently depleted by the succinate dehydrogenase complex. The metabolic engineering strategy identified included deletion of the primary succinate consuming reaction, Sdh3p, and interruption of glycolysis derived serine by deletion of 3-phosphoglycerate dehydrogenase, Ser3p/Ser33p. Pursuing these targets, a multi-gene deletion strain was constructed, and directed evolution with selection used to identify a succinate producing mutant. Physiological characterization coupled with integrated data analysis of transcriptome data in the metabolically engineered strain were used to identify 2nd-round metabolic engineering targets. The resulting strain represents a 30-fold improvement in succinate titer, and a 43-fold improvement in succinate yield on biomass, with only a 2.8-fold decrease in the specific growth rate compared to the reference strain. Intuitive genetic targets for either over-expression or interruption of succinate producing or consuming pathways, respectively, do not lead to increased succinate. Rather, we demonstrate how systems biology tools coupled with directed evolution and selection allows non-intuitive, rapid and substantial re-direction of carbon fluxes in S. cerevisiae, and hence show proof of concept that this is a potentially attractive cell factory for over-producing different platform chemicals

First observations of separated atmospheric nu_mu and bar{nu-mu} events in the MINOS detector

The complete 5.4 kton MINOS far detector has been taking data since the beginning of August 2003 at a depth of 2070 meters water-equivalent in the Soudan mine, Minnesota. This paper presents the first MINOS observations of nuµ and [overline nu ]µ charged-current atmospheric neutrino interactions based on an exposure of 418 days. The ratio of upward- to downward-going events in the data is compared to the Monte Carlo expectation in the absence of neutrino oscillations, giving Rup/downdata/Rup/downMC=0.62-0.14+0.19(stat.)±0.02(sys.). An extended maximum likelihood analysis of the observed L/E distributions excludes the null hypothesis of no neutrino oscillations at the 98% confidence level. Using the curvature of the observed muons in the 1.3 T MINOS magnetic field nuµ and [overline nu ]µ interactions are separated. The ratio of [overline nu ]µ to nuµ events in the data is compared to the Monte Carlo expectation assuming neutrinos and antineutrinos oscillate in the same manner, giving R[overline nu ][sub mu]/nu[sub mu]data/R[overline nu ][sub mu]/nu[sub mu]MC=0.96-0.27+0.38(stat.)±0.15(sys.), where the errors are the statistical and systematic uncertainties. Although the statistics are limited, this is the first direct observation of atmospheric neutrino interactions separately for nuµ and [overline nu ]µ

Phenotypic and Transcriptomic Response of Auxotrophic Mycobacterium avium Subsp. paratuberculosis leuD Mutant under Environmental Stress

Mycobacterium avium subsp. paratuberculosis (MAP) is the causative agent of severe gastroenteritis in cattle. To gain a better understanding of MAP virulence, we investigated the role of leuD gene in MAP metabolism and stress response. For this, we have constructed an auxotrophic strain of MAP by deleting the leuD gene using allelic exchange. The wildtype and mutant strains were then compared for metabolic phenotypic changes using Biolog phenotype microarrays. The responses of both strains to physiologically relevant stress conditions were assessed using DNA microarrays. Transcriptomic data was then analyzed in the context of cellular metabolic pathways and gene networks. Our results showed that deletion of leuD gene has a global effect on both MAP phenotypic and transcriptome response. At the metabolic level, the mutant strain lost the ability to utilize most of the carbon, nitrogen, sulphur, phosphorus and nutrient supplements as energy source. At the transcriptome level, more than 100 genes were differentially expressed in each of the stress condition tested. Systems level network analysis revealed that the differentially expressed genes were distributed throughout the gene network, thus explaining the global impact of leuD deletion in metabolic phenotype. Further, we find that leuD deletion impacted metabolic pathways associated with fatty acids. We verified this by experimentally estimating the total fatty acid content of both mutant and wildtype. The mutant strain had 30% less fatty acid content when compared to wildtype, thus supporting the results from transcriptional and computational analyses. Our results therefore reveal the intricate connection between the metabolism and virulence in MAP

Characterizing Complex Polysera Produced by Antigen-Specific Immunization through the Use of Affinity-Selected Mimotopes

BACKGROUND: Antigen-based (as opposed to whole organism) vaccines are actively being pursued for numerous indications. Even though different formulations may produce similar levels of total antigen-specific antibody, the composition of the antibody response can be quite distinct resulting in different levels of therapeutic activity. METHODOLOGY/PRINCIPAL FINDINGS: Using plasmid-based immunization against the proto-oncogene HER-2 as a model, we have demonstrated that affinity-selected epitope mimetics (mimotopes) can provide a defined signature of a polyclonal antibody response. Further, using novel computer algorithms that we have developed, these mimotopes can be used to predict epitope targets. CONCLUSIONS/SIGNIFICANCE: By combining our novel strategy with existing methods of epitope prediction based on physical properties of an individual protein, we believe that this method offers a robust method for characterizing the breadth of epitope-specificity within a specific polyserum. This strategy is useful as a tool for monitoring immunity following vaccination and can also be used to define relevant epitopes for the creation of novel vaccines

Proteomic Characterization of Cellular and Molecular Processes that Enable the Nanoarchaeum equitans-Ignicoccus hospitalis Relationship

Nanoarchaeum equitans, the only cultured representative of the Nanoarchaeota, is dependent on direct physical contact with its host, the hyperthermophile Ignicoccus hospitalis. The molecular mechanisms that enable this relationship are unknown. Using whole-cell proteomics, differences in the relative abundance of >75% of predicted protein-coding genes from both Archaea were measured to identify the specific response of I. hospitalis to the presence of N. equitans on its surface. A purified N. equitans sample was also analyzed for evidence of interspecies protein transfer. The depth of cellular proteome coverage achieved here is amongst the highest reported for any organism. Based on changes in the proteome under the specific conditions of this study, I. hospitalis reacts to N. equitans by curtailing genetic information processing (replication, transcription) in lieu of intensifying its energetic, protein processing and cellular membrane functions. We found no evidence of significant Ignicoccus biosynthetic enzymes being transported to N. equitans. These results suggest that, under laboratory conditions, N. equitans diverts some of its host's metabolism and cell cycle control to compensate for its own metabolic shortcomings, thus appearing to be entirely dependent on small, transferable metabolites and energetic precursors from I. hospitalis

Pathway Projector: Web-Based Zoomable Pathway Browser Using KEGG Atlas and Google Maps API

BACKGROUND: Biochemical pathways provide an essential context for understanding comprehensive experimental data and the systematic workings of a cell. Therefore, the availability of online pathway browsers will facilitate post-genomic research, just as genome browsers have contributed to genomics. Many pathway maps have been provided online as part of public pathway databases. Most of these maps, however, function as the gateway interface to a specific database, and the comprehensiveness of their represented entities, data mapping capabilities, and user interfaces are not always sufficient for generic usage. METHODOLOGY/PRINCIPAL FINDINGS: We have identified five central requirements for a pathway browser: (1) availability of large integrated maps showing genes, enzymes, and metabolites; (2) comprehensive search features and data access; (3) data mapping for transcriptomic, proteomic, and metabolomic experiments, as well as the ability to edit and annotate pathway maps; (4) easy exchange of pathway data; and (5) intuitive user experience without the requirement for installation and regular maintenance. According to these requirements, we have evaluated existing pathway databases and tools and implemented a web-based pathway browser named Pathway Projector as a solution. CONCLUSIONS/SIGNIFICANCE: Pathway Projector provides integrated pathway maps that are based upon the KEGG Atlas, with the addition of nodes for genes and enzymes, and is implemented as a scalable, zoomable map utilizing the Google Maps API. Users can search pathway-related data using keywords, molecular weights, nucleotide sequences, and amino acid sequences, or as possible routes between compounds. In addition, experimental data from transcriptomic, proteomic, and metabolomic analyses can be readily mapped. Pathway Projector is freely available for academic users at (http://www.g-language.org/PathwayProjector/)

Aberrations of TACC1 and TACC3 are associated with ovarian cancer

BACKGROUND: Dysregulation of the human Transforming Acidic Coiled Coil (TACC) genes is thought to be important in the development and progression of multiple myeloma, breast and gastric cancer. Recent, large-scale genomic analysis and Serial Analysis of Gene Expression data suggest that TACC1 and TACC3 may also be involved in the etiology of ovarian tumors from both familial and sporadic cases. Therefore, the aim of this study was to determine the occurrence of alterations of these TACCs in ovarian cancer. METHODS: Detection and scoring of TACC1 and TACC3 expression was performed by immunohistochemical analysis of the T-BO-1 tissue/tumor microarray slide from the Cooperative Human Tissue Network, Tissue Array Research Program (TARP) of the National Cancer Institute, National Institutes of Health, Bethesda, MD, USA. Tumors were categorized as either positive (greater than 10% of cells staining) or negative. Statistical analysis was performed using Fisher's exact test and p < 0.05 (single comparisons), and p < 0.02 (multiple comparisons) were considered to be significant. Transgenomics WAVE high performance liquid chromatography (dHPLC) was used to pre-screen the TACC3 gene in constitutional DNA from ovarian cancer patients and their unaffected relatives from 76 families from the Gilda Radner Familial Ovarian Cancer Registry. All variant patterns were then sequenced. RESULTS: This study demonstrated absence of at least one or both TACC proteins in 78.5% (51/65) of ovarian tumors tested, with TACC3 loss observed in 67.7% of tumors. The distribution pattern of expression of the two TACC proteins was different, with TACC3 loss being more common in serous papillary carcinoma compared with clear cell carcinomas, while TACC1 staining was less frequent in endometroid than in serous papillary tumor cores. In addition, we identified two constitutional mutations in the TACC3 gene in patients with ovarian cancer from the Gilda Radner Familial Ovarian Cancer Registry. These patients had previously tested negative for mutations in known ovarian cancer predisposing genes. CONCLUSION: When combined, our data suggest that aberrations of TACC genes, and TACC3 in particular, underlie a significant proportion of ovarian cancers. Thus, TACC3 could be a hitherto unknown endogenous factor that contributes to ovarian tumorigenesis