KDM6A showed a promotive activity for regulating the EMT in HK-2 cells
Bearing in mind the importance of epigenetic modification in the development of kidney disease, we firstly investigated the expression and significance of KDM6A in the HK-2 cells cultured in high glucose (HG) and mannitol (Mann) medium. First of all, we confirmed the change of HK-2 migration in the HG group, which is closely associated with epithelial-mesenchymal transition(28), HG-induced HK-2 cells showed a significant deteriorative migration relative to Mann groups, suggesting HK-2 cells displayed severe epithelial-mesenchymal transition when cultured in HG medium(Fig. 1A). The EMT maker, E-cadherin and TGFβ, were detected by immunoblotting and qPCR. The immunoblotting results showed an increased tendency of E-cadherin and TGFβ expression in HG-induced HK-2 cells comparing to Mann-induced group (Fig. 1B), consistent with the qPCR results (Fig. 1C). Interestingly, the expression of KDM6A in the HG group was also elevated both in transcription and translation level (Fig. 1B-C), indicating that upregulation of KDM6A may be associated with the EMT progress of HK-2 cell induced by HG.
To confirm the hypothesis, GSK-J4, a kind of KDM6A inhibitor, was utilized to treat HK-2 cells in HG medium with negative control (NC and GSK-J4). The migration detected by wound healing assay was decreased while the HK-2 cells were treated with GSK-J4 (Fig. D). The lower expression of E-cadherin and TGFβ was observed by western bolt and qPCR in the GSK-J4 groups, expounding the inhibition of KDM6A could receded the EMT degree in HG-induced HK-2 (Fig. E, F).
Those results illustrated that KDM6A might be able to propel the progression of epithelial-mesenchymal transition in renal tubular epithelial cells.
KDM6A exhibited a pivotal role in kidney injury of DN
Basing on the in vitro studies, we next investigated the in vivo role of KDM6A in diabetes kidney of DN using a STZ-induced mouse model (WT-NC groups and WT-STZ groups). We indeed found a higher level of KDM6A was expressed in kidney tissues of the 12-week STZ-induced diabetic mice relative to control groups, both measured by qPCR and Immunohistochemistry (IHC) (Fig. 2A-C). To further explore the influence of KDM6A on the progression of diabetes kidney injury, GSK-J4 was applied to treat diabetic mice with negative control (WT-STZ-NC and WT-STZ-G groups). The injury of diabetic kidney was evaluated through the morphology of renal proximal tubule and renal distal tubule which were associated with tubulointerstitial fibrosis. Obviously, the HE and PAS stain results in kidney of STZ-induced mice showed renal proximal tubular epithelial cells expansion and proximal tubule lumen reduction relative to normal mice, while the degree of proximal tubule injury in diabetic mice was dramatically decreased after treated with GKS-J4 (Fig. 2D, F). Similarly, severe fibrosis of renal distal tubule in diabetic mice was monitored and the progress of fibrosis was retarded in the presence of GKS-J4 (Fig. 2E, G). In those GSK-J4-induced or non-GSK-J4-induced diabetic mice, several typical signs of diabetic nephropathy, such as albumin to creatinine ratio (ACR), blood urea nitrogen (BUN), creatinine (Cr), proteinuria, was measured. All of those sighs were aggravated in STZ-induced mice and restored partly by GSK-J4 treatment (Fig. 2H-K), which were closely consistent with HE and PAS stain results. Those results indicated that the downregulation of KMD6A might alleviate the kidney injury of diabetic nephropathy mice.
Interestingly, the IHC and its quantitative results displayed a suppressed expression of E-cadherin, a maker of EMT and fibrosis process(14), and an elevated expression of TGF-β, another maker of EMT and fibrosis process(6), whereas both of them were altered by GSK-J4 (Fig. 2L-O). Treated diabetic mice with GSK-J4, the depressed E-cadherin was recovered and the increased TGF-β was downregulated, totally conforming to the previous results.
Together, the expression of KDM6A was improved in the kidney tissues of STZ-induced diabetic nephropathy mice and the inhibition of KDM6A resulted in a mitigative kidney injury, suggesting that KDM6A could promote the progress of diabetic nephropathy.
KDM6A was the upstream regulator of E-cadherin for repressing the its expression
Tons of reports have showed that repressing the epithelial cadherin (E-cadherin) could trigger the EMT and fibrosis directly (29, 30). Depending on the previous work, we assumed that KDM6A might be the upstream effector of E-cadherin and could inhibit the E-cadherin expression, resulting in EMT and even fibrosis of DN. To testify our assumption, HK-2 cells stably overexpressed E-cadherin (HK2-vehicle and HK2-E-cadherin groups) was constructed. Overexpressing E-cadherin could inhibit the EMT process and reduced the migration of HK-2 cells comparing to HK2-vehicle groups when cells were cultured in HG medium (Fig. 3A). However, as the increasing of E-cadherin in HK2-E-cadherin groups, the expression of KDM6A showed no difference in both transcription and translation level between HK2-vehicle and HK2-E-cadherin groups (Fig. 3B, C). It seemed that KDM6A expression might be not affected by overexpressing E-cadherin.
Furthermore, the silencing RNA, targeting E-cadherin (si-E-cadherin), was applied to knockdown the E-cadherin in HK-2 cells as well as negative control (si-NC) when cells were cultured in mannitol medium. In this condition, the migration of HK-2 was improved dramatically with the E-cadherin knockdown (Fig. 3D). Notably, the expression of KDM6A was either still in a stable state even if the EMT was aggravated by E-cadherin knockdown (Fig. 3E, F). All those results showed the expression of KDM6A would not be impacted only by variation of E-cadherin expression, suggesting that KDM6A was the upstream regulator of E-cadherin.
KDM6A was the target of miR-199b-3p during the renal injury
Considering the importance of miRNA in development of diabetic nephropathy, we wanted to determine whether miRNA participated in the progression of DN by targeting KDM6A mRNA. For mining potential miRNA candidates targeting KDM6A, we predicted possible miRNA by several databases, such as PITA, RNA22, miRmap, microT, miRanda, PicTar, and TargetScan, Finally, we found 8 miRNAs, hsa-miR-19a-3p, hsa-miR-19b-3p, hsa-miR-23a-3p, hsa-miR-199a-3p, hsa-miR-23b-3p, hsa-miR-199b-3p, hsa-miR-758-3p, hsa-miR-142-3p. Each of miRNAs were exhibited from at least 5 databases (Table 1). Among of them, hsa-miR-142-3p, hsa-miR-19a-3p, hsa-miR-19b-3p, hsa-miR-23a-3p and hsa-miR-23b-3p have been proved for binding to KDM6A in kidney (26, 31–33). Additionally, the fold-change of hsa-miR-199a-3p, hsa-miR-199b-3p, hsa-miR-758-3p between normal kidney and kidney with DN was analyzed by GEO DataSets. The results of GSE51674, a dataset acquired from the progression of kidney fibrosis in human Diabetic Nephropathy, showed noticeable change of hsa-miR-199a-3p, hsa-miR-199b-3p expression, but no difference in hsa-miR-758-3p expression (Fig. S1A, B, Table S1). The change of hsa-miR-199b-3p expression was remarkably outstanding compared with other 2 miRNAs.
Table 1
The predictive results of miRNA targeting KDM6A
miRNA ID | miRNA name | Databeses | Numbers of datebase |
PITA | RNA22 | miRmap | microT | miRanda | PicTar | TargetScan |
MIMAT0000073 | hsa-miR-19a-3p | √ | × | × | √ | √ | √ | √ | 5 |
MIMAT0000074 | hsa-miR-19b-3p | √ | × | × | √ | √ | √ | √ | 5 |
MIMAT0000078 | hsa-miR-23a-3p | √ | × | × | √ | √ | √ | √ | 5 |
MIMAT0000232 | hsa-miR-199a-3p | √ | × | √ | √ | × | √ | √ | 5 |
MIMAT0000418 | hsa-miR-23b-3p | √ | × | × | √ | √ | √ | √ | 5 |
MIMAT0000434 | hsa-miR-142-3p | √ | × | √ | √ | √ | √ | √ | 6 |
MIMAT0003879 | hsa-miR-758-3p | √ | × | √ | √ | √ | √ | √ | 6 |
MIMAT0004563 | hsa-miR-199b-3p | √ | × | √ | √ | × | √ | √ | 5 |
To confirm the results, qPCR assay was performed for testing those miRNA expression in HG-induced HK-2 cells and the results showed a similar trend with precious work (Fig. 4A). We also detected the miRNA expression in diabetes kidney using a STZ-induced mouse model. The results was compatible with the consequence in HG-treated HK-2 cell (Fig. 4B).
Considering the greatest change of miR-199b-3p, we assumed that miR-199b-3p could bind to KDM6A mRNA. According to the predicted alignment of miR-199b-3p binding to KDM6A mRNA (Fig. 4C), WT and mutant KDM6A 3’UTR luciferase plasmids as well as miR-199b-3p mimics and mi-NC were designed for dual-luciferase reporter assay in order to confirm the prediction. Comparing with cells co-transfected with WT KDM6A 3’UTR and mi-NC, cells co-transfected with WT KDM6A 3’UTR and miR-199b-3p mimics displayed a conspicuous reduction of luciferase activity. However, when the WT KDM6A 3’UTR was replaced by mutant KDM6A 3’UTR, no luciferase activity change could be observed regardless of co-transfecting with miR-199b-3p mimics or mi-NC (Fig. D).
Those results suggested that miR-199b-3p expression was changed in the EMT of renal tubular epithelial cells and could target to KDM6A, which means miR-199b-3p may play a regulatory role in tubulointerstitial fibrosis of DN through downregulate KDM6A.
MiR-199b-3p showed a protective effect against epithelial-mesenchymal transition of renal tubular epithelial cells through downregulating KDM6A
For further explore the role of miR-199b-3p in the process of EMT in renal tubular epithelial cells, we transfected the HK-2 cells with lentivirus containing miR-199b-3p mimics or miR-NC and selected the stable cell line (HK2-miRNA and HK2-miNC) and cultured them in HG medium. The expression of miR-199b-3p in both cell lines was tested (Fig. S2A). As the increasing expression of miR-199b-3p in HK2-miRNA, the expression of KDM6A, measured by western blot, was inhibited in HK2-miRNA cells (Fig. 5A) with a slightly lower transcription level of KDM6A detected by the qPCR (Fig. S2A). The expression of TGFβ and E-cadherin was also tested. Decreased expression of TGFβ was observed, whereas the expression of E-cadherin was upregulated relative to HK2-miNC cells (Fig. 5A). The similar tendency was also monitored in transcription level (Fig. S2B). Furthermore, the HK2-miNC cells in HG showed a more serious EMT than HK2-miRNA cells in HG with the deteriorative migration, suggesting that miR-199b-3p could halt HG-induce EMT of HK-2 cells dramatically (Fig. 5B). Those results indicated that miR-199b-3p might be able to downregulate the KDM6A and further protected HK-2 cell from the HG-induced injury.
In order to confirm the influence of miR-199b-3p on HK-2 cells through downregulating KDM6A, we transfected HK2-miRNA cells with a miR-199b-3p inhibitor (inhibitor) or negative control (inhibitor NC) to determine whether inhibition of miR-199b-3p could exacerbate HG-induced EMT. The qPCR results showed the successful inhibition of miR-199b-3p because of the inhibitor, accompanying with the restored expression of KDM6A (Fig S2. C). Consistently, the maker of EMT, TGFβ and E-cadherin, showed an increased and decreased tendency respectively, no matter in mRNA level or protein level (Fig. 5F, Fig. S2D). As we supposed, the expression of KDM6A in translation level was certainly restored due to the inhibitor existence (Fig. 5C). In the presence of inhibitor, migration of HK2-miRNA cells were increased, which indicated HK2-miRNA cells were preferred to mesenchymal-like performance compared with inhibitor NC group (Fig. 5D). Those results suggested that the restoration of KDM6A due to repression of miR-199b-3p could aggravate the EMT process.
Together, miR-199b-3p could prevent the degeneration of EMT in renal tubular epithelial cells through suppressing KDM6A expression.
MiR-199b-3p could inhibit the kidney injury caused by diabetic nephropathy by regulating KDM6A expression
In order to further confirm the function of miR-199b-3p in vivo, we induced the diabetic nephropathy in miR-199b-3p deficient mice (miR-199b-3p−/−, Fig. S3A) using STZ as well as negative control (miR-199b-3p−/−-NC and miR-199b-3p−/−-STZ). Several typical signs of DN were measured in different groups and all of them showed a similar tendency. The level of ACR, BUN, Cr and proteinuria in miR-199b-3p−/−-STZ groups were significantly greater than WT-STZ groups, indicating the protective function was removed due to the miR-199b-3p deficiency (Fig. 6A-D). Furthermore, the pathological features of STZ-induced diabetic kidney in WT and miR-199b-3p−/− mice were observed through HE and PAS stain. Also in the induction of STZ, a more serious renal proximal tubular epithelial cells expansion as well as a more severe proximal tubule lumen reduction was emerged in miR-199b-3p−/− mouse compared to WT mouse (Fig. 6E, G). Consistently, more severe fibrosis of renal distal tubule was monitored due to the miR-199b-3p deficiency (Fig. 6F, H). Comparing with WT-STZ groups, the IHC and qPCR results in miR-199b-3p−/−-STZ groups showed an increased expression of TGFβ and a decreased expression of E-cadherin, which was consistent with pathological results (Fig. 6I-J, L-M, Fig. S4B). However, the higher level of KDM6A was expressed in diabetic kidney tissues without miR-199b-3p existence, no matter in transcription or translation level (Fig. 6K, N, Fig. S4B).
Summarily, based on those in vivo results, we evidently found the absence of miR-199b-3p was closely associated with the exacerbation of kidney injury in STZ-induced diabetic mice by abolishing the repression of KDM6A, indicating that miR-199b-3p played a protective prole in kidney injury of diabetic nephropathy through downregulating KDM6A.