1. Nicotine exacerbates DN in mice
To examine the effects of Nicotine and hyperglycemia on kidney injuries, we first determined albuminuria and BUN in treated mice. The Results showed that Nicotine alone did not increase albuminuria or BUN, whereas hyperglycemia significantly increased both (Fig. 1A, B). Interestingly, we found that the combination of Nicotine and hyperglycemia further increased albuminuria and BUN compared with hyperglycemia alone (Fig. 1A, B). We then detected the expression of Kidney Injury Molecule-1 (KIM-1) and Neutrophil gelatinase-associated lipocalin (NGAL), two biomarkers for kidney injury, by real-time PCR. The results showed that Nicotine alone did not increase the expression of either gene, while its combination with hyperglycemia significantly increased the expression of both genes compared to hyperglycemia alone (Fig. 1C, D).
We performed PAS staining to quantify the kidney glycogen accumulation in control and experimental mice in response to Nicotine and hyperglycemia. Nicotine or hyperglycemia alone did not increase the PAS staining in kidney tissues; in contrast, the combination of Nicotine and hyperglycemia significantly increased the PAS staining in both glomerular and tubular tissue, suggesting Nicotine and hyperglycemia synergistically increased kidney glycogen content (Fig. 2). We also examined apoptosis caused by Nicotine and hyperglycemia by using TUNEL staining. We found the combination of Nicotine and hyperglycemia-induced apoptosis in a more significant number of kidney cells than in either treatment alone (Fig. 3).
Nephrin is a vital slit diaphragm (SD) molecule, a major constituent of the glomerular filtration barrier. A decreased nephrin expression is associated with podocyte injury and kidney dysfunction (Jourdan, et al. 2018; Poulaki, et al. 2020). We performed real-time PCR and immunofluorescence staining to examine the effects of Nicotine and hyperglycemia on nephrin expression. Both modalities showed that the combined treatment (Nicotine and hyperglycemia) decreased nephrin protein and mRNA expressions compared to the treatment with Nicotine or hyperglycemia alone (Fig. 4).
Taken together, these results demonstrated that Nicotine exacerbated DN in mice.
2. Transcriptome changes induced by nicotine-exacerbated DN.
To uncover the underlying molecular mechanisms for nicotine-exacerbated DN, we performed RNA-seq analysis to compare the transcriptomes among the treatments of control (Con), nicotine (Nic), diabetic mellitus (DM) group, and combination of Nicotine and diabetic mellitus (Nic + DM). We randomly selected three mice from each group and extracted total RNAs from their kidneys to construct sequencing libraries. We obtained about 2.95 G of clean data with an average mapping rate of 95.76% per sample to the reference genome (data not shown). A total of 20110 genes were identified in all samples (Table S1). To confirm the effects of different treatments on gene expression in mouse kidneys, we analyzed the gene expression between the experiment groups (Nic, DM, and Nic + DM) and the control group. Compared with the control group, 4171 differentially expressed genes (DEGs) were identified in the Nic + DM group, of which 2063 genes were down-regulated and 2180 were up-regulated. (Fig. 5A, Table S2). To explore the genes strongly reflected in the process of nicotine promoting DN damage, the parameters were selected as log2FC > = 5 (up-regulated) or log2FC<=-5 (down-regulated) but -log10 (p-value) > 5 at the same time. We found that five up-regulated genes (Cyp4a12b, Lcn2, Gbp10, Gbp8, Ubd, Grem1) and eight down-regulated genes (Gm6300, Plin1, Cfd, Car3, Kap, Nat8l, Lrrc15, and Cacna1i) were strongly responsive to the Nicotine promoting DN damage (Fig. 5A).
We also analyzed DEGs in the Nic group vs. Con group and the DM group vs. Con group, respectively. Compared with the control group, we identified that 1610 genes were down-regulated and 1632 genes were up-regulated in the Nic group, 755 genes were down-regulated, and 936 genes were up-regulated in the DM group (Figure S1, Table S3, S4). Here, we found that the down-regulated gene Gm6300 and up-regulated genes Ubd, Gm14391, and Art2a-ps strongly responded to hyperglycemia. These genes presented as a differential expression in the DM group vs. the Con group. The Gm6300 and Ubd genes were also obviously differentially expressed in the Nic + DM group vs. Con group. These two genes may be specific genes involved in regulating hyperglycemia damaging kidneys. Moreover, we found that some genes related to oxidative stress were differentially expressed. For example, CYBB and Ncf1, which promote ROS production, were up-regulated, while SOD1, which inhibits oxidative stress, was significantly down-regulated in the Nic + DM group vs. the Con group (Table S2).
3. Special gene group responded in nicotine-exacerbated DN mice.
We performed a Fuzzy C-means clustering analysis to figure out the genes playing important roles in Nicotine and hyperglycemia-induced kidney damage. Our results showed that 1082 genes clustered into four unique expression patterns in 4 groups of samples (Fig. 5C, Table S5). Interestingly, 259 genes were clustered into cluster 2, and the expression patterns of these genes presented an exciting phenomenon (Table S5). The overall expression trend of these genes is that the expression levels of the first three groups (Con, Nic, and DM) are relatively close, but there is a sudden up-regulation trend in the fourth group (Nic + DM) (Fig. 5C). The up-regulated genes in Nic + DM vs. Con and the fold change (FC) were much larger than those in DM vs. Con; candidate genes were selected for Fuzzy C means clustering analysis. These 259 genes of cluster 2 may be a potential hub genome that specifically increases and responds to nicotine-induced DN kidney damage nicotine to promote DN damage.
Further, we performed the Venn analysis to identify genes involved in Nicotine promoting kidney damage in DN. The selected 50 genes had the highest up-regulation levels, were in cluster 2, and were placed in both the Nic and the DM groups. Venn diagram showed seven common genes (Ubd, Saa1, Grem1, Gbp8, Lcn2, Mat1, Chil3) among the four treatment groups. These genes might play critical roles in the process of Nicotine and high glucose-induced kidney damage (Fig. 5B).
4. The Grem1 may be a vital gene in the process of Nicotine exacerbating mouse DN.
To identify genes that play a vital role in Nicotine-exacerbated DN in mice, we first screened the fold changes (FC) of the six common genes in the Nic + DM group compared to the Con group. Our results showed that the genes Ubd, Saa1, and Grem1were the top three most up-regulated genes (Fig. 6A). We also analyzed the expression trends of 6 common genes in the three comparison groups, including Nic + DM vs. Con, Nic + DM vs. Nic, and Nic + DM vs. DM. Surprisingly, only Grem1 of the top three genes (Ubd, Saa1, and Grem1) was up-regulated in the three comparison groups (adopted p value < 0.05, Table 2). Furthermore, our data showed that Grem1 (log2FC = 2.854042328) and Grem2 (log2FC = 2.508648596) genes were also up-regulated in the DM vs Con group (Table S4), but the difference in up-regulation was not large. However, Grem1 (log2FC = 5.731750645) was significantly up-regulated than Grem2 (log2FC = 3.655715778) in the DM + Nic vs Con group (Table S2).
Table 2
The fold changes of identified seven differentially expressed genes in three groups
|
|
Nic + DM vs. Con
|
Nic + DM vs. Nic
|
Nic + DM vs. DM
|
Gene id
|
Gene name
|
log2FC
|
p-value
|
log2FC
|
p-value
|
log2FC
|
p-value
|
24108
|
Ubd
|
6.4330023
|
2.14E-08
|
5.568934739
|
3.42102E-06
|
1.405613845
|
0.260412665
|
20208
|
Saa1
|
6.191799165
|
0.000046
|
3.271862017
|
0.009054783
|
1.43047504
|
0.159121216
|
23892
|
Grem1
|
5.731750645
|
0.000000226
|
5.001079469
|
4.89209E-07
|
2.883274271
|
0.002156905
|
76074
|
Gbp8
|
5.54548219
|
0.000000125
|
3.928375351
|
6.72766E-05
|
1.457002581
|
0.158508517
|
16819
|
Lcn2
|
5.40400181
|
5.93E-12
|
5.181678867
|
1.61987E-11
|
2.290049156
|
0.02696154
|
11720
|
Mat1a
|
5.043812868
|
0.003280271
|
2.448511208
|
0.094333698
|
1.361990523
|
0.47284171
|
12655
|
Chil3
|
4.947779857
|
0.0000217
|
3.158978799
|
0.000422072
|
2.538680922
|
0.012031274
|
Grem1 is a critical member of the TGF-beta signaling pathway(Marquez-Exposito, et al. 2020), so we suspect that the TGF-beta signaling pathway involves nicotine-promoting kidney damage in DN. We analyzed the DEGs in the Nic + DM vs. Con group to test this hypothesis. We found that 24 DEGs were mapped to the TGF-beta pathway, including 12 down-regulated and 12 up-regulated genes (Fig. 6B). It is worth noting that Grem1 was still the gene with the most significant fold change among the 24 DEGs (Fig. 6B). Notably, the Grem2 gene was the second most up-regulated gene mapped to the TGF-beta pathway. These results suggest that the 24 DEGs represented by Grem1 in the TGF-beta signaling pathway were critical in regulating Nicotine to promote kidney damage in DN. Accordingly, the TGF-beta pathway is likely to be one of the core pathways involved in the nicotine-exacerbated DN.
To confirm the expression of Grem1 in nicotine-exacerbated DN, Real-time PCR and Western blot were performed to examine its expression in mouse kidneys. Our results showed that Nicotine had no significant effect on Grem1 expression in mice with blood glucose levels less than 200 dg/ml; however, it significantly increased Grem1 expression in the kidneys of diabetic mice (Fig. 6C).
In addition, we further confirmed the expression of Grem1 in mouse kidneys through immunofluorescent staining. As displayed in Fig. 7, nicotine increased Grem1 expression in the glomerular and tubular cells under the condition of hyperglycemia. Considering the importance of podocyte for kidney function, we also examined Grem1 expression in podocyte by immunofluorescence staining. The results showed that some Grem1-expressing cells also expressed nephrin, a molecular marker for podocyte, indicating that nicotine could increase Grem1 expression in podocyte (Fig. 8).
5. The phosphorylation of the Smad signaling pathway, being downstream of the TGF-beta pathway, responds to the Nicotine and exacerbates hyperglycemia-induced Grem1 increasing
Grem1 plays its function through 3 signaling pathways: 1) activation of Smad2/3; 2) suppression of Smad1/5/8; 3) activation of VEGFR2. We also examined the effects of Nicotine and hyperglycemia on these downstream signaling pathways. Western blotting and IHC staining showed that the combination of Nicotine and hyperglycemia significantly increased the activation of Smad2/3 and decreased the activation of Smad1/5/8 (Fig. 9A-D) compared with Nicotine or hyperglycemia alone.
To further confirm the activity changes of Smads, we examined the mRNA of Id1 and Id4, the targeting genes for Smad1/5/8, and Snail, the targeting gene for Smad2/3. The results showed that the combination of Nicotine and hyperglycemia significantly decreased the mRNA levels of Id1 and Id4 but increased the level of Snail (Fig. 9E). At the same time, our RNA-seq data also showed that the Id1 and Id4 were down-regulated, and Snail was up-regulated in Nic + DM vs. Con group (Table S2). These data suggest that Nicotine promotes activation of Smad2/3 but suppresses Smad1/5/8 in hyperglycemia.
6. Knocking down Gram1 partially attenuated Nicotine and high glucose-induced podocyte injury
To examine the role of Grem1 in nicotine-exacerbated DN, we treated human podocytes with high glucose and Nicotine. The results showed that high glucose and nicotine co-treatment significantly increased Grem1 expression (Fig. 10A). Cleaved caspase 3, a biomarker of cell apoptosis, was also raised with Nicotine and high glucose; however, the expression of nephrin was significantly decreased (Fig. 10A).
To establish the causation between Grem1 expression and podocyte injury, we knocked down Grem1 expression with its specific siRNA. The results showed that knocking down Grem1 expression significantly reduced cleaved caspase 3 but partially restored nephrin expression (Fig. 10B). These results demonstrated that Grem1 played a vital role in nicotine-exacerbated podocyte injury.