Single sample scoring of molecular phenotypes.

BACKGROUND: Gene set scoring provides a useful approach for quantifying concordance between sample transcriptomes and selected molecular signatures. Most methods use information from all samples to score an individual sample, leading to unstable scores in small data sets and introducing biases from sample composition (e.g. varying numbers of samples for different cancer subtypes). To address these issues, we have developed a truly single sample scoring method, and associated R/Bioconductor package singscore ( https://bioconductor.org/packages/singscore ). RESULTS: We use multiple cancer data sets to compare singscore against widely-used methods, including GSVA, z-score, PLAGE, and ssGSEA. Our approach does not depend upon background samples and scores are thus stable regardless of the composition and number of samples being scored. In contrast, scores obtained by GSVA, z-score, PLAGE and ssGSEA can be unstable when less data are available (NS < 25). The singscore method performs as well as the best performing methods in terms of power, recall, false positive rate and computational time, and provides consistently high and balanced performance across all these criteria. To enhance the impact and utility of our method, we have also included a set of functions implementing visual analysis and diagnostics to support the exploration of molecular phenotypes in single samples and across populations of data. CONCLUSIONS: The singscore method described here functions independent of sample composition in gene expression data and thus it provides stable scores, which are particularly useful for small data sets or data integration. Singscore performs well across all performance criteria, and includes a suite of powerful visualization functions to assist in the interpretation of results. This method performs as well as or better than other scoring approaches in terms of its power to distinguish samples with distinct biology and its ability to call true differential gene sets between two conditions. These scores can be used for dimensional reduction of transcriptomic data and the phenotypic landscapes obtained by scoring samples against multiple molecular signatures may provide insights for sample stratification.

Identifying Key Regulators of Keratinization in Lung Squamous Cell Cancer Using Integrated TCGA Analysis.

Keratinization is one of lung squamous cell cancer's (LUSC) hallmark histopathology features. Epithelial cells produce keratin to protect their integrity from external harmful substances. In addition to their roles as cell protectors, recent studies have shown that keratins have important roles in regulating either normal cell or tumor cell functions. The objective of this study is to identify the genes and microRNAs (miRNAs) that act as key regulators of the keratinization process in LUSC. To address this goal, we classified LUSC samples from GDC-TCGA databases based on their keratinization molecular signatures. Then, we performed differential analyses of genes, methylation, and miRNA expression between high keratinization and low keratinization samples. By reconstruction and analysis of the differentially expressed genes (DEGs) network, we found that TP63 and SOX2 were the hub genes that were highly connected to other genes and displayed significant correlations with several keratin genes. Methylation analysis showed that the P63, P73, and P53 DNA-binding motif sites were significantly enriched for differentially methylated probes. We identified SNAI2, GRHL3, TP63, ZNF750, and FOXE1 as the top transcription factors associated with these binding sites. Finally, we identified 12 miRNAs that influence the keratinization process by using miRNA-mRNA correlation analysis.