Group additive regression models for genomic data analysis

One important problem in genomic research is to identify genomic features such as gene expression data or DNA single nucleotide polymorphisms (SNPs) that are related to clinical phenotypes. Often these genomic data can be naturally divided into biologically meaningful groups such as genes belonging to the same pathways or SNPs within genes. In this paper, we propose group additive regression models and a group gradient descent boosting procedure for identifying groups of genomic features that are related to clinical phenotypes. Our simulation results show that by dividing the variables into appropriate groups, we can obtain better identification of the group features that are related to the phenotypes. In addition, the prediction mean square errors are also smaller than the component-wise boosting procedure. We demonstrate the application of the methods to pathway-based analysis of microarray gene expression data of breast cancer. Results from analysis of a breast cancer microarray gene expression data set indicate that the pathways of metalloendopeptidases (MMPs) and MMP inhibitors, as well as cell proliferation, cell growth, and maintenance are important to breast cancer-specific survival.     

A Model-Based Joint Identification of Differentially Expressed Genes and Phenotype-Associated Genes

Over the last decade, many analytical methods and tools have been developed for microarray data. The detection of differentially expressed genes (DEGs) among different treatment groups is often a primary purpose of microarray data analysis. In addition, association studies investigating the relationship between genes and a phenotype of interest such as survival time are also popular in microarray data analysis. Phenotype association analysis provides a list of phenotype-associated genes (PAGs). However, it is sometimes necessary to identify genes that are both DEGs and PAGs. We consider the joint identification of DEGs and PAGs in microarray data analyses. The first approach we used was a naïve approach that detects DEGs and PAGs separately and then identifies the genes in an intersection of the list of PAGs and DEGs. The second approach we considered was a hierarchical approach that detects DEGs first and then chooses PAGs from among the DEGs or vice versa. In this study, we propose a new model-based approach for the joint identification of DEGs and PAGs. Unlike the previous two-step approaches, the proposed method identifies genes simultaneously that are DEGs and PAGs. This method uses standard regression models but adopts different null hypothesis from ordinary regression models, which allows us to perform joint identification in one-step. The proposed model-based methods were evaluated using experimental data and simulation studies. The proposed methods were used to analyze a microarray experiment in which the main interest lies in detecting genes that are both DEGs and PAGs, where DEGs are identified between two diet groups and PAGs are associated with four phenotypes reflecting the expression of leptin, adiponectin, insulin-like growth factor 1, and insulin. Model-based approaches provided a larger number of genes, which are both DEGs and PAGs, than other methods. Simulation studies showed that they have more power than other methods. Through analysis of data from experimental microarrays and simulation studies, the proposed model-based approach was shown to provide a more powerful result than the naïve approach and the hierarchical approach. Since our approach is model-based, it is very flexible and can easily handle different types of covariates.     

A Nonparametric Likelihood Ratio Test to Identify Differentially Expressed Genes from Microarray Data

Microarray experiments contribute significantly to the progress in disease treatment by enabling a precise and early diagnosis. One of the major objectives of microarray experiments is to identify differentially expressed genes under various conditions. The statistical methods currently used to analyse microarray data are inadequate, mainly due to the lack of understanding of the distribution of microarray data. We present a nonparametric likelihood ratio (NPLR) test to identify differentially expressed genes using microarray data. The NPLR test is highly robust against extreme values and does not assume the distribution of the parent population. Simulation studies show that the NPLR test is more powerful than some of the commonly used methods, such as the two-sample t-test, the Mann-Whitney U-test and significance analysis of microarrays (SAM). When applied to microarray data, we found that the NPLR test identifies more differentially expressed genes than its competitors. The asymptotic distribution of the NPLR test statistic and the p-value function is presented. The application of the NPLR method is shown, using both synthetic and real-life data. The biological significance of some of the genes detected only by the NPLR method is discussed.     

Molecular pattern discovery based on penalized matrix decomposition

