Multi-scalar data integration decoding risk genes for chronic kidney disease

BMC Nephrol. 2024 Oct 18;25(1):364. doi: 10.1186/s12882-024-03798-2.

Abstract

Background: Chronic Kidney Disease (CKD) impacts over 10% of the global population, and recent advancements in high-throughput analytical technologies are uncovering the complex physiology underlying this condition. By integrating Genome-Wide Association Studies (GWAS), RNA sequencing (RNA-seq/RNA array), and single-cell RNA sequencing (scRNA-seq) data, our study aimed to explore the genes and cell types relevant to CKD traits.

Methods: GWAS summary data for end-stage renal failure (ESRD) and decreased eGFR (CKD) with or without diabetes and (micro)proteinuria were obtained from the GWAS Catalog and the UK Biobank (UKB) database. Two gene Expression Omnibus (GEO) transcriptome datasets were used to establish glomerular and tubular gene expression differences between CKD patients and healthy individuals. Two scRNA-seq datasets were utilized to obtain the expression of key genes at the single-cell level. The expression profile, differentially expressed genes (DEGs), gene-gene interaction, and pathway enrichment were analysed for these CKD risk genes.

Results: A total of 779 distinct SNPs were identified from GWAS across different CKD traits, involving 681 genes. While many of these risk genes are specific to the CKD traits of renal failure, decreased eGFR, and (micro)proteinuria, they share common pathways, including extracellular matrix (ECM). ECM modeling was enriched in upregulated glomerular and tubular DEGs from CKD kidneys compared to healthy controls, with the expression of relevant collagen genes, such as COL1A2, prevalent in fibroblasts/myofibroblasts. Additionally, immune responses, including T cell differentiation, were dysregulated in CKD kidneys. The late podocyte signature gene THSD7A was enriched in podocytes but downregulated in CKD. We also highlighted that the regulated risk genes of CKD are mainly expressed in tubular cells and immune cells in the kidney.

Conclusions: Our integrated analysis highlight the genes, pathways, and relevant cell types associational with the pathogenesis of kidney traits, as a basis for further mechanistic studies to understand the pathogenesis of CKD.

Keywords: CKD; GWAS; Integrated analysis; RNA-seq; scRNA-seq.

MeSH terms

  • Genetic Predisposition to Disease
  • Genome-Wide Association Study*
  • Glomerular Filtration Rate
  • Humans
  • Kidney Failure, Chronic / genetics
  • Polymorphism, Single Nucleotide
  • Renal Insufficiency, Chronic* / genetics
  • Sequence Analysis, RNA
  • Single-Cell Analysis
  • Transcriptome