Chromatin accessibility analysis and architectural profiling of human kidneys reveal key cell types and a regulator of diabetic kidney disease

Diabetes is the leading cause of kidney disease that progresses to kidney failure. However, the key molecular and cellular pathways involved in diabetic kidney disease (DKD) pathogenesis are largely unknown. Here, we performed a comparative analysis of adult human kidneys by examining cell type-specific chromatin accessibility by single-nucleus ATAC-seq (snATAC-seq) and analyzing three-dimensional chromatin architecture via high-throughput chromosome conformation capture (Hi-C method) of paired samples. We mapped the cell type-specific and DKD-specific open chromatin landscape and found that genetic variants associated with kidney diseases were significantly enriched in the proximal tubule- (PT) and injured PT-specific open chromatin regions in samples from patients with DKD. BACH1 was identified as a core transcription factor of injured PT cells; its binding target genes were highly associated with fibrosis and inflammation, which were also key features of injured PT cells. Additionally, Hi-C analysis revealed global chromatin architectural changes in DKD, accompanied by changes in local open chromatin patterns. Combining the snATAC-seq and Hi-C data identified direct target genes of BACH1, and indicated that BACH1 binding regions showed increased chromatin contact frequency with promoters of their target genes in DKD. Thus, our multi-omics analysis revealed BACH1 target genes in injured PTs and highlighted the role of BACH1 as a novel regulator of tubular inflammation and fibrosis. Diabetes is the leading cause of kidney disease that progresses to kidney failure. However, the key molecular and cellular pathways involved in diabetic kidney disease (DKD) pathogenesis are largely unknown. Here, we performed a comparative analysis of adult human kidneys by examining cell type-specific chromatin accessibility by single-nucleus ATAC-seq (snATAC-seq) and analyzing three-dimensional chromatin architecture via high-throughput chromosome conformation capture (Hi-C method) of paired samples. We mapped the cell type-specific and DKD-specific open chromatin landscape and found that genetic variants associated with kidney diseases were significantly enriched in the proximal tubule- (PT) and injured PT-specific open chromatin regions in samples from patients with DKD. BACH1 was identified as a core transcription factor of injured PT cells; its binding target genes were highly associated with fibrosis and inflammation, which were also key features of injured PT cells. Additionally, Hi-C analysis revealed global chromatin architectural changes in DKD, accompanied by changes in local open chromatin patterns. Combining the snATAC-seq and Hi-C data identified direct target genes of BACH1, and indicated that BACH1 binding regions showed increased chromatin contact frequency with promoters of their target genes in DKD. Thus, our multi-omics analysis revealed BACH1 target genes in injured PTs and highlighted the role of BACH1 as a novel regulator of tubular inflammation and fibrosis.