Fig. 1: Heterogeneous response of NPC43 cells to lytic induction treatment.

A UMAP representation of NPC43 cells at three distinct treatment time points, highlighting distinct clustering patterns. B UMAP plots illustrating the cellular states at different treatment time points. C Heatmap depicting the correlation between pseudo-bulks of different cellular states at various treatment time points, based on the top 1000 variable genes. Heatmap was hierarchically ordered. D Violin plots showing the activation levels of the EGFR, MAPK, PI3K, and NF-κB pathways in NPC43 cells. Pathway activity scores were inferred using PROGENy. Wilcoxon signed-rank test was performed. E Dot plot illustrating the expression levels of markers associated with different cellular states, demonstrating their conservation across treatment time points. F Bar plots showing the percentile distribution of cellular states at different treatment time points. Composition analysis was performed with scCODA, setting FDR to 0.05, and using cycling cells as the reference cell type. The proportions of NR cells, ECM-related cells (2), and Differentiated cells showed credible changes. G Scatter plots depicting the activity scores of Gene Ontologies (GOs) differentially activated in cellular states. H Scatter plot highlighting the most specific Differentially Expressed genes (DEs) of NR cells at T48. I Dot plot showcasing the expression levels of surface markers of NR cells in different cellular states, indicating conservation across different treatment time points.