EDGE: extraction and analysis of differential gene expression

EDGE (Extraction of Differential Gene Expression) is an open source, point-and-click software program for the significance analysis of DNA microarray experiments. EDGE can perform both standard and time course differential expression analysis. The functions are based on newly developed statistical theory and methods. This document introduces the EDGE software package.     

生物信息学基因表达差异分析

基因的差异化表达由多种因素共同导致,并且与许多疾病的发生和发展有密切联系,对差异化表达的基因进行生物信息学以及生物统计学的分析对于研究细胞调节机制和疾病机理有着重要意义。目前,对差异化表达的基因有以下几种主流的研究方法：DNA微阵列（DNA microarray）,抑制性消减杂交（SSH）,基因表达连续性分析（SAGE）,代表性差异分析（RDA）,以及mRNA差异显示PCR（mRNA DDRT-PCR）。目前许多基因差异化表达数据是建立在时段（time series）基础上,因此对基于时间变化的基因差异化表达分析变得尤为重要。本文将对差异化表达基因的几种主流方法进行详细阐述,并介绍一种基于傅里叶函数的时段基因差异化表达分析。     

EDGE: extraction and analysis of differential gene expression

The authors would like to apologise for the     

Clustering gene expression patterns.

Recent advances in biotechnology allow researchers to measure expression levels for thousands of genes simultaneously, across different conditions and over time. Analysis of data produced by such experiments offers potential insight into gene function and regulatory mechanisms. A key step in the analysis of gene expression data is the detection of groups of genes that manifest similar expression patterns. The corresponding algorithmic problem is to cluster multicondition gene expression patterns. In this paper we describe a novel clustering algorithm that was developed for analysis of gene expression data. We define an appropriate stochastic error model on the input, and prove that under the conditions of the model, the algorithm recovers the cluster structure with high probability. The running time of the algorithm on an n-gene dataset is O[n2[log(n)]c]. We also present a practical heuristic based on the same algorithmic ideas. The heuristic was implemented and its performance is demonstrated on simulated data and on real gene expression data, with very promising results.     

A novel strategy for identifying differential gene expression: An improved method of differential display analysis.

We propose a novel alternative approach, an advanced method for recently developed strategies, for identifying differentially expressed genes. Firstly, double-stranded cDNAs were digested using Sau3AI and the 3'-end restriction fragments of the cDNA were ligated to a double-stranded adapter. Next, the restriction fragments were directly amplified using several combinations of adapter-specific primers and FITC-labeled oligo dT primers. The selected cDNA fragments were displayed on a polyacrylamide gel. Neither nested PCR nor purification of 3'-end fragments are necessary. We examined the validity of this approach by evaluating gene expression changes during granulocytic differentiation of HL-60 cells. This method can theoretically detect almost all gene expression changes more rapidly and through simpler manipulations than by any other approach.     

Outlier sums for differential gene expression analysis

We propose a method for detecting genes that, in a disease group, exhibit unusually high gene expression in some but not all samples. This can be particularly useful in cancer studies, where mutations that can amplify or turn off gene expression often occur in only a minority of samples. In real and simulated examples, the new method often exhibits lower false discovery rates than simple t-statistic thresholding. We also compare our approach to the recent cancer profile outlier analysis proposal of Tomlins and others (2005).     

Global analysis of gene expression by differential display

Differential display (DD) is one of the most commonly used approaches for identifying differentially expressed genes. However, there has been lack of an accurate guidance on how many DD polymerase chain reaction (PCR) primer combinations are needed to display most of the genes expressed in a eukaryotic cell. This study critically evaluated the gene coverage by DD as a function of the number of arbitrary primers, the number of 3′ bases of an arbitrary primer required to completely match an mRNA target sequence, the additional 5′ base match(s) of arbitrary primers in first-strand cDNA recognition, and the length of mRNA tails being analyzed. The resulting new DD mathematical model predicts that 80 to 160 arbitrary 13mers, when used in combinations with 3 one-base anchored oligo-dT primers, would allow any given mRNA within a eukaryotic cell to be detected with a 74% to 93% probability, respectively. The prediction was supported by both computer simulation of the DD process and experimental data from a comprehensive fluorescent DD screening for target genes of tumor-suppressor p53. Thus, this work provides a theoretical foundation upon which global analysis of gene expression by DD can be pursued.     

