Matrix factorization analyses of the immgenT scRNA-seq data

workflowr

• Summary
• Checks
• Past versions

Last updated: 2026-08-27

Checks: 2 0

Knit directory: immgenT-GP-analysis/analysis/

This reproducible R Markdown analysis was created with workflowr (version 1.7.2). The Checks tab describes the reproducibility checks that were applied when the results were created. The Past versions tab lists the development history.

---

R Markdown file: up-to-date

Great! Since the R Markdown file has been committed to the Git repository, you know the exact version of the code that produced these results.

Repository version: 659ff00

Great! You are using Git for version control. Tracking code development and connecting the code version to the results is critical for reproducibility.

The results in this page were generated with repository version 659ff00. See the Past versions tab to see a history of the changes made to the R Markdown and HTML files.

Note that you need to be careful to ensure that all relevant files for the analysis have been committed to Git prior to generating the results (you can use wflow_publish or wflow_git_commit). workflowr only checks the R Markdown file, but you know if there are other scripts or data files that it depends on. Below is the status of the Git repository when the results were generated:

Ignored files:
    Ignored:    .DS_Store
    Ignored:    .claude/
    Ignored:    analysis/.DS_Store
    Ignored:    analysis/.Rhistory
    Ignored:    analysis/assets/.DS_Store
    Ignored:    captions/
    Ignored:    code/.DS_Store
    Ignored:    code/other/topic_flashier_20250212.R
    Ignored:    code/other/topic_wrapper_20250215_alldata_backfit.sh
    Ignored:    data
    Ignored:    experiments/
    Ignored:    figures/.DS_Store
    Ignored:    figures/Previous/.DS_Store
    Ignored:    figures/Previous/bits/.DS_Store
    Ignored:    figures/Previous/bits/Figure 1/.DS_Store
    Ignored:    figures/Previous/bits/Figure 2/.DS_Store
    Ignored:    figures/Previous/bits/Figure 3/.DS_Store
    Ignored:    figures/Previous/bits/Figure 4/.DS_Store
    Ignored:    figures/Previous/bits/Figure 6/.DS_Store
    Ignored:    figures/Previous/bits/Figure 7/.DS_Store
    Ignored:    figures/Previous/bits/Figure S1/.DS_Store
    Ignored:    figures/Previous/bits/Figure S2/.DS_Store
    Ignored:    figures/Previous/bits/Figure S3/.DS_Store
    Ignored:    figures/Previous/bits/Figure S6/.DS_Store
    Ignored:    figures/Previous/bits/Figure S7/.DS_Store
    Ignored:    figures/final-selected/.DS_Store
    Ignored:    figures/final-selected/Figure 1/.DS_Store
    Ignored:    figures/final-selected/Figure 2/.DS_Store
    Ignored:    figures/final-selected/Figure 4/.DS_Store
    Ignored:    figures/final-selected/Figure S1/.DS_Store
    Ignored:    figures/final-selected/Figure S4/.DS_Store
    Ignored:    figures/templates_20260729/
    Ignored:    log/
    Ignored:    output/.DS_Store
    Ignored:    output/Figure2/
    Ignored:    output/Figure7b/7b_cell_metadata.csv.gz
    Ignored:    plan/
    Ignored:    tables/
    Ignored:    tmp/

Note that any generated files, e.g. HTML, png, CSS, etc., are not included in this status report because it is ok for generated content to have uncommitted changes.

---

These are the previous versions of the repository in which changes were made to the R Markdown (analysis/index.Rmd) and HTML (docs/index.html) files. If you’ve configured a remote Git repository (see ?wflow_git_remote), click on the hyperlinks in the table below to view the files as they were in that past version.

File	Version	Author	Date	Message
Rmd	659ff00	Ziang Zhang	2026-08-27	Extended Data Figure 8: per-lineage cluster structure plots
html	dd76d38	Ziang Zhang	2026-08-19	Build site: Figure 7 and index pages rebuilt
Rmd	801dbf7	Ziang Zhang	2026-08-19	Figure 7 pages: present Figure 7 as one figure spanning EBMF and RQVI
html	eeca07b	Ziang Zhang	2026-08-05	Keep pre-refactor provenance in panel comments off the published pages
Rmd	5651d0e	Ziang Zhang	2026-08-05	Extended Data tables: reorder to six, rebuild Table 1, drop internal notes
html	dd71400	Ziang Zhang	2026-07-31	Build site: Methods section on cell filtering by total GP membership
Rmd	3f88d67	Ziang Zhang	2026-07-31	Methods: document how L_pm_filtered / F_pm_filtered are obtained
html	6926d1b	Ziang Zhang	2026-07-31	Build site: Extended Data Table 1 annotations + tabbed top-15 version
Rmd	80f61c0	Ziang Zhang	2026-07-31	Extended Data Table 1: manual annotations + longer internal version
html	cbcec52	Ziang Zhang	2026-07-30	Build site: Extended Data Figure naming
Rmd	66aa029	Ziang Zhang	2026-07-30	Name the Extended Data figures as published on the site
html	ac650a0	Ziang Zhang	2026-07-30	Build site: Figure S5 (ex-S6a) and Figure S6 as a-f
Rmd	c9b020f	Ziang Zhang	2026-07-30	Split Figure S6’s protein-program heatmap out as Figure S5
html	732ac8d	Ziang Zhang	2026-07-29	Build site: caption alignment with captions_20260729_final.docx
Rmd	8d73953	Ziang Zhang	2026-07-28	Align captions with captions_20260729_final.docx
html	ae21d37	Ziang Zhang	2026-07-28	Build site: republish after the reorder commits
html	d538aa2	Ziang Zhang	2026-07-28	Build site: reordered Figures 6 / S6 / S3 and the new Figure 7b page
Rmd	4c07670	Ziang Zhang	2026-07-28	Reorder Figures 6, S6 and S3; make the ex-S5 figure Figure 7b
html	029b0ae	Ziang Zhang	2026-07-28	Build site.
Rmd	0f5b5da	Ziang Zhang	2026-07-28	Align all figure captions with captions_20260728_final.docx
html	d0ebed0	Ziang Zhang	2026-07-28	Build site: Extended Data Table 7 as the manual protein threshold list
Rmd	cf77f4e	Ziang Zhang	2026-07-28	Extended Data Table 7: swap GP gating alignment for the manual protein thresholds
html	3fc3789	Ziang Zhang	2026-07-27	Republish all 24 pages
html	1390a03	Ziang Zhang	2026-07-27	Republish all 24 pages
html	adaef21	Ziang Zhang	2026-07-27	Build site: panel fixes and PDF-derived assets
html	5b19858	Ziang Zhang	2026-07-27	Build site.
html	9862b6d	Ziang Zhang	2026-07-23	Build site.
Rmd	98d2924	Ziang Zhang	2026-07-23	Reword Figure S5 page for a publication audience
html	9398c72	Ziang Zhang	2026-07-23	Publish Figure S5 workflowr page
Rmd	b9f4f58	Ziang Zhang	2026-07-23	Add Figure S5: EBMF vs matched-RQVI level2-cluster comparison
Rmd	b197499	Ziang Zhang	2026-07-16	Restructure main figures (renumber + split Figure 1)
html	b197499	Ziang Zhang	2026-07-16	Restructure main figures (renumber + split Figure 1)
Rmd	f0c481d	Ziang Zhang	2026-07-14	Document CITE-seq fixed-loading EBMF workflow
html	f0c481d	Ziang Zhang	2026-07-14	Document CITE-seq fixed-loading EBMF workflow
Rmd	c344aaa	Ziang Zhang	2026-07-13	Add Figure S4 mean-loading heatmaps
html	c344aaa	Ziang Zhang	2026-07-13	Add Figure S4 mean-loading heatmaps
Rmd	cba8dd4	Ziang Zhang	2026-07-12	Add Methods page: Flashier matrix factorization
html	cba8dd4	Ziang Zhang	2026-07-12	Add Methods page: Flashier matrix factorization
html	c5cbc77	Ziang Zhang	2026-07-08	Build site.
Rmd	fb8c855	Ziang Zhang	2026-07-08	Add Extended Data Table 8: comprehensive gene signature matrix (long + wide)
html	92d1365	Ziang Zhang	2026-07-05	Build site.
Rmd	d07f797	Ziang Zhang	2026-07-05	Replace Table S1/S2/S3 with 7 fully-reproducible Extended Data Tables
html	92021bf	Ziang Zhang	2026-07-02	Build site.
html	827c89b	Ziang Zhang	2026-07-02	Build site.
Rmd	f9db962	Ziang Zhang	2026-07-02	Simplify layout: drop old code/script folders, rename
html	c6e5086	Ziang Zhang	2026-07-02	Build site.
html	5a79883	Ziang Zhang	2026-07-02	Build site.
Rmd	2b0e445	Ziang Zhang	2026-07-02	Fix GitHub source links to point at the new
html	cf1d0ac	Ziang Zhang	2026-07-02	Build site.
Rmd	06b2461	Ziang Zhang	2026-07-02	Initial commit: immgenT-GP-analysis
html	06b2461	Ziang Zhang	2026-07-02	Initial commit: immgenT-GP-analysis

