# SEPAR enables spatial metagene discovery and associated molecular pattern characterization in spatial transcriptomics and multi-omics datasets

**Authors:** Lei Zhang, Ying Zhu, Shuqin Zhang

PMC · DOI: 10.1038/s42003-025-09340-w · Communications Biology · 2025-12-10

## TL;DR

SEPAR is a new computational tool that helps analyze spatial transcriptomics data by identifying localized gene patterns and improving the understanding of spatial molecular interactions.

## Contribution

SEPAR introduces a novel unsupervised framework using spatial metagenes to uncover localized gene programs and spatial molecular interactions in SRT and multi-omics data.

## Key findings

- SEPAR identifies biologically meaningful gene ontologies and cell types linked to metagene patterns.
- SEPAR detects spatially variable genes with higher accuracy and enhances biological signals through gene refinement.
- SEPAR uncovers co-localized molecule correlations in spatial CITE-seq and gene-peak relationships in MISAR-seq.

## Abstract

Spatially resolved transcriptomics (SRT) profiles gene expressions at near- or sub-cellular resolution while preserving their spatial context, yet interpreting SRT data to understand spatial cellular and molecular organization remains challenging. Most existing computational methods focus on global spatial domains but overlook localized structures driven by specific gene subsets. Here, we introduce SEPAR, an unsupervised framework that leverages spatial metagenes to analyze SRT data by integrating gene activity and spatial neighborhood relationships. It enables multiple downstream analyses including: identifying metagene pattern-specific genes, detecting spatially variable genes (SVGs), delineating spatial domains, and refining expression signals. Evaluated on diverse datasets, SEPAR reveals biologically meaningful gene ontologies and cell types in gene sets linked to metagene patterns, identifies SVGs with higher accuracy, and enhances biological signals with gene refinement. In spatial multi-omics data, it uncovers co-localized molecule correlations in spatial CITE-seq and coordinated gene-peak relationships in MISAR-seq, offering insights into spatial molecular interactions.

SEPAR is a computational approach that uses spatial metagenes to reveal localized gene programs, identify spatially variable genes, delineate spatial domains, and refine gene expressions for spatial transcriptomics, and uncover coordinated interactions for spatial multi-omics.

## Full-text entities

