Gene expression profiling in osteoblast biology: bioinformatic tools

This review focuses on using microarray data on a clonal osteoblast cell model to demonstrate how various current and future bioinformatic tools can be used to understand, at a more global or comprehensible level, how cells grow and differentiate. In this example, BMP2 was used to stimulate growth and differentiation of osteoblast to a mineralized matrix. A discussion is included on various methods for clustering gene expression data, statistical evaluation of data, and various new tools that can be used to derive deeper insight into a particular biological problem. How these tools can be obtained is also discussed. New tools for the biologists to compare their datasets with others, as well as examples of future bioinformatic tools that can be used for developing gene networks and pathways for a given set of data are included and discussed.     

Use or abuse of bioinformatic tools: a response to Samach

In a recent paper, we described for the first time the effects of fruit on the expression of putative homologues of genes involved in flowering pathways. It was our aim to provide insight into the molecular mechanisms underlying alternate bearing in citrus. However, a bioinformatics-based critique of our and other related papers has been given by Samach in the preceding Viewpoint article in this issue of Annals of Botany. The use of certain bioinformatic tools in a context of structural rather than functional genomics can cast doubts about the veracity of a large amount of data published in recent years. In this response, the contentions raised by Samach are analysed, and rebuttals of his criticisms are presented.     

Mass spectrometry-based proteomics combined with bioinformatic tools for bacterial classification

Timely classification and identification of bacteria is of vital importance in many areas of public health. We present a mass spectrometry (MS)-based proteomics approach for bacterial classification. In this method, a bacterial proteome database is derived from all potential protein coding open reading frames (ORFs) found in 170 fully sequenced bacterial genomes. Amino acid sequences of tryptic peptides obtained by LC-ESI MS/MS analysis of the digest of bacterial cell extracts are assigned to individual bacterial proteomes in the database. Phylogenetic profiles of these peptides are used to create a matrix of sequence-to-bacterium assignments. These matrixes, viewed as specific assignment bitmaps, are analyzed using statistical tools to reveal the relatedness between a test bacterial sample and the microorganism database. It is shown that, if a sufficient amount of sequence information is obtained from the MS/MS experiments, a bacterial sample can be classified to a strain level by using this proteomics method, leading to its positive identification.     

Inverse pharmacology: Approaches and tools for introducing druggability into engineered proteins

A major feature of twenty-first century medical research is the development of therapeutic strategies that use ‘biologics’ (large molecules, usually engineered proteins) and living cells instead of, or as well as, the small molecules that were the basis of pharmacology in earlier eras. The high power of these techniques can bring correspondingly high risk, and therefore the need for the potential for external control. One way of exerting control on therapeutic proteins is to make them responsive to small molecules; in a clinical context, these small molecules themselves have to be safe. Conventional pharmacology has resulted in thousands of small molecules licensed for use in humans, and detailed structural data on their binding to their protein targets. In principle, these data can be used to facilitate the engineering of drug-responsive modules, taken from natural proteins, into synthetic proteins. This has been done for some years (for example, Cre-ERT2) but usually in a painstaking manner. Recently, we have developed the bioinformatic tool SynPharm to facilitate the design of drug-responsive proteins. In this review, we outline the history of the field, the design and use of the Synpharm tool, and describe our own experiences in engineering druggability into the Cpf1 effector of CRISPR gene editing.     

Next-generation teaching: a template for bringing genomic and bioinformatic tools into the classroom

The recent increase in accessibility and scale of genetic data available through next-generation sequencing (NGS) technology has transformed biological inquiry. As a direct result, the application and analysis of NGS data has quickly become an important skill for future scientists. However, the steep learning curve for applying NGS technology to biological questions, including the complexity of sample preparation for sequencing and the analysis of large data sets, are deterrents to the integration of NGS into undergraduate education. Here, we present a course-based undergraduate research experience (CURE) designed to aid in overcoming these limitations through NGS investigations of prokaryotic diversity. Specifically, we use 16S rRNA sequencing to explore patterns of diversity stemming from student-directed hypothesis development. This CURE addresses three learning objectives: (1) it provides a forum for experimental design hypothesis generation, (2) it introduces modern genomic tools through a hands-on experience generating an NGS data-set, and (3) it provides students with an introductory experience in bioinformatics.     

