Judging the quality of gene expression-based clustering methods using gene annotation

Insight into the role of stress response and toxic mechanism induced by Chloro-haloacetonitrile in vitro

Chloro-haloacetonitrile (Cl-HAN), belongs to a group of nitrogenous disinfection by-products (N-DBPs) found in surface water, and are known to pose a major risk to the safety of human drinking water. However, the exact biological toxicity mechanism and the extent of the stress response caused by Cl-HAN remain unclear, resulting in a lack of effective measures to control its presence. Thus, the quantitative toxicological genomics and bioinformatics methods were applied to explore the effects of three chloro-haloacetonitriles (Cl-HANs) on the transcription of fusion genes under varying concentrations of stress in E. coli over 2-hour period. The initial stress response and their toxic mechanism were analyzed. The study also identified the molecular toxicity endpoint, and the core genes that are responsible for the specific toxicity of different Cl-HANs. Cl-HANs exhibited concentration-dependent characteristics of toxic effects, and caused changes in gene expression related oxidative and membrane stress. The stress response results showed that dichloroacetonitrile (dCAN) still caused significant DNA damage under the lowest concentration stress. Chloroacetonitrile (CAN) and trichloroacetonitrile (tCAN) exhibited lower genetic toxicity levels at 513 μg/L and 10.7 μg/L, respectively. The toxic effects of tCAN were widespread. And there was a good correlation between the molecular endpoint (EC-TELI1.5) and the phenotypic endpoint (LD50) with rp=-0.8634 (P=0.0593). In all concentrations of stress in CAN, dCAN, and tCAN, the number of overexpressed genes shared was 15, 2, and 14, respectively. Furthermore, bioinformatics analysis demonstrated that Cl-HANs affected genes associated with general stress pathways, such as cell biochemistry and physical homeostasis, resulting in changes in biological processes. And for CAN-induced DNA damage, polA played a dominant role, while katG, oxyR, and ahpC were the core genes involved in oxidative stress induced by dCAN and tCAN, respectively. These findings provide valuable data for the toxic effect of Cl-HANs.

Type I Diabetes Mellitus Suppresses Experimental Skin Carcinogenesis

This study explores the previously uncharted territory of the effects of ultraviolet (UV) radiation on diabetic skin, compared to its well-documented impact on normal skin, particularly focusing on carcinogenesis and aging. Employing hairless SKH-hr2, Type 1 and 2 diabetic, and nondiabetic male mice, the research subjected these to UV radiation thrice weekly for eight months. The investigation included comprehensive assessments of photoaging and photocarcinogenesis in diabetic versus normal skin, measuring factors such as hydration, trans-epidermal water loss, elasticity, skin thickness, melanin, sebum content, stratum corneum exfoliation and body weight, alongside photo documentation. Additionally, oxidative stress and the presence of hydrophilic antioxidants (uric acid and glutathione) in the stratum corneum were evaluated. Histopathological examination post-sacrifice provided insights into the morphological changes. Findings reveal that under UV exposure, Type 1 diabetic skin showed heightened dehydration, thinning, and signs of accelerated aging. Remarkably, Type 1 diabetic mice did not develop squamous cell carcinoma or pigmented nevi, contrary to normal and Type 2 diabetic skin. This unexpected resistance to UV-induced skin cancers in Type 1 diabetic skin prompts a crucial need for further research to uncover the underlying mechanisms providing this resistance.

In vitro toxicity assessment of haloacetamides via a toxicogenomics assay

It is important to study the stress effects and mechanisms of haloacetamide (HAcAm) disinfection byproducts to reveal their health hazards. In this context, toxicological g was applied to evaluate the effects of four HAcAms, revealing the status of gene expression on Escherichia coli in different stress response types (oxidative, protein, membrane, general, DNA). This study revealed that the main toxic action modes of these HAcAms were general and membrane stresses by high-resolution, real-time gene expression profiling combined with clustering analysis. The results of time-gene evaluation showed that the presence of chloroacetamide (CAcAm) and bromoacetamide (BAcAm) generated more reactive oxygen species, thus activating oxidative stress. Trichloroacetamide (tCAcAm) induced altered expression of glutathione marker genes and membrane stress-related genes, and iodoacetamide (IAcAm) caused severe DNA damage by damaging DNA strands and individual nucleotides mainly through damage to nucleic acids and bases. Furthermore, quantitative structure-activity relationship (QSAR) modelling results indicated that the biological activities of HAcAms were related to their quantum chemical and topological properties.

A Novel Calibration Step in Gene Co-Expression Network Construction

High-throughput technologies such as DNA microarrays and RNA-sequencing are used to measure the expression levels of large numbers of genes simultaneously. To support the extraction of biological knowledge, individual gene expression levels are transformed to Gene Co-expression Networks (GCNs). In a GCN, nodes correspond to genes, and the weight of the connection between two nodes is a measure of similarity in the expression behavior of the two genes. In general, GCN construction and analysis includes three steps; 1) calculating a similarity value for each pair of genes 2) using these similarity values to construct a fully connected weighted network 3) finding clusters of genes in the network, commonly called modules. The specific implementation of these three steps can significantly impact the final output and the downstream biological analysis. GCN construction is a well-studied topic. Existing algorithms rely on relatively simple statistical and mathematical tools to implement these steps. Currently, software package WGCNA appears to be the most widely accepted standard. We hypothesize that the raw features provided by sequencing data can be leveraged to extract modules of higher quality. A novel preprocessing step of the gene expression data set is introduced that in effect calibrates the expression levels of individual genes, before computing pairwise similarities. Further, the similarity is computed as an inner-product of positive vectors. In experiments, this provides a significant improvement over WGCNA, as measured by aggregate p-values of the gene ontology term enrichment of the computed modules.

Differential and Common Signatures of miRNA Expression and Methylation in Childhood Central Nervous System Malignancies: An Experimental and Computational Approach

Epigenetic modifications are considered of utmost significance for tumor ontogenesis and progression. Especially, it has been found that miRNA expression, as well as DNA methylation plays a significant role in central nervous system tumors during childhood. A total of 49 resected brain tumors from children were used for further analysis. DNA methylation was identified with methylation-specific MLPA and, in particular, for the tumor suppressor genes CASP8, RASSF1, MGMT, MSH6, GATA5, ATM1, TP53, and CADM1. miRNAs were identified with microarray screening, as well as selected samples, were tested for their mRNA expression levels. CASP8, RASSF1 were the most frequently methylated genes in all tumor samples. Simultaneous methylation of genes manifested significant results with respect to tumor staging, tumor type, and the differentiation of tumor and control samples. There was no significant dependence observed with the methylation of one gene promoter, rather with the simultaneous presence of all detected methylated genes' promoters. miRNA expression was found to be correlated to gene methylation. Epigenetic regulation appears to be of major importance in tumor progression and pathophysiology, making it an imperative field of study.

Bioinformatic Analysis of Temporal and Spatial Proteome Alternations During Infections

