Hypergraph Representations of scRNA-seq Data for Improved Clustering with Random Walks

TL;DR

This paper introduces hypergraph-based methods for analyzing scRNA-seq data, capturing multi-way relationships and improving clustering performance over traditional network approaches, especially in weakly modular data.

Contribution

The paper proposes two novel hypergraph walk algorithms, DIPHW and CoMem-DIPHW, that outperform existing methods by leveraging higher-order relationships in scRNA-seq data.

Findings

Hypergraph methods outperform traditional networks in clustering accuracy.

CoMem-DIPHW effectively integrates multiple data representations.

Significant improvements in weakly modular datasets.

Abstract

Analysis of single-cell RNA sequencing data is often conducted through network projections such as coexpression networks, primarily due to the abundant availability of network analysis tools for downstream tasks. However, this approach has several limitations: loss of higher-order information, inefficient data representation caused by converting a sparse dataset to a fully connected network, and overestimation of coexpression due to zero-inflation. To address these limitations, we propose conceptualizing scRNA-seq expression data as hypergraphs, which are generalized graphs in which the hyperedges can connect more than two vertices. In the context of scRNA-seq data, the hypergraph nodes represent cells and the edges represent genes. Each hyperedge connects all cells where its corresponding gene is actively expressed and records the expression of the gene across different cells. This hypergraph conceptualization enables us to explore multi-way relationships beyond the pairwise interactions in coexpression networks without loss of information. We propose two novel clustering methods: (1) the Dual-Importance Preference Hypergraph Walk (DIPHW) and (2) the Coexpression and Memory-Integrated Dual-Importance Preference Hypergraph Walk (CoMem-DIPHW). They outperform established methods on both simulated and real scRNA-seq datasets. The improvement brought by our proposed methods is especially significant when data modularity is weak. Furthermore, CoMem-DIPHW incorporates the gene coexpression network, cell coexpression network, and the cell-gene expression hypergraph from the single-cell abundance counts data altogether for embedding computation. This approach accounts for both the local level information from single-cell level gene expression and the global level information from the pairwise similarity in the two coexpression networks.