TranscriptomeBrowser: a powerful and flexible toolbox to explore productively the transcriptional landscape of the Gene Expression Omnibus database

SciGeneX: enhancing transcriptional analysis through gene module detection in single-cell and spatial transcriptomics data

The standard pipeline to analyze single-cell RNA-seq or spatial transcriptomics data focuses on a gene-centric approach that overlooks the collective behavior of genes. However, understanding cell populations necessitates recognizing intricate combinations of activated and repressed pathways. Therefore, a broader view of gene behavior offers more accurate insights into cellular heterogeneity in single-cell or spatial transcriptomics data. Here, we describe SciGeneX (Single-cell informative Gene eXplorer), a R package implementing a neighborhood analysis and a graph partitioning method to generate co-expression gene modules. These modules, whether shared or restricted to cell populations, collectively reflect cellular heterogeneity. Their combinations are able to highlight specific cell populations, even rare ones. SciGeneX uncovers rare and novel cell populations that were not observed before in human thymus spatial transcriptomics data. We show that SciGeneX outperforms existing methods on both artificial and experimental datasets. Overall, SciGeneX will aid in unravelling cellular and molecular diversity in single-cell and spatial transcriptomics studies.

Multifactorial Modes of Action of Arsenic Trioxide in Cancer Cells as Analyzed by Classical and Network Pharmacology

Arsenic trioxide is a traditional remedy in Chinese Medicine since ages. Nowadays, it is clinically used to treat acute promyelocytic leukemia (APL) by targeting PML/RARA. However, the drug's activity is broader and the mechanisms of action in other tumor types remain unclear. In this study, we investigated molecular modes of action by classical and network pharmacological approaches. CEM/ADR5000 resistance leukemic cells were similar sensitive to As2O3 as their wild-type counterpart CCRF-CEM (resistance ratio: 1.88). Drug-resistant U87.MG ΔEGFR glioblastoma cells harboring mutated epidermal growth factor receptor were even more sensitive (collateral sensitive) than wild-type U87.MG cells (resistance ratio: 0.33). HCT-116 colon carcinoma p53-/- knockout cells were 7.16-fold resistant toward As2O3 compared to wild-type cells. Forty genes determining cellular responsiveness to As2O3 were identified by microarray and COMPARE analyses in 58 cell lines of the NCI panel. Hierarchical cluster analysis-based heat mapping revealed significant differences between As2O3 sensitive cell lines and resistant cell lines with p-value: 1.86 × 10-5. The genes were subjected to Galaxy Cistrome gene promoter transcription factor analysis to predict the binding of transcription factors. We have exemplarily chosen NF-kB and AP-1, and indeed As2O3 dose-dependently inhibited the promoter activity of these two transcription factors in reporter cell lines. Furthermore, the genes identified here and those published in the literature were assembled and subjected to Ingenuity Pathway Analysis for comprehensive network pharmacological approaches that included all known factors of resistance of tumor cells to As2O3. In addition to pathways related to the anticancer effects of As2O3, several neurological pathways were identified. As arsenic is well-known to exert neurotoxicity, these pathways might account for neurological side effects. In conclusion, the activity of As2O3 is not restricted to acute promyelocytic leukemia. In addition to PML/RARA, numerous other genes belonging to diverse functional classes may also contribute to its cytotoxicity. Network pharmacology is suited to unravel the multifactorial modes of action of As2O3.

E2F transcription factors associated with up-regulated genes in glioblastoma

BACKGROUND

Glioblastoma is considered to the most common and malignant brain tumor in adults. Patients have a median survival of approximately one year from diagnosis due to poor response to therapy.

OBJECTIVE

We applied bioinformatics approaches to predict transcription factors (TF) that are deregulated in glioblastoma in an attempt to point out molecular targets for therapy.

METHODS

Up-regulated genes in glioblastoma selected from public microarray data were submitted to two TF association analyses. Thereafter, the expression levels of TF obtained in the overlap of analyses were assessed by RT-qPCR carried out in seven glioblastoma cell lines (T98, U251, U138, U87, U343, M059J, and M059K).

RESULTS

E2F1 and E2F4 were highlighted in both TF analyses. However, only E2F1 was confirmed as significantly up-regulated in all glioblastoma cell lines in vitro.

CONCLUSION

E2F1 is a potential common regulator of differentially expressed genes in glioblastoma, despite the genetic heterogeneity of tumor cells.