Microbial pathogens have evolved numerous mechanisms to hijack host's systems, thus causing disease. This is mediated by alterations in the combined host-pathogen proteome in time and space. Mass spectrometry-based proteomics approaches have been developed and tailored to map disease progression. The result is complex multidimensional data that pose numerous analytic challenges for downstream interpretation. However, a systematic review of approaches for the downstream analysis of such data has been lacking in the field. In this review, we detail the steps of a typical temporal and spatial analysis, including data pre-processing steps (i.e., quality control, data normalization, the imputation of missing values, and dimensionality reduction), different statistical and machine learning approaches, validation, interpretation, and the extraction of biological information from mass spectrometry data. We also discuss current best practices for these steps based on a collection of independent studies to guide users in selecting the most suitable strategies for their dataset and analysis objectives. Moreover, we also compiled the list of commonly used R software packages for each step of the analysis. These could be easily integrated into one's analysis pipeline. Furthermore, we guide readers through various analysis steps by applying these workflows to mock and host-pathogen interaction data from public datasets. The workflows presented in this review will serve as an introduction for data analysis novices, while also helping established users update their data analysis pipelines. We conclude the review by discussing future directions and developments in temporal and spatial proteomics and data analysis approaches. Data analysis codes, prepared for this review are available from https://github.com/BabuLab-UofR/TempSpac, where guidelines and sample datasets are also offered for testing purposes.

Co-Deregulated miRNA Signatures in Childhood Central Nervous System Tumors: In Search for Common Tumor miRNA-Related Mechanics

Simple Summary

Childhood tumors of the central nervous system (CNS) constitute a grave disease and their diagnosis is difficult to be handled. To gain better knowledge of the tumor’s biology, it is essential to understand the underlying mechanisms of the disease. MicroRNAs (miRNAs) are small noncoding RNAs that are dysregulated in many types of CNS tumors and regulate their occurrence and development through specific signal pathways. However, different types of CNS tumors’ area are characterized by different deregulated miRNAs. Here, we hypothesized that CNS tumors could have commonly deregulated miRNAs, i.e., miRNAs that are simultaneously either upregulated or downregulated in all tumor types compared to the normal brain tissue, irrespectively of the tumor sub-type and/or diagnosis. The only criterion is that they are present in brain tumors. This approach could lead us to the discovery of miRNAs that could be used as pan-CNS tumoral therapeutic targets, if successful.

Abstract

Despite extensive experimentation on pediatric tumors of the central nervous system (CNS), related to both prognosis, diagnosis and treatment, the understanding of pathogenesis and etiology of the disease remains scarce. MicroRNAs are known to be involved in CNS tumor oncogenesis. We hypothesized that CNS tumors possess commonly deregulated miRNAs across different CNS tumor types. Aim: The current study aims to reveal the co-deregulated miRNAs across different types of pediatric CNS tumors. Materials: A total of 439 CNS tumor samples were collected from both in-house microarray experiments as well as data available in public databases. Diagnoses included medulloblastoma, astrocytoma, ependydoma, cortical dysplasia, glioblastoma, ATRT, germinoma, teratoma, yoc sac tumors, ocular tumors and retinoblastoma. Results: We found miRNAs that were globally up- or down-regulated in the majority of the CNS tumor samples. MiR-376B and miR-372 were co-upregulated, whereas miR-149, miR-214, miR-574, miR-595 and miR-765 among others, were co-downregulated across all CNS tumors. Receiver-operator curve analysis showed that miR-149, miR-214, miR-574, miR-595 and miR765 could distinguish between CNS tumors and normal brain tissue. Conclusions: Our approach could prove significant in the search for global miRNA targets for tumor diagnosis and therapy. To the best of our knowledge, there are no previous reports concerning the present approach.

A genome-wide portrait of pervasive drug contaminants

Using a validated yeast chemogenomic platform, we characterized the genome-wide effects of several pharmaceutical contaminants, including three N-nitrosamines (NDMA, NDEA and NMBA), two related compounds (DMF and 4NQO) and several of their metabolites. A collection of 4800 non-essential homozygous diploid yeast deletion strains were screened in parallel and the strain abundance was quantified by barcode sequencing. These data were used to rank deletion strains representing genes required for resistance to the compounds to delineate affected cellular pathways and to visualize the global cellular effects of these toxins in an easy-to-use searchable database. Our analysis of the N-nitrosamine screens uncovered genes (via their corresponding homozygous deletion mutants) involved in several evolutionarily conserved pathways, including: arginine biosynthesis, mitochondrial genome integrity, vacuolar protein sorting and DNA damage repair. To investigate why NDMA, NDEA and DMF caused fitness defects in strains lacking genes of the arginine pathway, we tested several N-nitrosamine metabolites (methylamine, ethylamine and formamide), and found they also affected arginine pathway mutants. Notably, each of these metabolites has the potential to produce ammonium ions during their biotransformation. We directly tested the role of ammonium ions in N-nitrosamine toxicity by treatment with ammonium sulfate and we found that ammonium sulfate also caused a growth defect in arginine pathway deletion strains. Formaldehyde, a metabolite produced from NDMA, methylamine and formamide, and which is known to cross-link free amines, perturbed deletion strains involved in chromatin remodeling and DNA repair pathways. Finally, co-administration of N-nitrosamines with ascorbic or ferulic acid did not relieve N-nitrosamine toxicity. In conclusion, we used parallel deletion mutant analysis to characterize the genes and pathways affected by exposure to N-nitrosamines and related compounds, and provide the data in an accessible, queryable database.

Differential Expression and PAH Degradation: What Burkholderia vietnamiensis G4 Can Tell Us?

Petroleum is the major energy matrix in the world whose refining generates chemical byproducts that may damage the environment. Among such waste, polycyclic aromatic hydrocarbons (PAH) are considered persistent pollutants. Sixteen of these are considered priority for remediation, and among them is benzo(a)pyrene. Amid remediation techniques, bioremediation stands out. The genus Burkholderia is amongst the microorganisms known for being capable of degrading persistent compounds; its strains are used as models to study such ability. High-throughput sequencing allows researchers to reach a wider knowledge about biodegradation by bacteria. Using transcripts and mRNA analysis, the genomic regions involved in this aptitude can be detected. To unravel these processes, we used the model B. vietnamiensis strain G4 in two experimental groups: one was exposed to benzo(a)pyrene and the other one (control) was not. Six transcriptomes were generated from each group aiming to compare gene expression and infer which genes are involved in degradation pathways. One hundred fifty-six genes were differentially expressed in the benzo(a)pyrene exposed group, from which 33% are involved in catalytic activity. Among these, the most significant genomic regions were phenylacetic acid degradation protein paaN, involved in the degradation of organic compounds to obtain energy; oxidoreductase FAD-binding subunit, related to the regulation of electrons within groups of dioxygenase enzymes with potential to cleave benzene rings; and dehydrogenase, described as accountable for phenol degradation. These data provide the basis for understanding the bioremediation of benzo(a)pyrene and the possible applications of this strain in polluted environments.

