Identifying candidate genes for wood formation in poplar based on microarray network analysis and graph theory

In poplar, genetic research on wood properties is very important for the improvement of wood quality. Studies of wood formation genes at each developmental stage using modern biotechnology have often been limited to several genes or gene families. Because of the complex regulatory network involved in the co-expression and interactions of thousands of genes, however, the genetic mechanisms of wood formation must be surveyed on a genome-wide scale. In this study, we identified wood formation-related genes using a differentially co-expressed (DCE) gene subset approach based on biological networks inferred from microarray data. Gene co-expression networks in leaf, root, and wood tissues were first constructed and topologically analyzed using microarray data collected from the Gene Expression Omnibus. The DCE gene modules in wood-forming tissue were then detected based on graph theory, which was followed by gene ontology (GO) enrichment analysis and GO annotation of probe sets. Finally, 72 probe sets were identified in the largest cohesive subgroup of the DCE gene network in wood tissue, with most of the probe sets associated with wood formation-related biological processes and GO cellular component categories. The approach described in this paper provides an effective strategy to identify wood formation genes in poplar and should contribute to the better understanding of the genetic and molecular mechanisms underlying wood properties in trees.     

LncRNA and mRNA integration network reconstruction reveals novel key regulators in esophageal squamous-cell carcinoma

Many experimental and computational studies have identified key protein coding genes in initiation and progression of esophageal squamous cell carcinoma (ESCC). However, the number of researches that tried to reveal the role of long non-coding RNAs (lncRNAs) in ESCC has been limited. LncRNAs are one of the important regulators of cancers which are transcribed dominantly in the genome and in various conditions. The main goal of this study was to use a systems biology approach to predict novel lncRNAs as well as protein coding genes associated with ESCC and assess their prognostic values. By using microarray expression data for mRNAs and lncRNAs from a large number of ESCC patients, we utilized “Weighted Gene Co-expression Network Analysis” (WGCNA) method to make a big coding-non-coding gene co-expression network, and discovered important functional modules. Gene set enrichment and pathway analysis revealed major biological processes and pathways involved in these modules. After selecting some protein coding genes involved in biological processes and pathways related to cancer, we used “LncTar”, a computational tool to predict potential interactions between these genes and lncRNAs. By combining interaction results with Pearson correlations, we introduced some novel lncRNAs with putative key regulatory roles in the network. Survival analysis with Kaplan-Meier estimator and Log-rank test statistic confirmed that most of the introduced genes are associated with poor prognosis in ESCC. Overall, our study reveals novel protein coding genes and lncRNAs associated with ESCC, along with their predicted interactions. Based on the promising results of survival analysis, these genes can be used as good estimators of patients' survival, or even can be analyzed further as new potential signatures or targets for the therapy of ESCC disease.     

Bayesian Network Expansion Identifies New ROS and Biofilm Regulators

Signaling and regulatory pathways that guide gene expression have only been partially defined for most organisms. However, given the increasing number of microarray measurements, it may be possible to reconstruct such pathways and uncover missing connections directly from experimental data. Using a compendium of microarray gene expression data obtained from Escherichia coli, we constructed a series of Bayesian network models for the reactive oxygen species (ROS) pathway as defined by EcoCyc. A consensus Bayesian network model was generated using those networks sharing the top recovered score. This microarray-based network only partially agreed with the known ROS pathway curated from the literature and databases. A top network was then expanded to predict genes that could enhance the Bayesian network model using an algorithm we termed ‘BN+1’. This expansion procedure predicted many stress-related genes (e.g., dusB and uspE), and their possible interactions with other ROS pathway genes. A term enrichment method discovered that biofilm-associated microarray data usually contained high expression levels of both uspE and gadX. The predicted involvement of gene uspE in the ROS pathway and interactions between uspE and gadX were confirmed experimentally using E. coli reporter strains. Genes gadX and uspE showed a feedback relationship in regulating each other''s expression. Both genes were verified to regulate biofilm formation through gene knockout experiments. These data suggest that the BN+1 expansion method can faithfully uncover hidden or unknown genes for a selected pathway with significant biological roles. The presently reported BN+1 expansion method is a generalized approach applicable to the characterization and expansion of other biological pathways and living systems.     

