Multi-level reproducibility of signature hubs in human interactome for breast cancer metastasis

Inferring Virus-Host relationship between HPV and its host Homo sapiens using protein interaction network

Human papilloma virus (HPV) is a serious threat to human life globally with over 100 genotypes including cancer causing high risk HPVs. Study on protein interaction maps of pathogens with their host is a recent trend in ‘omics’ era and has been practiced by researchers to find novel drug targets. In current study, we construct an integrated protein interaction map of HPV with its host human in Cytoscape and analyze it further by using various bioinformatics tools. We found out 2988 interactions between 12 HPV and 2061 human proteins among which we identified MYLK, CDK7, CDK1, CDK2, JAK1 and 6 other human proteins associated with multiple viral oncoproteins. The functional enrichment analysis of these top-notch key genes is performed using KEGG pathway and Gene Ontology analysis, which reveals that the gene set is enriched in cell cycle a crucial cellular process, and the second most important pathway in which the gene set is involved is viral carcinogenesis. Among the viral proteins, E7 has the highest number of associations in the network followed by E6, E2 and E5. We found out a group of genes which is not targeted by the existing drugs available for HPV infections. It can be concluded that the molecules found in this study could be potential targets and could be used by scientists in their drug design studies.

Dynamic modeling and network approaches for omics time course data: overview of computational approaches and applications

Inferring networks and dynamics of genes, proteins, cells and other biological entities from high-throughput biological omics data is a central and challenging issue in computational and systems biology. This is essential for understanding the complexity of human health, disease susceptibility and pathogenesis for Predictive, Preventive, Personalized and Participatory (P4) system and precision medicine. The delineation of the possible interactions of all genes/proteins in a genome/proteome is a task for which conventional experimental techniques are ill suited. Urgently needed are rapid and inexpensive computational and statistical methods that can identify interacting candidate disease genes or drug targets out of thousands that can be further investigated or validated by experimentations. Moreover, identifying biological dynamic systems, and simultaneously estimating the important kinetic structural and functional parameters, which may not be experimentally accessible could be important directions for drug-disease-gene network studies. In this article, we present an overview and comparison of recent developments of dynamic modeling and network approaches for time-course omics data, and their applications to various biological systems, health conditions and disease statuses. Moreover, various data reduction and analytical schemes ranging from mathematical to computational to statistical methods are compared including their merits, drawbacks and limitations. The most recent software, associated web resources and other potentials for the compared methods are also presented and discussed in detail.

Identification of transcription factor-miRNA-lncRNA feed-forward loops in breast cancer subtypes

Previous studies have demonstrated that transcription factor-miRNA-gene feed-forward loops (FFLs) played important roles in tumorigenesis. However, the lncRNA-involved FFLs have not been explored very well. Understanding the characteristics of lncRNA-involved FFLs in breast cancer subtypes may be a key question with clinical implications. In this study, we firstly constructed an integrated background regulatory network. Then, based on mRNA, miRNA, and lncRNA differential expression, we identified 147, 140, 284, 1031 dysregulated FFLs for luminal A, luminal B, HER2+ and basal-like subtype of breast cancer, respectively. Importantly, the known breast cancer-associated lncRNAs and miRNAs were enriched in the identified dysregulated FFLs. Through merging the dysregulated FFLs, we constructed the regulatory sub-network for each subtype. We found that all sub-networks were enriched in the well-known cancer-related pathways, such as cell cycle, pathways in cancer. Next, we also identified potential prognostic FFLs for subtypes of breast cancer, such as the hsa-miR-182-5p_JUN_XIST in basal-like subtype. Finally, we also discussed the potential application of inferring the candidate drugs for breast cancer treatment through modulating the lncRNA expression in the dysregulated FFLs. Collectively, this study elucidated the roles of lncRNA-involved FFLs in breast cancer subtypes, which could contribute to understanding breast cancer pathogenesis and improving the treatment.

Analyzing of Molecular Networks for Human Diseases and Drug Discovery