- **Genes:** Tnnc2 (troponin C2, fast) [NCBI Gene 21925] {aka Tncs}, Myh3 (myosin, heavy polypeptide 3, skeletal muscle, embryonic) [NCBI Gene 17883] {aka MyHC-emb, Myhs-e, Myhse}, VIP (vasoactive intestinal peptide) [NCBI Gene 7432] {aka PHM27}, Tnnt1 (troponin T1, skeletal, slow) [NCBI Gene 21955] {aka Tnt, sTnT, ssTnT}, Ttn (titin) [NCBI Gene 22138] {aka 1100001C23Rik, 2310036G12Rik, 2310057K23Rik, 2310074I15Rik, D330041I19Rik, D830007G01Rik}, Tnni1 (troponin I, skeletal, slow 1) [NCBI Gene 21952] {aka 2700018B22Rik, ssTnI}, Tfrc (transferrin receptor) [NCBI Gene 22042] {aka 2610028K12Rik, CD71, E430033M20Rik, Mtvr1, TFR, TFR1}, Tcf21 (transcription factor 21) [NCBI Gene 21412] {aka Pod-1, Pod1, bHLHa23, capsulin, epc, epicardin}, Hexb (hexosaminidase B) [NCBI Gene 15212], Kcnip2 (Kv channel-interacting protein 2) [NCBI Gene 80906] {aka KChIP2}, Myrf (myelin regulatory factor) [NCBI Gene 225908] {aka 6030439E18, Gm1804, Gm98, Mrf}, Acta1 (actin alpha 1, skeletal muscle) [NCBI Gene 11459] {aka Acta-2, Acts, Actsk-1}, Myh8 (myosin, heavy polypeptide 8, skeletal muscle, perinatal) [NCBI Gene 17885] {aka 4832426G23Rik, MHCp, MyHC-pn, Myhs-p, Myhsp}, Col1a2 (collagen, type I, alpha 2) [NCBI Gene 12843] {aka Col1a-2, Cola-2, Cola2, oim}, Des (desmin) [NCBI Gene 13346], Myf5 (myogenic factor 5) [NCBI Gene 17877] {aka B130010J22Rik, Myf-5, bHLHc2}, Neb (nebulin) [NCBI Gene 17996], Tbx21 (T-box 21) [NCBI Gene 57765] {aka TBT1, Tbet, Tblym}, Neurod1 (neurogenic differentiation 1) [NCBI Gene 18012] {aka BETA2, BHF-1, Nd1, Neurod, bHLHa3}, Xcl1 (chemokine (C motif) ligand 1) [NCBI Gene 16963] {aka ATAC, LTN, Lptn, SCM-1, SCM-1a, Scyc1}, Rorb (RAR-related orphan receptor beta) [NCBI Gene 225998] {aka Nr1f2, RZR-beta, RZRB, Rorbeta, hstp}, Tnni2 (troponin I, skeletal, fast 2) [NCBI Gene 21953], PVALB (parvalbumin) [NCBI Gene 5816] {aka D22S749}, Foxp2 (forkhead box P2) [NCBI Gene 114142] {aka 2810043D05Rik, CAG-16, D0Kist7}, Neurog2 (neurogenin 2) [NCBI Gene 11924] {aka Atoh4, Math4A, bHLHa8, ngn-2, ngn2}, Nt5e (5' nucleotidase, ecto) [NCBI Gene 23959] {aka 2210401F01Rik, 5'-NT, CD73, NT, Nt5, eNT}, Tcf12 (transcription factor 12) [NCBI Gene 21406] {aka A130037E08Rik, ALF1, HEB, HEBAlt, HTF-4, HTF4}, Cd74 (CD74 antigen (invariant polypeptide of major histocompatibility complex, class II antigen-associated)) [NCBI Gene 16149] {aka CLIP, DHLAG, HLADG, Ia-GAMMA, Ii}, Neurod2 (neurogenic differentiation 2) [NCBI Gene 18013] {aka Ndrf, bHLHa1}, Plp1 (proteolipid protein (myelin) 1) [NCBI Gene 18823] {aka DM20, Plp, jimpy, jp, msd, rsh}, Cux2 (cut-like homeobox 2) [NCBI Gene 13048] {aka 1700051K22Rik, Cutl2, Cux-2}, Tnfrsf9 (tumor necrosis factor receptor superfamily, member 9) [NCBI Gene 21942] {aka 4-1BB, A930040I11Rik, CDw137, Cd137, ILA, Ly63}, Airn (antisense Igf2r RNA) [NCBI Gene 104103] {aka 2810051F02Rik, 2810434M15Rik, Air, B930018I07Rik, D17Ertd663e, IGF2RAS}, Mybpc1 (myosin binding protein C, slow-type) [NCBI Gene 109272] {aka 8030451F13Rik}, LAMP5 (lysosome associated membrane protein 5) [NCBI Gene 24141] {aka BAD-LAMP, BADLAMP, C20orf103, LAMP-5, UNC-46}, Tnnt3 (troponin T3, skeletal, fast) [NCBI Gene 21957] {aka fTnT}, Ttr (transthyretin) [NCBI Gene 22139] {aka prealbumin}, Actc1 (actin, alpha, cardiac muscle 1) [NCBI Gene 11464] {aka Actc-1}, Ywhab (tyrosine 3-monooxygenase/tryptophan 5-monooxygenase activation protein, beta polypeptide) [NCBI Gene 54401] {aka 1300003C17Rik, 14-3-3b}, Myog (myogenin) [NCBI Gene 17928] {aka MYF4, bHLHc3, myo}, Neurod6 (neurogenic differentiation 6) [NCBI Gene 11922] {aka Atoh2, Math-2, Math2, Nex, Nex1m, bHLHa2}, Slc17a7 (solute carrier family 17 (sodium-dependent inorganic phosphate cotransporter), member 7) [NCBI Gene 72961] {aka 2900052E22Rik, Vglut1}, Myod1 (myogenic differentiation 1) [NCBI Gene 17927] {aka MYF3, MyoD, Myod-1, bHLHc1}, SNCG (synuclein gamma) [NCBI Gene 6623] {aka BCSG1, SR}, Foxj1 (forkhead box J1) [NCBI Gene 15223] {aka FKHL-13, HFH-4, Hfh4}, Tbata (thymus, brain and testes associated) [NCBI Gene 65971] {aka 1700021K02Rik, Spatial, Titest}, Pou2f2 (POU domain, class 2, transcription factor 2) [NCBI Gene 18987] {aka Oct-2, Oct2a, Oct2b, Oct2c, Oct2d, Otf-2}, SST (somatostatin) [NCBI Gene 6750] {aka SMST, SST1}, Atoh1 (atonal bHLH transcription factor 1) [NCBI Gene 11921] {aka Hath1, MATH-1, Math1, bHLHa14}, Gfap (glial fibrillary acidic protein) [NCBI Gene 14580], Cnr1 (cannabinoid receptor 1) [NCBI Gene 12801] {aka CB-R, CB1, CB1A, CB1B, CB1R}, Omp (olfactory marker protein) [NCBI Gene 18378]
- **Diseases:** cancer (MESH:D009369), SVGs (MESH:D008569), CRC (MESH:D015179)
- **Chemicals:** lipid (MESH:D008055), H&amp;E (MESH:D006371), CITE (-)
- **Species:** Mus musculus (house mouse, species) [taxon 10090], Homo sapiens (human, species) [taxon 9606]

## Full text

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

## Figures

7 figures with captions in the complete paper: https://tomesphere.com/paper/PMC12820152/full.md

## References

1 references — full list in the complete paper: https://tomesphere.com/paper/PMC12820152/full.md

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