---

This site documents the gene program (GP) factorization analysis of the ImmGen-T single-cell RNA-seq/CITE-seq data: ~200 GPs learned by empirical Bayes matrix factorization (EBMF), and how they relate to cell lineage, tissue, activation state, transcription factors, and surface protein expression.

Each page below covers one published figure: for every panel, the R code that produced it, the figure itself, and a short caption explaining what it shows.

Main figures

• Figure 1: Overview of the 200 gene programs in immgenT
• Figure 2: Active genes and cells in gene programs
• Figure 3: T cell lineage-related GPs
• Figure 4: GPs associated with T cell activation
• Figure 5: GP expression in tissues is driven by T cell composition
• Figure 6: Surface proteins associated with GP activity
• Figure 7b: Cluster-level activity of EBMF GPs and their matched RQVI GPs

Figure 7 (RQVI validates EBMF-derived gene programs with a non-linear deep learning framework) brings the EBMF programs documented here together with RQVI, a non-linear deep-learning factorization of the same cells. Panel b, the cluster-level comparison of the two factorizations, is documented on the page above; the RQVI panels are documented in the companion repository for the RQVI framework.

Extended Data figures

Extended Data Figures 1-8. Each page’s URL keeps the FigureS* form, but every figure is named the way it is published.

• Extended Data Figure 1: Reproducibility across IGTs and batch-effect evaluation
• Extended Data Figure 2: Tz- and CD8-specific GPs (GP30 and GP58)
• Extended Data Figure 3: Characterizing activation-associated GPs
• Extended Data Figure 4: Healthy non-thymocyte GP mean-loading heatmaps
• Extended Data Figure 5: Protein programs across GPs
• Extended Data Figure 6: CD69-associated GPs and further surface-protein gating strategies

Extended Data Figure 7 (RQVI scalability and representative matches to lineage-associated EBMF GPs) is documented in the same companion repository for the RQVI framework.

• Extended Data Figure 8: Cluster-level GP membership within each T cell lineage

Extended Data tables

• Extended Data Table 1: summary of GP characteristics and comments
• Extended Data Table 2: gene signature per GP
• Extended Data Table 3: GP AUC by lineage
• Extended Data Table 4: GP activity during activation
• Extended Data Table 5: GP AUC by tissue
• Extended Data Table 6: GP AUC by cluster

Methods

• Fitting EBMF using flashier (including how L_pm_filtered / F_pm_filtered are obtained by filtering cells on total GP membership)
• Fitting CITE-seq protein programs with fixed scRNA loadings

Reproducing these figures

Every page above is generated by a single R script in script/, which sources shared plotting/data-loading helpers from code/R/. To regenerate a figure yourself, run (from the repository root):

Rscript script/Figure1.R