Dataset Information

EMDomics: a robust and powerful method for the identification of genes differentially expressed between heterogeneous classes.

ABSTRACT:

Motivation

A major goal of biomedical research is to identify molecular features associated with a biological or clinical class of interest. Differential expression analysis has long been used for this purpose; however, conventional methods perform poorly when applied to data with high within class heterogeneity.

Results

To address this challenge, we developed EMDomics, a new method that uses the Earth mover's distance to measure the overall difference between the distributions of a gene's expression in two classes of samples and uses permutations to obtain q-values for each gene. We applied EMDomics to the challenging problem of identifying genes associated with drug resistance in ovarian cancer. We also used simulated data to evaluate the performance of EMDomics, in terms of sensitivity and specificity for identifying differentially expressed gene in classes with high within class heterogeneity. In both the simulated and real biological data, EMDomics outperformed competing approaches for the identification of differentially expressed genes, and EMDomics was significantly more powerful than conventional methods for the identification of drug resistance-associated gene sets. EMDomics represents a new approach for the identification of genes differentially expressed between heterogeneous classes and has utility in a wide range of complex biomedical conditions in which sample classes show within class heterogeneity.

Availability and implementation

The R package is available at http://www.bioconductor.org/packages/release/bioc/html/EMDomics.html.

SUBMITTER: Nabavi S

PROVIDER: S-EPMC4743632 | biostudies-literature | 2016 Feb

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

EMDomics: a robust and powerful method for the identification of genes differentially expressed between heterogeneous classes.

Nabavi Sheida S Schmolze Daniel D Maitituoheti Mayinuer M Malladi Sadhika S Beck Andrew H AH

Bioinformatics (Oxford, England) 20151029 4

<h4>Motivation</h4>A major goal of biomedical research is to identify molecular features associated with a biological or clinical class of interest. Differential expression analysis has long been used for this purpose; however, conventional methods perform poorly when applied to data with high within class heterogeneity.<h4>Results</h4>To address this challenge, we developed EMDomics, a new method that uses the Earth mover's distance to measure the overall difference between the distributions of ...[more]

PMID: 26515818

Similar Datasets

Project description:BackgroundMicroarrays have been widely used for the analysis of gene expression and several commercial platforms are available. The combined use of multiple platforms can overcome the inherent biases of each approach, and may represent an alternative that is complementary to RT-PCR for identification of the more robust changes in gene expression profiles. In this paper, we combined statistical and functional analysis for the cross platform validation of two oligonucleotide-based technologies, Affymetrix (AFFX) and Applied Biosystems (ABI), and for the identification of differentially expressed genes.ResultsIn this study, we analysed differentially expressed genes after treatment of an ovarian carcinoma cell line with a cell cycle inhibitor. Treated versus control RNA was analysed for expression of 16425 genes represented on both platforms. We assessed reproducibility between replicates for each platform using CAT plots, and we found it high for both, with better scores for AFFX. We then applied integrative correlation analysis to assess reproducibility of gene expression patterns across studies, bypassing the need for normalizing expression measurements across platforms. We identified 930 genes as differentially expressed on AFFX and 908 on ABI, with approximately 80% common to both platforms. Despite the different absolute values, the range of intensities of the differentially expressed genes detected by each platform was similar. ABI showed a slightly higher dynamic range in FC values, which might be associated with its detection system. 62/66 genes identified as differentially expressed by Microarray were confirmed by RT-PCR.ConclusionIn this study we present a cross-platform validation of two oligonucleotide-based technologies, AFFX and ABI. We found good reproducibility between replicates, and showed that both platforms can be used to select differentially expressed genes with substantial agreement. Pathway analysis of the affected functions identified themes well in agreement with those expected for a cell cycle inhibitor, suggesting that this procedure is appropriate to facilitate the identification of biologically relevant signatures associated with compound treatment. The high rate of confirmation found for both common and platform-specific genes suggests that the combination of platforms may overcome biases related to probe design and technical features, thereby accelerating the identification of trustworthy differentially expressed genes.

