Integrated analysis of DNA copy number and gene expression microarray data using gene sets

Background  

Genes that play an important role in tumorigenesis are expected to show association between DNA copy number and RNA expression. Optimal power to find such associations can only be achieved if analysing copy number and gene expression jointly. Furthermore, some copy number changes extend over larger chromosomal regions affecting the expression levels of multiple resident genes.     

Technology for microarray analysis of gene expression

The past year has demonstrated the versatility of microarrays for the analysis of whole model-organism genomes and has seen the development of chips to measure the expression of 40,000 human genes. Microarray technology has also become considerably more robust and sensitive. Technology enhancements include the use of noncontact printing methods, improved 2-color sample preparation, and statistically based software for data analysis.     

Robust cluster analysis of microarray gene expression data with the number of clusters determined biologically

MOTIVATION: The success of each method of cluster analysis depends on how well its underlying model describes the patterns of expression. Outlier-resistant and distribution-insensitive clustering of genes are robust against violations of model assumptions. RESULTS: A measure of dissimilarity that combines advantages of the Euclidean distance and the correlation coefficient is introduced. The measure can be made robust using a rank order correlation coefficient. A robust graphical method of summarizing the results of cluster analysis and a biological method of determining the number of clusters are also presented. These methods are applied to a public data set, showing that rank-based methods perform better than log-based methods. AVAILABILITY: Software is available from http://www.davidbickel.com.     

Redefinition of Affymetrix probe sets by sequence overlap with cDNA microarray probes reduces cross-platform inconsistencies in cancer-associated gene expression measurements

Background  

Comparison of data produced on different microarray platforms often shows surprising discordance. It is not clear whether this discrepancy is caused by noisy data or by improper probe matching between platforms. We investigated whether the significant level of inconsistency between results produced by alternative gene expression microarray platforms could be reduced by stringent sequence matching of microarray probes. We mapped the short oligo probes of the Affymetrix platform onto cDNA clones of the Stanford microarray platform. Affymetrix probes were reassigned to redefined probe sets if they mapped to the same cDNA clone sequence, regardless of the original manufacturer-defined grouping. The NCI-60 gene expression profiles produced by Affymetrix HuFL platform were recalculated using these redefined probe sets and compared to previously published cDNA measurements of the same panel of RNA samples.     

Statistical design and the analysis of gene expression microarray data

Gene expression microarrays are an innovative technology with enormous promise to help geneticists explore and understand the genome. Although the potential of this technology has been clearly demonstrated, many important and interesting statistical questions persist. We relate certain features of microarrays to other kinds of experimental data and argue that classical statistical techniques are appropriate and useful. We advocate greater attention to experimental design issues and a more prominent role for the ideas of statistical inference in microarray studies.     

Subtraction-coupled custom microarray analysis for gene discovery and gene expression studies in the CNS

The revolution in our knowledge about the genomes of organisms gives rise to the question, what do we do with this information? The development of techniques allowing high throughput analysis of RNA and protein expression, such as cDNA microarrays, provide for genome-wide analysis of gene expression. These analyses will help bridge the gap between systems and molecular neuroscience. This review discusses the advantages of using a subtractive hybridization technique, such as a representational difference analysis, to generate a custom cDNA microarray enriched for genes relevant to investigating complex, heterogeneous tissues such as those involved in the chemical senses. Real and hypothetical examples of these experiments are discussed. Benefits of this approach over traditional microarray techniques include having a more relevant clone set, the potential for gene discovery and the creation of a new tool to investigate similar systems. Potential pitfalls may include PCR artifacts and the need for sequencing. However, these disadvantages can be overcome so that the coupling of subtraction techniques to microarray screening can be a fruitful approach to a variety of experimental systems.     

Detection of DNA copy number abnormality by microarray expression analysis

