Abstract Background and Aims The burden of cardiovascular disease (CVD) in patients with chronic kidney disease (CKD) remains very high. Endothelial cell (EC) dysfunction is an important key player in the initiation and progression of both pathologies. More importantly, endothelial dysfunction occurs early in both disease courses and is potentially reversible, making it a promising target for effective preventive strategies. Recently, accumulating single-cell RNA sequencing-based evidence has highlighted EC heterogeneity in different pathophysiological contexts. However, while biological processes such as systemic inflammation, oxidative stress fibrosis are being defined as comm EC in CKD and CVD, a comprehensive understanding of at single-cell resolution is lacking. Here, we characterized renal EC subpopulations by single-nucleus RNA sequencing (snRNA-seq) to investigate alterations in CKD patients healthy controls and identify potential drivers of EC dysfunction in CKD. Method Droplet-based snRNA-seq (10X Genomics Chromium) was performed on cryopreserved kidney biopsy cores from 6 human patients with CKD (eGFR < 45 mL/min/1.73 m2) and 3 healthy controls. CKD patients included 2 patients with primary focal segmental glomerulosclerosis (FSGS) and 4 with secondary/maladaptive FSGS (2 patients with diabetes mellitus, 2 female patients in total); all patients had a Framingham 10-year cardiovascular risk score > 20, but no established CVD. Healthy control samples were pre-perfusion biopsies from patients undergoing living kidney donation (including 2 female patients). EC subclusters were annotated based on the expression of canonical EC marker genes. Differential gene expression analysis for patient samples vs. controls included all ECs and DEGs were filtered for log2FC > 0.25 or < −0.25 and adjusted P-value of < .05. Gene set enrichment analysis (GSEA) on the identified DEGs was performed using GO terms. Results A total of 4 081 ECs were present in the dataset, consisting of 5 EC subtypes in both the patient and control cohorts (Panel A). DEGs in CKD patients vs. controls included upregulated genes related to organogenesis and cell differentiation (MEF2C, MEIS2), vascular development (FLRT2), and angiogenesis regulators (VEGFC, PRCP, HDAC9) (Panel B). Notably, VEGFC and FLRT2, are involved in directing EC motility and vascular patterning. In accordance, GSEA showed enrichment for GO biological processes related to vascular endothelial growth factor signaling pathway, endothelial cell differentiation, endothelial development, regulation of angiogenesis, and endothelial cell migration (Panel C). This suggests a potential link between these genes and aberrant angiogenesis, analogous to previous findings in diabetic nephropathy. Some but not all of these changes were also reported for renal ECs of early diabetic nephropathy patients. Of note, the majority of patients studied in our cohort were non-diabetic, warranting further investigation into renal endothelial cell transcriptional alterations that are general to renal disease or specific to particular renal conditions such as diabetic nephropathy. Differential gene expression and GSEA analysis of the EC subclusters in CKD indicated involvement of the peritubular capillary and ascending vasa recta ECs in EC migration / angiogenesis-related responses, overlapping with the overall analysis described above. Glomerular capillary ECs showed upregulation of immune-related genes and processes, reflecting the inflammatory status of the kidney disease. Conclusion snRNA-seq analysis revealed heterogeneity and subtype-specific gene expression alterations of renal ECs from CKD patients compared to healthy controls. Differential functions include enhanced angiogenic (potentially compensatory) and immune-related responses. These findings suggest that specific biological functions may direct different roles of the renal EC subtypes in CKD and will be further validated for their cluster-specificity and for their therapeutic relevance.
Read full abstract