CX3CR1 is dysregulated in blood and brain from schizophrenia patients

The molecular mechanisms underlying schizophrenia remain largely unknown. Although schizophrenia is a mental disorder, there is increasing evidence to indicate that inflammatory processes driven by diverse environmental factors play a significant role in its development. With gene expression studies having been conducted across a variety of sample types, e.g., blood and postmortem brain, it is possible to investigate convergent signatures that may reveal interactions between the immune and nervous systems in schizophrenia pathophysiology. We conducted two meta-analyses of schizophrenia microarray gene expression data (N=474) and non-psychiatric control (N=485) data from postmortem brain and blood. Then, we assessed whether significantly dysregulated genes in schizophrenia could be shared between blood and brain. To validate our findings, we selected a top gene candidate and analyzed its expression by RT-qPCR in a cohort of schizophrenia subjects stabilized by atypical antipsychotic monotherapy (N=29) and matched controls (N=31). Meta-analyses highlighted inflammation as the major biological process associated with schizophrenia and that the chemokine receptor CX3CR1 was significantly down-regulated in schizophrenia. This differential expression was also confirmed in our validation cohort. Given both the recent data demonstrating selective CX3CR1 expression in subsets of neuroimmune cells, as well as behavioral and neuropathological observations of CX3CR1 deficiency in mouse models, our results of reduced CX3CR1 expression adds further support for a role played by monocyte/microglia in the neurodevelopment of schizophrenia.

Intrinsic disorder in PTEN and its interactome confers structural plasticity and functional versatility

IDPs, while structurally poor, are functionally rich by virtue of their flexibility and modularity. However, how mutations in IDPs elicit diseases, remain elusive. Herein, we have identified tumor suppressor PTEN as an intrinsically disordered protein (IDP) and elucidated the molecular principles by which its intrinsically disordered region (IDR) at the carboxyl-terminus (C-tail) executes its functions. Post-translational modifications, conserved eukaryotic linear motifs and molecular recognition features present in the C-tail IDR enhance PTEN's protein-protein interactions that are required for its myriad cellular functions. PTEN primary and secondary interactomes are also enriched in IDPs, most being cancer related, revealing that PTEN functions emanate from and are nucleated by the C-tail IDR, which form pliable network-hubs. Together, PTEN higher order functional networks operate via multiple IDP-IDP interactions facilitated by its C-tail IDR. Targeting PTEN IDR and its interaction hubs emerges as a new paradigm for treatment of PTEN related pathologies.

Integrative analysis of the transcriptome profiles observed in type 1, type 2 and gestational diabetes mellitus reveals the role of inflammation

BACKGROUND

Type 1 diabetes (T1D) is an autoimmune disease, while type 2 (T2D) and gestational diabetes (GDM) are considered metabolic disturbances. In a previous study evaluating the transcript profiling of peripheral mononuclear blood cells obtained from T1D, T2D and GDM patients we showed that the gene profile of T1D patients was closer to GDM than to T2D. To understand the influence of demographical, clinical, laboratory, pathogenetic and treatment features on the diabetes transcript profiling, we performed an analysis integrating these features with the gene expression profiles of the annotated genes included in databases containing information regarding GWAS and immune cell expression signatures.

METHODS

Samples from 56 (19 T1D, 20 T2D, and 17 GDM) patients were hybridized to whole genome one-color Agilent 4x44k microarrays. Non-informative genes were filtered by partitioning, and differentially expressed genes were obtained by rank product analysis. Functional analyses were carried out using the DAVID database, and module maps were constructed using the Genomica tool.

RESULTS

The functional analyses were able to discriminate between T1D and GDM patients based on genes involved in inflammation. Module maps of differentially expressed genes revealed that modulated genes: i) exhibited transcription profiles typical of macrophage and dendritic cells; ii) had been previously associated with diabetic complications by association and by meta-analysis studies, and iii) were influenced by disease duration, obesity, number of gestations, glucose serum levels and the use of medications, such as metformin.

CONCLUSION

This is the first module map study to show the influence of epidemiological, clinical, laboratory, immunopathogenic and treatment features on the transcription profiles of T1D, T2D and GDM patients.

The CO-Regulation Database (CORD): a tool to identify coordinately expressed genes

Background