A reliable and precise identification of the type of tumors is crucial to the effective treatment of cancer. With the rapid development of microarray technologies, tumor clustering based on gene expression data is becoming a powerful approach to cancer class discovery. In this paper, we apply the penalized matrix decomposition (PMD) to gene expression data to extract metasamples for clustering. The extracted metasamples capture the inherent structures of samples belong to the same class. At the same time, the PMD factors of a sample over the metasamples can be used as its class indicator in return. Compared with the conventional methods such as hierarchical clustering (HC), self-organizing maps (SOM), affinity propagation (AP) and nonnegative matrix factorization (NMF), the proposed method can identify the samples with complex classes. Moreover, the factor of PMD can be used as an index to determine the cluster number. The proposed method provides a reasonable explanation of the inconsistent classifications made by the conventional methods. In addition, it is able to discover the modules in gene expression data of conterminous developmental stages. Experiments on two representative problems show that the proposed PMD-based method is very promising to discover biological phenotypes.     

Differential expression of extracellular matrix proteins in senescent and young human fibroblasts: A comparative proteomics and microarray study

The extracellular matrix (ECM) provides an essential structural framework for cell attachment, proliferation, and differentiation, and undergoes progressive changes during senescence. To investigate changes in protein expression in the extracellular matrix between young and senescent fibroblasts, we compared proteomic data (LTQ-FT) with cDNA microarray results. The peptide counts from the proteomics analysis were used to evaluate the level of ECM protein expression by young cells and senescent cells, and ECM protein expression data were compared with the microarray data. After completing the comparative analysis, we grouped the genes into four categories. Class I included genes with increased expression levels in both analyses, while class IV contained genes with reduced expression in both analyses. Class II and Class III contained genes with an inconsistent expression pattern. Finally, we validated the comparative analysis results by examining the expression level of the specific gene from each category using Western blot analysis and semiquantitative RT-PCR. Our results demonstrate that comparative analysis can be used to identify differentially expressed genes.     

On-chip hybridization kinetics for optimization of gene expression experiments

DNA microarray technology is a powerful tool for getting an overview of gene expression in biological samples. Although the successful use of microarray-based expression analysis was demonstrated in a number of applications, the main problem with this approach is the fact that expression levels deduced from hybridization experiments do not necessarily correlate with RNA concentrations. Moreover oligonucleotide probes corresponding to the same gene can give different hybridization signals. Apart from cross-hybridizations and differential splicing, this could be due to secondary structures of probes or targets. In addition, for low-copy genes, hybridization equilibrium may be reached after hybridization times much longer than the one commonly used (overnight, i.e., 15 h). Thus, hybridization signals could depend on kinetic properties of the probe, which may vary between different oligonucleotide probes immobilized on the same microarray. To validate this hypothesis, on-chip hybridization kinetics and duplex thermostability analysis were performed using oligonucleotide microarrays containing 50-mer probes corresponding to 10 mouse genes. We demonstrate that differences in hybridization kinetics between the probes exist and can influence the interpretation of expression data. In addition, we show that using on-chip hybridization kinetics, quantification of targets is feasible using calibration curves.     

Pop’s Pipes: poplar gene expression data analysis pipelines

We developed multiple gene expression pipelines and assembled them into a web-based tool called Pop’s Pipes to facilitate preprocessing and analysis of substantial poplar gene expression data. The input data can be spatiotemporal microarray and RNA-seq data from comparable tissues, time points, or treatment-vs-control conditions. Pop’s Pipes can be used to identify differentially expressed genes between one or multiple paired tissues, time points, or treatment-vs-control conditions in a single in silico analysis. The differentially expressed genes (DEGs) obtained for each comparison will be automatically analyzed by Pop’s Pipes for identifying significantly enriched gene ontologies and interpro protein domains. Also, significantly changed metabolic pathways across all input data sets will be identified. We also integrated a pipeline into Pop's Pipes for constructing any of three type gene ontology trees when a short list of gene ontologies from biological processes, molecular functions, or cellular components is used as an input. The resulting information from Pop’s Pipes enables scrutiny to create spatiotemporal models and hypotheses to understand how poplar develops and functions. Pop’s Pipes can analyze a microarray or RNA-seq data set with 10 time points in 4–10 h, with each time point containing three replicates of treatments and three controls. Such a data set usually takes a bioinformatician a few months to a year to analyze. Pop’s Pipes can thus save users tremendous amounts of research time when large numbers of comparative data need to be analyzed.     

Microarray analysis of antimicrobial resistance genes in Salmonella enterica from preharvest poultry environment