Project description:Osteocytes represent the most abundant cellular component of mammalian bones with important functions in bone mass maintenance and remodeling. To elucidate the differential gene expression between osteoblasts and osteocytes we completed a comprehensive analysis of their gene profiles. Selective identification of these two mature populations was achieved by utilization of visual markers of bone lineage cells. We have utilized dual GFP reporter mice in which osteocytes are expressing GFP (topaz) directed by the DMP1 promoter, while osteoblasts are identified by expression of GFP (cyan) driven by 2.3 kb of the Col1a1 promoter. Histological analysis of 7-day-old neonatal calvaria confirmed the expression pattern of DMP1GFP in osteocytes and Col2.3 in osteoblasts and osteocytes. To isolate distinct populations of cells we utilized fluorescent activated cell sorting (FACS). Cell suspensions were subjected to RNA extraction, in vitro transcription and labeling of cDNA and gene expression was analyzed using the Illumina WG-6v1 BeadChip. Following normalization of raw data from four biological replicates, 3444 genes were called present in all three sorted cell populations: GFP negative, Col2.3cyan(+) (osteoblasts), and DMP1topaz(+) (preosteocytes and osteocytes). We present the genes that showed in excess of a 2-fold change for gene expression between DMP1topaz(+) and Col2.3cyan(+) cells. The selected genes were classified and grouped according to their associated gene ontology terms. Genes clustered to osteogenesis and skeletal development such as Bmp4, Bmp8a, Dmp1, Enpp1, Phex and Ank were highly expressed in DMP1topaz(+)cells. Most of the genes encoding extracellular matrix components and secreted proteins had lower expression in DMP1topaz(+) cells, while most of the genes encoding plasma membrane proteins were increased. Interestingly a large number of genes associated with muscle development and function and with neuronal phenotype were increased in DMP1topaz(+) cells, indicating some new aspects of osteocyte biology. Although a large number of genes differentially expressed in DMP1topaz(+) and Col2.3cyan(+) cells in our study have already been assigned to bone development and physiology, for most of them we still lack any substantial data. Therefore, isolation of osteocyte and osteoblast cell populations and their subsequent microarray analysis allowed us to identify a number or genes and pathways with potential roles in regulation of bone mass.

