Trajectory inference has radically enhanced single-cell RNA-seq research by enabling the study of dynamic changes in gene expression. Downstream of trajectory inference, it is vital to discover genes that are (i) associated with the lineages in the trajectory, or (ii) differentially expressed between lineages, to illuminate the underlying biological processes. Current data analysis procedures, however, either fail to exploit the continuous resolution provided by trajectory inference, or fail to pinpoint the exact types of differential expression. We introduce tradeSeq, a powerful generalized additive model framework based on the negative binomial distribution that allows flexible inference of both within-lineage and between-lineage differential expression. By incorporating observation-level weights, the model additionally allows to account for zero inflation. We evaluate the method on simulated datasets and on real datasets from droplet-based and full-length protocols, and show that it yields biological insights through a clear interpretation of the data.

A critical step in the analysis of large genome-wide gene expression datasets is the use of module detection methods to group genes into co-expression modules. Because of limitations of classical clustering methods, numerous alternative module detection methods have been proposed, which improve upon clustering by handling co-expression in only a subset of samples, modelling the regulatory network, and/or allowing overlap between modules. In this study we use known regulatory networks to do a comprehensive and robust evaluation of these different methods. Overall, decomposition methods outperform all other strategies, while we do not find a clear advantage of biclustering and network inference-based approaches on large gene expression datasets. Using our evaluation workflow, we also investigate several practical aspects of module detection, such as parameter estimation and the use of alternative similarity measures, and conclude with recommendations for the further development of these methods.

Heterogeneity between different macrophage populations has become a defining feature of this lineage. However, the conserved factors defining macrophages remain largely unknown. The transcription factor ZEB2 is best described for its role in epithelial to mesenchymal transition; however, its role within the immune system is only now being elucidated. We show here that Zeb2 expression is a conserved feature of macrophages. Using Clec4f-cre, Itgax-cre, and Fcgr1-cre mice to target five different macrophage populations, we found that loss of ZEB2 resulted in macrophage disappearance from the tissues, coupled with their subsequent replenishment from bone-marrow precursors in open niches. Mechanistically, we found that ZEB2 functioned to maintain the tissue-specific identities of macrophages. In Kupffer cells, ZEB2 achieved this by regulating expression of the transcription factor LXRα, removal of which recapitulated the loss of Kupffer cell identity and disappearance. Thus, ZEB2 expression is required in macrophages to preserve their tissue-specific identities.