Integrative ATAC-seq and RNA-seq analysis associated with diabetic nephropathy and identification of novel targets for treatment by dapagliflozin

Dapagliflozin (DAPA) are clinically effective in improving diabetic nephropathy (DN). However, whether and how chromatin accessibility changed by DN responds to DAPA treatment is unclear. Therefore, we performed ATAC-seq, RNA-seq, and weighted gene correlation network analysis to identify the chromatin accessibility, the messenger RNA (mRNA) expression, and the correlation between clinical phenotypes and mRNA expression using kidney from three mouse groups: db/m mice (Controls), db/db mice (case group), and those treated with DAPA (treatment group). RNA-Seq and ATAC-seq conjoint analysis revealed many overlapping pathways and networks suggesting that the transcriptional changes of DN and DAPA intervention largely occured dependently on chromatin remodeling. Specifically, the results showed that some key signal transduction pathways, such as immune dysfunction, glucolipid metabolism, oxidative stress and xenobiotic and endobiotic metabolism, were repeatedly enriched in the analysis of the RNA-seq data alone, as well as combined analysis with ATAC-seq data. Furthermore, we identified some candidate genes (UDP glucuronosyltransferase 1 family, Dock2, Tbc1d10c, etc.) and transcriptional regulators (KLF6 and GFI1) that might be associated with DN and DAPA restoration. These reversed genes and regulators confirmed that pathways related to immune response and metabolism pathways were critically involved in DN progression. We first reported the restoration action of DAPA on the transcriptome and chromatin accessibility landscape changes caused by DN. Furthermore, identified hub genes, transcriptional regulators, and pathways could be novel biomarkers or targeted for adjunctive therapy.