SBML export interface for the systems biology toolbox for MATLAB

In this application note, we present an Systems biology markup language (SBML) export interface for the Systems Biology Toolbox for MATLAB. This interface allows modelers to automatically convert models, represented in the toolbox's own format (SBmodels) to SBML files. Since SBmodels do not explicitly contain all the information that is required to generate SBML, the necessary information is gathered by parsing SBmodels. The export can be done in two different ways. First, it is possible to call the export from the command line, thereby directly converting a model to an SBML file. The second option is to inspect and edit the conversion results with the help of a graphical user interface and to subsequently export the model to SBML. Availability: The SBML export interface has been integrated into the Systems Biology Toolbox for MATLAB, which is open source and freely available from http://www.sbtoolbox.org. The website also contains a tutorial, extensive documentation and examples.     

jmzIdentML API: A Java interface to the mzIdentML standard for peptide and protein identification data

We present a Java application programming interface (API), jmzIdentML, for the Human Proteome Organisation (HUPO) Proteomics Standards Initiative (PSI) mzIdentML standard for peptide and protein identification data. The API combines the power of Java Architecture of XML Binding (JAXB) and an XPath-based random-access indexer to allow a fast and efficient mapping of extensible markup language (XML) elements to Java objects. The internal references in the mzIdentML files are resolved in an on-demand manner, where the whole file is accessed as a random-access swap file, and only the relevant piece of XMLis selected for mapping to its corresponding Java object. The APIis highly efficient in its memory usage and can handle files of arbitrary sizes. The APIfollows the official release of the mzIdentML (version 1.1) specifications and is available in the public domain under a permissive licence at http://www.code.google.com/p/jmzidentml/.     

Regulation of oxidative phosphorylation by adrenaline and insulin, effect of the insulin-dependent cytoplasmic factor

The insulin-dependent cytoplasmic factor-regulator (IDR) was shown to inhibit pyruvate and succinate oxidation by the rat liver mitochondria. In view of the fact that insulin increases and adrenaline decreases IDR activity in the liver cytoplasm it is suggested that oxidation of the substrate in mitochondria in vivo is regulated by changes in IDR content in cytoplasm. It was shown that adrenaline-activated oxidation apart from decreased activity of IDR in cytoplasm is induced by a different mechanism not related to IDR content.     

GOfetcher: a database with complex searching facility for gene ontology

MOTIVATION: An important contribution to the Gene Ontology (GO) project is to develop tools that facilitate the creation, maintenance and use of ontologies. Several tools have been created for communicating and using the GO project. However, a limitation with most of these tools is that they suffer from lack of a comprehensive search facility. We developed a web application, GOfetcher, with a very comprehensive search facility for the GO project and a variety of output formats for the results. GOfetcher has three different levels for searching the GO: 'Quick Search', 'Advanced Search' and 'Upload Files' for searching. The application includes a unique search option which generates gene information given a nucleotide or protein accession number which can then be used in generating GO information. The output data in GOfetcher can be saved into several different formats; including spreadsheet, comma-separated values and the extensible markup language (XML) format. The database is available at http://mcbc.usm.edu/gofetcher/.     

NCL: a C++ class library for interpreting data files in NEXUS format

The NEXUS Class Library (NCL) is a collection of C++ classes designed to simplify interpreting data files written in the NEXUS format used by many computer programs for phylogenetic analyses. The NEXUS format allows different programs to share the same data files, even though none of the programs can interpret all of the data stored therein. Because users are not required to reformat the data file for each program, use of the NEXUS format prevents cut-and-paste errors as well as the proliferation of copies of the original data file. The purpose of making the NCL available is to encourage the use of the NEXUS format by making it relatively easy for programmers to add the ability to interpret NEXUS files in newly developed software. AVAILABILITY: The NCL is freely available under the GNU General Public License from http://hydrodictyon.eeb.uconn.edu/ncl/ Supplementary information: Documentation for the NCL (general information and source code documentation) is available in HTML format at http://hydrodictyon.eeb.uconn.edu/ncl/     

cluML

