Microarray probe expression measures, data normalization and statistical validation

DNA microarray technology is a high-throughput method for gaining information on gene function. Microarray technology is based on deposition/synthesis, in an ordered manner, on a solid surface, of thousands of EST sequences/genes/oligonucleotides. Due to the high number of generated datapoints, computational tools are essential in microarray data analysis and mining to grasp knowledge from experimental results. In this review, we will focus on some of the methodologies actually available to define gene expression intensity measures, microarray data normalization, and statistical validation of differential expression.     

A generalized likelihood ratio test to identify differentially expressed genes from microarray data

MOTIVATION: Microarray technology emerges as a powerful tool in life science. One major application of microarray technology is to identify differentially expressed genes under various conditions. Currently, the statistical methods to analyze microarray data are generally unsatisfactory, mainly due to the lack of understanding of the distribution and error structure of microarray data. RESULTS: We develop a generalized likelihood ratio (GLR) test based on the two-component model proposed by Rocke and Durbin to identify differentially expressed genes from microarray data. Simulation studies show that the GLR test is more powerful than commonly used methods, like the fold-change method and the two-sample t-test. When applied to microarray data, the GLR test identifies more differentially expressed genes than the t-test, has a lower false discovery rate and shows more consistency over independently repeated experiments. AVAILABILITY: The approach is implemented in software called GLR, which is freely available for downloading at http://www.cc.utah.edu/~jw27c60     

Unbiased pattern detection in microarray data series

MOTIVATION: Following the advent of microarray technology in recent years, the challenge for biologists is to identify genes of interest from the thousands of genetic expression levels measured in each microarray experiment. In many cases the aim is to identify pattern in the data series generated by successive microarray measurements. RESULTS: Here we introduce a new method of detecting pattern in microarray data series which is independent of the nature of this pattern. Our approach provides a measure of the algorithmic compressibility of each data series. A series which is significantly compressible is much more likely to result from simple underlying mechanisms than series which are incompressible. Accordingly, the gene associated with a compressible series is more likely to be biologically significant. We test our method on microarray time series of yeast cell cycle and show that it blindly selects genes exhibiting the expected cyclic behaviour as well as detecting other forms of pattern. Our results successfully predict two independent non-microarray experimental studies.     

An inexpensive gel electrophoresis-based polymerase chain reaction method for quantifying mRNA levels

In order to engage their students in a core methodology of the new genomics era, an ever-increasing number of faculty at primarily undergraduate institutions are gaining access to microarray technology. Their students are conducting successful microarray experiments designed to address a variety of interesting questions. A next step in these teaching and research laboratory projects is often validation of the microarray data for individual selected genes. In the research community, this usually involves the use of real-time polymerase chain reaction (PCR), a technology that requires instrumentation and reagents that are prohibitively expensive for most undergraduate institutions. The results of a survey of faculty teaching undergraduates in classroom and research settings indicate a clear need for an alternative approach. We sought to develop an inexpensive and student-friendly gel electrophoresis-based PCR method for quantifying messenger RNA (mRNA) levels using undergraduate researchers as models for students in teaching and research laboratories. We compared the results for three selected genes measured by microarray analysis, real-time PCR, and the gel electrophoresis-based method. The data support the use of the gel electrophoresis-based method as an inexpensive, convenient, yet reliable alternative for quantifying mRNA levels in undergraduate laboratories.     

A simple and robust algorithm for microarray data clustering based on gene population-variance ratio metric

With the advent of the microarray technology, the field of life science has been greatly revolutionized, since this technique allows the simultaneous monitoring of the expression levels of thousands of genes in a particular organism. However, the statistical analysis of expression data has its own challenges, primarily because of the huge amount of data that is to be dealt with, and also because of the presence of noise, which is almost an inherent characteristic of microarray data. Clustering is one tool used to mine meaningful patterns from microarray data. In this paper, we present a novel method of clustering yeast microarray data, which is robust and yet simple to implement. It identifies the best clusters from a given dataset on the basis of the population of the clusters as well as the variance of the feature values of the members from the cluster-center. It has been found to yield satisfactory results even in the presence of noisy data.     