Molecular networks represent the interactions and relations of genes/proteins, and also encode molecular mechanisms of biological processes, development and diseases. Among the molecular networks, protein-protein Interaction Networks (PINs) have become effective platforms for uncovering the molecular mechanisms of diseases and drug discovery. PINs have been constructed for various organisms and utilized to solve many biological problems. In human, most proteins present their complex functions by interactions with other proteins, and the sum of these interactions represents the human protein interactome. Especially in the research on human disease and drugs, as an emerging tool, the PIN provides a platform to systematically explore the molecular complexities of specific diseases and the references for drug design. In this review, we summarized the commonly used approaches to aid disease research and drug discovery with PINs, including the network topological analysis, identification of novel pathways, drug targets and sub-network biomarkers for diseases. With the development of bioinformatic techniques and biological networks, PINs will play an increasingly important role in human disease research and drug discovery.

CytoGTA: A cytoscape plugin for identifying discriminative subnetwork markers using a game theoretic approach

In recent years, analyzing genome-wide expression profiles to find genetic markers has received much attention as a challenging field of research aiming at unveiling biological mechanisms behind complex disorders. The identification of reliable and reproducible markers has lately been achieved by integrating genome-scale functional relationships and transcriptome datasets, and a number of algorithms have been developed to support this strategy. In this paper, we present a promising and easily applicable tool to accomplish this goal, namely CytoGTA, which is a Cytoscape plug-in that relies on an optimistic game theoretic approach (GTA) for identifying subnetwork markers. Given transcriptomic data of two phenotype classes and interactome data, this plug-in offers discriminative markers for the two classes. The high performance of CytoGTA would not have been achieved if the strategy of GTA was not implemented in Cytoscape. This plug-in provides a simple-to-use platform, convenient for biological researchers to interactively work with and visualize the structure of subnetwork markers. CytoGTA is one of the few available Cytoscape plug-ins for marker identification, which shows superior performance to existing methods.

Computational dynamic approaches for temporal omics data with applications to systems medicine

Modeling and predicting biological dynamic systems and simultaneously estimating the kinetic structural and functional parameters are extremely important in systems and computational biology. This is key for understanding the complexity of the human health, drug response, disease susceptibility and pathogenesis for systems medicine. Temporal omics data used to measure the dynamic biological systems are essentials to discover complex biological interactions and clinical mechanism and causations. However, the delineation of the possible associations and causalities of genes, proteins, metabolites, cells and other biological entities from high throughput time course omics data is challenging for which conventional experimental techniques are not suited in the big omics era. In this paper, we present various recently developed dynamic trajectory and causal network approaches for temporal omics data, which are extremely useful for those researchers who want to start working in this challenging research area. Moreover, applications to various biological systems, health conditions and disease status, and examples that summarize the state-of-the art performances depending on different specific mining tasks are presented. We critically discuss the merits, drawbacks and limitations of the approaches, and the associated main challenges for the years ahead. The most recent computing tools and software to analyze specific problem type, associated platform resources, and other potentials for the dynamic trajectory and interaction methods are also presented and discussed in detail.

GTA: a game theoretic approach to identifying cancer subnetwork markers

The identification of genetic markers (e.g. genes, pathways and subnetworks) for cancer has been one of the most challenging research areas in recent years.

Integrative network-based approach identifies key genetic elements in breast invasive carcinoma

BACKGROUND

Breast cancer is a genetically heterogeneous type of cancer that belongs to the most prevalent types with a high mortality rate. Treatment and prognosis of breast cancer would profit largely from a correct classification and identification of genetic key drivers and major determinants driving the tumorigenesis process. In the light of the availability of tumor genomic and epigenomic data from different sources and experiments, new integrative approaches are needed to boost the probability of identifying such genetic key drivers. We present here an integrative network-based approach that is able to associate regulatory network interactions with the development of breast carcinoma by integrating information from gene expression, DNA methylation, miRNA expression, and somatic mutation datasets.

RESULTS