Multi-membership gene regulation in pathway based microarray analysis

Background

Gene expression analysis has been intensively researched for more than a decade. Recently, there has been elevated interest in the integration of microarray data analysis with other types of biological knowledge in a holistic analytical approach. We propose a methodology that can be facilitated for pathway based microarray data analysis, based on the observation that a substantial proportion of genes present in biochemical pathway databases are members of a number of distinct pathways. Our methodology aims towards establishing the state of individual pathways, by identifying those truly affected by the experimental conditions based on the behaviour of such genes. For that purpose it considers all the pathways in which a gene participates and the general census of gene expression per pathway.

Results

We utilise hill climbing, simulated annealing and a genetic algorithm to analyse the consistency of the produced results, through the application of fuzzy adjusted rand indexes and hamming distance. All algorithms produce highly consistent genes to pathways allocations, revealing the contribution of genes to pathway functionality, in agreement with current pathway state visualisation techniques, with the simulated annealing search proving slightly superior in terms of efficiency.

Conclusions

We show that the expression values of genes, which are members of a number of biochemical pathways or modules, are the net effect of the contribution of each gene to these biochemical processes. We show that by manipulating the pathway and module contribution of such genes to follow underlying trends we can interpret microarray results centred on the behaviour of these genes.     

Pathway analysis using random forests classification and regression

MOTIVATION: Although numerous methods have been developed to better capture biological information from microarray data, commonly used single gene-based methods neglect interactions among genes and leave room for other novel approaches. For example, most classification and regression methods for microarray data are based on the whole set of genes and have not made use of pathway information. Pathway-based analysis in microarray studies may lead to more informative and relevant knowledge for biological researchers. RESULTS: In this paper, we describe a pathway-based classification and regression method using Random Forests to analyze gene expression data. The proposed methods allow researchers to rank important pathways from externally available databases, discover important genes, find pathway-based outlying cases and make full use of a continuous outcome variable in the regression setting. We also compared Random Forests with other machine learning methods using several datasets and found that Random Forests classification error rates were either the lowest or the second-lowest. By combining pathway information and novel statistical methods, this procedure represents a promising computational strategy in dissecting pathways and can provide biological insight into the study of microarray data. AVAILABILITY: Source code written in R is available from http://bioinformatics.med.yale.edu/pathway-analysis/rf.htm.     

Highly interconnected genes in disease-specific networks are enriched for disease-associated polymorphisms

Background

Complex diseases are associated with altered interactions between thousands of genes. We developed a novel method to identify and prioritize disease genes, which was generally applicable to complex diseases.

Results

We identified modules of highly interconnected genes in disease-specific networks derived from integrating gene-expression and protein interaction data. We examined if those modules were enriched for disease-associated SNPs, and could be used to find novel genes for functional studies. First, we analyzed publicly available gene expression microarray and genome-wide association study (GWAS) data from 13, highly diverse, complex diseases. In each disease, highly interconnected genes formed modules, which were significantly enriched for genes harboring disease-associated SNPs. To test if such modules could be used to find novel genes for functional studies, we repeated the analyses using our own gene expression microarray and GWAS data from seasonal allergic rhinitis. We identified a novel gene, FGF2, whose relevance was supported by functional studies using combined small interfering RNA-mediated knock-down and gene expression microarrays. The modules in the 13 complex diseases analyzed here tended to overlap and were enriched for pathways related to oncological, metabolic and inflammatory diseases. This suggested that this union of the modules would be associated with a general increase in susceptibility for complex diseases. Indeed, we found that this union was enriched with GWAS genes for 145 other complex diseases.

Conclusions

Modules of highly interconnected complex disease genes were enriched for disease-associated SNPs, and could be used to find novel genes for functional studies.     

Microarray expression profiling of tissues from mice with uniparental duplications of Chromosomes 7 and 11 to identify imprinted genes