Multi-class cancer classification by total principal component regression (TPCR) using microarray gene expression data

DNA microarray technology provides a promising approach to the diagnosis and prognosis of tumors on a genome-wide scale by monitoring the expression levels of thousands of genes simultaneously. One problem arising from the use of microarray data is the difficulty to analyze the high-dimensional gene expression data, typically with thousands of variables (genes) and much fewer observations (samples), in which severe collinearity is often observed. This makes it difficult to apply directly the classical statistical methods to investigate microarray data. In this paper, total principal component regression (TPCR) was proposed to classify human tumors by extracting the latent variable structure underlying microarray data from the augmented subspace of both independent variables and dependent variables. One of the salient features of our method is that it takes into account not only the latent variable structure but also the errors in the microarray gene expression profiles (independent variables). The prediction performance of TPCR was evaluated by both leave-one-out and leave-half-out cross-validation using four well-known microarray datasets. The stabilities and reliabilities of the classification models were further assessed by re-randomization and permutation studies. A fast kernel algorithm was applied to decrease the computation time dramatically. (MATLAB source code is available upon request.).     

Microarrays in biology and medicine

The remarkable speed with which biotechnology has become critical to the practice of life sciences owes much to a series of technological revolutions. Microarray is the latest invention in this ongoing technological revolution. This technology holds the promise to revolutionize the future of biology and medicine unlike any other technology that preceded it. Development of microarray technology has significantly changed the way questions about diseases and/or biological phenomena are addressed. This is because microarrays facilitate monitoring the expression of thousands of genes or proteins in a single experiment. This enormous power of microarrays has enabled scientists to monitor thousands of genes and their products in a given living organism in one experiment, and to understand how these genes function in an orchestrated manner. Obtaining such a global view of life at the molecular level was impossible using conventional molecular biological techniques. However, despite all the progress made in developing this technology, microarray is yet to reach a point where all data are obtained, analyzed, and shared in a standardized fashion. The present article is a brief overview of microarray technologies and their applications with an emphasis on DNA microarray.     

Coupled analysis of gene expression and chromosomal location

Microarray technology can be used to assess simultaneously global changes in expression of mRNA or genomic DNA copy number among thousands of genes in different biological states. In many cases, it is desirable to determine if altered patterns of gene expression correlate with chromosomal abnormalities or assess expression of genes that are contiguous in the genome. We describe a method, differential gene locus mapping (DIGMAP), which aligns the known chromosomal location of a gene to its expression value deduced by microarray analysis. The method partitions microarray data into subsets by chromosomal location for each gene interrogated by an array. Microarray data in an individual subset can then be clustered by physical location of genes at a subchromosomal level based upon ordered alignment in genome sequence. A graphical display is generated by representing each genomic locus with a colored cell that quantitatively reflects its differential expression value. The clustered patterns can be viewed and compared based on their expression signatures as defined by differential values between control and experimental samples. In this study, DIGMAP was tested using previously published studies of breast cancer analyzed by comparative genomic hybridization (CGH) and prostate cancer gene expression profiles assessed by cDNA microarray experiments. Analysis of the breast cancer CGH data demonstrated the ability of DIGMAP to deduce gene amplifications and deletions. Application of the DIGMAP method to the prostate data revealed several carcinoma-related loci, including one at 16q13 with marked differential expression encompassing 19 known genes including 9 encoding metallothionein proteins. We conclude that DIGMAP is a powerful computational tool enabling the coupled analysis of microarray data with genome location.     

肿瘤微阵列数据统计分析概述