Aims:  To detect antimicrobial resistance genes in Salmonella isolates from turkey flocks using the microarray technology.
Methods and Results:  A 775 gene probe oligonucleotide microarray was used to detect antimicrobial resistance genes in 34 isolates. All tetracycline-resistant Salmonella harboured tet(A) , tet(C) or tet(R) , with the exception of one Salmonella serotype Heidelberg isolate. The sul1 gene was detected in 11 of 16 sulfisoxazole-resistant isolates. The aadA , aadA1 , aadA2 , strA or strB genes were found in aminoglycoside-resistant isolates of Salm. Heidelberg, Salmonella serotype Senftenberg and untypeable Salmonella . The prevalence of mobile genetic elements, such as class I integron and transposon genes, in drug-resistant Salmonella isolates suggested that these elements may contribute to the dissemination of antimicrobial resistance genes in the preharvest poultry environment. Hierarchical clustering analysis demonstrated a close relationship between drug-resistant phenotypes and the corresponding antimicrobial resistance gene profiles.
Conclusions:  Salmonella serotypes isolated from the poultry environment carry multiple genes that can render them resistant to several antimicrobials used in poultry and humans.
Significance and Impact of the Study:  Multiple antimicrobial resistance genes in environmental Salmonella isolates could be identified efficiently by microarray analysis. Hierarchical clustering analysis of the data was also found to be a useful tool for analysing emerging patterns of drug resistance.     

The Baumgartner-Weiss-Schindler test for the detection of differentially expressed genes in replicated microarray experiments

MOTIVATION: An important application of microarray experiments is to identify differentially expressed genes. Because microarray data are often not distributed according to a normal distribution nonparametric methods were suggested for their statistical analysis. Here, the Baumgartner-Weiss-Schindler test, a novel and powerful test based on ranks, is investigated and compared with the parametric t-test as well as with two other nonparametric tests (Wilcoxon rank sum test, Fisher-Pitman permutation test) recently recommended for the analysis of gene expression data. RESULTS: Simulation studies show that an exact permutation test based on the Baumgartner-Weiss-Schindler statistic B is preferable to the other three tests. It is less conservative than the Wilcoxon test and more powerful, in particular in case of asymmetric or heavily tailed distributions. When the underlying distribution is symmetric the differences in power between the tests are relatively small. Thus, the Baumgartner-Weiss-Schindler is recommended for the usual situation that the underlying distribution is a priori unknown. AVAILABILITY: SAS code available on request from the authors.     

Nonparametric pathway-based regression models for analysis of genomic data

High-throughout genomic data provide an opportunity for identifying pathways and genes that are related to various clinical phenotypes. Besides these genomic data, another valuable source of data is the biological knowledge about genes and pathways that might be related to the phenotypes of many complex diseases. Databases of such knowledge are often called the metadata. In microarray data analysis, such metadata are currently explored in post hoc ways by gene set enrichment analysis but have hardly been utilized in the modeling step. We propose to develop and evaluate a pathway-based gradient descent boosting procedure for nonparametric pathways-based regression (NPR) analysis to efficiently integrate genomic data and metadata. Such NPR models consider multiple pathways simultaneously and allow complex interactions among genes within the pathways and can be applied to identify pathways and genes that are related to variations of the phenotypes. These methods also provide an alternative to mediating the problem of a large number of potential interactions by limiting analysis to biologically plausible interactions between genes in related pathways. Our simulation studies indicate that the proposed boosting procedure can indeed identify relevant pathways. Application to a gene expression data set on breast cancer distant metastasis identified that Wnt, apoptosis, and cell cycle-regulated pathways are more likely related to the risk of distant metastasis among lymph-node-negative breast cancer patients. Results from analysis of other two breast cancer gene expression data sets indicate that the pathways of Metalloendopeptidases (MMPs) and MMP inhibitors, as well as cell proliferation, cell growth, and maintenance are important to breast cancer relapse and survival. We also observed that by incorporating the pathway information, we achieved better prediction for cancer recurrence.     

FiGS: a filter-based gene selection workbench for microarray data

Background  

The selection of genes that discriminate disease classes from microarray data is widely used for the identification of diagnostic biomarkers. Although various gene selection methods are currently available and some of them have shown excellent performance, no single method can retain the best performance for all types of microarray datasets. It is desirable to use a comparative approach to find the best gene selection result after rigorous test of different methodological strategies for a given microarray dataset.     