Project description:Soybean is a major source of protein and oil and a primary feedstock for biodiesel production. Research on soybean seed composition and yield has revealed that protein, oil and yield are controlled quantitatively and quantitative trait loci (QTL) have been identified for each of these traits. However, very limited information is available regarding the genetic mechanisms controlling seed composition and yield. To help address this deficiency, we used Affymetrix Soybean GeneChips® to identify genes that are differentially expressed between developing seeds of the Minsoy and Archer soybean cultivars, which differ in seed weight, yield, protein content and oil content. A total of 700 probe sets were found to be expressed at significantly different (defined as having an adjusted p-value below or equal to 0.05 and an at least 2-fold difference) levels between the two cultivars at one or more of the three developmental stages and in at least one of the two years assayed. Comparison of data from soybeans collected in two different years revealed that 97 probe sets were expressed at significantly different levels in both years. Functional annotations were assigned to 78% of these 97 probe sets based on the SoyBase Affymetrix™ GeneChip® Soybean Genome Array Annotation. Genes involved in receptor binding/activity and protein binding are overrepresented among the group of 97 probe sets that were differentially expressed in both years assayed. Probe sets involved in growth/development, signal transduction, transcription, defense/stress response and protein and lipid metabolism were also identified among the 97 probe sets and their possible implications in the regulation of agronomic traits are discussed. As the Minsoy and Archer soybean cultivars differ with respect to seed size, yield, protein content and lipid content, some of the differentially expressed probe sets identified in this study may thus play important roles in controlling these traits. Others of these probe sets may be involved in regulation of general seed development or metabolism. All microarray data and expression values after GCRMA are available at the Gene Expression Omnibus (GEO) at NCBI (http://www.ncbi.nlm.nih.gov/geo), under accession number GSE21598.

Project description:Myogenesis is a central step in prenatal myofiber formation, postnatal myofiber hypertrophy, and muscle damage repair in adulthood. RNA-Seq technology has greatly helped reveal the molecular mechanism of myogenesis, but batch effects in different experiments inevitably lead to misinterpretation of differentially expressed genes (DEGs). We previously applied the robust rank aggregation (RRA) method to effectively circumvent batch effects across multiple RNA-Seq datasets from 3T3-L1 cells. Here, we also used the RRA method to integrate nine RNA-Seq datasets from C2C12 cells and obtained 3140 robust DEGs between myoblasts and myotubes, which were then validated with array expression profiles and H3K27ac signals. The upregulated robust DEGs were highly enriched in gene ontology (GO) terms related to muscle cell differentiation and development. Considering that the cooperative binding of transcription factors (TFs) to enhancers to regulate downstream gene expression is a classical epigenetic mechanism, differentially expressed TFs (DETFs) were screened, and potential novel myogenic factors (MAF, BCL6, and ESR1) with high connection degree in protein-protein interaction (PPI) network were presented. Moreover, KLF5 cooperatively binds with the three key myogenic factors (MYOD, MYOG, and MEF2D) in C2C12 cells. Motif analysis speculates that the binding of MYOD and MYOG is KLF5-independent, while MEF2D is KLF5-dependent. It was revealed that KLF5-binding sites could be exploited to filter redundant MYOD-, MYOG-, and MEF2D-binding sites to focus on key enhancers for myogenesis. Further functional annotation of KLF5-binding sites suggested that KLF5 may regulate myogenesis through the PI3K-AKt signaling pathway, Rap1 signaling pathway, and the Hippo signaling pathway. In general, our study provides a wealth of untapped candidate targets for myogenesis and contributes new insights into the core regulatory mechanisms of myogenesis relying on KLF5-binding signal.

Project description:Microarrays have been widely used for the analysis of gene expression and several commercial platforms are available. The combined use of multiple platforms can overcome the inherent biases of each approach, and may represent an alternative that is complementary to RT-PCR for identification of the more robust changes in gene expression profiles. In this paper, we combined statistical and functional analysis for the cross platform validation of two oligonucleotide-based technologies, Affymetrix (AFFX) and Applied Biosystems (ABI), and for the identification of differentially expressed genes. In this study, we analysed differentially expressed genes after treatment of an ovarian carcinoma cell line with a cell cycle inhibitor. Treated versus control RNA was analysed for expression of 16425 genes represented on both platforms. We assessed reproducibility between replicates for each platform using CAT plots, and we found it high for both, with better scores for AFFX. We then applied integrative correlation analysis to assess reproducibility of gene expression patterns across studies, bypassing the need for normalizing expression measurements across platforms. We identified 930 genes as differentially expressed on AFFX and 908 on ABI, with ~80% common to both platforms. Despite the different absolute values, the range of intensities of the differentially expressed genes detected by each platform was similar. ABI showed a slightly higher dynamic range in FC values, which might be associated with its detection system. 62/66 genes identified as differentially expressed by Microarray were confirmed by RT-PCR. In this study we present a cross-platform validation of two oligonucleotide-based technologies, AFFX and ABI. We found good reproducibility between replicates, and showed that both platforms can be used to select differentially expressed genes with substantial agreement. Pathway analysis of the affected functions identified themes well in agreement with those expected for a cell cycle inhibitor, suggesting that this procedure is appropriate to facilitate the identification of biologically relevant signatures associated with compound treatment. The high rate of confirmation found for both common and platform-specific genes suggests that the combination of platforms may overcome biases related to probe design and technical features, thereby accelerating the identification of trustworthy differentially expressed genes. accepted for publication on BMC Bioinformatics Keywords: platform comparison

Dataset Information

EMDomics: a robust and powerful method for the identification of genes differentially expressed between heterogeneous classes.

Motivation

Results

Availability and implementation

Publications

EMDomics: a robust and powerful method for the identification of genes differentially expressed between heterogeneous classes.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets