Abstract:Background: Gene signatures are important to represent the molecular changes in the disease genomes or the cells in specific conditions, and have been often used to separate samples into different groups for better research or clinical treatment. While many methods and applications have been available in literature, there still lack powerful ones that can take account of the complex data and detect the most informative signatures. Methods: In this article, we present a new framework for identifying gene signatures using Pareto-optimal cluster size identification for RNA-seq data. We first performed pre-filtering steps and normalization, then utilized the empirical Bayes test in Limma package to identify the differentially expressed genes (DEGs). Next, we used a multi-objective optimization technique, "Multi-objective optimization for collecting cluster alternatives" (MOCCA in R package) on these DEGs to find Pareto-optimal cluster size, and then applied k-means clustering to the RNA-seq data based on the optimal cluster size. The best cluster was obtained through computing the average Spearman's Correlation Score among all the genes in pair-wise manner belonging to the module. The best cluster is treated as the signature for the respective disease or cellular condition. Results: We applied our framework to a cervical cancer RNA-seq dataset, which included 253 squamous cell carcinoma (SCC) samples and 22 adenocarcinoma (ADENO) samples. We identified a total of 582 DEGs by Limma analysis of SCC versus ADENO samples. Among them, 260 are up-regulated genes and 322 are down-regulated genes. Using MOCCA, we obtained seven Pareto-optimal clusters. The best cluster has a total of 35 DEGs consisting of all-upregulated genes. For validation, we ran PAMR (prediction analysis for microarrays) classifier on the selected best cluster, and assessed the classification performance. Our evaluation, measured by sensitivity, specificity, precision, and accuracy, showed high confidence. Conclusions: Our framework identified a multi-objective based cluster that is treated as a signature that can classify the disease and control group of samples with higher classification performance (accuracy 0.935) for the corresponding disease. Our method is useful to find signature for any RNA-seq or microarray data.

Molecular Signatures from Gene Expression Data

Increasing stability and interpretability of gene expression signatures

Extraction and Comparison of Gene Expression Patterns from 2d Rna in Situ Hybridization Images

Identification of Gene Signatures for COAD Using Feature Selection and Bayesian Network Approaches

DNA-based molecular classifiers for the profiling of gene expression signatures

A simple and robust method for connecting small-molecule drugs using gene-expression signatures

Finding disease specific alterations in the co-expression of genes

Gene signatures for cancer research: A 25-year retrospective and future avenues

Computational Dissection of Tumor Expression Profiles.

Disease-related Gene Expression Analysis Using an Ensemble Statistical Test Method

Molecular Signature of Cancer at Gene Level or Pathway Level? Case Studies of Colorectal Cancer and Prostate Cancer Microarray Data.

Identification of gene signatures from RNA-seq data using Pareto-optimal cluster algorithm

The Iterative Signature Algorithm for the analysis of large scale gene expression data

ConSIG: Consistent Discovery of Molecular Signature from OMIC Data

Gesearch: An Interactive Gui Tool For Identifying Gene Expression Signature

Motif-guided sparse decomposition of gene expression data for regulatory module identification

A Common Gene Expression Signature Analysis Method for Multiple Types of Cancer

Prognostic relevance of gene-expression signatures

Data driven refinement of gene expression signatures for enrichment analysis

Towards Mechanism Classifiers: Expression-Anchored Gene Ontology Signature Predicts Clinical Outcome in Lung Adenocarcinoma Patients.

Identification of significant features in DNA microarray data