Factorized embeddings learns rich and biologically meaningful embedding spaces using factorized tensor decomposition

MOTIVATION

The recent development of sequencing technologies revolutionized our understanding of the inner workings of the cell as well as the way disease is treated. A single RNA sequencing (RNA-Seq) experiment, however, measures tens of thousands of parameters simultaneously. While the results are information rich, data analysis provides a challenge. Dimensionality reduction methods help with this task by extracting patterns from the data by compressing it into compact vector representations.

RESULTS

We present the factorized embeddings (FE) model, a self-supervised deep learning algorithm that learns simultaneously, by tensor factorization, gene and sample representation spaces. We ran the model on RNA-Seq data from two large-scale cohorts and observed that the sample representation captures information on single gene and global gene expression patterns. Moreover, we found that the gene representation space was organized such that tissue-specific genes, highly correlated genes as well as genes participating in the same GO terms were grouped. Finally, we compared the vector representation of samples learned by the FE model to other similar models on 49 regression tasks. We report that the representations trained with FE rank first or second in all of the tasks, surpassing, sometimes by a considerable margin, other representations.

AVAILABILITY AND IMPLEMENTATION

A toy example in the form of a Jupyter Notebook as well as the code and trained embeddings for this project can be found at: https://github.com/TrofimovAssya/FactorizedEmbeddings.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Lag penalized weighted correlation for time series clustering

BACKGROUND

The similarity or distance measure used for clustering can generate intuitive and interpretable clusters when it is tailored to the unique characteristics of the data. In time series datasets generated with high-throughput biological assays, measurements such as gene expression levels or protein phosphorylation intensities are collected sequentially over time, and the similarity score should capture this special temporal structure.

RESULTS

We propose a clustering similarity measure called Lag Penalized Weighted Correlation (LPWC) to group pairs of time series that exhibit closely-related behaviors over time, even if the timing is not perfectly synchronized. LPWC aligns time series profiles to identify common temporal patterns. It down-weights aligned profiles based on the length of the temporal lags that are introduced. We demonstrate the advantages of LPWC versus existing time series and general clustering algorithms. In a simulated dataset based on the biologically-motivated impulse model, LPWC is the only method to recover the true clusters for almost all simulated genes. LPWC also identifies clusters with distinct temporal patterns in our yeast osmotic stress response and axolotl limb regeneration case studies.

CONCLUSIONS

LPWC achieves both of its time series clustering goals. It groups time series with correlated changes over time, even if those patterns occur earlier or later in some of the time series. In addition, it refrains from introducing large shifts in time when searching for temporal patterns by applying a lag penalty. The LPWC R package is available at https://github.com/gitter-lab/LPWC and CRAN under a MIT license.

Multiple organ dysfunction after trauma

Background

The nature of multiple organ dysfunction syndrome (MODS) after traumatic injury is evolving as resuscitation practices advance and more patients survive their injuries to reach critical care. The aim of this study was to characterize contemporary MODS subtypes in trauma critical care at a population level.

Methods

Adult patients admitted to major trauma centre critical care units were enrolled in this 4‐week point‐prevalence study. MODS was defined by a daily total Sequential Organ Failure Assessment (SOFA) score of more than 5. Hierarchical clustering of SOFA scores over time was used to identify MODS subtypes.

Results

Some 440 patients were enrolled, of whom 245 (55·7 per cent) developed MODS. MODS carried a high mortality rate (22·0 per cent versus 0·5 per cent in those without MODS; P < 0·001) and 24·0 per cent of deaths occurred within the first 48 h after injury. Three patterns of MODS were identified, all present on admission. Cluster 1 MODS resolved early with a median time to recovery of 4 days and a mortality rate of 14·4 per cent. Cluster 2 had a delayed recovery (median 13 days) and a mortality rate of 35 per cent. Cluster 3 had a prolonged recovery (median 25 days) and high associated mortality rate of 46 per cent. Multivariable analysis revealed distinct clinical associations for each form of MODS; 24‐hour crystalloid administration was associated strongly with cluster 1 (P = 0·009), traumatic brain injury with cluster 2 (P = 0·002) and admission shock severity with cluster 3 (P = 0·003).

Conclusion

Contemporary MODS has at least three distinct types based on patterns of severity and recovery. Further characterization of MODS subtypes and their underlying pathophysiology may lead to future opportunities for early stratification and targeted interventions.

Measuring similarity between gene interaction profiles

Background

Gene and protein interaction data are often represented as interaction networks, where nodes stand for genes or gene products and each edge stands for a relationship between a pair of gene nodes. Commonly, that relationship within a pair is specified by high similarity between profiles (vectors) of experimentally defined interactions of each of the two genes with all other genes in the genome; only gene pairs that interact with similar sets of genes are linked by an edge in the network. The tight groups of genes/gene products that work together in a cell can be discovered by the analysis of those complex networks.

Results

We show that the choice of the similarity measure between pairs of gene vectors impacts the properties of networks and of gene modules detected within them. We re-analyzed well-studied data on yeast genetic interactions, constructed four genetic networks using four different similarity measures, and detected gene modules in each network using the same algorithm. The four networks induced different numbers of putative functional gene modules, and each similarity measure induced some unique modules. In an example of a putative functional connection suggested by comparing genetic interaction vectors, we predict a link between SUN-domain proteins and protein glycosylation in the endoplasmic reticulum.

Conclusions

The discovery of molecular modules in genetic networks is sensitive to the way of measuring similarity between profiles of gene interactions in a cell. In the absence of a formal way to choose the “best” measure, it is advisable to explore the measures with different mathematical properties, which may identify different sets of connections between genes.

Electronic supplementary material

The online version of this article (10.1186/s12859-019-3024-x) contains supplementary material, which is available to authorized users.

Computational methods to identify bimodal gene expression and facilitate personalized treatment in cancer patients

Standard methods for detecting cancer-associated genes rely on comparison of sample means between cancer patients and healthy controls. While such methods have successfully identified several oncogenes and tumor suppressor genes, they neglect to account for heterogeneity within the cancer population. Genetic mutations, translocations, and amplifications are often inconsistent across tumors, and instead they often affect smaller subsets of patients. This concept gives rise to the idea of bimodally expressed genes, or genes that display two modes of expression within one population. Analysis of bimodal gene expression has been explored via a variety of techniques including test statistics and clustering. In this review, we summarize the methodologies used to quantify bimodal gene expression and address the utility of these genes in patient stratification and specialized therapeutics in breast and lung cancer. Finally we discuss the limitations and future directions for bimodal genes in the era of high-throughput sequencing and personalized medicine.

Visual Pattern-Driven Exploration of Big Data