Meta-analysis of gene expression array databases has the potential to reveal information about gene function. The identification of gene-gene interactions may be inferred from gene expression information but such meta-analysis is often limited to a single microarray platform. To address this limitation, we developed a gene-centered approach to analyze differential expression across thousands of gene expression experiments and created the CO-Regulation Database (CORD) to determine which genes are correlated with a queried gene.

Results

Using the GEO and ArrayExpress database, we analyzed over 120,000 group by group experiments from gene microarrays to determine the correlating genes for over 30,000 different genes or hypothesized genes. CORD output data is presented for sample queries with focus on genes with well-known interaction networks including p16 (CDKN2A), vimentin (VIM), MyoD (MYOD1). CDKN2A, VIM, and MYOD1 all displayed gene correlations consistent with known interacting genes.

Conclusions

We developed a facile, web-enabled program to determine gene-gene correlations across different gene expression microarray platforms. Using well-characterized genes, we illustrate how CORD's identification of co-expressed genes contributes to a better understanding a gene's potential function. The website is found at http://cord-db.org.

Navigating the multilayered organization of eukaryotic signaling: a new trend in data integration

The ever-increasing capacity of biological molecular data acquisition outpaces our ability to understand the meaningful relationships between molecules in a cell. Multiple databases were developed to store and organize these molecular data. However, emerging fundamental questions about concerted functions of these molecules in hierarchical cellular networks are poorly addressed. Here we review recent advances in the development of publically available databases that help us analyze the signal integration and processing by multilayered networks that specify biological responses in model organisms and human cells

Alterations of immune response of Non-Small Cell Lung Cancer with Azacytidine

Innovative therapies are needed for advanced Non-Small Cell Lung Cancer (NSCLC). We have undertaken a genomics based, hypothesis driving, approach to query an emerging potential that epigenetic therapy may sensitize to immune checkpoint therapy targeting PD-L1/PD-1 interaction. NSCLC cell lines were treated with the DNA hypomethylating agent azacytidine (AZA - Vidaza) and genes and pathways altered were mapped by genome-wide expression and DNA methylation analyses. AZA-induced pathways were analyzed in The Cancer Genome Atlas (TCGA) project by mapping the derived gene signatures in hundreds of lung adeno (LUAD) and squamous cell carcinoma (LUSC) samples. AZA up-regulates genes and pathways related to both innate and adaptive immunity and genes related to immune evasion in a several NSCLC lines. DNA hypermethylation and low expression of IRF7, an interferon transcription factor, tracks with this signature particularly in LUSC. In concert with these events, AZA up-regulates PD-L1 transcripts and protein, a key ligand-mediator of immune tolerance. Analysis of TCGA samples demonstrates that a significant proportion of primary NSCLC have low expression of AZA-induced immune genes, including PD-L1. We hypothesize that epigenetic therapy combined with blockade of immune checkpoints - in particular the PD-1/PD-L1 pathway - may augment response of NSCLC by shifting the balance between immune activation and immune inhibition, particularly in a subset of NSCLC with low expression of these pathways. Our studies define a biomarker strategy for response in a recently initiated trial to examine the potential of epigenetic therapy to sensitize patients with NSCLC to PD-1 immune checkpoint blockade.

Expansion of autoreactive unresponsive CD21-/low B cells in Sjögren's syndrome-associated lymphoproliferation

OBJECTIVE

Primary Sjögren's syndrome (SS) is an autoimmune disease associated with a high risk of developing non-Hodgkin's lymphoma. This study was undertaken to determine the nature of B cells driving lymphoproliferation in primary SS.

METHODS

B cell subsets and function were analyzed in peripheral blood from 66 adult patients with primary SS (including 14 patients with B cell lymphoproliferative disease [LPD]) and 30 healthy donors, using flow cytometry, calcium mobilization, and gene array analysis. The reactivity of recombinant antibodies isolated from single B cells from patients with primary SS and LPD was tested using an enzyme-linked immunosorbent assay.

RESULTS

We observed an expansion of an unusual CD21-/low B cell population that correlated with lymphoproliferation in patients with primary SS. A majority of CD21-/low B cells from patients with primary SS expressed autoreactive antibodies, which recognized nuclear and cytoplasmic structures. These B cells belonged to the memory compartment, since their Ig genes were mutated. They were unable to induce calcium flux, become activated, or proliferate in response to B cell receptor and/or CD40 triggering, suggesting that these autoreactive B cells may be anergic. However, CD21-/low B cells from patients with primary SS remained responsive to Toll-like receptor (TLR) stimulation. Molecules specifically expressed in CD21-/low B cells that are likely to induce their unresponsive stage were detected in gene array analyses.

