# scRDEN: single-cell dynamic gene rank differential expression network and robust trajectory inference

**Authors:** Han Zhang, Wei Zhang, Xiaoying Zheng, Yuanyuan Li

PMC · DOI: 10.1038/s41598-025-01969-1 · Scientific Reports · 2025-05-15

## TL;DR

The paper introduces scRDEN, a method for analyzing gene expression in single cells to better understand cell differentiation and regulatory mechanisms.

## Contribution

scRDEN introduces a novel framework for robust trajectory inference and dynamic gene rank differential expression network analysis in single-cell RNA sequencing data.

## Key findings

- scRDEN successfully captures stable cell subpopulations and identifies marker genes in human and mouse cell differentiation datasets.
- The method reveals non-monotonic trends in network properties, suggesting mechanisms of cell differentiation into stable functions.
- scRDEN performs exceptionally well on large-scale, multi-branched datasets like mouse dentate gyrus data.

## Abstract

The remarkable advancement of single-cell RNA sequencing (scRNA-seq) technology has empowered researchers to probe gene expression at the single-cell level with unprecedented precision. To gain a profound understanding of the heterogeneity inherent in cell fate determination, a central challenge lies in the comprehensive analysis of the dynamic regulatory alterations that underlie transcriptional differences and the accurate inference of the differentiation trajectory. Here, we propose the method scRDEN, a robust framework that infers important cell sub-populations and differential expression networks of multiple genes along the differentiation directions of each branch by converting the unstable gene expression values in cells into relatively stable gene-gene interactions (global features) and extracting the order of differential expression (network features), and further integrating the expression features of different dimension reduction methods. When applied to five published scRNA-seq datasets from human and mouse cell differentiation, scRDEN not only successfully captures the stable cell subpopulations with potential marker genes, measures the transcriptional differences of gene pairs to identify the rank differential expression network along the differentiation direction of each branch. In addition, in multiple gene rank differential expression networks, the rank expression directly related to transcription factors/marker genes shows a significant strengthening and weakening trend along with their expression changes, and the distribution of diversity and cluster coefficient show a non-monotonic change trend, including the cases of increasing first and then decreasing or decreasing first and then increasing. This may correspond to the mechanism of cells gradually differentiating into stable functions. It is particularly noteworthy that scRDEN method yielded exceptional results when applied to the large-scale, multi-branched, double-batch mouse dentate gyrus data. This outstanding performance provides novel and valuable insights into large-scale, multi-batch trajectory inference and the study of transcriptional mechanism regulation during the processes of differentiation and development.

## Linked entities

- **Species:** Homo sapiens (taxon 9606), Mus musculus (taxon 10090)

## Full-text entities