Pattern extraction algorithms are enabling insights into the ever-growing amount of today's datasets by translating reoccurring data properties into compact representations. Yet, a practical problem arises: With increasing data volumes and complexity also the number of patterns increases, leaving the analyst with a vast result space. Current algorithmic and especially visualization approaches often fail to answer central overview questions essential for a comprehensive understanding of pattern distributions and support, their quality, and relevance to the analysis task. To address these challenges, we contribute a visual analytics pipeline targeted on the pattern-driven exploration of result spaces in a semi-automatic fashion. Specifically, we combine image feature analysis and unsupervised learning to partition the pattern space into interpretable, coherent chunks, which should be given priority in a subsequent in-depth analysis. In our analysis scenarios, no ground-truth is given. Thus, we employ and evaluate novel quality metrics derived from the distance distributions of our image feature vectors and the derived cluster model to guide the feature selection process. We visualize our results interactively, allowing the user to drill down from overview to detail into the pattern space and demonstrate our techniques in two case studies on Earth observation and biomedical genomic data.

Identification of gene-specific DNA methylation signature for Colorectal Cancer

BACKGROUND

Colorectal Cancer (CC), a common disease causing approximately million deaths annually, has been the third most frequent type of malignancy. We aimed to identify gene-specific DNA methylation signature to function as prognostic and predictive markers for CC patient survival.

METHODS

Expression profiles of gene-specific DNA methylation and the corresponding clinical information of 201 CC patients were downloaded from The Cancer Genome Atlas (TCGA) dataset and differentially expressed gene-specific DNA methylation was identified after tumor subtype classification. A risk score model was further built by analyzing the expression data of these gene-specific DNA methylations from the training dataset of CC patients.

RESULTS

Totally, 214 gene-specific DNA methylations were found to be expressed significantly between different subtypes of CC, including 150 up-regulated and 64 down-regulated ones. Up-regulated gene-specific DNA methylation accounted for 70.1% and the down-regulated gene-specific DNA methylation accounted for 29.9%. Hereinto, six gene-specific DNA methylations were obtained, including methy_vimentin and methy_ TFPI2, which were found significantly correlated with overall survival status of patients with CC.

CONCLUSIONS

With the six gene-specific DNA methylation signatures, patients in the training set were divided into low-risk and high- risk groups. What's more, gene-specific DNA methylation target genes were highly associated with protein phosphorylation, which indicated that further research on phosphorylation of target gene-coding protein might provide new sight on the treatment of CC.

Inferring Unknown Biological Function by Integration of GO Annotations and Gene Expression Data

Characterizing genes with semantic information is an important process regarding the description of gene products. In spite that complete genomes of many organisms have been already sequenced, the biological functions of all of their genes are still unknown. Since experimentally studying the functions of those genes, one by one, would be unfeasible, new computational methods for gene functions inference are needed. We present here a novel computational approach for inferring biological function for a set of genes with previously unknown function, given a set of genes with well-known information. This approach is based on the premise that genes with similar behaviour should be grouped together. This is known as the guilt-by-association principle. Thus, it is possible to take advantage of clustering techniques to obtain groups of unknown genes that are co-clustered with genes that have well-known semantic information (GO annotations). Meaningful knowledge to infer unknown semantic information can therefore be provided by these well-known genes. We provide a method to explore the potential function of new genes according to those currently annotated. The results obtained indicate that the proposed approach could be a useful and effective tool when used by biologists to guide the inference of biological functions for recently discovered genes. Our work sets an important landmark in the field of identifying unknown gene functions through clustering, using an external source of biological input. A simple web interface to this proposal can be found at http://fich.unl.edu.ar/sinc/webdemo/gamma-am/.

Expression-based clustering of CAZyme-encoding genes of Aspergillus niger

Background

The Aspergillus niger genome contains a large repertoire of genes encoding carbohydrate active enzymes (CAZymes) that are targeted to plant polysaccharide degradation enabling A. niger to grow on a wide range of plant biomass substrates. Which genes need to be activated in certain environmental conditions depends on the composition of the available substrate. Previous studies have demonstrated the involvement of a number of transcriptional regulators in plant biomass degradation and have identified sets of target genes for each regulator. In this study, a broad transcriptional analysis was performed of the A. niger genes encoding (putative) plant polysaccharide degrading enzymes. Microarray data focusing on the initial response of A. niger to the presence of plant biomass related carbon sources were analyzed of a wild-type strain N402 that was grown on a large range of carbon sources and of the regulatory mutant strains ΔxlnR, ΔaraR, ΔamyR, ΔrhaR and ΔgalX that were grown on their specific inducing compounds.

Results

The cluster analysis of the expression data revealed several groups of co-regulated genes, which goes beyond the traditionally described co-regulated gene sets. Additional putative target genes of the selected regulators were identified, based on their expression profile. Notably, in several cases the expression profile puts questions on the function assignment of uncharacterized genes that was based on homology searches, highlighting the need for more extensive biochemical studies into the substrate specificity of enzymes encoded by these non-characterized genes. The data also revealed sets of genes that were upregulated in the regulatory mutants, suggesting interaction between the regulatory systems and a therefore even more complex overall regulatory network than has been reported so far.

Conclusions

Expression profiling on a large number of substrates provides better insight in the complex regulatory systems that drive the conversion of plant biomass by fungi. In addition, the data provides additional evidence in favor of and against the similarity-based functions assigned to uncharacterized genes.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4164-x) contains supplementary material, which is available to authorized users.

Incorporating gene ontology into fuzzy relational clustering of microarray gene expression data

The product of gene expression works together in the cell for each living organism in order to achieve different biological processes. Many proteins are involved in different roles depending on the environment of the organism for the functioning of the cell. In this paper, we propose gene ontology (GO) annotations based semi-supervised clustering algorithm called GO fuzzy relational clustering (GO-FRC) where one gene is allowed to be assigned to multiple clusters which are the most biologically relevant behavior of genes. In the clustering process, GO-FRC utilizes useful biological knowledge which is available in the form of a gene ontology, as a prior knowledge along with the gene expression data. The prior knowledge helps to improve the coherence of the groups concerning the knowledge field. The proposed GO-FRC has been tested on the two yeast (Saccharomyces cerevisiae) expression profiles datasets (Eisen and Dream5 yeast datasets) and compared with other state-of-the-art clustering algorithms. Experimental results imply that GO-FRC is able to produce more biologically relevant clusters with the use of the small amount of GO annotations.

Data-analysis strategies for image-based cell profiling

Image-based cell profiling is a high-throughput strategy for the quantification of phenotypic differences among a variety of cell populations. It paves the way to studying biological systems on a large scale by using chemical and genetic perturbations. The general workflow for this technology involves image acquisition with high-throughput microscopy systems and subsequent image processing and analysis. Here, we introduce the steps required to create high-quality image-based (i.e., morphological) profiles from a collection of microscopy images. We recommend techniques that have proven useful in each stage of the data analysis process, on the basis of the experience of 20 laboratories worldwide that are refining their image-based cell-profiling methodologies in pursuit of biological discovery. The recommended techniques cover alternatives that may suit various biological goals, experimental designs, and laboratories' preferences.