CONCLUSION

Patients with primary SS who display high frequencies of autoreactive and unresponsive CD21-/low B cells are susceptible to developing lymphoproliferation. These cells remain in peripheral blood controlled by functional anergy instead of being eliminated, and chronic antigenic stimulation through TLR stimulation may create a favorable environment for breaking tolerance and activating these cells.

TranscriptomeBrowser 3.0: introducing a new compendium of molecular interactions and a new visualization tool for the study of gene regulatory networks

Background

Deciphering gene regulatory networks by in silico approaches is a crucial step in the study of the molecular perturbations that occur in diseases. The development of regulatory maps is a tedious process requiring the comprehensive integration of various evidences scattered over biological databases. Thus, the research community would greatly benefit from having a unified database storing known and predicted molecular interactions. Furthermore, given the intrinsic complexity of the data, the development of new tools offering integrated and meaningful visualizations of molecular interactions is necessary to help users drawing new hypotheses without being overwhelmed by the density of the subsequent graph.

Results

We extend the previously developed TranscriptomeBrowser database with a set of tables containing 1,594,978 human and mouse molecular interactions. The database includes: (i) predicted regulatory interactions (computed by scanning vertebrate alignments with a set of 1,213 position weight matrices), (ii) potential regulatory interactions inferred from systematic analysis of ChIP-seq experiments, (iii) regulatory interactions curated from the literature, (iv) predicted post-transcriptional regulation by micro-RNA, (v) protein kinase-substrate interactions and (vi) physical protein-protein interactions. In order to easily retrieve and efficiently analyze these interactions, we developed In-teractomeBrowser, a graph-based knowledge browser that comes as a plug-in for Transcriptome-Browser. The first objective of InteractomeBrowser is to provide a user-friendly tool to get new insight into any gene list by providing a context-specific display of putative regulatory and physical interactions. To achieve this, InteractomeBrowser relies on a "cell compartments-based layout" that makes use of a subset of the Gene Ontology to map gene products onto relevant cell compartments. This layout is particularly powerful for visual integration of heterogeneous biological information and is a productive avenue in generating new hypotheses. The second objective of InteractomeBrowser is to fill the gap between interaction databases and dynamic modeling. It is thus compatible with the network analysis software Cytoscape and with the Gene Interaction Network simulation software (GINsim). We provide examples underlying the benefits of this visualization tool for large gene set analysis related to thymocyte differentiation.

Conclusions

The InteractomeBrowser plugin is a powerful tool to get quick access to a knowledge database that includes both predicted and validated molecular interactions. InteractomeBrowser is available through the TranscriptomeBrowser framework and can be found at: http://tagc.univ-mrs.fr/tbrowser/. Our database is updated on a regular basis.

Comprehensive literature review and statistical considerations for microarray meta-analysis

With the rapid advances of various high-throughput technologies, generation of ‘-omics’ data is commonplace in almost every biomedical field. Effective data management and analytical approaches are essential to fully decipher the biological knowledge contained in the tremendous amount of experimental data. Meta-analysis, a set of statistical tools for combining multiple studies of a related hypothesis, has become popular in genomic research. Here, we perform a systematic search from PubMed and manual collection to obtain 620 genomic meta-analysis papers, of which 333 microarray meta-analysis papers are summarized as the basis of this paper and the other 249 GWAS meta-analysis papers are discussed in the next companion paper. The review in the present paper focuses on various biological purposes of microarray meta-analysis, databases and software and related statistical procedures. Statistical considerations of such an analysis are further scrutinized and illustrated by a case study. Finally, several open questions are listed and discussed.

Using MCL to extract clusters from networks

MCL is a general purpose cluster algorithm for both weighted and unweighted networks. The algorithm utilises network topology as well as edge weights, is highly scalable and has been applied in a wide variety of bioinformatic methods. In this chapter, we give protocols and case studies for clustering of networks derived from, respectively, protein sequence similarities and gene expression profile correlations.

Immune response and mitochondrial metabolism are commonly deregulated in DMD and aging skeletal muscle

