Reconstructing the Developmental Trajectory of Adipocytes in Human Adipose Tissue Using Single-Cell RNA Sequencing

Humasak Tommy Argo Simanjuntak; Weny S. M Sitinjak

arxiv: 2606.27657 · v1 · pith:3MMVYA7Lnew · submitted 2026-06-26 · 🧬 q-bio.GN · cs.AI

Reconstructing the Developmental Trajectory of Adipocytes in Human Adipose Tissue Using Single-Cell RNA Sequencing

Weny S. M Sitinjak , Humasak Tommy Argo Simanjuntak This is my paper

Pith reviewed 2026-06-29 02:23 UTC · model grok-4.3

classification 🧬 q-bio.GN cs.AI

keywords single-cell RNA sequencingadipocyte differentiationadipose tissueIGF signalingFGF signalingtrajectory inferencedepot-specific differencesextracellular matrix remodeling

0 comments

The pith

Single-cell RNA sequencing of human adipose tissue reconstructs adipocyte development as a sequence of 15 cell clusters connected by 7 transitional states, with IGF and FGF pathways as the dominant signals.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper applies single-cell RNA sequencing to samples of human adipose tissue to trace how adipocytes form from progenitor cells. It identifies 15 transcriptionally distinct clusters, seven of which represent transitional states along a differentiation path. Sixteen signaling pathways are detected as active during this process, with insulin-like growth factor and fibroblast growth factor networks showing the strongest and most consistent activity. The analysis also finds that visceral fat undergoes extra extracellular-matrix remodeling not seen in subcutaneous fat, and that the two pathways operate in different spatial niches. These observations are presented as the first comprehensive map of human adipocyte development and as a source of potential intervention points for metabolic disease.

Core claim

Single-cell RNA sequencing of human adipose tissue identifies 15 transcriptionally distinct cell clusters that include seven transitional states, thereby reconstructing the developmental trajectory of adipocytes. Sixteen signaling pathways mediate communication between adipocytes and progenitors, with IGF and FGF pathways displaying the highest activity across stages. Depot-specific differences appear, including additional extracellular matrix remodeling in visceral adipocytes that is absent during subcutaneous differentiation. Spatial mapping places IGF signaling mainly in perivascular niches and FGF activity in zones of mature adipocytes.

What carries the argument

Trajectory inference applied to single-cell RNA sequencing data that groups cells into 15 clusters and 7 transitional states while scoring 16 active signaling pathways.

If this is right

IGF and FGF pathways remain active at every stage and therefore constitute candidate targets for promoting healthy fat expansion or limiting excess accumulation.
Visceral adipocytes perform additional extracellular-matrix remodeling that is not required in subcutaneous differentiation.
IGF signaling concentrates in perivascular locations while FGF signaling concentrates around mature adipocytes.
The identified networks supply a reference for designing interventions that modulate adipose tissue behavior in metabolic disorders.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

If the trajectory holds, drugs that selectively modulate IGF or FGF signaling could be tested for effects on fat depot expansion versus inflammation.
Depot-specific remodeling differences imply that therapies may need to be tailored to visceral versus subcutaneous tissue rather than applied uniformly.
Spatial separation of the two pathways suggests that local delivery or niche-targeted agents could improve specificity over systemic treatment.

Load-bearing premise

The 15 transcriptionally distinct clusters and 7 transitional states reflect a genuine developmental sequence rather than technical artifacts or unrelated cell populations.

What would settle it

A time-course experiment or independent marker validation that shows the 15 clusters do not form the inferred differentiation order would falsify the trajectory reconstruction.

read the original abstract

Obesity is a global health crisis associated with metabolic disorders such as type 2 diabetes and cardiovascular disease. This study employed single-cell RNA sequencing to reconstruct the developmental trajectory of human adipocytes from adipose tissue samples. Our analysis identified 15 transcriptionally distinct cell clusters, including 7 transitional states, revealing the dynamic process of adipocyte differentiation. We detected 16 functionally active signaling pathways mediating cellular communication between adipocytes and their progenitors. Among these, insulin-like growth factor (IGF) and fibroblast growth factor (FGF) pathways emerged as the most prominent networks, showing consistent activity across differentiation stages (p<0.05). The study revealed depot-specific differences, with visceral adipocytes undergoing additional extracellular matrix remodeling absent in subcutaneous differentiation. Spatial analysis further showed that IGF signaling was particularly active in perivascular niches, while FGF activity dominated in mature adipocyte zones. These results provide the first comprehensive map of human adipocyte development, highlighting IGF and FGF pathways as potential therapeutic targets. The identified signaling networks offer new insights for developing interventions to promote healthy adipose expansion or inhibit pathological fat accumulation. This work advances our fundamental understanding of adipose tissue biology while providing clinically relevant data for metabolic disorder treatments.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

Standard scRNA-seq trajectory analysis on human adipocytes, but the developmental map claims rest on unvalidated cluster ordering.

read the letter

The paper runs single-cell RNA sequencing on human adipose tissue and pulls out 15 transcriptionally distinct clusters with 7 transitional states, then flags IGF and FGF pathways as active during differentiation plus some depot differences such as extra extracellular matrix remodeling in visceral samples.

It does a reasonable job extending prior work by using actual human depot samples and adding a spatial angle that ties IGF activity to perivascular areas. That part gives a concrete picture of where signaling happens in the tissue.

The soft spot is the trajectory inference itself. The abstract presents the 15 clusters and 7 states as a developmental sequence and reports pathway activity across stages, but supplies no checks against known marker timing, no lineage-tracing comparison, and no discussion of batch or technical structure. Without those steps the ordering could easily be an artifact rather than a real differentiation path, which undercuts the "first comprehensive map" claim.

This is for labs already working on adipose biology who want to see human single-cell data layered on top of existing mouse or bulk studies. A reader focused on metabolic disease pathways might pull the IGF/FGF observations for follow-up, but anyone needing a reliable trajectory will have to wait for the methods and validation sections.

I would bring it to a reading group to talk through the trajectory method choices. I would not cite it until the validation is clearer. It deserves peer review so the full methods can be examined.

Referee Report

2 major / 2 minor

Summary. This manuscript uses single-cell RNA sequencing on human adipose tissue samples to reconstruct adipocyte developmental trajectories. It reports identification of 15 transcriptionally distinct clusters (including 7 transitional states), detection of 16 signaling pathways active during differentiation (with IGF and FGF pathways highlighted as most prominent at p<0.05), depot-specific differences (visceral adipocytes showing additional ECM remodeling not seen in subcutaneous), and spatial localization of signaling (IGF in perivascular niches, FGF in mature zones). The central claim is that these data constitute the first comprehensive map of human adipocyte development and identify IGF/FGF as therapeutic targets for metabolic disorders.

Significance. If the trajectory reconstruction and pathway activity claims hold after rigorous validation, the work would provide a useful resource for adipose biology by mapping differentiation stages and highlighting depot heterogeneity and spatial signaling niches. This could inform studies on healthy vs. pathological adipose expansion and suggest candidate pathways for intervention in obesity-related disease.

major comments (2)

[Abstract/Results] Abstract and Results (trajectory inference): The central claim that the 15 clusters and 7 transitional states constitute a true developmental trajectory (rather than static subpopulations or inference artifacts) is load-bearing for all downstream conclusions about differentiation stages, pathway activity, and the 'first comprehensive map.' No orthogonal validation is described, such as alignment with established markers (e.g., PPARG, CEBPA, FABP4 dynamics), comparison to bulk time-course data, or independent lineage tracing.
[Abstract/Methods] Abstract and Methods (signaling pathways): The statement that 16 pathways were detected as 'functionally active' with IGF/FGF 'most prominent' (p<0.05) across stages underpins the therapeutic-target conclusion, yet the abstract supplies no description of the cell-cell communication tool, statistical test, multiple-testing correction, or how activity was quantified across the inferred trajectory.

minor comments (2)

[Abstract] The abstract refers to 'spatial analysis' showing niche-specific signaling but provides no details on the spatial transcriptomics platform, resolution, or how niches were defined.
[Abstract] Depot-specific ECM remodeling is asserted as 'absent in subcutaneous differentiation' without quantification or statistical comparison between depots.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their thoughtful review and constructive feedback on our manuscript. We address each of the major comments below and will revise the manuscript accordingly to strengthen the presentation of our findings.

read point-by-point responses

Referee: [Abstract/Results] Abstract and Results (trajectory inference): The central claim that the 15 clusters and 7 transitional states constitute a true developmental trajectory (rather than static subpopulations or inference artifacts) is load-bearing for all downstream conclusions about differentiation stages, pathway activity, and the 'first comprehensive map.' No orthogonal validation is described, such as alignment with established markers (e.g., PPARG, CEBPA, FABP4 dynamics), comparison to bulk time-course data, or independent lineage tracing.