Time-resolved transcriptome analysis and lipid pathway reconstruction of the oleaginous green microalga Monoraphidium neglectum reveal a model for triacylglycerol and lipid hyperaccumulation

BACKGROUND

Oleaginous microalgae are promising production hosts for the sustainable generation of lipid-based bioproducts and as bioenergy carriers such as biodiesel. Transcriptomics of the lipid accumulation phase, triggered efficiently by nitrogen starvation, is a valuable approach for the identification of gene targets for metabolic engineering.

RESULTS

An explorative analysis of the detailed transcriptional response to different stages of nitrogen availability was performed in the oleaginous green alga Monoraphidium neglectum. Transcript data were correlated with metabolic data for cellular contents of starch and of different lipid fractions. A pronounced transcriptional down-regulation of photosynthesis became apparent in response to nitrogen starvation, whereas glucose catabolism was found to be up-regulated. An in-depth reconstruction and analysis of the pathways for glycerolipid, central carbon, and starch metabolism revealed that distinct transcriptional changes were generally found only for specific steps within a metabolic pathway. In addition to pathway analyses, the transcript data were also used to refine the current genome annotation. The transcriptome data were integrated into a database and complemented with data for other microalgae which were also subjected to nitrogen starvation. It is available at https://tdbmn.cebitec.uni-bielefeld.de.

CONCLUSIONS

Based on the transcriptional responses to different stages of nitrogen availability, a model for triacylglycerol and lipid hyperaccumulation is proposed, which involves transcriptional induction of thioesterases, differential regulation of lipases, and a re-routing of the central carbon metabolism. Over-expression of distinct thioesterases was identified to be a potential strategy to increase the oleaginous phenotype of M. neglectum, and furthermore specific lipases were identified as potential targets for future metabolic engineering approaches.

Ocean warming alters cellular metabolism and induces mortality in fish early life stages: A proteomic approach

Climate change has pervasive effects on marine ecosystems, altering biodiversity patterns, abundance and distribution of species, biological interactions, phenology, and organisms' physiology, performance and fitness. Fish early life stages have narrow thermal windows and are thus more vulnerable to further changes in water temperature. The aim of this study was to address the sensitivity and underlying molecular changes of larvae of a key fisheries species, the sea bream Sparus aurata, towards ocean warming. Larvae were exposed to three temperatures: 18°C (control), 24°C (warm) and 30°C (heat wave) for seven days. At the end of the assay, i) survival curves were plotted for each temperature treatment and ii) entire larvae were collected for proteomic analysis via 2D gel electrophoresis, image analysis and mass spectrometry. Survival decreased with increasing temperature, with no larvae surviving at 30°C. Therefore, proteomic analysis was only carried out for 18°C and 24°C. Larvae up-regulated protein folding and degradation, cytoskeletal re-organization, transcriptional regulation and the growth hormone while mostly down-regulating cargo transporting and porphyrin metabolism upon exposure to heat stress. No changes were detected in proteins related to energetic metabolism suggesting that larval fish may not have the energetic plasticity needed to sustain cellular protection in the long-term. These results indicate that despite proteome modulation, S. aurata larvae do not seem able to fully acclimate to higher temperatures as shown by the low survival rates. Consequently, elevated temperatures seem to have bottleneck effects during fish early life stages, and future ocean warming can potentially compromise recruitment's success of key fisheries species.

Toward reliable biomarker signatures in the age of liquid biopsies - how to standardize the small RNA-Seq workflow

Small RNA-Seq has emerged as a powerful tool in transcriptomics, gene expression profiling and biomarker discovery. Sequencing cell-free nucleic acids, particularly microRNA (miRNA), from liquid biopsies additionally provides exciting possibilities for molecular diagnostics, and might help establish disease-specific biomarker signatures. The complexity of the small RNA-Seq workflow, however, bears challenges and biases that researchers need to be aware of in order to generate high-quality data. Rigorous standardization and extensive validation are required to guarantee reliability, reproducibility and comparability of research findings. Hypotheses based on flawed experimental conditions can be inconsistent and even misleading. Comparable to the well-established MIQE guidelines for qPCR experiments, this work aims at establishing guidelines for experimental design and pre-analytical sample processing, standardization of library preparation and sequencing reactions, as well as facilitating data analysis. We highlight bottlenecks in small RNA-Seq experiments, point out the importance of stringent quality control and validation, and provide a primer for differential expression analysis and biomarker discovery. Following our recommendations will encourage better sequencing practice, increase experimental transparency and lead to more reproducible small RNA-Seq results. This will ultimately enhance the validity of biomarker signatures, and allow reliable and robust clinical predictions.

multiClust: An R-package for Identifying Biologically Relevant Clusters in Cancer Transcriptome Profiles

Clustering is carried out to identify patterns in transcriptomics profiles to determine clinically relevant subgroups of patients. Feature (gene) selection is a critical and an integral part of the process. Currently, there are many feature selection and clustering methods to identify the relevant genes and perform clustering of samples. However, choosing an appropriate methodology is difficult. In addition, extensive feature selection methods have not been supported by the available packages. Hence, we developed an integrative R-package called multiClust that allows researchers to experiment with the choice of combination of methods for gene selection and clustering with ease. Using multiClust, we identified the best performing clustering methodology in the context of clinical outcome. Our observations demonstrate that simple methods such as variance-based ranking perform well on the majority of data sets, provided that the appropriate number of genes is selected. However, different gene ranking and selection methods remain relevant as no methodology works for all studies.

Adaptional evolution of trichome in Caragana korshinskii to natural drought stress on the Loess Plateau, China

Caragana korshinskii is commonly employed to improve drought ecosystems on the Loess Plateau, although the molecular mechanism at work is poorly understood, particularly in terms of the plant's ability to tolerate drought stress. Water is the most severe limiting factor for plant growth on the Loess Plateau. The trichome is known to play an efficient role in reducing water loss through decreasing the rate of transpiration, so in this study, we focused on the trichome‐related gene expression of ecological adaptation in C. korshinskii under low precipitation conditions. In order to explore the responses of trichomes to drought, we selected two experimental sites from wet to dry along the Loess Plateau latitude gradient for observation. Micro‐phenomena through which trichomes grew denser and larger under reduced precipitation were observed using a scanning electron microscope; de novo transcriptomes and quantitative PCR were then used to explore and verify gene expression patterns of C. korshinskii trichomes. Results showed that GIS2,TTG1, and GL2 were upregulated (as key positive‐regulated genes on trichome development), while CPC was downregulated (negative‐regulated gene). Taken together, our data indicate that downstream genes of gibberellin and cytokinin signaling pathways, alongside several cytoskeleton‐related genes, contribute to modulating trichome development to enhance transpiration resistance ability and increase the resistance to drought stress in C. korshinskii.

Computational gene expression profiling under salt stress reveals patterns of co-expression