Microarray analysis allows the screening of thousands of identifiable genes in a single experiment. The challenge of this approach is to combine the new technology with established genetic tools to associate genes with specific biological function. In this study we have designed a screen to identify imprinted genes from mice with uniparental duplications of proximal Chromosomes (Chrs) 7 and 11, using microarray analysis. By comparing the expression patterns in embryonic and newborn tissues of maternally versus paternally inherited proximal Chrs 7 and 11, we have correctly identified four out of five known imprinted genes represented on a microarray. We have additionally identified two novel imprinted candidate genes as well as a differentially expressed clone that is a potential downstream target. Interpretation of the microarray data requires careful preparation of age- and strain-matched samples and attention to detail in tissue dissection technique. Received: 15 March 2001 / Accepted: 13 June 2001     

Inference of disease-related molecular logic from systems-based microarray analysis

Computational analysis of gene expression data from microarrays has been useful for medical diagnosis and prognosis. The ability to analyze such data at the level of biological modules, rather than individual genes, has been recognized as important for improving our understanding of disease-related pathways. It has proved difficult, however, to infer pathways from microarray data by deriving modules of multiple synergistically interrelated genes, rather than individual genes. Here we propose a systems-based approach called Entropy Minimization and Boolean Parsimony (EMBP) that identifies, directly from gene expression data, modules of genes that are jointly associated with disease. Furthermore, the technique provides insight into the underlying biomolecular logic by inferring a logic function connecting the joint expression levels in a gene module with the outcome of disease. Coupled with biological knowledge, this information can be useful for identifying disease-related pathways, suggesting potential therapeutic approaches for interfering with the functions of such pathways. We present an example providing such gene modules associated with prostate cancer from publicly available gene expression data, and we successfully validate the results on additional independently derived data. Our results indicate a link between prostate cancer and cellular damage from oxidative stress combined with inhibition of apoptotic mechanisms normally triggered by such damage.     

Toxicogenomic approach for assessing toxicant-related disease

The problems of identifying environmental factors involved in the etiology of human disease and performing safety and risk assessments of drugs and chemicals have long been formidable issues. Three principal components for predicting potential human health risks are: (1) the diverse structure and properties of thousands of chemicals and other stressors in the environment; (2) the time and dose parameters that define the relationship between exposure and disease; and (3) the genetic diversity of organisms used as surrogates to determine adverse chemical effects. The global techniques evolving from successful genomics efforts are providing new exciting tools with which to address these intractable problems of environmental health and toxicology. In order to exploit the scientific opportunities, the National Institute of Environmental Health Sciences has created the National Center for Toxicogenomics (NCT). The primary mission of the NCT is to use gene expression technology, proteomics and metabolite profiling to create a reference knowledge base that will allow scientists to understand mechanisms of toxicity and to be able to predict the potential toxicity of new chemical entities and drugs. A principal scientific objective underpinning the use of microarray analysis of chemical exposures is to demonstrate the utility of signature profiling of the action of drugs or chemicals and to utilize microarray methodologies to determine biomarkers of exposure and potential adverse effects. The initial approach of the NCT is to utilize proof-of-principle experiments in an effort to "phenotypically anchor" the altered patterns of gene expression to conventional parameters of toxicity and to define dose and time relationships in which the expression of such signature genes may precede the development of overt toxicity. The microarray approach is used in conjunction with proteomic techniques to identify specific proteins that may serve as signature biomarkers. The longer-range goal of these efforts is to develop a reference relational database of chemical effects in biological systems (CEBS) that can be used to define common mechanisms of toxicity, chemical and drug actions, to define cellular pathways of response, injury and, ultimately, disease. In order to implement this strategy, the NCT has created a consortium of research organizations and private sector companies to actively collaborative in populating the database with high quality primary data. The evolution of discrete databases to a knowledge base of toxicogenomics will be accomplished through establishing relational interfaces with other sources of information on the structure and activity of chemicals such as that of the National Toxicology Program (NTP) and with databases annotating gene identity, sequence, and function.     

Prediction of regulatory pathways using mRNA expression and protein interaction data: application to identification of galactose regulatory pathway