Our results showed strong association between regulatory elements from different data sources in terms of the mutual regulatory influence and genomic proximity. By analyzing different types of regulatory interactions, TF-gene, miRNA-mRNA, and proximity analysis of somatic variants, we identified 106 genes, 68 miRNAs, and 9 mutations that are candidate drivers of oncogenic processes in breast cancer. Moreover, we unraveled regulatory interactions among these key drivers and the other elements in the breast cancer network. Intriguingly, about one third of the identified driver genes are targeted by known anti-cancer drugs and the majority of the identified key miRNAs are implicated in cancerogenesis of multiple organs. Also, the identified driver mutations likely cause damaging effects on protein functions. The constructed gene network and the identified key drivers were compared to well-established network-based methods.

CONCLUSION

The integrated molecular analysis enabled by the presented network-based approach substantially expands our knowledge base of prospective genomic drivers of genes, miRNAs, and mutations. For a good part of the identified key drivers there exists solid evidence for involvement in the development of breast carcinomas. Our approach also unraveled the complex regulatory interactions comprising the identified key drivers. These genomic drivers could be further investigated in the wet lab as potential candidates for new drug targets. This integrative approach can be applied in a similar fashion to other cancer types, complex diseases, or for studying cellular differentiation processes.

Molecular interaction networks for the analysis of human disease: utility, limitations, and considerations

High-throughput '-omics' data can be combined with large-scale molecular interaction networks, for example, protein-protein interaction networks, to provide a unique framework for the investigation of human molecular biology. Interest in these integrative '-omics' methods is growing rapidly because of their potential to understand complexity and association with disease; such approaches have a focus on associations between phenotype and "network-type." The potential of this research is enticing, yet there remain a series of important considerations. Here, we discuss interaction data selection, data quality, the relative merits of using data from large high-throughput studies versus a meta-database of smaller literature-curated studies, and possible issues of sociological or inspection bias in interaction data. Other work underway, especially international consortia to establish data formats, quality standards and address data redundancy, and the improvements these efforts are making to the field, is also evaluated. We present options for researchers intending to use large-scale molecular interaction networks as a functional context for protein or gene expression data, including microRNAs, especially in the context of human disease.

VAN: an R package for identifying biologically perturbed networks via differential variability analysis

BACKGROUND

Large-scale molecular interaction networks are dynamic in nature and are of special interest in the analysis of complex diseases, which are characterized by network-level perturbations rather than changes in individual genes/proteins. The methods developed for the identification of differentially expressed genes or gene sets are not suitable for network-level analyses. Consequently, bioinformatics approaches that enable a joint analysis of high-throughput transcriptomics datasets and large-scale molecular interaction networks for identifying perturbed networks are gaining popularity. Typically, these approaches require the sequential application of multiple bioinformatics techniques - ID mapping, network analysis, and network visualization. Here, we present the Variability Analysis in Networks (VAN) software package: a collection of R functions to streamline this bioinformatics analysis.

FINDINGS

VAN determines whether there are network-level perturbations across biological states of interest. It first identifies hubs (densely connected proteins/microRNAs) in a network and then uses them to extract network modules (comprising of a hub and all its interaction partners). The function identifySignificantHubs identifies dysregulated modules (i.e. modules with changes in expression correlation between a hub and its interaction partners) using a single expression and network dataset. The function summarizeHubData identifies dysregulated modules based on a meta-analysis of multiple expression and/or network datasets. VAN also converts protein identifiers present in a MITAB-formatted interaction network to gene identifiers (UniProt identifier to Entrez identifier or gene symbol using the function generatePpiMap) and generates microRNA-gene interaction networks using TargetScan and Microcosm databases (generateMicroRnaMap). The function obtainCancerInfo is used to identify hubs (corresponding to significantly perturbed modules) that are already causally associated with cancer(s) in the Cancer Gene Census database. Additionally, VAN supports the visualization of changes to network modules in R and Cytoscape (visualizeNetwork and obtainPairSubset, respectively). We demonstrate the utility of VAN using a gene expression data from metastatic melanoma and a protein-protein interaction network from the Human Protein Reference Database.