cluML is a new markup language for microarray data clustering and cluster validity assessment. The XML-based format has been designed to address some of the limitations observed in traditional formats, such as inability to store multiple clustering (including biclustering) and validation results within a dataset. cluML is an effective tool to support biomedical knowledge representation in gene expression data analysis. Although cluML was developed for DNA microarray analysis applications, it can be effectively used for the representation of clustering and for the validation of other biomedical and physical data that has no limitations.     

Adaptive informatics for multifactorial and high-content biological data

Whereas genomic data are universally machine-readable, data from imaging, multiplex biochemistry, flow cytometry and other cell- and tissue-based assays usually reside in loosely organized files of poorly documented provenance. This arises because the relational databases used in genomic research are difficult to adapt to rapidly evolving experimental designs, data formats and analytic algorithms. Here we describe an adaptive approach to managing experimental data based on semantically typed data hypercubes (SDCubes) that combine hierarchical data format 5 (HDF5) and extensible markup language (XML) file types. We demonstrate the application of SDCube-based storage using ImageRail, a software package for high-throughput microscopy. Experimental design and its day-to-day evolution, not rigid standards, determine how ImageRail data are organized in SDCubes. We applied ImageRail to collect and analyze drug dose-response landscapes in human cell lines at single-cell resolution.     

A prototype object database for mitochondrial DNA variation

Surveys of biochemical and molecular genetic variation in natural populations have generated a wealth of data, but this valuable resource has not been adequately preserved. We hope to prevent further loss by establishing a community database for population genetic surveys. We explored the feasibility of a population genetics database by developing a prototype for animal mitochondrial DNA (mtDNA) surveys. This prototype includes the specification of a format for data files that are to be submitted to the database, an open-source object database that encapsulates data with methods to display and analyze data, and a website where data can be retrieved in either its original form or extensible markup language (XML). Data from more than 50 published surveys of mtDNA variation were retrieved from the literature and entered into the database. We hope that the population genetics community will support this project by contributing both data and expertise.     

A large disordered region confers a wide spanning volume to vertebrate Suppressor of Fused as shown in a trans-species solution study

Hedgehog (Hh) pathway inhibition by the conserved protein Suppressor of Fused (SuFu) is crucial to vertebrate development. By constrast, SuFu loss-of-function mutant has little effect in drosophila.Previous publications showed that the crystal structures of human and drosophila SuFu consist of two ordered domains that are capable of breathing motions upon ligand binding. However, the crystal structure of human SuFu does not give information about twenty N-terminal residues (IDR1) and an eighty-residue-long region predicted as disordered (IDR2) in the C-terminus, whose function is important for the pathway repression. These two intrinsically disordered regions (IDRs) are species-dependent.To obtain information about the IDR regions, we studied full-length SuFu’s structure in solution, both with circular dichroism and small angle X-ray scattering, comparing drosophila, zebrafish and human species, to better understand this considerable difference. Our studies show that, in spite of similar crystal structures restricted to ordered domains, drosophila and vertebrate SuFu have very different structures in solution. The IDR2 of vertebrates spans a large area, thus enabling it to reach for partners and be accessible for post-translational modifications. Furthermore, we show that the IDR2 region is highly conserved within phyla but varies in length and sequence, with insects having a shorter disordered region while that of vertebrates is broad and mobile. This major variation may explain the different phenotypes observed upon SuFu removal.     

Java editor for biological pathways

SUMMARY: A visual Java-based tool for drawing and annotating biological pathways was developed. This tool integrates the possibilities of charting elements with different attributes (size, color, labels), drawing connections between elements in distinct characteristics (color, structure, width, arrows), as well as adding links to molecular biology databases, promoter sequences, information on the function of the genes or gene products, and references. It is easy to use and system independent. The result of the editing process is a PNG (portable network graphics) file for the images and XML (extended markup language) file for the appropriate links.     

SPREAD: spatial phylogenetic reconstruction of evolutionary dynamics