Authors: We agree that orthogonal validation is important to support the inferred trajectory. In the revised manuscript, we will add a new section in Results showing the expression dynamics of canonical adipocyte differentiation markers (PPARG, CEBPA, FABP4) along the inferred pseudotime trajectory, demonstrating progressive upregulation consistent with known biology. We will also include a comparison of our trajectory to published bulk RNA-seq datasets from in vitro adipocyte differentiation time courses to confirm the ordering of stages. While independent lineage tracing is not possible in human tissue samples, these additions will provide stronger support for the developmental trajectory interpretation. revision: yes
Referee: [Abstract/Methods] Abstract and Methods (signaling pathways): The statement that 16 pathways were detected as 'functionally active' with IGF/FGF 'most prominent' (p<0.05) across stages underpins the therapeutic-target conclusion, yet the abstract supplies no description of the cell-cell communication tool, statistical test, multiple-testing correction, or how activity was quantified across the inferred trajectory.

Authors: We acknowledge that the abstract and methods lack sufficient detail on the signaling pathway analysis. In the revised version, we will expand the Methods section to describe the cell-cell communication tool employed (CellChat), the statistical framework for detecting active pathways, the multiple-testing correction applied (Benjamini-Hochberg), and the quantification of pathway activity scores across trajectory stages. The abstract will be updated to briefly mention the use of CellChat for inferring ligand-receptor interactions and the significance threshold used. revision: yes

Circularity Check

0 steps flagged

No significant circularity; empirical data analysis only

full rationale

The paper reports results from applying standard scRNA-seq clustering and trajectory inference to adipose tissue samples. No mathematical derivations, equations, or parameter-fitting steps are described that reduce predictions to inputs by construction. No self-citations, ansatzes, or uniqueness theorems are invoked as load-bearing elements in the provided text. Claims about 15 clusters, 7 states, and IGF/FGF pathways rest on data processing rather than self-referential definitions or fitted-input renamings. This is the expected non-finding for a purely empirical bioinformatics study.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

Only the abstract is available; no explicit free parameters, axioms, or invented entities are described in the text. The p<0.05 threshold is a standard statistical choice, not a fitted parameter unique to this work.

pith-pipeline@v0.9.1-grok · 5752 in / 1204 out tokens · 45136 ms · 2026-06-29T02:23:36.707729+00:00 · methodology

discussion (0)

Reference graph

Works this paper leans on

2 extracted references · 1 canonical work pages

[1]

The workflow begins with raw data processing and progresses through sequential analytical stages to uncover the molecular programs governing adipogenesis

Cell-Cell Communication Analysis: Ligand-receptor interactions (CellChatDB.human v1.6.0) quantify signaling networks (e.g., IGF1-IGF1R, FGF pathways) between adipocytes (ADIPOQ+/PLIN1+) and ASPCs (PDGFRA+/DCN+) [9]. The workflow begins with raw data processing and progresses through sequential analytical stages to uncover the molecular programs governing ...
[2]

harmony_1

Feature Selection by top genes To reduce dimensionality, 2,000 highly variable genes (HVGs) were selected from Figure 2b, which plots standardized variance (y-axis) against average expression (x-axis, 1e−04 to 1e+01). The majority of genes (29,833 non-variable) cluster in the low-variance region (red points), while 2,000 variable genes (black-labeled, e.g...

work page doi:10.1172/jci45887.2094 2025

[1] [1]

The workflow begins with raw data processing and progresses through sequential analytical stages to uncover the molecular programs governing adipogenesis

Cell-Cell Communication Analysis: Ligand-receptor interactions (CellChatDB.human v1.6.0) quantify signaling networks (e.g., IGF1-IGF1R, FGF pathways) between adipocytes (ADIPOQ+/PLIN1+) and ASPCs (PDGFRA+/DCN+) [9]. The workflow begins with raw data processing and progresses through sequential analytical stages to uncover the molecular programs governing ...

[2] [2]

harmony_1

Feature Selection by top genes To reduce dimensionality, 2,000 highly variable genes (HVGs) were selected from Figure 2b, which plots standardized variance (y-axis) against average expression (x-axis, 1e−04 to 1e+01). The majority of genes (29,833 non-variable) cluster in the low-variance region (red points), while 2,000 variable genes (black-labeled, e.g...

work page doi:10.1172/jci45887.2094 2025