CONCLUSIONS

Our package provides a comprehensive and user-friendly platform for the integrative analysis of -omics data to identify disease-associated network modules. This bioinformatics approach, which is essentially focused on the question of explaining phenotype with a 'network type' and in particular, how regulation is changing among different states of interest, is relevant to many questions including those related to network perturbations across developmental timelines.

Topology of molecular interaction networks

Molecular interactions are often represented as network models which have become the common language of many areas of biology. Graphs serve as convenient mathematical representations of network models and have themselves become objects of study. Their topology has been intensively researched over the last decade after evidence was found that they share underlying design principles with many other types of networks.Initial studies suggested that molecular interaction network topology is related to biological function and evolution. However, further whole-network analyses did not lead to a unified view on what this relation may look like, with conclusions highly dependent on the type of molecular interactions considered and the metrics used to study them. It is unclear whether global network topology drives function, as suggested by some researchers, or whether it is simply a byproduct of evolution or even an artefact of representing complex molecular interaction networks as graphs.Nevertheless, network biology has progressed significantly over the last years. We review the literature, focusing on two major developments. First, realizing that molecular interaction networks can be naturally decomposed into subsystems (such as modules and pathways), topology is increasingly studied locally rather than globally. Second, there is a move from a descriptive approach to a predictive one: rather than correlating biological network topology to generic properties such as robustness, it is used to predict specific functions or phenotypes.Taken together, this change in focus from globally descriptive to locally predictive points to new avenues of research. In particular, multi-scale approaches are developments promising to drive the study of molecular interaction networks further.

Genes dysregulated to different extent or oppositely in estrogen receptor-positive and estrogen receptor-negative breast cancers

Background

Directly comparing gene expression profiles of estrogen receptor-positive (ER+) and estrogen receptor-negative (ER−) breast cancers cannot determine whether differentially expressed genes between these two subtypes result from dysregulated expression in ER+ cancer or ER− cancer versus normal controls, and thus would miss critical information for elucidating the transcriptomic difference between the two subtypes.

Principal Findings

Using microarray datasets from TCGA, we classified the genes dysregulated in both ER+ and ER− cancers versus normal controls into two classes: (i) genes dysregulated in the same direction but to a different extent, and (ii) genes dysregulated to opposite directions, and then validated the two classes in RNA-sequencing datasets of independent cohorts. We showed that the genes dysregulated to a larger extent in ER+ cancers than in ER− cancers enriched in glycerophospholipid and polysaccharide metabolic processes, while the genes dysregulated to a larger extent in ER− cancers than in ER+ cancers enriched in cell proliferation. Phosphorylase kinase and enzymes of glycosylphosphatidylinositol (GPI) anchor biosynthesis were upregulated to a larger extent in ER+ cancers than in ER− cancers, whereas glycogen synthase and phospholipase A2 were downregulated to a larger extent in ER+ cancers than in ER− cancers. We also found that the genes oppositely dysregulated in the two subtypes significantly enriched with known cancer genes and tended to closely collaborate with the cancer genes. Furthermore, we showed the possibility that these oppositely dysregulated genes could contribute to carcinogenesis of ER+ and ER− cancers through rewiring different subpathways.

Conclusions

GPI-anchor biosynthesis and glycogenolysis were elevated and hydrolysis of phospholipids was depleted to a larger extent in ER+ cancers than in ER− cancers. Our findings indicate that the genes oppositely dysregulated in the two subtypes are potential cancer genes which could contribute to carcinogenesis of both ER+ and ER− cancers through rewiring different subpathways.

Disturbed protein-protein interaction networks in metastatic melanoma are associated with worse prognosis and increased functional mutation burden

