Fast, Flexible Analysis of Differences in Cellular Composition with Crumblr

Overview

Journal bioRxiv

Date 2025 Feb 20

PMID 39975411

Authors

Gabriel E Hoffman

Panos Roussos

Affiliations

Soon will be listed here.

Abstract

Changes in cell type composition play an important role in human health and disease. Recent advances in single-cell technology have enabled the measurement of cell type composition at increasing cell lineage resolution across large cohorts of individuals. Yet this raises new challenges for statistical analysis of these compositional data to identify changes in cell type frequency. We introduce crumblr (DiseaseNeurogenomics.github.io/crumblr), a scalable statistical method for analyzing count ratio data using precision-weighted linear mixed models incorporating random effects for complex study designs. Uniquely, crumblr performs statistical testing at multiple levels of the cell lineage hierarchy using a multivariate approach to increase power over tests of one cell type. In simulations, crumblr increases power compared to existing methods while controlling the false positive rate. We demonstrate the application of crumblr to published single-cell RNA-seq datasets for aging, tuberculosis infection in T cells, bone metastases from prostate cancer, and SARS-CoV-2 infection.

References

Phipson B, Sim C, Porrello E, Hewitt A, Powell J, Oshlack A . propeller: testing for differences in cell type proportions in single cell data. Bioinformatics. 2022; 38(20):4720-4726. PMC: 9563678. DOI: 10.1093/bioinformatics/btac582. View

Buttner M, Ostner J, Muller C, Theis F, Schubert B . scCODA is a Bayesian model for compositional single-cell data analysis. Nat Commun. 2021; 12(1):6876. PMC: 8616929. DOI: 10.1038/s41467-021-27150-6. View

Love M, Huber W, Anders S . Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014; 15(12):550. PMC: 4302049. DOI: 10.1186/s13059-014-0550-8. View

Mogilenko D, Shchukina I, Artyomov M . Immune ageing at single-cell resolution. Nat Rev Immunol. 2021; 22(8):484-498. PMC: 8609266. DOI: 10.1038/s41577-021-00646-4. View

Yazar S, Alquicira-Hernandez J, Wing K, Senabouth A, Gordon M, Andersen S . Single-cell eQTL mapping identifies cell type-specific genetic control of autoimmune disease. Science. 2022; 376(6589):eabf3041. DOI: 10.1126/science.abf3041. View

Law C, Chen Y, Shi W, Smyth G . voom: Precision weights unlock linear model analysis tools for RNA-seq read counts. Genome Biol. 2014; 15(2):R29. PMC: 4053721. DOI: 10.1186/gb-2014-15-2-r29. View

. A blood atlas of COVID-19 defines hallmarks of disease severity and specificity. Cell. 2022; 185(5):916-938.e58. PMC: 8776501. DOI: 10.1016/j.cell.2022.01.012. View

Hoffman G, Schadt E . variancePartition: interpreting drivers of variation in complex gene expression studies. BMC Bioinformatics. 2016; 17(1):483. PMC: 5123296. DOI: 10.1186/s12859-016-1323-z. View

Braun E, Danan-Gotthold M, Borm L, Lee K, Vinsland E, Lonnerberg P . Comprehensive cell atlas of the first-trimester developing human brain. Science. 2023; 382(6667):eadf1226. DOI: 10.1126/science.adf1226. View

10.

Emani P, Liu J, Clarke D, Jensen M, Warrell J, Gupta C . Single-cell genomics and regulatory networks for 388 human brains. Science. 2024; 384(6698):eadi5199. PMC: 11365579. DOI: 10.1126/science.adi5199. View

11.

Lin D, Sullivan P . Meta-analysis of genome-wide association studies with overlapping subjects. Am J Hum Genet. 2009; 85(6):862-72. PMC: 2790578. DOI: 10.1016/j.ajhg.2009.11.001. View

12.

Rood J, Maartens A, Hupalowska A, Teichmann S, Regev A . Impact of the Human Cell Atlas on medicine. Nat Med. 2022; 28(12):2486-2496. DOI: 10.1038/s41591-022-02104-7. View

13.

Filippov I, Schauser L, Peterson P . An integrated single-cell atlas of blood immune cells in aging. NPJ Aging. 2024; 10(1):59. PMC: 11606963. DOI: 10.1038/s41514-024-00185-x. View

14.

Terekhova M, Swain A, Bohacova P, Aladyeva E, Arthur L, Laha A . Single-cell atlas of healthy human blood unveils age-related loss of NKG2CGZMBCD8 memory T cells and accumulation of type 2 memory T cells. Immunity. 2023; 56(12):2836-2854.e9. DOI: 10.1016/j.immuni.2023.10.013. View

15.

Ritchie M, Phipson B, Wu D, Hu Y, Law C, Shi W . limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015; 43(7):e47. PMC: 4402510. DOI: 10.1093/nar/gkv007. View

16.

Kfoury Y, Baryawno N, Severe N, Mei S, Gustafsson K, Hirz T . Human prostate cancer bone metastases have an actionable immunosuppressive microenvironment. Cancer Cell. 2021; 39(11):1464-1478.e8. PMC: 8578470. DOI: 10.1016/j.ccell.2021.09.005. View

17.

Hoffman G, Roussos P . Dream: powerful differential expression analysis for repeated measures designs. Bioinformatics. 2020; 37(2):192-201. PMC: 8055218. DOI: 10.1093/bioinformatics/btaa687. View

18.

Mittelbrunn M, Kroemer G . Hallmarks of T cell aging. Nat Immunol. 2021; 22(6):687-698. DOI: 10.1038/s41590-021-00927-z. View

19.

Han B, Duong D, Sul J, de Bakker P, Eskin E, Raychaudhuri S . A general framework for meta-analyzing dependent studies with overlapping subjects in association mapping. Hum Mol Genet. 2016; 25(9):1857-66. PMC: 4986332. DOI: 10.1093/hmg/ddw049. View

20.

Halabi S, Kelly W, Ma H, Zhou H, Solomon N, Fizazi K . Meta-Analysis Evaluating the Impact of Site of Metastasis on Overall Survival in Men With Castration-Resistant Prostate Cancer. J Clin Oncol. 2016; 34(14):1652-9. PMC: 4872320. DOI: 10.1200/JCO.2015.65.7270. View