To conduct the cluster differential abundance test, edgeR (3.38.1)71 was used. The input for edgeR consisted of cluster counts across all conditions and replicates, which were subsequently subjected to “estimateDisp”, “glmQLFit”, and “glmQLFTest“72. Significance was determined based on a false discovery rate threshold of <0.05 for each test.
