Intercellular communication analysis of the human retinal pigment epithelial and choroidal cells predicts pathways associated with aging, cellular senescence and age-related macular degeneration

The retinal pigment epithelium (RPE) and the choroid are ocular tissues with fundamental roles in supporting neuroretinal function. The pathogenesis of age-related macular degeneration (AMD), a leading cause of irreversible blindness for which aging is the highest risk factor is closely linked with progressive impairment of various functions of these tissues. Cellular senescence, marked by cell cycle arrest and secretion of proinflammatory factors, is known to be associated with aging and has been proposed as a potential driver of AMD. Here, we investigated the role played by intercellular communication in the RPE/choroid within the context of aging, senescence and AMD. We inferred cell–cell interactions in the RPE/choroid by applying CellChat and scDiffCom on a publicly available scRNA-seq dataset from three human donors with and without AMD. We identified age-regulated ligand and receptor genes by using limma on a separate publicly available bulk microarray dataset providing RPE/choroid samples at multiple time points. Cellular senescence was investigated by assigning a score to each cell and each sample of these scRNA-seq and microarray datasets, respectively, based on the expression of key signature genes determined by a previous senescence meta-analysis. We identified VEGF-, BMP-and tenascin-mediated pathways supporting some of the strongest cell–cell interactions between RPE cells, fibroblasts and choroidal endothelial cells and as strong intercellular communication pathways related to both aging and senescence. Their signaling strength was enhanced between subpopulations of cells having high senescence scores. Predominant ligands of these pathways were upregulated with age whereas predominant receptors were downregulated. Globally, we also observed that cells from AMD samples presented slightly bigger senescence scores than normal cells and that the senescence score positively correlated with age in bulk samples (R = 0.26, value of p