TLRs: differential adapter utilization by toll-like receptors mediates TLR-specific patterns of gene expression

In response to microbial or environmental "danger" signals, represented by structural motifs not normally expressed by cells, Toll-like receptors mediate intracellular signaling that leads to inflammatory gene expression. In response to agonists, TLR aggregation enables the recruitment and/or activation of TLR-specific adapter molecules. To date, four adapter proteins have been identified: MyD88, TIRAP/Mal, TRIF/TICAM-1, and TIRP/TRAM/TICAM-2. The interaction of the different TLRs with distinct combinations of adapter molecules creates a platform to which additional kinases, transacting factors, and possibly other molecules are recruited, events that lead, ultimately, to gene expression. Given the rapidity with which such interactions have been described, we have attempted to summarize our current understanding of the adapters that are so essential for TLR signaling and provide a working model for future studies.     

Random forests-based differential analysis of gene sets for gene expression data

In DNA microarray studies, gene-set analysis (GSA) has become the focus of gene expression data analysis. GSA utilizes the gene expression profiles of functionally related gene sets in Gene Ontology (GO) categories or priori-defined biological classes to assess the significance of gene sets associated with clinical outcomes or phenotypes. Many statistical approaches have been proposed to determine whether such functionally related gene sets express differentially (enrichment and/or deletion) in variations of phenotypes. However, little attention has been given to the discriminatory power of gene sets and classification of patients.     

Quantitative set analysis for gene expression: a method to quantify gene set differential expression including gene-gene correlations

Enrichment analysis of gene sets is a popular approach that provides a functional interpretation of genome-wide expression data. Existing tests are affected by inter-gene correlations, resulting in a high Type I error. The most widely used test, Gene Set Enrichment Analysis, relies on computationally intensive permutations of sample labels to generate a null distribution that preserves gene–gene correlations. A more recent approach, CAMERA, attempts to correct for these correlations by estimating a variance inflation factor directly from the data. Although these methods generate P-values for detecting gene set activity, they are unable to produce confidence intervals or allow for post hoc comparisons. We have developed a new computational framework for Quantitative Set Analysis of Gene Expression (QuSAGE). QuSAGE accounts for inter-gene correlations, improves the estimation of the variance inflation factor and, rather than evaluating the deviation from a null hypothesis with a P-value, it quantifies gene-set activity with a complete probability density function. From this probability density function, P-values and confidence intervals can be extracted and post hoc analysis can be carried out while maintaining statistical traceability. Compared with Gene Set Enrichment Analysis and CAMERA, QuSAGE exhibits better sensitivity and specificity on real data profiling the response to interferon therapy (in chronic Hepatitis C virus patients) and Influenza A virus infection. QuSAGE is available as an R package, which includes the core functions for the method as well as functions to plot and visualize the results.     

MOST: detecting cancer differential gene expression

We propose a new statistics for the detection of differentially expressed genes when the genes are activated only in a subset of the samples. Statistics designed for this unconventional circumstance has proved to be valuable for most cancer studies, where oncogenes are activated for a small number of disease samples. Previous efforts made in this direction include cancer outlier profile analysis (Tomlins and others, 2005), outlier sum (Tibshirani and Hastie, 2007), and outlier robust t-statistics (Wu, 2007). We propose a new statistics called maximum ordered subset t-statistics (MOST) which seems to be natural when the number of activated samples is unknown. We compare MOST to other statistics and find that the proposed method often has more power then its competitors.