For disseminated melanoma, new prognostic biomarkers and therapeutic targets are urgently needed. The organization of protein-protein interaction networks was assessed via the transcriptomes of four independent studies of metastatic melanoma and related to clinical outcome and MAP-kinase pathway mutations (BRAF/NRAS). We also examined patient outcome-related differences in a predicted network of microRNAs and their targets. The 32 hub genes with the most reproducible survival-related disturbances in co-expression with their protein partner genes included oncogenes and tumor suppressors, previously known correlates of prognosis, and other proteins not previously associated with melanoma outcome. Notably, this network-based gene set could classify patients according to clinical outcomes with 67-80% accuracy among cohorts. Reproducibly disturbed networks were also more likely to have a higher functional mutation burden than would be expected by chance. The disturbed regions of networks are therefore markers of clinically relevant, selectable tumor evolution in melanoma which may carry driver mutations.

Extracting a few functionally reproducible biomarkers to build robust subnetwork-based classifiers for the diagnosis of cancer

In microarray-based case-control studies of a disease, people often attempt to identify a few diagnostic or prognostic markers amongst the most significant differentially expressed (DE) genes. However, the reproducibility of DE genes identified in different studies for a disease is typically very low. To tackle the problem, we could evaluate the reproducibility of DE genes across studies and define robust markers for disease diagnosis using disease-associated protein-protein interaction (PPI) subnetwork. Using datasets for four cancer types, we found that the most significant DE genes in cancer exhibit consistent up- or down-regulation in different datasets. For each cancer type, the 5 (or 10) most significant DE genes separately extracted from different datasets tend to be significantly coexpressed and closely connected in the PPI subnetwork, thereby indicating that they are highly reproducible at the PPI level. Consequently, we were able to build robust subnetwork-based classifiers for cancer diagnosis.

Systems biology for molecular life sciences and its impact in biomedicine

Modern systems biology is already contributing to a radical transformation of molecular life sciences and biomedicine, and it is expected to have a real impact in the clinical setting in the next years. In this review, the emergence of systems biology is contextualized with a historic overview, and its present state is depicted. The present and expected future contribution of systems biology to the development of molecular medicine is underscored. Concerning the present situation, this review includes a reflection on the "inflation" of biological data and the urgent need for tools and procedures to make hidden information emerge. Descriptions of the impact of networks and models and the available resources and tools for applying them in systems biology approaches to molecular medicine are provided as well. The actual current impact of systems biology in molecular medicine is illustrated, reviewing two cases, namely, those of systems pharmacology and cancer systems biology. Finally, some of the expected contributions of systems biology to the immediate future of molecular medicine are commented.

Differential network entropy reveals cancer system hallmarks

The cellular phenotype is described by a complex network of molecular interactions. Elucidating network properties that distinguish disease from the healthy cellular state is therefore of critical importance for gaining systems-level insights into disease mechanisms and ultimately for developing improved therapies. By integrating gene expression data with a protein interaction network we here demonstrate that cancer cells are characterised by an increase in network entropy. In addition, we formally demonstrate that gene expression differences between normal and cancer tissue are anticorrelated with local network entropy changes, thus providing a systemic link between gene expression changes at the nodes and their local correlation patterns. In particular, we find that genes which drive cell-proliferation in cancer cells and which often encode oncogenes are associated with reductions in network entropy. These findings may have potential implications for identifying novel drug targets.

A Steiner tree-based method for biomarker discovery and classification in breast cancer metastasis

Background

Metastatic breast cancer is a leading cause of cancer-related deaths in women worldwide. DNA microarray has become an important tool to help identify biomarker genes for improving the prognosis of breast cancer. Recently, it was shown that pathway-level relationships between genes can be incorporated to build more robust classification models and to obtain more useful biological insight from such models. Due to the unavailability of complete pathways, protein-protein interaction (PPI) network is becoming more popular to researcher and opens a new way to investigate the developmental process of breast cancer.

Methods

In this study, a network-based method is proposed to combine microarray gene expression profiles and PPI network for biomarker discovery for breast cancer metastasis. The key idea in our approach is to identify a small number of genes to connect differentially expressed genes into a single component in a PPI network; these intermediate genes contain important information about the pathways involved in metastasis and have a high probability of being biomarkers.

Results

