Wynand Alkema

BioVenn – a web application for the comparison and visualization of biological lists using area-proportional Venn diagrams

BackgroundIn many genomics projects, numerous lists containing biological identifiers are produced. Often it is useful to see the overlap between different lists, enabling researchers to quickly observe similarities and differences between the data sets they are analyzing. One of the most popular methods to visualize the overlap and differences between data sets is the Venn diagram: a diagram consisting of two or more circles in which each circle corresponds to a data set, and the overlap between the circles corresponds to the overlap between the data sets. Venn diagrams are especially useful when they are area-proportional i.e. the sizes of the circles and the overlaps correspond to the sizes of the data sets. Currently there are no programs available that can create area-proportional Venn diagrams connected to a wide range of biological databases.ResultsWe designed a web application named BioVenn to summarize the overlap between two or three lists of identifiers, using area-proportional Venn diagrams. The user only needs to input these lists of identifiers in the textboxes and push the submit button. Parameters like colors and text size can be adjusted easily through the web interface. The position of the text can be adjusted by drag-and-drop principle. The output Venn diagram can be shown as an SVG or PNG image embedded in the web application, or as a standalone SVG or PNG image. The latter option is useful for batch queries. Besides the Venn diagram, BioVenn outputs lists of identifiers for each of the resulting subsets. If an identifier is recognized as belonging to one of the supported biological databases, the output is linked to that database. Finally, BioVenn can map Affymetrix and EntrezGene identifiers to Ensembl genes.ConclusionBioVenn is an easy-to-use web application to generate area-proportional Venn diagrams from lists of biological identifiers. It supports a wide range of identifiers from the most used biological databases currently available. Its implementation on the World Wide Web makes it available for use on any computer with internet connection, independent of operating system and without the need to install programs locally. BioVenn is freely accessible at http://www.cmbi.ru.nl/cdd/biovenn/.

Literature mining for the discovery of hidden connections between drugs, genes and diseases

The scientific literature represents a rich source for retrieval of knowledge on associations between biomedical concepts such as genes, diseases and cellular processes. A commonly used method to establish relationships between biomedical concepts from literature is co-occurrence. Apart from its use in knowledge retrieval, the co-occurrence method is also well-suited to discover new, hidden relationships between biomedical concepts following a simple ABC-principle, in which A and C have no direct relationship, but are connected via shared B-intermediates. In this paper we describe CoPub Discovery, a tool that mines the literature for new relationships between biomedical concepts. Statistical analysis using ROC curves showed that CoPub Discovery performed well over a wide range of settings and keyword thesauri. We subsequently used CoPub Discovery to search for new relationships between genes, drugs, pathways and diseases. Several of the newly found relationships were validated using independent literature sources. In addition, new predicted relationships between compounds and cell proliferation were validated and confirmed experimentally in an in vitro cell proliferation assay. The results show that CoPub Discovery is able to identify novel associations between genes, drugs, pathways and diseases that have a high probability of being biologically valid. This makes CoPub Discovery a useful tool to unravel the mechanisms behind disease, to find novel drug targets, or to find novel applications for existing drugs.

Prednisolone-induced differential gene expression in mouse liver carrying wild type or a dimerization-defective glucocorticoid receptor

BackgroundGlucocorticoids (GCs) control expression of a large number of genes via binding to the GC receptor (GR). Transcription may be regulated either by binding of the GR dimer to DNA regulatory elements or by protein-protein interactions of GR monomers with other transcription factors. Although the type of regulation for a number of individual target genes is known, the relative contribution of both mechanisms to the regulation of the entire transcriptional program remains elusive. To study the importance of GR dimerization in the regulation of gene expression, we performed gene expression profiling of livers of prednisolone-treated wild type (WT) and mice that have lost the ability to form GR dimers (GRdim).ResultsThe GR target genes identified in WT mice were predominantly related to glucose metabolism, the cell cycle, apoptosis and inflammation. In GRdim mice, the level of prednisolone-induced gene expression was significantly reduced compared to WT, but not completely absent. Interestingly, for a set of genes, involved in cell cycle and apoptosis processes and strongly related to Foxo3a and p53, induction by prednisolone was completely abolished in GRdim mice. In contrast, glucose metabolism-related genes were still modestly upregulated in GRdim mice upon prednisolone treatment. Finally, we identified several novel GC-inducible genes from which Fam107a, a putative histone acetyltransferase complex interacting protein, was most strongly dependent on GR dimerization.ConclusionsThis study on prednisolone-induced effects in livers of WT and GRdim mice identified a number of interesting candidate genes and pathways regulated by GR dimers and sheds new light onto the complex transcriptional regulation of liver function by GCs.

Application of text mining in the biomedical domain

In recent years the amount of experimental data that is produced in biomedical research and the number of papers that are being published in this field have grown rapidly. In order to keep up to date with developments in their field of interest and to interpret the outcome of experiments in light of all available literature, researchers turn more and more to the use of automated literature mining. As a consequence, text mining tools have evolved considerably in number and quality and nowadays can be used to address a variety of research questions ranging from de novo drug target discovery to enhanced biological interpretation of the results from high throughput experiments. In this paper we introduce the most important techniques that are used for a text mining and give an overview of the text mining tools that are currently being used and the type of problems they are typically applied for.

Molecular pathway profiling of T lymphocyte signal transduction pathways; Th1 and Th2 genomic fingerprints are defined by TCR and CD28-mediated signaling

BackgroundT lymphocytes are orchestrators of adaptive immunity. Naïve T cells may differentiate into Th1, Th2, Th17 or iTreg phenotypes, depending on environmental co-stimulatory signals. To identify genes and pathways involved in differentiation of Jurkat T cells towards Th1 and Th2 subtypes we performed comprehensive transcriptome analyses of Jurkat T cells stimulated with various stimuli and pathway inhibitors. Results from these experiments were validated in a human experimental setting using whole blood and purified CD4+ Tcells.ResultsCalcium-dependent activation of T cells using CD3/CD28 and PMA/CD3 stimulation induced a Th1 expression profile reflected by increased expression of T-bet, RUNX3, IL-2, and IFNγ, whereas calcium-independent activation via PMA/CD28 induced a Th2 expression profile which included GATA3, RXRA, CCL1 and Itk. Knock down with siRNA and gene expression profiling in the presence of selective kinase inhibitors showed that proximal kinases Lck and PKCθ are crucial signaling hubs during T helper cell activation, revealing a clear role for Lck in Th1 development and for PKCθ in both Th1 and Th2 development. Medial signaling via MAPkinases appeared to be less important in these pathways, since specific inhibitors of these kinases displayed a minor effect on gene expression. Translation towards a primary, whole blood setting and purified human CD4+ T cells revealed that PMA/CD3 stimulation induced a more pronounced Th1 specific, Lck and PKCθ dependent IFNγ production, whereas PMA/CD28 induced Th2 specific IL-5 and IL-13 production, independent of Lck activation. PMA/CD3-mediated skewing towards a Th1 phenotype was also reflected in mRNA expression of the master transcription factor Tbet, whereas PMA/CD28-mediated stimulation enhanced GATA3 mRNA expression in primary human CD4+ Tcells.ConclusionsThis study identifies stimulatory pathways and gene expression profiles for in vitro skewing of T helper cell activation. PMA/CD3 stimulation enhances a Th1-like response in an Lck and PKCθ dependent fashion, whereas PMA/CD28 stimulation results in a Th2-like phenotype independent of the proximal TCR-tyrosine kinase Lck. This approach offers a robust and fast translational in vitro system for skewed T helper cell responses in Jurkat T cells, primary human CD4+ Tcells and in a more complex matrix such as human whole blood.

Literature-based compound profiling: application to toxicogenomics

INTRODUCTIONnTo reduce continuously increasing costs in drug development, adverse effects of drugs need to be detected as early as possible in the process. In recent years, compound-induced gene expression profiling methodologies have been developed to assess compound toxicity, including Gene Ontology term and pathway over-representation analyses. The objective of this study was to introduce an additional approach, in which literature information is used for compound profiling to evaluate compound toxicity and mode of toxicity.nnnMETHODSnGene annotations were built by text mining in Medline abstracts for retrieval of co-publications between genes, pathology terms, biological processes and pathways. This literature information was used to generate compound-specific keyword fingerprints, representing over-represented keywords calculated in a set of regulated genes after compound administration. To see whether keyword fingerprints can be used for assessment of compound toxicity, we analyzed microarray data sets of rat liver treated with 11 hepatotoxicants.nnnRESULTSnAnalysis of keyword fingerprints of two genotoxic carcinogens, two nongenotoxic carcinogens, two peroxisome proliferators and two randomly generated gene sets, showed that each compound produced a specific keyword fingerprint that correlated with the experimentally observed histopathological events induced by the individual compounds. By contrast, the random sets produced a flat aspecific keyword profile, indicating that the fingerprints induced by the compounds reflect biological events rather than random noise. A more detailed analysis of the keyword profiles of diethylhexylphthalate, dimethylnitrosamine and methapyrilene (MPy) showed that the differences in the keyword fingerprints of these three compounds are based upon known distinct modes of action. Visualization of MPy-linked keywords and MPy-induced genes in a literature network enabled us to construct a mode of toxicity proposal for MPy, which is in agreement with known effects of MPy in literature.nnnCONCLUSIONnCompound keyword fingerprinting based on information retrieved from literature is a powerful approach for compound profiling, allowing evaluation of compound toxicity and analysis of the mode of action.

Org 214007-0: A Novel Non-Steroidal Selective Glucocorticoid Receptor Modulator with Full Anti-Inflammatory Properties and Improved Therapeutic Index

Glucocorticoids (GCs) such as prednisolone are potent immunosuppressive drugs but suffer from severe adverse effects, including the induction of insulin resistance. Therefore, development of so-called Selective Glucocorticoid Receptor Modulators (SGRM) is highly desirable. Here we describe a non-steroidal Glucocorticoid Receptor (GR)-selective compound (Org 214007-0) with a binding affinity to GR similar to that of prednisolone. Structural modelling of the GR-Org 214007-0 binding site shows disturbance of the loop between helix 11 and helix 12 of GR, confirmed by partial recruitment of the TIF2-3 peptide. Using various cell lines and primary human cells, we show here that Org 214007-0 acts as a partial GC agonist, since it repressed inflammatory genes and was less effective in induction of metabolic genes. More importantly, in vivo studies in mice indicated that Org 214007-0 retained full efficacy in acute inflammation models as well as in a chronic collagen-induced arthritis (CIA) model. Gene expression profiling of muscle tissue derived from arthritic mice showed a partial activity of Org 214007-0 at an equi-efficacious dosage of prednisolone, with an increased ratio in repression versus induction of genes. Finally, in mice Org 214007-0 did not induce elevated fasting glucose nor the shift in glucose/glycogen balance in the liver seen with an equi-efficacious dose of prednisolone. All together, our data demonstrate that Org 214007-0 is a novel SGRMs with an improved therapeutic index compared to prednisolone. This class of SGRMs can contribute to effective anti-inflammatory therapy with a lower risk for metabolic side effects.

ss-TEA: Entropy based identification of receptor specific ligand binding residues from a multiple sequence alignment of class A GPCRs.

BackgroundG-protein coupled receptors (GPCRs) are involved in many different physiological processes and their function can be modulated by small molecules which bind in the transmembrane (TM) domain. Because of their structural and sequence conservation, the TM domains are often used in bioinformatics approaches to first create a multiple sequence alignment (MSA) and subsequently identify ligand binding positions. So far methods have been developed to predict the common ligand binding residue positions for class A GPCRs.ResultsHere we present 1) ss-TEA, a method to identify specific ligand binding residue positions for any receptor, predicated on high quality sequence information. 2) The largest MSA of class A non olfactory GPCRs in the public domain consisting of 13324 sequences covering most of the species homologues of the human set of GPCRs. A set of ligand binding residue positions extracted from literature of 10 different receptors shows that our method has the best ligand binding residue prediction for 9 of these 10 receptors compared to another state-of-the-art method.ConclusionsThe combination of the large multi species alignment and the newly introduced residue selection method ss-TEA can be used to rapidly identify subfamily specific ligand binding residues. This approach can aid the design of site directed mutagenesis experiments, explain receptor function and improve modelling. The method is also available online via GPCRDB at http://www.gpcr.org/7tm/.

Prednisolone-induced changes in gene-expression profiles in healthy volunteers.

BACKGROUNDnPrednisolone and other glucocorticoids (GCs) are potent anti-inflammatory and immunosuppressive drugs. However, prolonged use at a medium or high dose is hampered by side effects of which the metabolic side effects are most evident. Relatively little is known about their effect on gene-expression in vivo, the effect on cell subpopulations and the relation to the efficacy and side effects of GCs.nnnAIMnTo identify and compare prednisolone-induced gene signatures in CD4⁺ T lymphocytes and CD14⁺ monocytes derived from healthy volunteers and to link these signatures to underlying biological pathways involved in metabolic adverse effects.nnnMATERIALS & METHODSnWhole-genome expression profiling was performed on CD4⁺ T lymphocytes and CD14⁺ monocytes derived from healthy volunteers treated with prednisolone. Text-mining analyses was used to link genes to pathways involved in metabolic adverse events.nnnRESULTSnInduction of gene-expression was much stronger in CD4⁺ T lymphocytes than in CD14⁺ monocytes with respect to fold changes, but the number of truly cell-specific genes where a strong prednisolone effect in one cell type was accompanied by a total lack of prednisolone effect in the other cell type, was relatively low. Subsequently, a large set of genes was identified with a strong link to metabolic processes, for some of which the association with GCs is novel.nnnCONCLUSIONnThe identified gene signatures provide new starting points for further study into GC-induced transcriptional regulation in vivo and the mechanisms underlying GC-mediated metabolic side effects.