New bioinformatic tools for analysis of nucleotide modifications in eukaryotic rRNA

This report presents a valuable new bioinformatics package for research on rRNA nucleotide modifications in the ribosome, especially those created by small nucleolar RNA:protein complexes (snoRNPs). The interactive service, which is not available elsewhere, enables a user to visualize the positions of pseudouridines, 2'-O-methylations, and base methylations in three-dimensional space in the ribosome and also in linear and secondary structure formats of ribosomal RNA. Our tools provide additional perspective on where the modifications occur relative to functional regions within the rRNA and relative to other nearby modifications. This package of new tools is presented as a major enhancement of an existing but significantly upgraded yeast snoRNA database available publicly at http://people.biochem.umass.edu/sfournier/fournierlab/snornadb/. The other key features of the enhanced database include details of the base pairing of snoRNAs with target RNAs, genomic organization of the yeast snoRNA genes, and information on corresponding snoRNAs and modifications in other model organisms.     

ALLERDB database and integrated bioinformatic tools for assessment of allergenicity and allergic cross-reactivity

Databases and computational tools are increasingly important in the study of allergies, particularly in the assessment of allergenicity and allergic cross-reactivity. ALLERDB database contains sequences of allergens and information on reported cross-reactivity between allergens. It focuses on analysis of allergenicity and allergic cross-reactivity of clinically relevant protein allergens. The official IUIS allergen data were extracted from the IUIS Allergen Nomenclature Sub-Committee website, and their sequence information from the public databases, and reference publications. The analysis tools assist allergen data analysis and retrieval, and include keyword searching, BLAST, prediction of allergenicity, modification of BLAST that displays cross-reactive allergens, and graphics representation of cross-reactivity data. ALLERDB is new brand of allergen databases with a rich set of tools for sequence comparison, pattern identification, and visualization of results. It is accessible at http://research.i2r.a-star.edu.sg/Templar/DB/Allergen.     

Integration of QTL and bioinformatic tools to identify candidate genes for triglycerides in mice

To identify genetic loci influencing lipid levels, we performed quantitative trait loci (QTL) analysis between inbred mouse strains MRL/MpJ and SM/J, measuring triglyceride levels at 8 weeks of age in F2 mice fed a chow diet. We identified one significant QTL on chromosome (Chr) 15 and three suggestive QTL on Chrs 2, 7, and 17. We also carried out microarray analysis on the livers of parental strains of 282 F2 mice and used these data to find cis-regulated expression QTL. We then narrowed the list of candidate genes under significant QTL using a "toolbox" of bioinformatic resources, including haplotype analysis; parental strain comparison for gene expression differences and nonsynonymous coding single nucleotide polymorphisms (SNP); cis-regulated eQTL in livers of F2 mice; correlation between gene expression and phenotype; and conditioning of expression on the phenotype. We suggest Slc25a7 as a candidate gene for the Chr 7 QTL and, based on expression differences, five genes (Polr3 h, Cyp2d22, Cyp2d26, Tspo, and Ttll12) as candidate genes for Chr 15 QTL. This study shows how bioinformatics can be used effectively to reduce candidate gene lists for QTL related to complex traits.     

RNA structure: bioinformatic analysis

The range of functions ascribed to RNA molecules has grown considerably during recent years. Consequently, the analysis and comparison of RNA sequences have become recurrent tasks in molecular biology. Because the biological function of an RNA is expressed more by its folded architecture than by its sequence, original computational tools adapted to the multifaceted RNA functions have to be developed. Such tools, recently published, enable a user to solve classical problems related to RNA research: constructing 'structural' multiple alignments, inferring complete structures and structural motifs from RNA alignments, or searching structural homology in genomic databases.     

Gaining insights into cancer biology through exploration of the cancer secretome using proteomic and bioinformatic tools

Introduction: Tumor-associated proteins released by cancer cells and by tumor stroma cells, referred as ‘cancer secretome’, represent a valuable resource for discovery of potential cancer biomarkers. The last decade was marked by a great increase in number of studies focused on various aspects of cancer secretome including, composition and identification of components externalized by malignant cells and by the components of tumor microenvironment.

