Single-Cell Transcriptomics Reveals Endothelial Plasticity During Diabetic Atherogenesis

Atherosclerosis is the leading cause of cardiovascular diseases, which is also the primary cause of mortality among diabetic patients. Endothelial cell (EC) dysfunction is a critical early step in the development of atherosclerosis and aggravated in the presence of concurrent diabetes. Although the heterogeneity of the organ-specific ECs has been systematically analyzed at the single-cell level in healthy conditions, their transcriptomic changes in diabetic atherosclerosis remain largely unexplored. Here, we carried out a single-cell RNA sequencing (scRNA-seq) study using EC-enriched single cells from mouse heart and aorta after 12 weeks feeding of a standard chow or a diabetogenic high-fat diet with cholesterol. We identified eight EC clusters, three of which expressed mesenchymal markers, indicative of an endothelial-to-mesenchymal transition (EndMT). Analyses of the marker genes, pathways, and biological functions revealed that ECs are highly heterogeneous and plastic both in normal and atherosclerotic conditions. The metabolic transcriptomic analysis further confirmed that EndMT-derived fibroblast-like cells are prominent in atherosclerosis, with diminished fatty acid oxidation and enhanced biological functions, including regulation of extracellular-matrix organization and apoptosis. In summary, our data characterized the phenotypic and metabolic heterogeneity of ECs in diabetes-associated atherogenesis at the single-cell level and paves the way for a deeper understanding of endothelial cell biology and EC-related cardiovascular diseases.

Single-Cell Transcriptome Analysis Reveals Six Subpopulations Reflecting Distinct Cellular Fates in Senescent Mouse Embryonic Fibroblasts

Replicative senescence is a hallmark of aging, which also contributes to individual aging. Mouse embryonic fibroblasts (MEFs) provide a convenient replicative senescence model. However, the heterogeneity of single MEFs during cellular senescence has remained unclear. Here, we conducted single-cell RNA sequencing on senescent MEFs. Principal component analysis showed obvious heterogeneity among these MEFs such that they could be divided into six subpopulations. Three types of gene expression analysis revealed distinct expression features of these six subpopulations. Trajectory analysis revealed three distinct lineages during MEF senescence. In the main lineage, some senescence-associated secretory phenotypes were upregulated in a subset of cells from senescent clusters, which could not be distinguished in a previous bulk study. In the other two lineages, a possibility of escape from cell cycle arrest and coupling between translation-related genes and ATP synthesis-related genes were also discovered. Additionally, we found co-expression of transcription factor HOXD8 coding gene and its potential target genes in the main lineage. Overexpression of Hoxd8 led to senescence-associated phenotypes, suggesting HOXD8 is a new regulator of MEF senescence. Together, our single-cell sequencing on senescent MEFs largely expanded the knowledge of a basic cell model for aging research.

Intratumoral Heterogeneity: More Than Just Mutations

Most human tumors are composed of genetically and phenotypically heterogeneous cancer cell populations, which poses a major challenge for the clinical management of cancer patients. Advances of single-cell technologies have allowed the profiling of tumors at unprecedented depth, which, in combination with newly developed computational tools, enable the dissection of tumor evolution with increasing precision. However, our understanding of mechanisms that regulate intratumoral heterogeneity and our ability to modulate it has been lagging behind. Recent data demonstrate that epigenetic regulators, including histone demethylases, may control the cell-to-cell variability of transcriptomes and chromatin profiles and they may modulate therapeutic responses via this function. Thus, the therapeutic targeting of epigenetic enzymes may be used to decrease intratumoral cellular heterogeneity and treatment resistance, when used in combination with other types of agents.

KDM5 Histone Demethylase Activity Links Cellular Transcriptomic Heterogeneity to Therapeutic Resistance

Members of the KDM5 histone H3 lysine 4 demethylase family are associated with therapeutic resistance, including endocrine resistance in breast cancer, but the underlying mechanism is poorly defined. Here we show that genetic deletion of KDM5A/B or inhibition of KDM5 activity increases sensitivity to anti-estrogens by modulating estrogen receptor (ER) signaling and by decreasing cellular transcriptomic heterogeneity. Higher KDM5B expression levels are associated with higher transcriptomic heterogeneity and poor prognosis in ER⁺ breast tumors. Single-cell RNA sequencing, cellular barcoding, and mathematical modeling demonstrate that endocrine resistance is due to selection for pre-existing genetically distinct cells, while KDM5 inhibitor resistance is acquired. Our findings highlight the importance of cellular phenotypic heterogeneity in therapeutic resistance and identify KDM5A/B as key regulators of this process.