Duchenne Muscular Dystrophy (DMD) is a complex process involving multiple pathways downstream of the primary genetic insult leading to fatal muscle degeneration. Aging muscle is a multifactorial neuromuscular process characterized by impaired muscle regeneration leading to progressive atrophy. We hypothesized that these chronic atrophying situations may share specific myogenic adaptative responses at transcriptional level according to tissue remodeling. Muscle biopsies from four young DMD and four AGED subjects were referred to a group of seven muscle biopsies from young subjects without any neuromuscular disorder and explored through a dedicated expression microarray. We identified 528 differentially expressed genes (out of 2,745 analyzed), of which 328 could be validated by an exhaustive meta-analysis of public microarray datasets referring to DMD and Aging in skeletal muscle. Among the 328 validated co-expressed genes, 50% had the same expression profile in both groups and corresponded to immune/fibrosis responses and mitochondrial metabolism. Generalizing these observed meta-signatures with large compendia of public datasets reinforced our results as they could be also identified in other pathological processes and in diverse physiological conditions. Focusing on the common gene signatures in these two atrophying conditions, we observed enrichment in motifs for candidate transcription factors that may coordinate either the immune/fibrosis responses (ETS1, IRF1, NF1) or the mitochondrial metabolism (ESRRA). Deregulation in their expression could be responsible, at least in part, for the same transcriptome changes initiating the chronic muscle atrophy. This study suggests that distinct pathophysiological processes may share common gene responses and pathways related to specific transcription factors.

NCBI GEO: archive for functional genomics data sets--10 years on

A decade ago, the Gene Expression Omnibus (GEO) database was established at the National Center for Biotechnology Information (NCBI). The original objective of GEO was to serve as a public repository for high-throughput gene expression data generated mostly by microarray technology. However, the research community quickly applied microarrays to non-gene-expression studies, including examination of genome copy number variation and genome-wide profiling of DNA-binding proteins. Because the GEO database was designed with a flexible structure, it was possible to quickly adapt the repository to store these data types. More recently, as the microarray community switches to next-generation sequencing technologies, GEO has again adapted to host these data sets. Today, GEO stores over 20,000 microarray- and sequence-based functional genomics studies, and continues to handle the majority of direct high-throughput data submissions from the research community. Multiple mechanisms are provided to help users effectively search, browse, download and visualize the data at the level of individual genes or entire studies. This paper describes recent database enhancements, including new search and data representation tools, as well as a brief review of how the community uses GEO data. GEO is freely accessible at http://www.ncbi.nlm.nih.gov/geo/.

Sex-related differences in gene expression following Coxiella burnetii infection in mice: potential role of circadian rhythm

Background

Q fever, a zoonosis due to Coxiella burnetii infection, exhibits sexual dimorphism; men are affected more frequently and severely than women for a given exposure. Here we explore whether the severity of C. burnetii infection in mice is related to differences in male and female gene expression profiles.

Methodology/Principal Findings

Mice were infected with C. burnetii for 24 hours, and gene expression was measured in liver cells using microarrays. Multiclass analysis identified 2,777 probes for which expression was specifically modulated by C. burnetti infection. Only 14% of the modulated genes were sex-independent, and the remaining 86% were differentially expressed in males and females. Castration of males and females showed that sex hormones were responsible for more than 60% of the observed gene modulation, and this reduction was most pronounced in males. Using functional annotation of modulated genes, we identified four clusters enriched in males that were related to cell-cell adhesion, signal transduction, defensins and cytokine/Jak-Stat pathways. Up-regulation of the IL-10 and Stat-3 genes may account for the high susceptibility of men with Q fever to C. burnetii infection and autoantibody production. Two clusters were identified in females, including the circadian rhythm pathway, which consists of positive (Clock, Arntl) and negative (Per) limbs of a feedback loop. We found that Clock and Arntl were down-modulated whereas Per was up-regulated; these changes may be associated with efficient bacterial elimination in females but not in males, in which an exacerbated host response would be prominent.

Conclusion

This large-scale study revealed for the first time that circadian rhythm plays a major role in the anti-infectious response of mice, and it provides a new basis for elucidating the role of sexual dimorphism in human infections.

A novel tissue-specific meta-analysis approach for gene expression predictions, initiated with a mammalian gene expression testis database

BACKGROUND

In the recent years, there has been a rise in gene expression profiling reports. Unfortunately, it has not been possible to make maximum use of available gene expression data. Many databases and programs can be used to derive the possible expression patterns of mammalian genes, based on existing data. However, these available resources have limitations. For example, it is not possible to obtain a list of genes that are expressed in certain conditions. To overcome such limitations, we have taken up a new strategy to predict gene expression patterns using available information, for one tissue at a time.

