Abstract: Bulk RNA sequencing (RNA-seq) averages signals across heterogeneous cell populations. Computational deconvolution methods aim to infer cell type composition and cell type-specific gene expression from bulk data, but their performance in equine samples has not been evaluated. In this study, we assessed the ability of computational deconvolution to recover cellular composition and differential expression signals in bronchoalveolar lavage fluid (BALF) from horses with severe equine asthma (SEA) and controls (CTL). Cryopreserved BALF samples from six SEA and five CTL horses previously analyzed by scRNA-seq were used to generate bulk RNA-seq data. The matched scRNA-seq dataset served as the reference for deconvolution. Performance was evaluated by comparing deconvolution raw and mRNA-corrected estimates with scRNA-seq cell proportions. Differential expression between SEA and CTL was analyzed on bulk RNA-seq, deconvoluted expression profiles, and scRNA-seq pseudobulk data. Deconvolution primarily captured mRNA-derived cell type proportions rather than true cell counts: agreement with scRNA-seq cell counts was moderate (r = 0.62; 95% CI 0.45-0.75) but improved after mRNA content correction (r = 0.83; 95% CI 0.74-0.89). Comparison with mRNA-weighted scRNA-seq proportions showed near-perfect concordance (r = 0.98; 95% CI 0.97-0.99). Cell type-specific performance varied, with stronger correlations for B cells and dendritic cells and weaker performance for neutrophils, T cells and monocytes/macrophages. Recovery of cell type-specific differential expression was inconsistent, frequently showing cross-lineage signal spillover. Although both approaches detected a Th17 signature in SEA, most deconvolution-derived differentially expressed genes overlapped with conventional bulk RNA-seq results. Deconvolution of bulk RNA-seq did not reliably estimate cell counts or provide substantial biological insight beyond conventional bulk analysis, highlighting the value of scRNA-seq for resolving cell type-specific disease mechanisms in equine asthma.
The Equine Research Bank provides access to a large database of publicly available scientific literature. Inclusion in the Research Bank does not imply endorsement of study methods or findings by Mad Barn.
This research summary has been generated with artificial intelligence and may contain errors and omissions. Refer to the original study to confirm details provided. Submit correction.
Overview
This study evaluates computational methods to estimate cell type composition and gene expression from bulk RNA sequencing data of horse lung fluid, comparing these estimates to matched single-cell RNA sequencing data.
The research focuses on horses with severe equine asthma and healthy controls, assessing how well bulk data deconvolution reveals cell types and disease-associated gene expression patterns.
Background
Bulk RNA-seq vs. Single-cell RNA-seq: Bulk RNA-seq measures gene expression averaged across all cells in a sample, losing information about individual cell types.
Computational Deconvolution: Algorithms attempt to separate bulk RNA-seq data into contributions from different cell types to infer cellular composition and cell type-specific gene expression.
Equine BALF: Bronchoalveolar lavage fluid obtained from horse lungs contains a mixture of immune and lung cells relevant for studying equine asthma.
Equine Asthma Context: Severe equine asthma (SEA) is an inflammatory lung disease in horses that shares some features with human asthma, making cell type-specific expression analysis important.
Methods
Sample Collection: BALF samples were collected from 6 horses with SEA and 5 control horses.
Data Generation: The BALF samples were cryopreserved and subjected to both bulk RNA-seq and matched single-cell RNA-seq (scRNA-seq) for detailed cellular profiling.
Deconvolution Algorithms: Computational methods were applied to the bulk RNA-seq data using the scRNA-seq dataset as a reference to estimate:
Cell type proportions
Cell type-specific gene expression profiles
Evaluation of Performance: The estimated cell proportions from deconvolution were compared to actual cell proportions observed in scRNA-seq both as raw cell counts and after correcting for mRNA content per cell type.
Differential Expression Analysis: Compared bulk RNA-seq, deconvoluted expression data, and scRNA-seq pseudobulk data to identify genes differing between SEA and controls.
Key Findings
Cell Proportion Estimation:
Deconvolution better estimated the mRNA content contribution of cell types rather than actual cell numbers.
Moderate correlation (r = 0.62) was found between raw deconvolution estimates and true scRNA-seq cell counts.
Correlation improved substantially to r = 0.83 after correcting for differences in mRNA content per cell type.
Near-perfect correlation (r = 0.98) was observed when comparing deconvolution outputs to mRNA-weighted cell proportions from scRNA-seq.
Variation by Cell Type:
Stronger performance in estimating proportions for B cells and dendritic cells.
Poorer estimates for neutrophils, T cells, and monocytes/macrophages, likely because of variable mRNA content or overlapping gene expression profiles.
Differential Expression and Biological Insights:
Cell-type specific differential expression signals were inconsistently recovered from deconvoluted data.
Cross-lineage signal spillover occurred, where gene expression changes attributed to one cell type appeared in others incorrectly.
Both scRNA-seq and deconvolution detected a Th17 immune response signature characteristic of SEA.
Most differentially expressed genes identified by deconvolution overlapped with those from traditional bulk RNA-seq, indicating limited added value.
Conclusions and Implications
Deconvolution of bulk RNA-seq, while useful for estimating mRNA-weighted cell type proportions, does not reliably estimate actual cell counts in equine BALF.
It also fails to consistently recover accurate cell type-specific gene expression patterns and disease-associated changes.
Single-cell RNA-seq remains a more powerful and reliable approach for understanding cell type-specific disease mechanisms in equine asthma.
This highlights limitations in applying current computational deconvolution methods to complex immune cell mixtures in veterinary samples and supports prioritizing single-cell technologies for future equine respiratory disease research.
Cite This Article
APA
Jagannathan V, Leeb T, Gerber V, Sage SE.
(2026).
Cell-Type Deconvolution of Equine BALF RNA-Seq: A Critical Comparison with Matched Single-Cell Data.
Genes (Basel), 17(7), 773.
https://doi.org/10.3390/genes17070773
Luca B.A., Steen C.B., Matusiak M., Azizi A., Varma S., Zhu C., Przybyl J., Espín-Pérez A., Diehn M., Alizadeh A.A.. Atlas of Clinically Distinct Cell States and Ecosystems across Human Solid Tumors. Cell 2021;184:5482–5496.e28.
Chu T., Wang Z., Pe’er D., Danko C.G.. Cell Type and Gene Expression Deconvolution with BayesPrism Enables Bayesian Integrative Analysis across Bulk and Single-Cell RNA Sequencing in Oncology. Nat. Cancer 2022;3:505–517.
Conning-Rowland M., Cheng C.W., Brown O., Giannoudi M., Levelt E., Roberts L.D., Griffin K.J., Cubbon R.M.. Application of CIBERSORTx and BayesPrism to Deconvolution of Bulk RNA-Seq Data from Human Myocardium and Skeletal Muscle. Heliyon 2025;11:e42499.
Liao Y., Smyth G.K., Shi W.. featureCounts: An Efficient General Purpose Program for Assigning Sequence Reads to Genomic Features. Bioinformatics 2014;30:923–930.
Racle J, de Jonge K, Baumgaertner P, Speiser DE, Gfeller D. Simultaneous Enumeration of Cancer and Immune Cell Types from Bulk Tumor Gene Expression Data. eLife 2017;6:e26476.
Finotello F, Mayer C, Plattner C, Laschober G, Rieder D, Hackl H, Krogsdam A, Loncova Z, Posch W, Wilflingseder D. Molecular and Pharmacological Modulators of the Tumor Immune Contexture Revealed by Deconvolution of RNA-Seq Data. Genome Med 2019;11:34.
Zaitsev K, Bambouskova M, Swain A, Artyomov MN. Complete Deconvolution of Cellular Mixtures Based on Linearity of Transcriptional Signatures. Nat Commun 2019;10:2209.
Newman AM, Steen CB, Liu CL, Gentles AJ, Chaudhuri AA, Scherer F, Khodadoust MS, Esfahani MS, Luca BA, Steiner D. Determining Cell Type Abundance and Expression from Bulk Tissues with Digital Cytometry. Nat Biotechnol 2019;37:773–782.
He L, Davila-Velderrain J, Sumida TS, Hafler DA, Kellis M, Kulminski AM. NEBULA Is a Fast Negative Binomial Mixed Model for Differential or Co-Expression Analysis of Large-Scale Multi-Subject Single-Cell Data. Commun Biol 2021;4:629.
Alfei F, Kanev K, Hofmann M, Wu M, Ghoneim HE, Roelli P, Utzschneider DT, von Hoesslin M, Cullen JG, Fan Y. TOX Reinforces the Phenotype and Longevity of Exhausted T Cells in Chronic Viral Infection. Nature 2019;571:265–269.
Wing JB, Ise W, Kurosaki T, Sakaguchi S. Regulatory T Cells Control Antigen-Specific Expansion of Tfh Cell Number and Humoral Immune Responses via the Coreceptor CTLA-4. Immunity 2014;41:1013–1025.
Zotos D, Coquet JM, Zhang Y, Light A, D’Costa K, Kallies A, Corcoran LM, Godfrey DI, Toellner K-M, Smyth MJ. IL-21 Regulates Germinal Center B Cell Differentiation and Proliferation through a B Cell–Intrinsic Mechanism. J Exp Med 2010;207:365–378.
Joshi S, Pantalena L-C, Liu XK, Gaffen SL, Liu H, Rohowsky-Kochan C, Ichiyama K, Yoshimura A, Steinman L, Christakos S. 1,25-Dihydroxyvitamin D 3 Ameliorates Th17 Autoimmunity via Transcriptional Modulation of Interleukin-17A. Mol Cell Biol 2011;31:3653–3669.
de la Fuente H, Cruz-Adalia A, Martinez del Hoyo G, Cibrián-Vera D, Bonay P, Pérez-Hernández D, Vázquez J, Navarro P, Gutierrez-Gallego R, Ramirez-Huesca M. The Leukocyte Activation Receptor CD69 Controls T Cell Differentiation through Its Interaction with Galectin-1. Mol Cell Biol 2014;34:2479–2487.
Miyazaki H, Wannakul T, Yang S, Yang D, Karasawa A, Shishido A, Cao R, Yamamoto Y, Kagawa Y, Kobayashi S. FABP7 in Hepatic Macrophages Promotes Fibroblast Activation and CD4+ T-Cell Migration by Regulating M2 Polarization During Liver Fibrosis. J Immunol Res 2025;2025:6987981.
Shimizu Y, Dobashi K. CC-Chemokine CCL15 Expression and Possible Implications for the Pathogenesis of IgE-Related Severe Asthma. Mediat. Inflamm. 2012;2012:475253.
Prikk K, Maisi P, Pirilä E, Reintam M-A, Salo T, Sorsa T, Sepper R. Airway Obstruction Correlates with Collagenase-2 (MMP-8) Expression and Activation in Bronchial Asthma. Lab. Investig. 2002;82:1535–1545.
Skerlavaj B, Scocchi M, Gennaro R, Risso A, Zanetti M. Structural and Functional Analysis of Horse Cathelicidin Peptides. Antimicrob. Agents Chemother. 2001;45:715–722.
Han Y, Ji X, Rao J, Zhang Y, Chen X, Gong F. TNFAIP3, PADI4, and CXCR4 as Biomarkers of Neutrophil Extracellular Traps in Severe Asthma. Int. Immunopharmacol. 2026;170:116118.
Jiang S, Park D.W, Tadie J-M, Gregoire M, Deshane J, Pittet J.F, Abraham E, Zmijewski J.W. Human Resistin Promotes Neutrophil Proinflammatory Activation and Neutrophil Extracellular Trap Formation and Increases Severity of Acute Lung Injury. J. Immunol. 2014;192:4795–4803.
Quoc Q.L, Choi Y, Thi Bich T.C, Yang E-M, Shin Y.S, Park H-S. S100A9 in Adult Asthmatic Patients: A Biomarker for Neutrophilic Asthma. Exp. Mol. Med. 2021;53:1170–1179.
Sprenkeler E.G.G, Zandstra J, van Kleef N.D, Goetschalckx I, Verstegen B, Aarts C.E.M, Janssen H, Tool A.T.J, van Mierlo G, van Bruggen R. S100A8/A9 Is a Marker for the Release of Neutrophil Extracellular Traps and Induces Neutrophil Activation. Cells 2022;11:236.
Cremades-Jimeno L, de Pedro M.Á, López-Ramos M, Sastre J, Mínguez P, Fernández I.M, Baos S, Cárdaba B. Prioritizing Molecular Biomarkers in Asthma and Respiratory Allergy Using Systems Biology. Front. Immunol. 2021;12:1004.
Liu L, Zhang X, Liu Y, Zhang L, Zheng J, Wang J, Hansbro P.M, Wang L, Wang G, Hsu A.C.Y. Chitinase-like Protein YKL-40 Correlates with Inflammatory Phenotypes, Anti-Asthma Responsiveness and Future Exacerbations. Respir. Res. 2019;20:95.
Taketomi Y, Sugiki T, Saito T, Ishii S-I, Hisada M, Suzuki-Nishimura T, Uchida M.K, Moon T-C, Chang H-W, Natori Y. Identification of NDRG1 as an Early Inducible Gene during In Vitro Maturation of Cultured Mast Cells. Biochem. Biophys. Res. Commun. 2003;306:339–346.
Lu S, Yang J, Yan L, Liu J, Wang J.J, Jain R, Yu J. Transcriptome Size Matters for Single-Cell RNA-Seq Normalization and Bulk Deconvolution. Nat. Commun. 2025;16:1246.
Avila Cobos F, Alquicira-Hernandez J, Powell J.E, Mestdagh P, De Preter K. Benchmarking of Cell Type Deconvolution Pipelines for Transcriptomics Data. Nat. Commun. 2020;11:5650.
Dietrich A, Merotto L, Pelz K, Eder B, Zackl C, Reinisch K, Edenhofer F, Marini F, Sturm G, List M. Omnideconv: A Unifying Framework for Using and Benchmarking Single-Cell-Informed Deconvolution of Bulk RNA-Seq Data. Genome Biol. 2026;27:6.
Meng G, Tang W, Huang E, Li Z, Feng H. A Comprehensive Assessment of Cell Type-Specific Differential Expression Methods in Bulk Data. Brief. Bioinform. 2023;24:bbac516.
Madissoon E, Wilbrey-Clark A, Miragaia R.J, Saeb-Parsy K, Mahbubani K.T, Georgakopoulos N, Harding P, Polanski K, Huang N, Nowicki-Osuch K. ScRNA-Seq Assessment of the Human Lung, Spleen, and Esophagus Tissue Stability after Cold Preservation. Genome Biol. 2020;21:1.
Guillaumet-Adkins A, Rodríguez-Esteban G, Mereu E, Mendez-Lago M, Jaitin D.A, Villanueva A, Vidal A, Martinez-Marti A, Felip E, Vivancos A. Single-Cell Transcriptome Conservation in Cryopreserved Cells and Tissues. Genome Biol. 2017;18:45.
Denisenko E, Guo B.B, Jones M, Hou R, de Kock L, Lassmann T, Poppe D, Clément O, Simmons R.K, Lister R. Systematic Assessment of Tissue Dissociation and Storage Biases in Single-Cell and Single-Nucleus RNA-Seq Workflows. Genome Biol. 2020;21:130.
Browne D.J, Miller C.M, Doolan D.L. Technical Pitfalls When Collecting, Cryopreserving, Thawing, and Stimulating Human T-Cells. Front. Immunol. 2024;15:1382192.