Multiplexed single-cell imaging reveals diverging subpopulations with distinct senescence phenotypes during long-term senescence induction.

Cellular senescence is a phenotypic state that contributes to the progression of age-related disease through secretion of pro-inflammatory factors known as the senescence-associated secretory phenotype (SASP). Understanding the process by which healthy cells become senescent and develop SASP factors is critical for improving the identification of senescent cells and, ultimately, understanding tissue dysfunction. Here, we reveal how the duration of cellular stress modulates the SASP in distinct subpopulations of senescent cells. We used multiplex, single-cell imaging to build a proteomic map of senescence induction in human epithelial cells induced to senescence over the course of 31 days. We map how the expression of SASP proteins increases alongside other known senescence markers such as p53, p21, and p16. The aggregated population of cells responded to etoposide with an accumulation of stress response factors over the first 11 days, followed by a plateau in most proteins. At the single-cell level, however, we identified two distinct senescence cell populations, one defined primarily by larger nuclear area and the second by higher protein concentrations. Trajectory inference suggested that cells took one of two discrete molecular paths from unperturbed healthy cells, through a common transitional subpopulation, and ending at the discrete terminal senescence phenotypes. Our results underscore the importance of using single-cell proteomics to identify the mechanistic pathways governing the transition from senescence induction to a mature state of senescence characterized by the SASP.