We applied this approach on two breast cancer microarray datasets, and for both cases we identified significant numbers of well-known biomarker genes for breast cancer metastasis. Those selected genes are significantly enriched with biological processes and pathways related to cancer carcinogenic process, and, importantly, have much higher stability across different datasets than in previous studies. Furthermore, our selected genes significantly increased cross-data classification accuracy of breast cancer metastasis.

Conclusions

The randomized Steiner tree based approach described in this study is a new way to discover biomarker genes for breast cancer, and improves the prediction accuracy of metastasis. Though the analysis is limited here only to breast cancer, it can be easily applied to other diseases.

Network analysis of gene lists for finding reproducible prognostic breast cancer gene signatures

Many genome-scale studies in molecular biology deliver results in the form of a ranked list of gene names, accordingly to some scoring method. There is always the question how many top-ranked genes to consider for further analysis, for example, in order creating a diagnostic or predictive gene signature for a disease. This question is usually approached from a statistical point of view, without considering any biological properties of top-ranked genes or how they are related to each other functionally. Here we suggest a new method for selecting a number of genes in a ranked gene list such that this set forms the Optimally Functionally Enriched Network (OFTEN), formed by known physical interactions between genes or their products. The method allows associating a network with the gene list, providing easier interpretation of the results and classifying the genes or proteins accordingly to their position in the resulting network. We demonstrate the method on four breast cancer datasets and show that 1) the resulting gene signatures are more reproducible from one dataset to another compared to standard statistical procedures and 2) the overlap of these signatures has significant prognostic potential. The method is implemented in BiNoM Cytoscape plugin (http://binom.curie.fr).

Reproducibility and concordance of differential DNA methylation and gene expression in cancer

Background

Hundreds of genes with differential DNA methylation of promoters have been identified for various cancers. However, the reproducibility of differential DNA methylation discoveries for cancer and the relationship between DNA methylation and aberrant gene expression have not been systematically analysed.

Methodology/Principal Findings

Using array data for seven types of cancers, we first evaluated the effects of experimental batches on differential DNA methylation detection. Second, we compared the directions of DNA methylation changes detected from different datasets for the same cancer. Third, we evaluated the concordance between methylation and gene expression changes. Finally, we compared DNA methylation changes in different cancers. For a given cancer, the directions of methylation and expression changes detected from different datasets, excluding potential batch effects, were highly consistent. In different cancers, DNA hypermethylation was highly inversely correlated with the down-regulation of gene expression, whereas hypomethylation was only weakly correlated with the up-regulation of genes. Finally, we found that genes commonly hypomethylated in different cancers primarily performed functions associated with chronic inflammation, such as ‘keratinization’, ‘chemotaxis’ and ‘immune response’.

Conclusions

Batch effects could greatly affect the discovery of DNA methylation biomarkers. For a particular cancer, both differential DNA methylation and gene expression can be reproducibly detected from different studies with no batch effects. While DNA hypermethylation is significantly linked to gene down-regulation, hypomethylation is only weakly correlated with gene up-regulation and is likely to be linked to chronic inflammation.

Protein interactions: mapping interactome networks to support drug target discovery and selection

Proteins are biomolecular structures that build the microscopic working machinery of any living system. Proteins within the cells and biological systems do not act alone, but rather team up into macromolecular structures enclosing intricate physicochemical dynamic connections to undertake biological functions. A critical step towards unraveling the complex molecular relationships in living systems is the mapping of protein-to-protein physical "interactions". The complete map of protein interactions that can occur in a living organism is called the "interactome". Achieving an adequate atlas of all the protein interactions within a living system should allow to build its interaction network and to identity the "central nodes" that can be critical for the function, the homeostasis, and the movement of such system. Focusing on human studies, the data about the human interactome are most relevant for current biomedical research, because it is clear that the location of the proteins in the interactome network will allow to evaluate their centrality and to redefine the potential value of each protein as a drug target. This chapter presents our current knowledge on the human protein-protein interactome and explains how such knowledge can help us to select adequate targets for drugs.

Targeting oncogenic protein-protein interactions by diversity oriented synthesis and combinatorial chemistry approaches

We are currently witnessing a decline in the development of efficient new anticancer drugs, despite the salient efforts made on all fronts of cancer drug discovery. This trend presumably relates to the substantial heterogeneity and the inherent biological complexity of cancer, which hinder drug development success. Protein-protein interactions (PPIs) are key players in numerous cellular processes and aberrant interruption of this complex network provides a basis for various disease states, including cancer. Thus, it is now believed that cancer drug discovery, in addition to the design of single-targeted bioactive compounds, should also incorporate diversity-oriented synthesis (DOS) and other combinatorial strategies in order to exploit the ability of multi-functional scaffolds to modulate multiple protein-protein interactions (biological hubs). Throughout the review, we highlight the chemistry driven approaches to access diversity space for the discovery of small molecules that disrupt oncogenic PPIs, namely the p53-Mdm2, Bcl-2/Bcl-xL-BH3, Myc-Max, and p53-Mdmx/Mdm2 interactions.

Extensive increase of microarray signals in cancers calls for novel normalization assumptions

When using microarray data for studying a complex disease such as cancer, it is a common practice to normalize data to force all arrays to have the same distribution of probe intensities regardless of the biological groups of samples. The assumption underlying such normalization is that in a disease the majority of genes are not differentially expressed genes (DE genes) and the numbers of up- and down-regulated genes are roughly equal. However, accumulated evidences suggest gene expressions could be widely altered in cancer, so we need to evaluate the sensitivities of biological discoveries to violation of the normalization assumption. Here, we analyzed 7 large Affymetrix datasets of pair-matched normal and cancer samples for cancers collected in the NCBI GEO database. We showed that in 6 of these 7 datasets, the medians of perfect match (PM) probe intensities increased in cancer state and the increases were significant in three datasets, suggesting the assumption that all arrays have the same median probe intensities regardless of the biological groups of samples might be misleading. Then, we evaluated the effects of three currently most widely used normalization algorithms (RMA, MAS5.0 and dChip) on the selection of DE genes by comparing them with LVS which relies less on the above-mentioned assumption. The results showed using RMA, MAS5.0 and dChip may produce lots of false results of down-regulated DE genes while missing many up-regulated DE genes. At least for cancer study, normalizing all arrays to have the same distribution of probe intensities regardless of the biological groups of samples might be misleading. Thus, most current normalizations based on unreliable assumptions may distort biological differences between normal and cancer samples. The LVS algorithm might perform relatively well due to that it relies less on the above-mentioned assumption. Also, our results indicate that genes may be widely up-regulated in most human cancer.

Evaluating the consistency of differential expression of microRNA detected in human cancers

Differential expression of microRNA (miRNA) is involved in many human diseases and could potentially be used as a biomarker for disease diagnosis, prognosis, and therapy. However, inconsistency has often been found among differentially expressed miRNAs identified in various studies when using miRNA arrays for a particular disease such as a cancer. Before broadly applying miRNA arrays in a clinical setting, it is critical to evaluate inconsistent discoveries in a rational way. Thus, using data sets from 2 types of cancers, our study shows that the differentially expressed miRNAs detected from multiple experiments for each cancer exhibit stable regulation direction. This result also indicates that miRNA arrays could be used to reliably capture the signals of the regulation direction of differentially expressed miRNAs in cancer. We then assumed that 2 differentially expressed miRNAs with the same regulation direction in a particular cancer play similar functional roles if they regulate the same set of cancer-associated genes. On the basis of this hypothesis, we proposed a score to assess the functional consistency between differentially expressed miRNAs separately extracted from multiple studies for a particular cancer. We showed although lists of differentially expressed miRNAs identified from different studies for each cancer were highly variable, they were rather consistent at the level of function. Thus, the detection of differentially expressed miRNAs in various experiments for a certain disease tends to be functionally reproducible and capture functionally related differential expression of miRNAs in the disease.