- **Genes:** Fhl2 (four and a half LIM domains 2) [NCBI Gene 14200] {aka FHL-2, SLIM-3, SLIM3}, Cst3 (cystatin C) [NCBI Gene 13010] {aka CysC}, Eomes (eomesodermin) [NCBI Gene 13813] {aka TBR-2, Tbr2}, Tnnc2 (troponin C2, fast) [NCBI Gene 21925] {aka Tncs}, Hmga2 (high mobility group AT-hook 2) [NCBI Gene 15364] {aka 9430083A20Rik, HMGI-C, Hmgic, pg, pygmy}, Fxyd6 (FXYD domain-containing ion transport regulator 6) [NCBI Gene 59095] {aka 0610030I18Rik, Php}, Fth1 (ferritin heavy polypeptide 1) [NCBI Gene 14319] {aka FHC, Fth, HFt, MFH}, Snap25 (synaptosomal-associated protein 25) [NCBI Gene 20614] {aka Bdr, GENA70, SNAP-25, SUP, sp}, Dner (delta/notch-like EGF repeat containing) [NCBI Gene 227325] {aka A930026D19Rik, BET, Bret}, Tagln2 (transgelin 2) [NCBI Gene 21346] {aka 2700094C18Rik, SM22beta, Sm22B, Sm22a}, Timp1 (tissue inhibitor of metalloproteinase 1) [NCBI Gene 21857] {aka Clgi, EPA, TIMP-1, TPA-S1, Timp}, Ube2c (ubiquitin-conjugating enzyme E2C) [NCBI Gene 68612] {aka 1110015A16Rik, D2Ertd695e}, Fabp5 (fatty acid binding protein 5, epidermal) [NCBI Gene 16592] {aka E-FABP, Fabpe, Klbp, PA-FABP, mal1}, Scgb3a2 (secretoglobin, family 3A, member 2) [NCBI Gene 117158] {aka LuLeu1, Pnsp1, UGRP1, Utgrp1}, Ascl1 (achaete-scute family bHLH transcription factor 1) [NCBI Gene 17172] {aka ASH1, Mash1, bHLHa46}, S100a6 (S100 calcium binding protein A6 (calcyclin)) [NCBI Gene 20200] {aka 2A9, 5B10, CALCYCLIN, Cacy, PRA}, Fabp7 (fatty acid binding protein 7, brain) [NCBI Gene 12140] {aka B-FABP, BFABP, Blbp, MRG}, Snca (synuclein, alpha) [NCBI Gene 20617] {aka NACP, alpha-Syn, alphaSYN}, C1ql2 (complement component 1, q subcomponent-like 2) [NCBI Gene 226359] {aka Adii, CTRP10}, Dbi (diazepam binding inhibitor) [NCBI Gene 13167] {aka ACBD1, Acbp, EP, endozepine}, Cspg4 (chondroitin sulfate proteoglycan 4) [NCBI Gene 121021] {aka 4732461B14Rik, AN2, Cspg4a, NG2}, Sftpc (surfactant associated protein C) [NCBI Gene 20389] {aka Bricd6, SP-C, SP5, SPC, Sftp-2, Sftp2}, Tnnt1 (troponin T1, skeletal, slow) [NCBI Gene 21955] {aka Tnt, sTnT, ssTnT}, Actb (actin, beta) [NCBI Gene 11461] {aka Actx, E430023M04Rik, beta-actin}, Birc5 (baculoviral IAP repeat-containing 5) [NCBI Gene 11799] {aka AAC-11, Api4, TIAP, survivin40}, Tubb3 (tubulin, beta 3 class III) [NCBI Gene 22152] {aka 3200002H15Rik, M(beta)3, M(beta)6}, Cadm1 (cell adhesion molecule 1) [NCBI Gene 54725] {aka 2900073G06Rik, 3100001I08Rik, Bl2, Igsf4, Igsf4a, Necl2}, Cldn11 (claudin 11) [NCBI Gene 18417] {aka Claudin-11, Claudin11, Osp, Otm}, Ftl1 (ferritin light polypeptide 1) [NCBI Gene 14325] {aka Ftl, Ftl-1, L-ferritin}, Ccnd2 (cyclin D2) [NCBI Gene 12444] {aka 2600016F06Rik, Vin-1, Vin1, cD2}, Cpe (carboxypeptidase E) [NCBI Gene 12876] {aka CPH, Cph-1, Cph1, NF-alpha1, fat}, Coro2b (coronin, actin binding protein, 2B) [NCBI Gene 235431] {aka CLIPINC, E130012P22Rik}, Actg1 (actin, gamma, cytoplasmic 1) [NCBI Gene 11465] {aka Actg, Actl, E51}
- **Diseases:** scRDEN (MESH:D030342), cancer (MESH:D009369)
- **Chemicals:** fatty acid (MESH:D005227), oxygen (MESH:D010100), calcium (MESH:D002118), RA (MESH:D011883), cytochalasin B (MESH:D003571), GDP-mannose (MESH:D006155)
- **Species:** Homo sapiens (human, species) [taxon 9606], Mus musculus (house mouse, species) [taxon 10090]
- **Cell lines:** MEF — Mus musculus (Mouse), Finite cell line (CVCL_9115), C6 — Rattus norvegicus (Rat), Rat malignant glioma, Cancer cell line (CVCL_0194), Neuron — Mus musculus (Mouse), Hybrid cell line (CVCL_U508), S2 — Drosophila melanogaster (Fruit fly), Spontaneously immortalized cell line (CVCL_Z232), fibroblasts — Mus musculus (Mouse), Spontaneously immortalized cell line (CVCL_0594)

## Full text

_Full body text omitted from this summary view._ Fetch the complete paper as Markdown: https://tomesphere.com/paper/PMC12081924/full.md

## Figures

5 figures with captions in the complete paper: https://tomesphere.com/paper/PMC12081924/full.md

## References

2 references — full list in the complete paper: https://tomesphere.com/paper/PMC12081924/full.md

---
Source: https://tomesphere.com/paper/PMC12081924