Search

Single-cell RNA-sequencing (scRNA-seq) has enabled the molecular profiling of thousands to millions of cells simultaneously in biologically heterogenous samples. Currently, common practice in scRNA-seq is to determine cell type labels through unsupervised clustering and the examination of cluster-specific genes. However, even small differences in analysis and parameter choice can greatly alter clustering solutions and thus impose great influence on which cell types are identified. Existing methods largely focus on determining the optimal number of robust clusters, which is not favorable for identifying cells of extremely low abundance due to their subtle contributions towards overall patterns of gene expression. Here we present a carefully designed framework, SCISSORS, which accurately profiles subclusters within major cluster(s) for the identification of rare cell types in scRNA-seq data. SCISSORS employs silhouette scoring for the estimation of heterogeneity of clusters and reveals rare cells in heterogenous clusters by implementing a multi-step, semi-supervised reclustering process. Additionally, SCISSORS provides a method for the identification of marker genes of rare cells, which may be used for further study. SCISSORS is wrapped around the popular Seurat R package and can be easily integrated into existing Seurat pipelines. SCISSORS, including source code and vignettes for two example datasets, is freely available at https://github.com/jrleary/SCISSORS.

Numerous aspects of cellular signaling are regulated by the kinome – the network of over 500 protein kinases that guides and modulates information transfer throughout the cell. The key role played by both individual kinases and assemblies of kinases organized into functional subnetworks leads to kinome dysregulation being a key driver of many diseases, particularly cancer. In the case of pancreatic ductal adenocarcinoma (PDAC), a variety of kinases and associated signaling pathways have been identified for their key role in the establishment of disease as well as its progression. However, the identification of additional relevant therapeutic targets has been slow and is further confounded by interactions between the tumor and the surrounding tumor microenvironment. Fundamentally, it is an open question as to the degree to which knowledge of the state of the kinome at the protein level is able to provide insight into the downstream phenotype of the cell.In this work, we attempt to link the state of the kinome, or kinotype, with cell viability in representative PDAC tumor and stroma cell lines. Through the application of both regression and classification models to independent kinome perturbation and kinase inhibitor cell screen data, we find that the inferred kinotype of a cell has a significant and predictive relationship with cell viability. While regression models perform poorly, we find that classification approaches are able to predict drug viability effects. We further find that models are able to identify a set of kinases whose behavior in response to perturbation drive the majority of viability responses in these cell lines. Using the models to predict new compounds with cell viability effects and not in the initial data set, we conducted a validation screen that confirmed the accuracy of the models. These results suggest that characterizing the state of the protein kinome provides significant opportunity for better understanding signaling behavior and downstream cell phenotypes, as well as providing insight into the broader design of potential therapy design for PDAC.