Areas covered: Here, we provide an overview of achievements in the proteomic analysis of the cancer secretome, elicited through the tumor-associated interstitial fluid recovered from malignant tissues ex vivo or the protein component of conditioned media obtained from cultured cancer cells in vitro. We summarize various bioinformatic tools and approaches and critically appraise their outcomes, focusing on problems and challenges that arise when applied for the analysis of cancer secretomic databases.

Expert commentary: Recent achievements in the omics- analysis of structural and metabolic aspects of altered cancer secretome contribute greatly to the various hallmarks of cancer including the identification of clinically significant biomarkers and potential targets for therapeutic intervention.     

Characterization of the core and accessory genomes of Pseudomonas aeruginosa using bioinformatic tools Spine and AGEnt

Background

Pseudomonas aeruginosa is an important opportunistic pathogen responsible for many infections in hospitalized and immunocompromised patients. Previous reports estimated that approximately 10% of its 6.6 Mbp genome varies from strain to strain and is therefore referred to as “accessory genome”. Elements within the accessory genome of P. aeruginosa have been associated with differences in virulence and antibiotic resistance. As whole genome sequencing of bacterial strains becomes more widespread and cost-effective, methods to quickly and reliably identify accessory genomic elements in newly sequenced P. aeruginosa genomes will be needed.

Results

We developed a bioinformatic method for identifying the accessory genome of P. aeruginosa. First, the core genome was determined based on sequence conserved among the completed genomes of twelve reference strains using Spine, a software program developed for this purpose. The core genome was 5.84 Mbp in size and contained 5,316 coding sequences. We then developed an in silico genome subtraction program named AGEnt to filter out core genomic sequences from P. aeruginosa whole genomes to identify accessory genomic sequences of these reference strains. This analysis determined that the accessory genome of P. aeruginosa ranged from 6.9-18.0% of the total genome, was enriched for genes associated with mobile elements, and was comprised of a majority of genes with unknown or unclear function. Using these genomes, we showed that AGEnt performed well compared to other publically available programs designed to detect accessory genomic elements. We then demonstrated the utility of the AGEnt program by applying it to the draft genomes of two previously unsequenced P. aeruginosa strains, PA99 and PA103.

Conclusions

The P. aeruginosa genome is rich in accessory genetic material. The AGEnt program accurately identified the accessory genomes of newly sequenced P. aeruginosa strains, even when draft genomes were used. As P. aeruginosa genomes become available at an increasingly rapid pace, this program will be useful in cataloging the expanding accessory genome of this bacterium and in discerning correlations between phenotype and accessory genome makeup. The combination of Spine and AGEnt should be useful in defining the accessory genomes of other bacterial species as well.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-737) contains supplementary material, which is available to authorized users.     

The GMOD Drupal bioinformatic server framework