微阵列技术已广泛应用于生物学和医学研究领域,如肿瘤的诊断和分型、预测和治疗,理解肿瘤的发生机制、生物学通路和基因网络。统计学方法在这一科学挑战中的地位至关重要。我们综述了微阵列实验数据分析的统计学方法最新发展,主要描述了微阵列数据的标准化、差异表达基因的统计学检验及微阵列技术在肿瘤治疗中的应用,重点介绍了时间序列微阵列数据分析方法和基因调控网络在肿瘤研究中的最新发展。     

用生物信息学方法处理基因芯片结果

基因芯片是一种能够同时检测大量基因在同一组织中表达情况的有力工具.利用前期工作筛选的2210个鼻咽癌差异表达基因和Biocarta信号通路资源库,构建了一个基于信号通路的基因相互作用网络.通过统计学分析,进一步筛选出一批对该基因相互作用网络具有重大影响的基因(特别是RAN、CEL、RELA).随后,采用RT-PCR方法检测候选基因在鼻咽癌活检组织中的表达,发现RAN和CEL基因在高达80%的鼻咽癌组织中高表达.进一步将网络分析结果和ArrayXPath软件分析的结果比较,共计有40%(32/80)基因结果吻合,这验证了网络分析方法的有效性和可行性.最终探索建立了新的分析基因芯片的方法.     

Accurate and robust gene selection for disease classification using a simple statistic

Discrimination of disease patients based on gene expression data is a crucial problem in clinical area. An important issue to solve this problem is to find a discriminative subset of genes from thousands of genes on a microarray or DNA chip. Aiming at finding informative genes for disease classification on microarray, we present a gene selection method based on the forward variable (gene) selection method (FSM) and show, using typical public microarray datasets, that our method can extract a small set of genes being crucial for discriminating different classes with a very high accuracy almost closed to perfect classification.     

Theoretically modeling microarray with the chemical equilibrium and thermodynamics

To date, the idea that microarray may shed the light on cellular processes by identifying groups of genes that appear to be co-expressed seems to remain a dream. This is partly because that there are some blank (meaning the knowledge is unavailable) or even erroneous areas in the fundamental theory in this field. This paper attempts to present the digest of microarray hybridization system with chemical thermodynamics, theoretically clarifying some misunderstandings and looking for answers to some critical questions around this technology, such as the mechanisms and conditions of quantitative measuring by hybridization reaction, the reasons of inconsistency of the data and the analysis results and the solutions, how to analyze the data, etc. A theoretical model for the next generation of microarray is proposed. We believe that this model is universal, laying the foundation for microarray technology from array design through the data analysis.     

Quantitative assessment of DNA microarrays--comparison with Northern blot analyses

DNA microarray is a powerful technology that provides the expression profile of thousands of genes. However, less attention has been paid to its quantitative aspect. In this study, we constructed a small-scale DNA microarray that contains 84 genes and characterized its quantitative aspect. Analyses with this microarray showed that 17 genes were induced, whereas 8 genes were suppressed at least twofold during the differentiation of mouse embryonic stem cells. When repeated with the same combination of fluorescent dyes for probe labeling, the microarray produced consistent data (correlation coefficient = 0.991). In contrast, data were less consistent when repeated with the reverse combination of dyes (correlation coefficient = 0.945). The effect of dye combination was particularly evident in several genes. Total RNA (15 microg) and poly(A) RNA (0.5 microg) showed comparable sensitivity and produced essentially identical data (correlation coefficient = 0.983). The sensitivity of the DNA microarrays was slightly inferior to that of Northern blot analyses. In most genes, data obtained with the two methods were consistent. However, in 4 of 46 genes compared, DNA microarrays failed to detect the expression changes that were revealed by Northern blot. These data demonstrated that DNA microarrays provide quantitative data comparable to Northern blot in general, but a few issues must be considered when analyzing data.     

Robust prostate cancer marker genes emerge from direct integration of inter-study microarray data