Plants respond differently to environmental conditions. Among various abiotic stresses, salt stress is a condition where excess salt in soil causes inhibition of plant growth. To understand the response of plants to the stress conditions, identification of the responsible genes is required. Clustering is a data mining technique used to group the genes with similar expression. The genes of a cluster show similar expression and function. We applied clustering algorithms on gene expression data of Solanum tuberosum showing differential expression in Capsicum annuum under salt stress. The clusters, which were common in multiple algorithms were taken further for analysis. Principal component analysis (PCA) further validated the findings of other cluster algorithms by visualizing their clusters in three-dimensional space. Functional annotation results revealed that most of the genes were involved in stress related responses. Our findings suggest that these algorithms may be helpful in the prediction of the function of co-expressed genes.

Identification of certain cancer-mediating genes using Gaussian fuzzy cluster validity index

In this article, we have used an index, called Gaussian fuzzy index (GFI), recently developed by the authors, based on the notion of fuzzy set theory, for validating the clusters obtained by a clustering algorithm applied on cancer gene expression data. GFI is then used for the identification of genes that have altered quite significantly from normal state to carcinogenic state with respect to their mRNA expression patterns. The effectiveness of the methodology has been demonstrated on three gene expression cancer datasets dealing with human lung, colon and leukemia. The performance of GFI is compared with 19 exiting cluster validity indices. The results are appropriately validated biologically and statistically. In this context, we have used biochemical pathways, p-value statistics of GO attributes, t-test and zscore for the validation of the results. It has been reported that GFI is capable of identifying high-quality enriched clusters of genes, and thereby is able to select more cancer-mediating genes.

Spectral gene set enrichment (SGSE)

Background

Gene set testing is typically performed in a supervised context to quantify the association between groups of genes and a clinical phenotype. In many cases, however, a gene set-based interpretation of genomic data is desired in the absence of a phenotype variable. Although methods exist for unsupervised gene set testing, they predominantly compute enrichment relative to clusters of the genomic variables with performance strongly dependent on the clustering algorithm and number of clusters.

Results

We propose a novel method, spectral gene set enrichment (SGSE), for unsupervised competitive testing of the association between gene sets and empirical data sources. SGSE first computes the statistical association between gene sets and principal components (PCs) using our principal component gene set enrichment (PCGSE) method. The overall statistical association between each gene set and the spectral structure of the data is then computed by combining the PC-level p-values using the weighted Z-method with weights set to the PC variance scaled by Tracy-Widom test p-values. Using simulated data, we show that the SGSE algorithm can accurately recover spectral features from noisy data. To illustrate the utility of our method on real data, we demonstrate the superior performance of the SGSE method relative to standard cluster-based techniques for testing the association between MSigDB gene sets and the variance structure of microarray gene expression data.

Conclusions

Unsupervised gene set testing can provide important information about the biological signal held in high-dimensional genomic data sets. Because it uses the association between gene sets and samples PCs to generate a measure of unsupervised enrichment, the SGSE method is independent of cluster or network creation algorithms and, most importantly, is able to utilize the statistical significance of PC eigenvalues to ignore elements of the data most likely to represent noise.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-015-0490-7) contains supplementary material, which is available to authorized users.

Transcriptional profiling of the spleen in progressive visceral leishmaniasis reveals mixed expression of type 1 and type 2 cytokine-responsive genes

BACKGROUND

The Syrian golden hamster (Mesocricetus aureus) has been used as a model to study infections caused by a number of human pathogens. Studies of immunopathogenesis in hamster infection models are challenging because of the limited availability of reagents needed to define cellular and molecular determinants.

RESULTS

We sequenced a hamster cDNA library and developed a first-generation custom cDNA microarray that included 5131 unique cDNAs enriched for immune response genes. We used this microarray to interrogate the hamster spleen response to Leishmania donovani, an intracellular protozoan that causes visceral leishmaniasis. The hamster model of visceral leishmaniasis is of particular interest because it recapitulates clinical and immunopathological features of human disease, including cachexia, massive splenomegaly, pancytopenia, immunosuppression, and ultimately death. In the microarray a differentially expressed transcript was identified as having at least a 2-fold change in expression between uninfected and infected groups and a False Discovery Rate of <5%. Following a relatively silent early phase of infection (at 7 and 14 days post-infection only 8 and 24 genes, respectively, were differentially expressed), there was dramatic upregulation of inflammatory and immune-related genes in the spleen (708 differentially expressed genes were evident at 28 days post-infection). The differentially expressed transcripts included genes involved in inflammation, immunity, and immune cell trafficking. Of particular interest there was concomitant upregulation of the IFN-γ and interleukin (IL)-4 signaling pathways, with increased expression of a battery of IFN-γ- and IL-4-responsive genes. The latter included genes characteristic of alternatively activated macrophages.

CONCLUSIONS

Transcriptional profiling was accomplished in the Syrian golden hamster, for which a fully annotated genome is not available. In the hamster model of visceral leishmaniasis, a robust and functional IFN-γ response did not restrain parasite load and progression of disease. This supports the accumulating evidence that macrophages are ineffectively activated to kill the parasite. The concomitant expression of IL-4/IL-13 and their downstream target genes, some of which were characteristic of alternative macrophage activation, are likely to contribute to this. Further dissection of mechanisms that lead to polarization of macrophages toward a permissive state is needed to fully understand the pathogenesis of visceral leishmaniasis.

Shifting-and-Scaling Correlation Based Biclustering Algorithm

The existence of various types of correlations among the expressions of a group of biologically significant genes poses challenges in developing effective methods of gene expression data analysis. The initial focus of computational biologists was to work with only absolute and shifting correlations. However, researchers have found that the ability to handle shifting-and-scaling correlation enables them to extract more biologically relevant and interesting patterns from gene microarray data. In this paper, we introduce an effective shifting-and-scaling correlation measure named Shifting and Scaling Similarity (SSSim), which can detect highly correlated gene pairs in any gene expression data. We also introduce a technique named Intensive Correlation Search (ICS) biclustering algorithm, which uses SSSim to extract biologically significant biclusters from a gene expression data set. The technique performs satisfactorily with a number of benchmarked gene expression data sets when evaluated in terms of functional categories in Gene Ontology database.

Noise Resistant Generalized Parametric Validity Index of Clustering for Gene Expression Data

Validity indices have been investigated for decades. However, since there is no study of noise-resistance performance of these indices in the literature, there is no guideline for determining the best clustering in noisy data sets, especially microarray data sets. In this paper, we propose a generalized parametric validity (GPV) index which employs two tunable parameters α and β to control the proportions of objects being considered to calculate the dissimilarities. The greatest advantage of the proposed GPV index is its noise-resistance ability, which results from the flexibility of tuning the parameters. Several rules are set to guide the selection of parameter values. To illustrate the noise-resistance performance of the proposed index, we evaluate the GPV index for assessing five clustering algorithms in two gene expression data simulation models with different noise levels and compare the ability of determining the number of clusters with eight existing indices. We also test the GPV in three groups of real gene expression data sets. The experimental results suggest that the proposed GPV index has superior noise-resistance ability and provides fairly accurate judgements.