We propose a novel technique that constructs gene regulatory networks from DNA microarray data and gene-protein databases and then applies Mason rule to systematically search for the most dominant regulators of the network. The algorithm then recommends the identified dominant regulator genes as the best candidates for future knock-out experiments. Actively choosing the genes for knock-out experiments allows optimal perturbation of the pathway and therefore produces the most informative DNA microarray data for pathway identification purposes. This approach is more practically advantageous in analysis of large pathways where the time and cost of DNA microarray data experiments can be reduced using the proposed optimal experiment design. The proposed method was successfully tested on the galactose regulatory network.     

KEGGanim: pathway animations for high-throughput data

MOTIVATION: Gene expression analysis with microarrays has become one of the most widely used high-throughput methods for gathering genome-wide functional data. Emerging -omics fields such as proteomics and interactomics introduce new information sources. With the rise of systems biology, researchers need to concentrate on entire complex pathways that guide individual genes and related processes. Bioinformatics methods are needed to link the existing knowledge about pathways with the growing amounts of experimental data. RESULTS: We present KEGGanim, a novel web-based tool for visualizing experimental data in biological pathways. KEGGanim produces animations and images of KEGG pathways using public or user uploaded high-throughput data. Pathway members are coloured according to experimental measurements, and animated over experimental conditions. KEGGanim visualization highlights dynamic changes over conditions and allows the user to observe important modules and key genes that influence the pathway. The simple user interface of KEGGanim provides options for filtering genes and experimental conditions. KEGGanim may be used with public or private data for 14 organisms with a large collection of public microarray data readily available. Most common gene and protein identifiers and microarray probesets are accepted for visualization input. AVAILABILITY: http://biit.cs.ut.ee/KEGGanim/.     

Average rank-based score to measure deregulation of molecular pathway gene sets

Background

Deregulation of biological pathways has been shown to be involved in the turmorigenesis of a variety of cancers. The co-regulation of pathways in tumor and normal tissues has not been studied in a systematic manner.

Results

In this study we propose a novel statistic named AR-score (average rank based score) to measure pathway activities based on microarray gene expression profiles. We calculate and compare the AR-scores of pathways in microarray datasets containing expression profiles for a wide range of cancer types as well as the corresponding normal tissues. We find that many pathways undergo significant activity changes in tumors with respect to normal tissues. AR-scores for a small subset of pathways are capable of distinguishing tumor from normal tissues or classifying tumor subtypes. In normal tissues many pathways are highly correlated in their activities, whereas their correlations reduce significantly in tumors and cancer cell lines. The co-expression of genes in the same pathways was also significantly perturbed in tumors.

Conclusions

The co-regulation of genes in the same pathways and co-regulation of different pathways are significantly perturbed in tumors versus normal tissues. Our method provides a useful tool for better understanding the mechanistic changes in tumors, which can also be used for exploring other biological problems.     

基于全局角度网络策略的复杂疾病风险通路识别

复杂疾病的发生发展与机体内生物学通路的功能紊乱有密切联系,从高通量数据出发,利用计算机辅助方法来研究疾病与通路间的关系具有重要意义.本文提出了一个新的基于网络的全局性通路识别方法.该方法利用蛋白质互作信息和通路的基因集组成信息构建复杂的蛋白质-通路网.然后,基于表达谱数据,通过随机游走算法从全局层面优化疾病风险通路.最终,通过扰动方式识别统计学显著的风险通路.将该网络运用于结肠直肠癌风险通路识别,识别出15个与结肠直肠癌发生与发展过程显著相关的通路.通过与其他通路识别方法(超几何检验,SPIA)相比较,该方法能够更有效识别出疾病相关的风险通路.     

An ensemble method for gene discovery based on DNA microarray data

The advent of DNA microarray technology has offered the promise of casting new insights onto deciphering secrets of life by monitoring activities of thousands of genes simultaneously. Current analyses of microarray data focus on precise classification of biological types, for example, tumor versus normal tissues. A further scientific challenging task is to extract disease-relevant genes from the bewildering amounts of raw data, which is one of the most critical themes in the post-genomic era, but it is generally ignored due to lack of an efficient approach. In this paper, we present a novel ensemble method for gene extraction that can be tailored to fulfill multiple biological tasks including (i) precise classification of biological types; (ii) disease gene mining; and (iii) target-driven gene networking. We also give a numerical application for(i) and (ii) using a public microarrary data set and set aside a separate paper to address (iii).