Abstract:Abstract Background Gene expression is regulated by transcription factors, cofactors, and epigenetic mechanisms. Coexpressed genes indicate similar functional categories and gene networks. Detecting gene-gene coexpression is important for understanding the underlying mechanisms of cellular function and human diseases. A common practice of identifying coexpressed genes is to test the correlation of expression in a set of genes. In single-cell RNA-seq data, an important challenge is the abundance of zero values, so-called “dropout”, which results in biased estimation of gene-gene correlations for downstream analyses. In recent years, efforts have been made to recover coexpressed genes in scRNA-seq data. Here, our goal is to detect coexpressed gene pairs to reduce the “dropout” effect in scRNA-seq data using a novel graph-based k-partitioning method by merging transcriptomically similar cells. Results We observed that the number of zero values was reduced among the merged transcriptomically similar cell clusters. Motivated by this observation, we leveraged a graph-based algorithm and develop an R package, scCorr, to recover the missing gene-gene correlation in scRNA-seq data that enables the reliable acquisition of cluster-based gene-gene correlations in three independent scRNA-seq datasets. The graphically partitioned cell clusters did not change the local cell community. For example, in scRNA-seq data from peripheral blood mononuclear cells (PBMCs), the gene-gene correlation estimated by scCorr outperformed the correlation estimated by the nonclustering method. Among 85 correlated gene pairs in a set of 100 clusters, scCorr detected 71 gene pairs, while the nonclustering method detected only 4 pairs of a dataset from PBMCs. The performance of scCorr was comparable to those of three previously published methods. As an example of downstream analysis using scCorr, we show that scCorr accurately identified a known cell type (i.e., CD4+ T cells) in PBMCs with a receiver operating characteristic area under the curve of 0.96. Conclusions Our results demonstrate that scCorr is a robust and reliable graph-based method for identifying correlated gene pairs, which is fundamental to network construction, gene-gene interaction, and cellular omic analyses. scCorr can be quickly and easily implemented to minimize zero values in scRNA-seq analysis and is freely available at https://github.com/CBIIT-CGBB/scCorr .

Class-specific Correlations of Gene Expressions: Identification and Their Effects on Clustering Analyses.

Using Matrix of Thresholding Partial Correlation Coefficients to Infer Regulatory Network

Clustering gene expression data based on predicted differential effects of GV interaction.

Bi-clustering interpretation and prediction of correlation between gene expression and protein abundance

A New Strategy of Cooperativity of Biclustering and Hierarchical Clustering: a Case of Analyzing Yeast Genomic Microarray Datasets.

A novel graph-based k-partitioning approach improves the detection of gene-gene correlations by single-cell RNA sequencing

Diverse correlation structures in gene expression data and their utility in improving statistical inference

Association Study between Gene Expression and Multiple Relevant Phenotypes with Cluster Analysis

Single Cell Clustering Based on Cell-Pair Differentiability Correlation and Variance Analysis

Finding Correlated Biclusters from Gene Expression Data

Inferring gene-gene interactions and functional modules using sparse canonical correlation analysis

Negative correlation based gene markers identification in integrative gene expression data

Using Formal Concept Analysis to Identify Negative Correlations in Gene Expression Data

Distance correlation application to gene co-expression network analysis

Fuzzy soft subspace clustering method for gene co-expression network analysis

Corsig: A General Framework For Estimating Statistical Significance Of Correlation And Its Application To Gene Co-Expression Analysis

Quantitative trait associated microarray gene expression data analysis

Biclustering Of Linear Patterns In Gene Expression Data

Symmetry and Complexity in Gene Association Networks Using the Generalized Correlation Coefficient

Detecting Boolean Asymmetric Relationships with a Loop Counting Technique and its Implications for Analyzing Heterogeneity within Gene Expression Datasets

Defining and Identifying Cell Sub-Crosstalk Pairs for Characterizing Cell–cell Communication Patterns