MOTIVATION: Next-generation sequencing technologies have led to the widespread use of -omic applications. As a result, there is now a pronounced bioinformatic bottleneck. The general model organism database (GMOD) tool kit (http://gmod.org) has produced a number of resources aimed at addressing this issue. It lacks, however, a robust online solution that can deploy heterogeneous data and software within a Web content management system (CMS). RESULTS: We present a bioinformatic framework for the Drupal CMS. It consists of three modules. First, GMOD-DBSF is an application programming interface module for the Drupal CMS that simplifies the programming of bioinformatic Drupal modules. Second, the Drupal Bioinformatic Software Bench (biosoftware_bench) allows for a rapid and secure deployment of bioinformatic software. An innovative graphical user interface (GUI) guides both use and administration of the software, including the secure provision of pre-publication datasets. Third, we present genes4all_experiment, which exemplifies how our work supports the wider research community. Conclusion: Given the infrastructure presented here, the Drupal CMS may become a powerful new tool set for bioinformaticians. The GMOD-DBSF base module is an expandable community resource that decreases development time of Drupal modules for bioinformatics. The biosoftware_bench module can already enhance biologists' ability to mine their own data. The genes4all_experiment module has already been responsible for archiving of more than 150 studies of RNAi from Lepidoptera, which were previously unpublished. Availability and implementation: Implemented in PHP and Perl. Freely available under the GNU Public License 2 or later from http://gmod-dbsf.googlecode.com.     

Database and tools for metabolic network analysis

Metabolic network analysis has attracted much attention in the area of systems biology. It has a profound role in understanding the key features of organism metabolic networks and has been successfully applied in several fields of systems biology, including in silico gene knockouts, production yield improvement using engineered microbial strains, drug target identification, and phenotype prediction. A variety of metabolic network databases and tools have been developed in order to assist research in these fields. Databases that comprise biochemical data are normally integrated with the use of metabolic network analysis tools in order to give a more comprehensive result. This paper reviews and compares eight databases as well as twenty one recent tools. The aim of this review is to study the different types of tools in terms of the features and usability, as well as the databases in terms of the scope and data provided. These tools can be categorised into three main types: standalone tools; toolbox-based tools; and web-based tools. Furthermore, comparisons of the databases as well as the tools are also provided to help software developers and users gain a clearer insight and a better understanding of metabolic network analysis. Additionally, this review also helps to provide useful information that can be used as guidance in choosing tools and databases for a particular research interest.     

Correlating observed odds ratios from lung cancer case-control studies to SNP functional scores predicted by bioinformatic tools

Bioinformatic tools are widely utilized to predict functional single nucleotide polymorphisms (SNPs) for genotyping in molecular epidemiological studies. However, the extent to which these approaches are mirrored by epidemiological findings has not been fully explored. In this study, we first surveyed SNPs examined in case-control studies of lung cancer, the most extensively studied cancer type. We then computed SNP functional scores using four popular bioinformatics tools: SIFT, PolyPhen, SNPs3D, and PMut, and determined their predictive potential using the odds ratios (ORs) reported. Spearman's correlation coefficient (r) for the association with SNP score from SIFT, PolyPhen, SNPs3D, and PMut, and the summary ORs were r=-0.36 (p=0.007), r=0.25 (p=0.068), r=-0.20 (p=0.205), and r=-0.12 (p=0.370), respectively. By creating a combined score using information from all four tools we were able to achieve a correlation coefficient of r=0.51 (p<0.001). These results indicate that scores of predicted functionality could explain a certain fraction of the lung cancer risk detected in genetic association studies and more accurate predictions may be obtained by combining information from a variety of tools. Our findings suggest that bioinformatic tools are useful in predicting SNP functionality and may facilitate future genetic epidemiological studies.     

Clarifying the signal network of salvianolic acid B using proteomic assay and bioinformatic analysis

Salvianolic acid B (SB) is a natural compound with protective effect against ischemia-reperfusion heart injury. However, the signal network of SB including both direct target proteins and downstream signal-related proteins has not been clarified. In the present study, epidermal growth factor receptor (EGFR) was predicted to be the most possible direct protein target of SB by INVDOCK, a ligand-protein inverse-docking algorithm. Possible signal-related proteins of SB in H9C2 cells, including both under normal condition and under ischemia-reperfusion injury, were searched using 2-DE analysis. Totally, 14 signal-related proteins were found. Finally, signal network from EGFR to the signal-related proteins was established using bioinformatic analysis. Interestingly, 9 of the 14 signal-related proteins could be included in a network together with EGFR through direct interaction or only one intermediate partner. The signal cascade from EGFR to heat shock protein 27 (HSP27) and mitofilin (IMMT, inner membrane mitochondrial protein) might be the most important cascade. The signal network was certified by measuring the binding affinity of SB to EGFR in vitro, the effect of SB on internalization and phosphorylation of EGFR, the effect of SB on viability and proliferation of H9C2 cells, and the expression of inner membrane mitochondrial protein in the presence of EGFR inhibitor AG 1478.     

Genetic code: introducing pyrrolysine

Monomethylamine methyltransferase of the archaebacterium Methanosarcina barkeri contains a novel amino acid, pyrrolysine, encoded by the termination codon UAG. Initial studies suggest that pyrrolysine may be co-translationally inserted during protein synthesis, probably by a mechanism analogous to that operating during selenocysteine incorporation.