MAT-Cell: A Multi-Agent Tree-Structured Reasoning Framework for Batch-Level Single-Cell Annotation
Organizations: Westlake University, Hangzhou, China · Shenzhen University of Advanced Technology, Shenzhen, China · Center for Artificial Intelligence and Robotics, Hong Kong Institute of Science and Innovation, Chinese Academy of Sciences, Hong Kong, China · University Key Laboratory of Information and Communication Security Backup and Recovery, Yunnan Minzu University, Kunming, China
Abstract
Automated single-cell annotation is difficult when the most abundant genes are not the most discriminative ones, or when a target state is poorly covered by a fixed reference atlas. GPTCelltype-style one-shot prompting allows large language models (LLMs) to produce plausible labels from generic expression signals, while reference-based annotators can force unfamiliar states into the nearest known category. We propose MAT-Cell, a prompt-driven framework for batch-level single-cell annotation that separates evidence grounding from label decision. MAT-Cell first uses Reverse Verification Query (RVQ) to combine tissue context, observed differentially expressed genes, and LLM-elicited biological priors into structured candidate-specific premises. Verifier agents then convert these premises into explicit premise-to-claim reasoning trees, and bounded multi-round debate compares,challenges, and revises the resulting claims before consensus or final adjudication.The returned Syllogistic Derivation Tree (SDT) provides an auditable debate trace rather than a formal proof of the annotation. In open-candidate benchmarks across five datasets, a locally deployed Qwen3-30B model with MAT-Cell achieves 75.5% average accuracy, compared with 64.2% for the strongest evaluated CoT baseline and 51.9% for the strongest evaluated scPilot variant. In oracle-candidate bench-marks across three species,MAT-Cell remains competitive across backbones, and local inference substantially reduces monetary cost for batch annotation. Code is available at: https://anonymous.4open.science/r/MATCell-4067
Figures & tables
| Method | Model | PBMC3K | Liver | Retina | Brain | Heart | Avg Acc. |
| CellTypist (Sci’22) | Non-LLM | 46.4 | 56.3 | 38.8 | 24.2 | 69.0 | 46.9 |
| CellMarker2.0 (NAR’23) | Non-LLM | 30.4 | 25.0 | 63.2 | 62.5 | 26.7 | 41.6 |
| CoT (NeurIPS’22) | Gemini 2.5 Pro | 62.5 4.4 | 78.2 0.8 | 49.5 1.2 | 66.9 1.8 | 64.0 2.8 | 64.2 |
| GPT-o1 | 67.5 5.2 | 70.4 1.3 | 49.0 1.4 | 66.9 4.2 | 61.3 1.8 | 63.0 | |
| Qwen3-14B | 57.8 3.1 | 64.8 8.5 | 66.3 7.3 | 64.0 6.0 | 51.1 10.5 | 60.8 | |
| GPT-4o | 68.8 0.0 | 80.4 2.5 | 44.2 2.2 | 45.5 2.0 | 58.0 3.8 | 59.4 |
| Method | Model | Human | Mouse | Monkey | Avg Acc. | ||||||
| Both | DEG | HEG | Both | DEG | HEG | Both | DEG | HEG | |||
| scCATCH (iScience’20) | Non-LLM | 33.9 1.0 | 30.2 1.2 | 28.6 1.4 | 37.4 0.9 | 36.7 0.5 | 30.3 1.0 | 41.8 0.7 | 42.2 0.8 | 20.0 0.8 | 33.5 |
| Cell-ID (Nat Biotech’21) | Non-LLM | 43.1 0.4 | 43.2 0.3 | 21.8 1.1 | 28.8 1.3 | 28.9 1.1 | 16.0 0.7 | 51.7 0.5 | 49.8 1.1 | 30.3 0.7 | 34.8 |
| Cell-o1 (arXiv’25) | Qwen2.5-7B | 42.9 1.2 | 40.9 1.5 | 24.3 1.1 | 39.0 1.6 | 39.4 1.7 | 23.2 0.8 | 69.5 1.8 | 68.5 1.9 | 50.3 1.2 | 44.2 |
| CoT (NeurIPS’22) | Qwen3-14B | 48.7 0.5 | 61.8 1.7 | 21.6 0.4 | 47.9 0.8 | 53.4 1.5 | 27.4 0.9 | 67.2 0.8 | 72.0 0.6 | 40.4 1.0 | 48.9 |
| Qwen3-30B | 55.5 0.3 | 62.7 0.5 | 21.3 0.4 | 47.6 2.2 | 52.1 0.7 | 25.6 0.4 | 64.9 1.1 | 71.6 1.1 | 35.7 0.2 | 48.6 | |
| ID | Config | H | M | Mk | Gap | RI |
| M0 | GPTCelltype+HEG | 16.5% | 17.8% | 33.1% | 16.6% | 73% |
| M1 | GPTCelltype+DEG | 45.0% | 38.7% | 64.4% | 25.7% | 78% |
| M2 | CoT+HEG | 21.3% | 25.6% | 35.7% | 14.4% | 75% |
| M3 | CoT+DEG | 62.7% | 52.1% | 71.6% | 19.5% | 85% |
| M4 | MAT-Cell+HEG | 28.2% | 49.9% | 39.5% | 21.7% | 80% |
| M5 | MAT-Cell+DEG | 76.4% | 82.5% | 75.9% | 6.6% | 97% |
| Metric | CoT | CoT | GPTCelltype | MAT-Cell |
| Model | GPT-4.1 | Gemini 2.5 Flash | GPT-4.1 | Qwen3-30B |
| Deployment | API | API | API | Local A100 |
| Accuracy | 68.1% | 62.9% | 64.9% | 82.5% |
| Batch size | 1 | 1 | 1 | 20 |
| Time/query | 8.0s | 6.2s | 1.93s | 9.17s |
| Cost/query | 0.412 ¥ | 0.356 ¥ | 0.1100 ¥ | 0.0079 ¥ |
Appendix figures & tables16 assets
Supplementary material from the paper’s appendix.
Appendix
| Analytical Dimension | Metric Description | Value |
| Part I: Overall Performance Gains | Cell-level Avg. Partial Accuracy (GPTCelltype baseline) | 0.4278 |
| Cell-level Avg. Partial Accuracy (MAT-Cell) | 0.7268 | |
| Part II: Debate Convergence | Pre-debate (First Round) Full Agreement Rate (Mean) | 0.7632 |
| Post-debate (Last Round) Full Agreement Rate (Mean) | 0.9052 | |
| Pearson (Pre vs. Post Agreement, per batch) | 0.7538 | |
| Part III: Decision Agent (DA) Efficacy | DA Fallback Frequency (Disputed Batches) | 11.50% ( ) |
| Cell type | Count |
| CD4+ T Cells | 1,157 |
| Classical Monocytes | 479 |
| B Cells/Dendritic Cells | 341 |
| Effector/Activated T Cells | 297 |
| Effector/Activated T Cells/NK Cells | 158 |
| Non-Classical Monocytes | 157 |
| Property | Human | Mouse | Monkey |
| Number of batches | 2,400 | 2,400 | 2,400 |
| Total cells | 27,588 | 27,583 | 25,871 |
| Number of cell types | 75 | 123 | 121 |
| Unique top genes | 5,583 | 6,941 | 7,757 |
| Unique DEG genes | 1,434 | 2,432 | 2,485 |
| Cells per batch | 7–15 | 7–15 | 7–15 |
| Cell type | Count |
| Oligodendrocyte | 1,593 |
| L2/3–6 intratelencephalic projecting glutamatergic neuron | 1,560 |
| Astrocyte | 1,544 |
| Oligodendrocyte precursor cell | 1,508 |
| L2/3 intratelencephalic projecting glutamatergic neuron | 1,497 |
| Pvalb GABAergic cortical interneuron | 1,488 |
| Dataset | Top-25 genes |
| Brain | CNTNAP2, DSCAM, DPP10, ROBO2, KCNIP4, GRIP1, ZNF385D, EPIC1, CA10, FSTL4, ARL15, HTR1F, FOXP2, FSTL5, MYO16, PTCHD4, CLSTN2, CPNE4, NRG1, DTNA, CBLN2, CDH9, SLIT3, SLIT2, UNC13C |
| Heart | RGS6, ANK3, LINC02147, CLIC5, HIGD1B, PIK3R5, SLC38A11, MLIP, TRDN-AS1, MYBPC3, XIRP2, MYH6, ACACB, PRKAG2, ACTA1, SH3RF2, PPP1R3C, ENSG00000230490, MYL7, LINC02552, ENSG00000258231, ENSG00000271959, HECW2, FRMD5, G0S2 |
| Liver | H3f3a, Ubb, Tmsb10, Hba-a1, Gpx1, Hba-a2, Hbb-bt, S100a8, S100a9, Apoe, Igfbp7, Sparc, S100a6, Cd24a, Pglyrp1, Gm5483, BC100530, Stfa1, Anxa1, Ifitm6, Serpina3k, Mup20, Gnmt, Cd7, Ccl4 |
| PBMC3K | B2M, RPL13, MALAT1, RPL21, TPT1, RPL10A, ACTB, RPL8, H3F3B, RPS3A, RPS5, EEF1D, RPS27A, FTH1, MT-CO2, CD74, FTL, OAZ1, CD37, CD79A, FCGR3A, LST1, COTL1, GNLY, GZMB |
| Retina | FTH1, RHO, APOE, RBP3, WIF1, GLUL, CLU, PTGDS, TF, MFGE8, MPP4, FRZB, CRABP1, SPP1, ENO1, DKK3, RLBP1, CA2, GPX3, CRYAB, CADPS, NEAT1, NR2E3, YPEL2, CNGA1 |
| Human | MALAT1, ACTB, ACTG1, B2M, TMSB4X, FTH1, GAPDH, RPL13, RPS27, RPL41, RPL21, EEF1A1, TPT1, RPS3A, RPL32, RPL3, RPS2, RPS18, RPS6, RPS12, RPL10, RPL34, RPS27A, RPL13A, RPL11 |
| Rule ID | Markers (IF) | Cell Type (THEN) |
| R1 | CD4, IL7R, TCF7 | CD4+ T-Cell |
| R2 | CD8A, CD8B, GZMA | CD8+ T-Cell |
| R3 | CD14, LYZ, FCGR3A | Monocyte |
| R4 | PPBP, PF4, TUBB1 | Platelet |
| Parameter | Value |
| Number of Agent Types | 4 |
| #ExpertAgent (EA) Instances | 1 |
| #SolveAgent (SA) Instances | 1 |
| Council Scale ( , #VerifierAgent/VA) | 5 |
| #DecisionAgent (DA) Instances | 1 |
| Top-K genes (Both view) | 25 |
| Metric | Value |
| Consensus Accuracy | |
| Partial Cell Accuracy | |
| Pairwise Disagreement Rate | |
| Prediction Correlation | |
| Cohen’s Kappa | |
| Mutual Information (nats) |
| Top-K | Human (DEG) | Monkey (DEG) | Mouse (DEG) |
| Top-5 | |||
| Top-10 | |||
| Top-15 | |||
| Top-25 | |||
| Top-50 |
| Seed | Human (DEG) | Monkey (DEG) | Mouse (DEG) | Deviation from Mean |
| Seed-1 | H:-0.21%, Mk:-0.37%, M:-1.52% | |||
| Seed-2 | H:+1.77%, Mk:-0.19%, M:-0.34% | |||
| Seed-3 | H:+0.87%, Mk:+0.34%, M:+0.01% | |||
| Seed-4 | H:-2.43%, Mk:-0.84%, M:+1.85% | |||
| Mean Std | Stable |
| Method | Human (both) | Monkey (both) | Mouse (both) |
| RVQ (5 markers) | |||
| RVQ (10 markers) | |||
| RVQ (15 markers) | |||
| RVQ (20 markers) |
| Method | Human (DEG) | Mouse (DEG) | Monkey (DEG) | |||
| Acc. | 95% CI | Acc. | 95% CI | Acc. | 95% CI | |
| Cell-o1 | [0.403, 0.414] | [0.388, 0.401] | [0.678, 0.692] | |||
| GPTCelltype Qwen-3-30B | [0.445, 0.455] | [0.384, 0.390] | [0.641, 0.647] | |||
| GPTCelltype DeepSeek-V3 | [0.627, 0.636] | [0.563, 0.571] | [0.468, 0.474] | |||
| GPTCelltype Gemini-2.5 | [0.705, 0.713] | [0.653, 0.665] | [0.856, 0.862] | |||
| GPTCelltype GPT-4.1 | [0.731, 0.736] | [0.644, 0.654] | [0.863, 0.866] | |||
| Comparison | Mean Diff | t-statistic | df | p-value |
| MAT-Cell vs. GPT-4.1 | 29 | ** | ||
| MAT-Cell vs. Gemini-2.5 | 29 | *** | ||
| MAT-Cell vs. DeepSeek-V3 | 29 | *** | ||
| MAT-Cell vs. GPTCelltype-Qwen-3-30B | 29 | *** |
| Comparison | Mean Diff | t-statistic | df | p-value |
| MAT-Cell vs. GPT-4.1 | 29 | ** | ||
| MAT-Cell vs. Gemini-2.5 | 29 | *** | ||
| MAT-Cell vs. DeepSeek-V3 | 29 | *** | ||
| MAT-Cell vs. GPTCelltype-Qwen | 29 | *** |
| Comparison | Human (d) | Mouse (d) | Monkey (d) |
| GPT-4.1 | 3.21 (very large) | 17.01 (very large) | 1.26 (large) |
| Gemini-2.5-Flash | 5.29 (very large) | 14.04 (very large) | 0.53 (medium) |
| DeepSeek-V3 | 11.48 (very large) | 31.18 (very large) | 6.07 (very large) |
| GPTCelltype-Qwen-3-30B | 26.01 (very large) | 59.38 (very large) | 8.86 (very large) |
LLM-Assisted Annotation
Large language model-assisted annotation uses AI to label, classify, and structure data, often with human review to improve scale and consistency.
See all 202 papers