A Simple but Powerful Heuristic Method for Accelerating k-Means Clustering of Large-Scale Data in Life Science

K-means clustering has been widely used to gain insight into biological systems from large-scale life science data. To quantify the similarities among biological data sets, Pearson correlation distance and standardized Euclidean distance are used most frequently; however, optimization methods have been largely unexplored. These two distance measurements are equivalent in the sense that they yield the same k-means clustering result for identical sets of k initial centroids. Thus, an efficient algorithm used for one is applicable to the other. Several optimization methods are available for the Euclidean distance and can be used for processing the standardized Euclidean distance; however, they are not customized for this context. We instead approached the problem by studying the properties of the Pearson correlation distance, and we invented a simple but powerful heuristic method for markedly pruning unnecessary computation while retaining the final solution. Tests using real biological data sets with 50-60K vectors of dimensions 10-2001 (~400 MB in size) demonstrated marked reduction in computation time for k = 10-500 in comparison with other state-of-the-art pruning methods such as Elkan's and Hamerly's algorithms. The BoostKCP software is available at http://mlab.cb.k.u-tokyo.ac.jp/~ichikawa/boostKCP/.

R. S. WebTool, a web server for random sampling-based significance evaluation of pairwise distances

Pairwise comparison of data vectors represents a large part of computational biology, especially with the continuous increase in genome-wide approaches yielding more information from more biological samples simultaneously. Gene clustering for function prediction as well as analyses of signalling pathways and the time-dependent dynamics of a system are common biological approaches that often rely on large dataset comparison. Different metrics can be used to evaluate the similarity between entities to be compared, such as correlation coefficients and distances. While the latter offers a more flexible way of measuring potential biological relationships between datasets, the significance of any given distance is highly dependent on the dataset and cannot be easily determined. Monte Carlo methods are robust approaches for evaluating the significance of distance values by multiple random permutations of the dataset followed by distance calculation. We have developed R. S. WebTool (http://rswebtool.kwaklab.org), a user-friendly online server for random sampling-based evaluation of distance significances that features an array of visualization and analysis tools to help non-bioinformaticist users extract significant relationships from random noise in distance-based dataset analyses.

Measure the Semantic Similarity of GO Terms Using Aggregate Information Content

The rapid development of gene ontology (GO) and huge amount of biomedical data annotated by GO terms necessitate computation of semantic similarity of GO terms and, in turn, measurement of functional similarity of genes based on their annotations. In this paper we propose a novel and efficient method to measure the semantic similarity of GO terms. The proposed method addresses the limitations in existing GO term similarity measurement techniques; it computes the semantic content of a GO term by considering the information content of all of its ancestor terms in the graph. The aggregate information content (AIC) of all ancestor terms of a GO term implicitly reflects the GO term's location in the GO graph and also represents how human beings use this GO term and all its ancestor terms to annotate genes. We show that semantic similarity of GO terms obtained by our method closely matches the human perception. Extensive experimental studies show that this novel method also outperforms all existing methods in terms of the correlation with gene expression data. We have developed web services for measuring semantic similarity of GO terms and functional similarity of genes using the proposed AIC method and other popular methods. These web services are available at http://bioinformatics.clemson.edu/G-SESAME.

Discovering up-regulated VEGF-C expression in swine umbilical vein endothelial cells by classical swine fever virus Shimen

Infection of domestic swine with the highly virulent Shimen strain of classical swine fever virus causes hemorrhagic lymphadenitis and diffuse hemorrhaging in infected swine. We analyzed patterns of gene expression for CSFV Shimen in swine umbilical vein endothelial cells (SUVECs). Transcription of the vascular endothelial growth factor (VEGF) C gene (VEGF-C) and translation of the corresponding protein were significantly up-regulated in SUVECs. Our findings suggest that VEGF-C is involved in mechanisms of acute infection caused by virulent strains of CSFV.

Pathway-based analysis of the hidden genetic heterogeneities in cancers

Many cancers apparently showing similar phenotypes are actually distinct at the molecular level, leading to very different responses to the same treatment. It has been recently demonstrated that pathway-based approaches are robust and reliable for genetic analysis of cancers. Nevertheless, it remains unclear whether such function-based approaches are useful in deciphering molecular heterogeneities in cancers. Therefore, we aimed to test this possibility in the present study. First, we used a NCI60 dataset to validate the ability of pathways to correctly partition samples. Next, we applied the proposed method to identify the hidden subtypes in diffuse large B-cell lymphoma (DLBCL). Finally, the clinical significance of the identified subtypes was verified using survival analysis. For the NCI60 dataset, we achieved highly accurate partitions that best fit the clinical cancer phenotypes. Subsequently, for a DLBCL dataset, we identified three hidden subtypes that showed very different 10-year overall survival rates (90%, 46% and 20%) and were highly significantly (P=0.008) correlated with the clinical survival rate. This study demonstrated that the pathway-based approach is promising for unveiling genetic heterogeneities in complex human diseases.

On the selection of appropriate distances for gene expression data clustering

BACKGROUND

Clustering is crucial for gene expression data analysis. As an unsupervised exploratory procedure its results can help researchers to gain insights and formulate new hypothesis about biological data from microarrays. Given different settings of microarray experiments, clustering proves itself as a versatile exploratory tool. It can help to unveil new cancer subtypes or to identify groups of genes that respond similarly to a specific experimental condition. In order to obtain useful clustering results, however, different parameters of the clustering procedure must be properly tuned. Besides the selection of the clustering method itself, determining which distance is going to be employed between data objects is probably one of the most difficult decisions.

RESULTS AND CONCLUSIONS

We analyze how different distances and clustering methods interact regarding their ability to cluster gene expression, i.e., microarray data. We study 15 distances along with four common clustering methods from the literature on a total of 52 gene expression microarray datasets. Distances are evaluated on a number of different scenarios including clustering of cancer tissues and genes from short time-series expression data, the two main clustering applications in gene expression. Our results support that the selection of an appropriate distance depends on the scenario in hand. Moreover, in each scenario, given the very same clustering method, significant differences in quality may arise from the selection of distinct distance measures. In fact, the selection of an appropriate distance measure can make the difference between meaningful and poor clustering outcomes, even for a suitable clustering method.

Cluster and Principal Component Analysis of Human Glioblastoma Multiforme (GBM) Tumor Proteome

BACKGROUND

Glioblastoma Multiforme (GBM) or grade IV astrocytoma is the most common and lethal adult malignant brain tumor. Several of the molecular alterations detected in gliomas may have diagnostic and/or prognostic implications. Proteomics has been widely applied in various areas of science, ranging from the deciphering of molecular pathogen nests of discuses.

METHODS

In this study proteins were extracted from the tumor and normal brain tissues and then the protein purity was evaluated by Bradford test and spectrophotometry. In this study, proteins were separated by 2-Dimensional Gel (2DG) electrophoresis method and the spots were then analyzed and compared using statistical data and specific software. Protein clustering analysis was performed on the list of proteins deemed significantly altered in glioblastoma tumors (t-test and one-way ANOVA; P< 0.05).

RESULTS

The 2D gel showed totally 876 spots. We reported, 172 spots were exhibited differently in expression level (fold > 2) for glioblastoma. On each analytical 2D gel, an average of 876 spots was observed. In this study, 188 spots exhibited up regulation of expression level, whereas the remaining 232 spots were decreased in glioblastoma tumor relative to normal tissue. Results demonstrate that functional clustering (up and down regulated) and Principal Component Analysis (PCA) has considerable merits in aiding the interpretation of proteomic data.

CONCLUSION

2D gel electrophoresis is the core of proteomics which permitted the separation of thousands of proteins. High resolution 2DE can resolve up to 5,000 proteins simultaneously. Using cluster analysis, we can also form groups of related variables, similar to what is practiced in factor analysis.

A gene ontology inferred from molecular networks

Ontologies have proven very useful for capturing knowledge as a hierarchy of terms and their interrelationships. In biology a major challenge has been to construct ontologies of gene function given incomplete biological knowledge and inconsistencies in how this knowledge is manually curated. Here we show that large networks of gene and protein interactions in Saccharomyces cerevisiae can be used to infer an ontology whose coverage and power are equivalent to those of the manually curated Gene Ontology (GO). The network-extracted ontology (NeXO) contains 4,123 biological terms and 5,766 term-term relations, capturing 58% of known cellular components. We also explore robust NeXO terms and term relations that were initially not cataloged in GO, a number of which have now been added based on our analysis. Using quantitative genetic interaction profiling and chemogenomics, we find further support for many of the uncharacterized terms identified by NeXO, including multisubunit structures related to protein trafficking or mitochondrial function. This work enables a shift from using ontologies to evaluate data to using data to construct and evaluate ontologies.

Integrative analysis of cancer-related signaling pathways

Identification and classification of cancer types and subtypes is a major issue in current cancer research. Whole genome expression profiling of cancer tissues is often the basis for such subtype classifications of tumors and different signatures for individual cancer types have been described. However, the search for best performing discriminatory gene-expression signatures covering more than one cancer type remains a relevant topic in cancer research as such a signature would help understanding the common changes in signaling networks in these disease types. In this work, we explore the idea of a top down approach for sample stratification based on a module-based network of cancer relevant signaling pathways. For assembly of this network, we consider several of the most established cancer pathways. We evaluate our sample stratification approach using expression data of human breast and ovarian cancer signatures. We show that our approach performs equally well to previously reported methods besides providing the advantage to classify different cancer types. Furthermore, it allows to identify common changes in network module activity of those cancer samples.

Comparative Analysis of Cluster Validity Indices in Identifying Some Possible Genes Mediating Certain Cancers

In this article, we compare the performance of 19 cluster validity indices, in identifying some possible genes mediating certain cancers, based on gene expression data. For the purpose of this comparison, we have developed a method. The proposed method involves cluster generation, selection of the best k-value or c-values, cluster identification, identifying the altered gene cluster, scoring an altered gene cluster and determining the best k-value or c-value exploring through biological repositories. The effectiveness of the method has been demonstrated on three gene expression data sets dealing with human lung cancer, colon cancer, and leukemia. Here, we have used three clustering algorithms, i.e., k-means, PAM and fuzzy c-means. We have used biochemical pathways related to these cancers and p-value statistics for validating the study.

Paradigm of tunable clustering using Binarization of Consensus Partition Matrices (Bi-CoPaM) for gene discovery

Clustering analysis has a growing role in the study of co-expressed genes for gene discovery. Conventional binary and fuzzy clustering do not embrace the biological reality that some genes may be irrelevant for a problem and not be assigned to a cluster, while other genes may participate in several biological functions and should simultaneously belong to multiple clusters. Also, these algorithms cannot generate tight clusters that focus on their cores or wide clusters that overlap and contain all possibly relevant genes. In this paper, a new clustering paradigm is proposed. In this paradigm, all three eventualities of a gene being exclusively assigned to a single cluster, being assigned to multiple clusters, and being not assigned to any cluster are possible. These possibilities are realised through the primary novelty of the introduction of tunable binarization techniques. Results from multiple clustering experiments are aggregated to generate one fuzzy consensus partition matrix (CoPaM), which is then binarized to obtain the final binary partitions. This is referred to as Binarization of Consensus Partition Matrices (Bi-CoPaM). The method has been tested with a set of synthetic datasets and a set of five real yeast cell-cycle datasets. The results demonstrate its validity in generating relevant tight, wide, and complementary clusters that can meet requirements of different gene discovery studies.

Quality assessment and data handling methods for Affymetrix Gene 1.0 ST arrays with variable RNA integrity

Background

RNA and microarray quality assessment form an integral part of gene expression analysis and, although methods such as the RNA integrity number (RIN) algorithm reliably asses RNA integrity, the relevance of RNA integrity in gene expression analysis as well as analysis methods to accommodate the possible effects of degradation requires further investigation. We investigated the relationship between RNA integrity and array quality on the commonly used Affymetrix Gene 1.0 ST array platform using reliable within-array and between-array quality assessment measures. The possibility of a transcript specific bias in the apparent effect of RNA degradation on the measured gene expression signal was evaluated after either excluding quality-flagged arrays or compensation for RNA degradation at different steps in the analysis.

Results

Using probe-level and inter-array quality metrics to assess 34 Gene 1.0 ST array datasets derived from historical, paired tumour and normal primary colorectal cancer samples, 7 arrays (20.6%), with a mean sample RIN of 3.2 (SD = 0.42), were flagged during array quality assessment while 10 arrays from samples with RINs < 7 passed quality assessment, including one sample with a RIN < 3. We detected a transcript length bias in RNA degradation in only 5.8% of annotated transcript clusters (p-value 0.05, FC ≥ |2|), with longer and shorter than average transcripts under- and overrepresented in quality-flagged samples respectively. Applying compensatory measures for RNA degradation performed at least as well as excluding quality-flagged arrays, as judged by hierarchical clustering, gene expression analysis and Ingenuity Pathway Analysis; importantly, use of these compensatory measures had the significant benefit of enabling lower quality array data from irreplaceable clinical samples to be retained in downstream analyses.

Conclusions

Here, we demonstrate an effective array-quality assessment strategy, which will allow the user to recognize lower quality arrays that can be included in the analysis once appropriate measures are applied to account for known or unknown sources of variation, such as array quality- and batch- effects, by implementing ComBat or Surrogate Variable Analysis. This approach of quality control and analysis will be especially useful for clinical samples with variable and low RNA qualities, with RIN scores ≥ 2.