Pathway analysis of microarray data via regression.

Pathway analysis of microarray data evaluates gene expression profiles of a priori defined biological pathways in association with a phenotype of interest. We propose a unified pathway-analysis method that can be used for diverse phenotypes including binary, multiclass, continuous, count, rate, and censored survival phenotypes. The proposed method also allows covariate adjustments and correlation in the phenotype variable that is encountered in longitudinal, cluster-sampled, and paired designs. These are accomplished by combining the regression-based test statistic for each individual gene in a pathway of interest into a pathway-level test statistic. Applications of the proposed method are illustrated with two real pathway-analysis examples: one evaluating relapse-associated gene expression involving a matched-pair binary phenotype in children with acute lymphoblastic leukemia; and the other investigating gene expression in breast cancer tissues in relation to patients' survival (a censored survival phenotype). Implementations for various phenotypes are available in R. Additionally, an Excel Add-in for a user-friendly interface is currently being developed.     

Clustering methods for microarray gene expression data

Within the field of genomics, microarray technologies have become a powerful technique for simultaneously monitoring the expression patterns of thousands of genes under different sets of conditions. A main task now is to propose analytical methods to identify groups of genes that manifest similar expression patterns and are activated by similar conditions. The corresponding analysis problem is to cluster multi-condition gene expression data. The purpose of this paper is to present a general view of clustering techniques used in microarray gene expression data analysis.     

Gene set enrichment analysis for multiple continuous phenotypes

Background

Gene set analysis (GSA) methods test the association of sets of genes with phenotypes in gene expression microarray studies. While GSA methods on a single binary or categorical phenotype abounds, little attention has been paid to the case of a continuous phenotype, and there is no method to accommodate correlated multiple continuous phenotypes.

Result

We propose here an extension of the linear combination test (LCT) to its new version for multiple continuous phenotypes, incorporating correlations among gene expressions of functionally related gene sets, as well as correlations among multiple phenotypes. Further, we extend our new method to its nonlinear version, referred as nonlinear combination test (NLCT), to test potential nonlinear association of gene sets with multiple phenotypes. Simulation study and a real microarray example demonstrate the practical aspects of the proposed methods.

Conclusion

The proposed approaches are effective in controlling type I errors and powerful in testing associations between gene-sets and multiple continuous phenotypes. They are both computationally effective. Naively (univariately) analyzing a group of multiple correlated phenotypes could be dangerous. R-codes to perform LCT and NLCT for multiple continuous phenotypes are available at http://www.ualberta.ca/~yyasui/homepage.html.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2105-15-260) contains supplementary material, which is available to authorized users.     

Statistical tests for differential expression in cDNA microarray experiments

Extracting biological information from microarray data requires appropriate statistical methods. The simplest statistical method for detecting differential expression is the t test, which can be used to compare two conditions when there is replication of samples. With more than two conditions, analysis of variance (ANOVA) can be used, and the mixed ANOVA model is a general and powerful approach for microarray experiments with multiple factors and/or several sources of variation.     

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     

Identifying Gene Interaction Enrichment for Gene Expression Data

Gene set analysis allows the inclusion of knowledge from established gene sets, such as gene pathways, and potentially improves the power of detecting differentially expressed genes. However, conventional methods of gene set analysis focus on gene marginal effects in a gene set, and ignore gene interactions which may contribute to complex human diseases. In this study, we propose a method of gene interaction enrichment analysis, which incorporates knowledge of predefined gene sets (e.g. gene pathways) to identify enriched gene interaction effects on a phenotype of interest. In our proposed method, we also discuss the reduction of irrelevant genes and the extraction of a core set of gene interactions for an identified gene set, which contribute to the statistical variation of a phenotype of interest. The utility of our method is demonstrated through analyses on two publicly available microarray datasets. The results show that our method can identify gene sets that show strong gene interaction enrichments. The enriched gene interactions identified by our method may provide clues to new gene regulation mechanisms related to the studied phenotypes. In summary, our method offers a powerful tool for researchers to exhaustively examine the large numbers of gene interactions associated with complex human diseases, and can be a useful complement to classical gene set analyses which only considers single genes in a gene set.