Score Matching for Differential Abundance Testing of Compositional High‐Throughput Sequencing Data

Fuente: Wiley Open Access
Gespeichert in:
Bibliographische Detailangaben
Hauptverfasser: Johannes Ostner, Hongzhe Li, Christian L. Müller
Format: Artículo Open Access
Veröffentlicht: Wiley 2026
Schlagworte:
Online-Zugang:
Tags: Tag hinzufügen
Keine Tags, Fügen Sie den ersten Tag hinzu!
_version_ 1867004291565223936
author Johannes Ostner
Hongzhe Li
Christian L. Müller
author_facet Johannes Ostner
Hongzhe Li
Christian L. Müller
Johannes Ostner
Hongzhe Li
Christian L. Müller
collection Wiley Open Access
contents Score Matching for Differential Abundance Testing of Compositional High‐Throughput Sequencing Data Johannes Ostner Hongzhe Li Christian L. Müller Statistics in Medicine ABSTRACT The class of a‐b power interaction models, proposed by [1], provides a general framework for modeling sparse compositional data with pairwise feature interactions. This class includes many distributions as special cases and enables modeling of zero entries through power transformations, making it particularly suitable for modern high‐throughput sequencing data with excess zeros, including single‐cell RNA‐Seq and microbial amplicon data. Here, we present an extension of this class of models that allows inclusion of covariate information, thus enabling accurate characterization of covariate dependencies in heterogeneous populations. Combining this model with a tailored differential abundance (DA) test leads to a novel DA testing scheme, cosmoDA , that can reduce the false positive detection rate caused by correlated features. cosmoDA uses penalized generalized score matching for parsimonious model fitting. We show on simulated benchmarks that cosmoDA can accurately estimate feature interactions in the presence of population heterogeneity and significantly reduces the false discovery rate when testing for differential abundance of correlated features. Using single‐cell and amplicon data, we illustrate cosmoDA 's ability to estimate data‐adaptive Box–Cox‐type data transformations and assess the impact of zero replacement and power transformations on downstream differential abundance results. cosmoDA is available at https://github.com/bio‐datascience/cosmoDA . 10.1002/sim.70534 http://creativecommons.org/licenses/by/4.0/
doi_str_mv 10.1002/sim.70534
format Artículo Open Access
id wiley_oa_10_1002_sim_70534
institution Wiley Open Access
license_str_mv http://creativecommons.org/licenses/by/4.0/
publishDate 2026
publisher Wiley
record_format wiley_oa
spellingShingle Score Matching for Differential Abundance Testing of Compositional High‐Throughput Sequencing Data
Johannes Ostner
Hongzhe Li
Christian L. Müller
Statistics in Medicine
Score Matching for Differential Abundance Testing of Compositional High‐Throughput Sequencing Data Johannes Ostner Hongzhe Li Christian L. Müller Statistics in Medicine ABSTRACT The class of a‐b power interaction models, proposed by [1], provides a general framework for modeling sparse compositional data with pairwise feature interactions. This class includes many distributions as special cases and enables modeling of zero entries through power transformations, making it particularly suitable for modern high‐throughput sequencing data with excess zeros, including single‐cell RNA‐Seq and microbial amplicon data. Here, we present an extension of this class of models that allows inclusion of covariate information, thus enabling accurate characterization of covariate dependencies in heterogeneous populations. Combining this model with a tailored differential abundance (DA) test leads to a novel DA testing scheme, cosmoDA , that can reduce the false positive detection rate caused by correlated features. cosmoDA uses penalized generalized score matching for parsimonious model fitting. We show on simulated benchmarks that cosmoDA can accurately estimate feature interactions in the presence of population heterogeneity and significantly reduces the false discovery rate when testing for differential abundance of correlated features. Using single‐cell and amplicon data, we illustrate cosmoDA 's ability to estimate data‐adaptive Box–Cox‐type data transformations and assess the impact of zero replacement and power transformations on downstream differential abundance results. cosmoDA is available at https://github.com/bio‐datascience/cosmoDA . 10.1002/sim.70534 http://creativecommons.org/licenses/by/4.0/
title Score Matching for Differential Abundance Testing of Compositional High‐Throughput Sequencing Data
topic Statistics in Medicine
url https://onlinelibrary.wiley.com/doi/10.1002/sim.70534