Gene copy-number abnormalities (CNAs) are characteristic of solid tumors and are found in association with developmental abnormalities and/or mental retardation. The ultimate impact of CNAs is exerted by the altered expression of encoded genes. We have utilized high-density oligonucleotide arrays from Affymetrix to identify DNA CNAs via their impact on mRNA expression levels. In these studies, we have used three different trisomic cell lines (trisomy 9, trisomy 18, trisomy 21) as models of CNAs and have compared mRNA expression in those trisomic cells with that observed in diploid cell lines of matched tissue origin. Our data clearly show that genes from CNA chromosome regions are substantially over-represented (P<0.000001 by chi-square analysis) in the differentially expressed subset from comparisons of all three trisomic cell lines with normal matching cells. In addition, we have been able to detect the origin of the duplication by a statistical scan for over-expressed genes. These data show that microarray detection of differential mRNA expression can be used to identify significant DNA CNAs.     

Optimization based tumor classification from microarray gene expression data

Background

An important use of data obtained from microarray measurements is the classification of tumor types with respect to genes that are either up or down regulated in specific cancer types. A number of algorithms have been proposed to obtain such classifications. These algorithms usually require parameter optimization to obtain accurate results depending on the type of data. Additionally, it is highly critical to find an optimal set of markers among those up or down regulated genes that can be clinically utilized to build assays for the diagnosis or to follow progression of specific cancer types. In this paper, we employ a mixed integer programming based classification algorithm named hyper-box enclosure method (HBE) for the classification of some cancer types with a minimal set of predictor genes. This optimization based method which is a user friendly and efficient classifier may allow the clinicians to diagnose and follow progression of certain cancer types.

Methodology/Principal Findings

We apply HBE algorithm to some well known data sets such as leukemia, prostate cancer, diffuse large B-cell lymphoma (DLBCL), small round blue cell tumors (SRBCT) to find some predictor genes that can be utilized for diagnosis and prognosis in a robust manner with a high accuracy. Our approach does not require any modification or parameter optimization for each data set. Additionally, information gain attribute evaluator, relief attribute evaluator and correlation-based feature selection methods are employed for the gene selection. The results are compared with those from other studies and biological roles of selected genes in corresponding cancer type are described.

Conclusions/Significance

The performance of our algorithm overall was better than the other algorithms reported in the literature and classifiers found in WEKA data-mining package. Since it does not require a parameter optimization and it performs consistently very high prediction rate on different type of data sets, HBE method is an effective and consistent tool for cancer type prediction with a small number of gene markers.     

Weighted analysis of microarray gene expression using maximum-likelihood

MOTIVATION: The numerical values of gene expression measured using microarrays are usually presented to the biological end-user as summary statistics of spot pixel data, such as the spot mean, median and mode. Much of the subsequent data analysis reported in the literature, however, uses only one of these spot statistics. This results in sub-optimal estimates of gene expression levels and a need for improvement in quantitative spot variation surveillance. RESULTS: This paper develops a maximum-likelihood method for estimating gene expression using spot mean, variance and pixel number values available from typical microarray scanners. It employs a hierarchical model of variation between and within microarray spots. The hierarchical maximum-likelihood estimate (MLE) is shown to be a more efficient estimator of the mean than the 'conventional' estimate using solely the spot mean values (i.e. without spot variance data). Furthermore, under the assumptions of our model, the spot mean and spot variance are shown to be sufficient statistics that do not require the use of all pixel data.The hierarchical MLE method is applied to data from both Monte Carlo (MC) simulations and a two-channel dye-swapped spotted microarray experiment. The MC simulations show that the hierarchical MLE method leads to improved detection of differential gene expression particularly when 'outlier' spots are present on the arrays. Compared with the conventional method, the MLE method applied to data from the microarray experiment leads to an increase in the number of differentially expressed genes detected for low cut-off P-values of interest.     

Differential analysis of DNA microarray gene expression data