RESULTS

The first step of this approach involved manual collection of maximum data derived from large-scale (genome-wide) gene expression studies, pertaining to mammalian testis. These data have been compiled into a Mammalian Gene Expression Testis-database (MGEx-Tdb). This process resulted in a richer collection of gene expression data compared to other databases/resources, for multiple testicular conditions. The gene-lists collected this way in turn were exploited to derive a 'consensus' expression status for each gene, across studies. The expression information obtained from the newly developed database mostly agreed with results from multiple small-scale studies on selected genes. A comparative analysis showed that MGEx-Tdb can retrieve the gene expression information more efficiently than other commonly used databases. It has the ability to provide a clear expression status (transcribed or dormant) for most genes, in the testis tissue, under several specific physiological/experimental conditions and/or cell-types.

CONCLUSIONS

Manual compilation of gene expression data, which can be a painstaking process, followed by a consensus expression status determination for specific locations and conditions, can be a reliable way of making use of the existing data to predict gene expression patterns. MGEx-Tdb provides expression information for 14 different combinations of specific locations and conditions in humans (25,158 genes), 79 in mice (22,919 genes) and 23 in rats (14,108 genes). It is also the first system that can predict expression of genes with a 'reliability-score', which is calculated based on the extent of agreements and contradictions across gene-sets/studies. This new platform is publicly available at the following web address: http://resource.ibab.ac.in/MGEx-Tdb/.

Protein Bioinformatics Infrastructure for the Integration and Analysis of Multiple High-Throughput "omics" Data

High-throughput “omics” technologies bring new opportunities for biological and biomedical researchers to ask complex questions and gain new scientific insights. However, the voluminous, complex, and context-dependent data being maintained in heterogeneous and distributed environments plus the lack of well-defined data standard and standardized nomenclature imposes a major challenge which requires advanced computational methods and bioinformatics infrastructures for integration, mining, visualization, and comparative analysis to facilitate data-driven hypothesis generation and biological knowledge discovery. In this paper, we present the challenges in high-throughput “omics” data integration and analysis, introduce a protein-centric approach for systems integration of large and heterogeneous high-throughput “omics” data including microarray, mass spectrometry, protein sequence, protein structure, and protein interaction data, and use scientific case study to illustrate how one can use varied “omics” data from different laboratories to make useful connections that could lead to new biological knowledge.

Genomics Portals: integrative web-platform for mining genomics data

Background

A large amount of experimental data generated by modern high-throughput technologies is available through various public repositories. Our knowledge about molecular interaction networks, functional biological pathways and transcriptional regulatory modules is rapidly expanding, and is being organized in lists of functionally related genes. Jointly, these two sources of information hold a tremendous potential for gaining new insights into functioning of living systems.

Results

Genomics Portals platform integrates access to an extensive knowledge base and a large database of human, mouse, and rat genomics data with basic analytical visualization tools. It provides the context for analyzing and interpreting new experimental data and the tool for effective mining of a large number of publicly available genomics datasets stored in the back-end databases. The uniqueness of this platform lies in the volume and the diversity of genomics data that can be accessed and analyzed (gene expression, ChIP-chip, ChIP-seq, epigenomics, computationally predicted binding sites, etc), and the integration with an extensive knowledge base that can be used in such analysis.

Conclusion

The integrated access to primary genomics data, functional knowledge and analytical tools makes Genomics Portals platform a unique tool for interpreting results of new genomics experiments and for mining the vast amount of data stored in the Genomics Portals backend databases. Genomics Portals can be accessed and used freely at http://GenomicsPortals.org.

SimCT: a generic tool to visualize ontology-based relationships for biological objects

We present a web-based service, SimCT, which allows to graphically display the relationships between biological objects (e.g. genes or proteins) based on their annotations to a biomedical ontology. The result is presented as a tree of these objects, which can be viewed and explored through a specific java applet designed to highlight relevant features. Unlike the numerous tools that search for overrepresented terms, SimCT draws a simplified representation of biological terms present in the set of objects, and can be applied to any ontology for which annotation data is available. Being web-based, it does not require prior installation, and provides an intuitive, easy-to-use service.

AVAILABILITY

http://tagc.univ-mrs.fr/SimCT.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.