MOTIVATION: DNA microarray data analysis has been used previously to identify marker genes which discriminate cancer from normal samples. However, due to the limited sample size of each study, there are few common markers among different studies of the same cancer. With the rapid accumulation of microarray data, it is of great interest to integrate inter-study microarray data to increase sample size, which could lead to the discovery of more reliable markers. RESULTS: We present a novel, simple method of integrating different microarray datasets to identify marker genes and apply the method to prostate cancer datasets. In this study, by applying a new statistical method, referred to as the top-scoring pair (TSP) classifier, we have identified a pair of robust marker genes (HPN and STAT6) by integrating microarray datasets from three different prostate cancer studies. Cross-platform validation shows that the TSP classifier built from the marker gene pair, which simply compares relative expression values, achieves high accuracy, sensitivity and specificity on independent datasets generated using various array platforms. Our findings suggest a new model for the discovery of marker genes from accumulated microarray data and demonstrate how the great wealth of microarray data can be exploited to increase the power of statistical analysis. CONTACT: leixu@jhu.edu.     

Positional artifacts in microarrays: experimental verification and construction of COP, an automated detection tool

Microarray technology is currently one of the most widely-used technologies in biology. Many studies focus on inferring the function of an unknown gene from its co-expressed genes. Here, we are able to show that there are two types of positional artifacts in microarray data introducing spurious correlations between genes. First, we find that genes that are close on the microarray chips tend to have higher correlations between their expression profiles. We call this the 'chip artifact'. Our calculations suggest that the carry-over during the printing process is one of the major sources of this type of artifact, which is later confirmed by our experiments. Based on our experiments, the measured intensity of a microarray spot contains 0.1% (for fully-hybridized spots) to 93% (for un-hybridized ones) of noise resulting from this artifact. Secondly, we, for the first time, show that genes that are close on the microtiter plates in microarray experiments also tend to have higher correlations. We call this the 'plate artifact'. Both types of artifacts exist with different severity in all cDNA microarray experiments that we analyzed. Therefore, we develop an automated web tool-COP (COrrelations by Positional artifacts) to detect these artifacts in microarray experiments. COP has been integrated with the microarray data normalization tool, ExpressYourself, which is available at http://bioinfo.mbb.yale.edu/ExpressYourself/. Together, the two can eliminate most of the common noises in microarray data.     

A Novel Strategy for Gene Selection of Microarray Data Based on Gene-to-Class Sensitivity Information

To obtain predictive genes with lower redundancy and better interpretability, a hybrid gene selection method encoding prior information is proposed in this paper. To begin with, the prior information referred to as gene-to-class sensitivity (GCS) of all genes from microarray data is exploited by a single hidden layered feedforward neural network (SLFN). Then, to select more representative and lower redundant genes, all genes are grouped into some clusters by K-means method, and some low sensitive genes are filtered out according to their GCS values. Finally, a modified binary particle swarm optimization (BPSO) encoding the GCS information is proposed to perform further gene selection from the remainder genes. For considering the GCS information, the proposed method selects those genes highly correlated to sample classes. Thus, the low redundant gene subsets obtained by the proposed method also contribute to improve classification accuracy on microarray data. The experiments results on some open microarray data verify the effectiveness and efficiency of the proposed approach.     

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     

A robust method for estimating gene expression states using Affymetrix microarray probe level data

Background  

Microarray technology is a high-throughput method for measuring the expression levels of thousand of genes simultaneously. The observed intensities combine a non-specific binding, which is a major disadvantage with microarray data. The Affymetrix GeneChip assigned a mismatch (MM) probe with the intention of measuring non-specific binding, but various opinions exist regarding usefulness of MM measures. It should be noted that not all observed intensities are associated with expressed genes and many of those are associated with unexpressed genes, of which measured values express mere noise due to non-specific binding, cross-hybridization, or stray signals. The implicit assumption that all genes are expressed leads to poor performance of microarray data analyses. We assume two functional states of a gene - expressed or unexpressed - and propose a robust method to estimate gene expression states using an order relationship between PM and MM measures.