Here, we review briefly the sources of experimental and biological variance that affect the interpretation of high-dimensional DNA microarray experiments. We discuss methods using a regularized t-test based on a Bayesian statistical framework that allow the identification of differentially regulated genes with a higher level of confidence than a simple t-test when only a few experimental replicates are available. We also describe a computational method for calculating the global false-positive and false-negative levels inherent in a DNA microarray data set. This method provides a probability of differential expression for each gene based on experiment-wide false-positive and -negative levels driven by experimental error and biological variance.     

A rapid and inexpensive labeling method for microarray gene expression analysis

Background

Polyethyleneimine (PEI), a cationic polymer, is one of the successful and widely used vectors for non-viral gene transfection in vitro. However, its in vivo application was greatly limited due to its high cytotoxicity and short duration of gene expression. To improve its biocompatibility and transfection efficiency, PEI has been modified with PEG, folic acid, and chloroquine in order to improve biocompatibility and enhance targeting.

Results

Poly(ε-caprolactone)-Pluronic-Poly(ε-caprolactone) (PCFC) was synthesized by ring-opening polymerization, and PCFC-g-PEI was obtained by Michael addition reaction with GMA-PCFC-GMA and polyethyleneimine (PEI, 25 kD). The prepared PCFC-g-PEI was characterized by ¹H-NMR, SEC-MALLS. Meanwhile, DNA condensation, DNase I protection, the particle size and zeta potential of PCFC-g-PEI/DNA complexes were also determined. According to the results of flow cytometry and MTT assay, the synthesized PCFC-g-PEI, with considerable transfection efficiency, had obviously lower cytotoxicity against 293 T and A549 cell lines compared with that of PEI 25 kD.

Conclusion

The cytotoxicity and in vitro transfection study indicated that PCFC-g-PEI copolymer prepared in this paper was a novel gene delivery system with lower cytotoxicity and considerable transfection efficiency compared with commercial PEI (25 kD).     

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.     

A rapid and inexpensive labeling method for microarray gene expression analysis

Background  

Global gene expression profiling by DNA microarrays is an invaluable tool in biological research. However, existing labeling methods are time consuming and costly and therefore often limit the scale of microarray experiments and sample throughput. Here we introduce a new, fast, inexpensive method for direct random-primed fluorescent labeling of eukaryotic cDNA for gene expression analysis and compare the results obtained on the NimbleGen microarray platform with two other widely-used labeling methods, namely the NimbleGen-recommended double-stranded cDNA protocol and the indirect (aminoallyl) method.     

Quantitative trait associated microarray gene expression data analysis

Selection on phenotypes may cause genetic change. To understand the relationship between phenotype and gene expression from an evolutionary viewpoint, it is important to study the concordance between gene expression and profiles of phenotypes. In this study, we use a novel method of clustering to identify genes whose expression profiles are related to a quantitative phenotype. Cluster analysis of gene expression data aims at classifying genes into several different groups based on the similarity of their expression profiles across multiple conditions. The hope is that genes that are classified into the same clusters may share underlying regulatory elements or may be a part of the same metabolic pathways. Current methods for examining the association between phenotype and gene expression are limited to linear association measured by the correlation between individual gene expression values and phenotype. Genes may be associated with the phenotype in a nonlinear fashion. In addition, groups of genes that share a particular pattern in their relationship to phenotype may be of evolutionary interest. In this study, we develop a method to group genes based on orthogonal polynomials under a multivariate Gaussian mixture model. The effect of each expressed gene on the phenotype is partitioned into a cluster mean and a random deviation from the mean. Genes can also be clustered based on a time series. Parameters are estimated using the expectation-maximization algorithm and implemented in SAS. The method is verified with simulated data and demonstrated with experimental data from 2 studies, one clusters with respect to severity of disease in Alzheimer's patients and another clusters data for a rat fracture healing study over time. We find significant evidence of nonlinear associations in both studies and successfully describe these patterns with our method. We give detailed instructions and provide a working program that allows others to directly implement this method in their own analyses.