SUMMARY: SPREAD is a user-friendly, cross-platform application to analyze and visualize Bayesian phylogeographic reconstructions incorporating spatial-temporal diffusion. The software maps phylogenies annotated with both discrete and continuous spatial information and can export high-dimensional posterior summaries to keyhole markup language (KML) for animation of the spatial diffusion through time in virtual globe software. In addition, SPREAD implements Bayes factor calculation to evaluate the support for hypotheses of historical diffusion among pairs of discrete locations based on Bayesian stochastic search variable selection estimates. SPREAD takes advantage of multicore architectures to process large joint posterior distributions of phylogenies and their spatial diffusion and produces visualizations as compelling and interpretable statistical summaries for the different spatial projections. AVAILABILITY: SPREAD is licensed under the GNU Lesser GPL and its source code is freely available as a GitHub repository: https://github.com/phylogeography/SPREAD CONTACT: filip.bielejec@rega.kuleuven.be.     

Protein kinases phosphorylate long disordered regions in intrinsically disordered proteins

Phosphorylation is a major post‐translational modification that plays a central role in signaling pathways. Protein kinases phosphorylate substrates (phosphoproteins) by adding phosphate at Ser/Thr or Tyr residues (phosphosites). A large amount of data identifying and describing phosphosites in phosphoproteins has been reported but the specificity of phosphorylation is not fully resolved. In this report, data of kinase‐substrate pairs identified by the Kinase‐Interacting Substrate Screening (KISS) method were used to analyze phosphosites in intrinsically disordered regions (IDRs) of intrinsically disordered proteins. We compared phosphorylated and nonphosphorylated IDRs and found that the phosphorylated IDRs were significantly longer than nonphosphorylated IDRs. The phosphorylated IDR is often the longest IDR (71%) in a phosphoprotein when only a single phosphosite exists in the IDR, and when the phosphoprotein has multiple phosphosites in an IDR(s), the phosphosites are primarily localized in a single IDR (78%) and this IDR is usually the longest one (81%). We constructed a stochastic model of phosphorylation to estimate the effect of IDR length. The model that accounted for IDR length produced more realistic results when compared with a model that excluded the IDR length. We propose that the IDR length is a significant determinant for locating kinase phosphorylation sites in phosphoproteins.     

Importing MAGE-ML format microarray data into BioConductor

The microarray gene expression markup language (MAGE-ML) is a widely used XML (eXtensible Markup Language) standard for describing and exchanging information about microarray experiments. It can describe microarray designs, microarray experiment designs, gene expression data and data analysis results. We describe RMAGEML, a new Bioconductor package that provides a link between cDNA microarray data stored in MAGE-ML format and the Bioconductor framework for preprocessing, visualization and analysis of microarray experiments. AVAILABILITY: http://www.bioconductor.org. Open Source.     

The jmzQuantML programming interface and validator for the mzQuantML data standard

