Comparative Analysis of Association Networks Using Single-Cell RNA Sequencing Data Reveals Perturbation-Relevant Gene Signatures

Single-cell RNA sequencing (scRNA-seq) data has elevated our understanding of systemic perturbations to organismal physiology at the individual cell level. However, despite the rich information content of scRNA-seq data, the relevance of genes to a perturbation is still commonly assessed through differential expression analysis. This approach provides a one-dimensional perspective of the transcriptomic landscape, risking the oversight of tightly controlled genes characterized by modest changes in expression but with profound downstream effects. We present GENIX (Gene Expression Network Importance eXamination), a novel platform for constructing gene association networks, equipped with an innovative network-based comparative model to uncover condition-relevant genes. To demonstrate the effectiveness of GENIX, we analyze influenza vaccine-induced immune responses in peripheral blood mononuclear cells (PBMCs) collected from recovered COVID-19 patients, shedding light on the mechanistic underpinnings of gender differences. Our methodology offers a promising avenue to identify genes relevant to perturbation responses in biological systems, expanding the scope of response signature discovery beyond differential gene expression analysis. HIGHLIGHTSO_LIConventional methods used to identify perturbation-relevant genes in scRNA-seq data rely on differential expression analysis, susceptible to overlooking essential genes. C_LIO_LIGENIX leverages cell-type-specific inferred gene association networks to identify condition-relevant genes and gene programs, irrespective of their specific expression alterations. C_LIO_LIGENIX provides insight into the gene-regulatory response to the influenza vaccine in naive and recovered COVID-19 patients, expanding on previously observed gender-specific differences. C_LI GRAPHICAL ABSTRACT O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=115 SRC="FIGDIR/small/556872v1_ufig1.gif" ALT="Figure 1"> View larger version (27K): org.highwire.dtl.DTLVardef@1837d3eorg.highwire.dtl.DTLVardef@1937860org.highwire.dtl.DTLVardef@c40114org.highwire.dtl.DTLVardef@22d3b9_HPS_FORMAT_FIGEXP M_FIG C_FIG