A Literature-Derived Knowledge Graph Augments the Interpretation of Single Cell RNA-seq Datasets

• Main Authors: Deeksha Doddahonnaiah, Patrick J. Lenehan, Travis K. Hughes, David Zemmour, Enrique Garcia-Rivera, A. J. Venkatakrishnan, Ramakrishna Chilaka, Apoorv Khare, Akhil Kasaraneni, Abhinav Garg, Akash Anand, Rakesh Barve, Viswanathan Thiagarajan, Venky Soundararajan
• Format: Article
• Language: English
• Published: MDPI AG 2021-06-01
• Series: Genes
• Subjects: single cell genomics; natural language processing
• Online Access: https://www.mdpi.com/2073-4425/12/6/898

Technology to generate single cell RNA-sequencing (scRNA-seq) datasets and tools to annotate them have advanced rapidly in the past several years. Such tools generally rely on existing transcriptomic datasets or curated databases of cell type defining genes, while the application of scalable natural language processing (NLP) methods to enhance analysis workflows has not been adequately explored. Here we deployed an NLP framework to objectively quantify associations between a comprehensive set of over 20,000 human protein-coding genes and over 500 cell type terms across over 26 million biomedical documents. The resultant gene-cell type associations (GCAs) are significantly stronger between a curated set of matched cell type-marker pairs than the complementary set of mismatched pairs (Mann Whitney <i>p</i> = 6.15 × 10⁻⁷⁶, r = 0.24; cohen’s D = 2.6). Building on this, we developed an augmented annotation algorithm (single cell Annotation via Literature Encoding, or scALE) that leverages GCAs to categorize cell clusters identified in scRNA-seq datasets, and we tested its ability to predict the cellular identity of 133 clusters from nine datasets of human breast, colon, heart, joint, ovary, prostate, skin, and small intestine tissues. With the optimized settings, the true cellular identity matched the top prediction in 59% of tested clusters and was present among the top five predictions for 91% of clusters. scALE slightly outperformed an existing method for reference data driven automated cluster annotation, and we demonstrate that integration of scALE can meaningfully improve the annotations derived from such methods. Further, contextualization of differential expression analyses with these GCAs highlights poorly characterized markers of well-studied cell types, such as CLIC6 and DNASE1L3 in retinal pigment epithelial cells and endothelial cells, respectively. Taken together, this study illustrates for the first time how the systematic application of a literature-derived knowledge graph can expedite and enhance the annotation and interpretation of scRNA-seq data.