The mzQuantML standard from the HUPO Proteomics Standards Initiative has recently been released, capturing quantitative data about peptides and proteins, following analysis of MS data. We present a Java application programming interface (API) for mzQuantML called jmzQuantML. The API provides robust bridges between Java classes and elements in mzQuantML files and allows random access to any part of the file. The API provides read and write capabilities, and is designed to be embedded in other software packages, enabling mzQuantML support to be added to proteomics software tools ( http://code.google.com/p/jmzquantml/ ). The mzQuantML standard is designed around a multilevel validation system to ensure that files are structurally and semantically correct for different proteomics quantitative techniques. In this article, we also describe a Java software tool ( http://code.google.com/p/mzquantml‐validator/ ) for validating mzQuantML files, which is a formal part of the data standard.     

Localization of the thyroid peroxidase autoantibody immunodominant region to a junctional region containing portions of the domains homologous to complement control protein and myeloperoxidase

Thyroid peroxidase (TPO) autoantibody epitopes are largely restricted to an immunodominant region (IDR) on the extracellular region of the native molecule. Localization of the IDR has been a longstanding and difficult goal. The TPO extracellular region comprises a large myeloperoxidase-like domain, linked to the plasma membrane by two smaller domains with homology to complement control protein (CCP) and epidermal growth factor (EGF), respectively. Recent studies have focused on the CCP- and EGF-like domains as the putative location of the TPO autoantibody IDR. To address this issue, we attempted to express on the surface of transfected cells native TPO in which the CCP- and EGF-like domains were deleted, either together or individually. We used a quartet of human monoclonal autoantibodies that define the TPO IDR, as well as polyclonal TPO autoantibodies in patients' sera, to detect these mutated TPO molecules by flow cytometry. The combined CCP/EGF-like domain deletion did not produce a signal with TPO autoantibodies but did not traffic to the cell surface. In contrast, both monoclonal and polyclonal autoantibodies recognized TPO with the juxtamembrane EGF-like domain deleted equally as well as the wild-type TPO on the cell surface. TPO with the CCP-like domain deleted expressed normally on the cell surface, as determined using the polyclonal mouse antiserum. Nevertheless, this modified TPO molecule was recognized very poorly by both the human monoclonal autoantibodies and the polyclonal autoantibodies in patients' sera. In conclusion, we have clearly excluded the juxtamembrane EGF-like domain as being part of the IDR. In contrast, a component of the CCP-like domain does contribute to the IDR. These data, together with findings from other studies, localize the TPO autoantibody IDR to the junction of the CCP-like domain and the much larger myeloperoxidase-like domain on TPO.     

The PSI semantic validator: A framework to check MIAPE compliance of proteomics data

The Human Proteome Organization's Proteomics Standards Initiative (PSI) promotes the development of exchange standards to improve data integration and interoperability. PSI specifies the suitable level of detail required when reporting a proteomics experiment (via the Minimum Information About a Proteomics Experiment), and provides extensible markup language (XML) exchange formats and dedicated controlled vocabularies (CVs) that must be combined to generate a standard compliant document. The framework presented here tackles the issue of checking that experimental data reported using a specific format, CVs and public bio‐ontologies (e.g. Gene Ontology, NCBI taxonomy) are compliant with the Minimum Information About a Proteomics Experiment recommendations. The semantic validator not only checks the XML syntax but it also enforces rules regarding the use of an ontology class or CV terms by checking that the terms exist in the resource and that they are used in the correct location of a document. Moreover, this framework is extremely fast, even on sizable data files, and flexible, as it can be adapted to any standard by customizing the parameters it requires: an XML Schema Definition, one or more CVs or ontologies, and a mapping file describing in a formal way how the semantic resources and the format are interrelated. As such, the validator provides a general solution to the common problem in data exchange: how to validate the correct usage of a data standard beyond simple XML Schema Definition validation. The framework source code and its various applications can be found at http://psidev.info/validator .     

Tabular strategies for metadata in ecology,evolution, and the environmental sciences

Data support knowledge development and theory advances in ecology and evolution. We are increasingly reusing data within our teams and projects and through the global, openly archived datasets of others. Metadata can be challenging to write and interpret, but it is always crucial for reuse. The value metadata cannot be overstated—even as a relatively independent research object because it describes the work that has been done in a structured format. We advance a new perspective and classify methods for metadata curation and development with tables. Tables with templates can be effectively used to capture all components of an experiment or project in a single, easy‐to‐read file familiar to most scientists. If coupled with the R programming language, metadata from tables can then be rapidly and reproducibly converted to publication formats including extensible markup language files suitable for data repositories. Tables can also be used to summarize existing metadata and store metadata across many datasets. A case study is provided and the added benefits of tables for metadata, a priori, are developed to ensure a more streamlined publishing process for many data repositories used in ecology, evolution, and the environmental sciences. In ecology and evolution, researchers are often highly tabular thinkers from experimental data collection in the lab and/or field, and representations of metadata as a table will provide novel research and reuse insights.     

GB‐to‐TNT: facilitating creation of matrices from GenBank and diagnosis of results in TNT

This paper presents a pipeline, implemented in an open‐source program called GB→TNT (GenBank‐to‐TNT), for creating large molecular matrices, starting from GenBank files and finishing with TNT matrices which incorporate taxonomic information in the terminal names. GB→TNT is designed to retrieve a defined genomic region from a bulk of sequences included in a GenBank file. The user defines the genomic region to be retrieved and several filters (genome, length of the sequence, taxonomic group, etc.); each genomic region represents a different data block in the final TNT matrix. GB→TNT first generates Fasta files from the input GenBank files, then creates an alignment for each of those (by calling an alignment program), and finally merges all the aligned files into a single TNT matrix. The new version of TNT can make use of the taxonomic information contained in the terminal names, allowing easy diagnosis of results, evaluation of fit between the trees and the taxonomy, and automatic labelling or colouring of tree branches with the taxonomic groups they represent. © The Willi Hennig Society 2012.